The pathway of ligand entry from the membrane bilayer to a lipid G protein-coupled receptor

Computational Methods for Modeling Lipid-Mediated Active Pharmaceutical Ingredient Delivery

Lipid-mediated delivery of active pharmaceutical ingredients (API) opened new possibilities in advanced therapies. By encapsulating an API into a lipid nanocarrier (LNC), one can safely deliver APIs not soluble in water, those with otherwise strong adverse effects, or very fragile ones such as nucleic acids. However, for the rational design of LNCs, a detailed understanding of the composition-structure-function relationships is missing. This review presents currently available computational methods for LNC investigation, screening, and design. The state-of-the-art physics-based approaches are described, with the focus on molecular dynamics simulations in all-atom and coarse-grained resolution. Their strengths and weaknesses are discussed, highlighting the aspects necessary for obtaining reliable results in the simulations. Furthermore, a machine learning, i.e., data-based learning, approach to the design of lipid-mediated API delivery is introduced. The data produced by the experimental and theoretical approaches provide valuable insights. Processing these data can help optimize the design of LNCs for better performance. In the final section of this Review, state-of-the-art of computer simulations of LNCs are reviewed, specifically addressing the compatibility of experimental and computational insights.

Exploring Hypertension: The Role of AT1 Receptors, Sartans, and Lipid Bilayers

The rational design of AT1 receptor antagonists represents a pivotal approach in the development of therapeutic agents targeting cardiovascular pathophysiology. Sartans, a class of compounds engineered to inhibit the binding and activation of Angiotensin II on the AT1 receptor, have demonstrated significant clinical efficacy. This review explores the multifaceted role of sartans in mitigating hypertension and related complications. We highlight the integration of crystallography, computational simulations, and NMR spectroscopy to elucidate sartan-AT1 receptor interactions, providing a foundation for the next-generation antagonist design. The review also delves into the challenges posed by the high lipophilicity and suboptimal bioavailability of sartans, emphasizing advancements in nanotechnology and novel drug delivery systems. Additionally, we discuss the impact of lipid bilayers on the AT1 receptor conformation and drug binding, underscoring the importance of the lipidic environment in receptor-drug interactions. We suggest that optimizing drug design to account for these factors could enhance the therapeutic potential of AT1 receptor antagonists, paving the way for improved cardiovascular health outcomes.

Homodimerization of CB2 cannabinoid receptor triggered by a bivalent ligand enhances cellular signaling

G protein-coupled receptors (GPCRs) exist within a landscape of interconvertible conformational states and in dynamic equilibrium between monomers and higher-order oligomers, both influenced by ligand binding. Here, we show that a homobivalent ligand formed by equal chromenopyrazole moieties as pharmacophores, connected by 14 methylene units, can modulate the dynamics of the cannabinoid CB2 receptor (CB2R) homodimerization by simultaneously binding both protomers of the CB2R-CB2R homodimer. Computational and pharmacological experiments showed that one of the ligand pharmacophores binds to the orthosteric site of one protomer, and the other pharmacophore to a membrane-oriented pocket between transmembranes 1 and 7 of the partner protomer. This results in unique pharmacological properties, including increased potency in Gi-mediated signaling and enhanced recruitment of β-arrestin. Thus, by modulating dimerization dynamics, it may be possible to fine-tune CB2R activity, potentially leading to improved therapeutic outcomes.

Metadynamics Study of Lipid-Mediated Antibacterial Toxin Binding to the EmrE Multiefflux Protein

EmrE is a bacterial efflux protein in the small multidrug-resistant (SMR) family present in Escherichia coli. Due to its small size, 110 residues in each dimer subunit, it is an ideal model system to study ligand-protein-membrane interactions. Here in our work, we have calculated the free energy landscape of benzyltrimetylammonium (BTMA) and tetraphenyl phosphonium (TPP) binding to EmrE using the enhanced sampling method-multiple walker metadynamics. We estimate that the free energy of BTMA binding to EmrE is -21.2 ± 3.3 kJ/mol and for TPP is -43.6 ± 3.8 kJ/mol. BTMA passes through two metastable states to reach the binding pocket, while TPP has a more complex binding landscape with four metastable states and one main binding site. Our simulations show that the ligands interact with the membrane lipids at a distance 1 nm away from the binding site which forms a broad local minimum, consistent for both BTMA and TPP. This site can be an alternate entry point for ligands to partition from the membrane into the protein, especially for bulky and/or branched ligands. We also observed the membrane lipid and C-terminal 110HisA form salt-bridge interactions with the helix-1 residue 22LysB. Our free energy estimates and clusters are in close agreement with experimental data and give us an atomistic view of the ligand-protein-lipid interactions. Understanding the binding pathway of these ligands can guide us in future design of ligands that can alter or halt the function of EmrE.

Cannabidiol at Nanomolar Concentrations Negatively Affects Signaling through the Adenosine A2A Receptor

Cannabidiol (CBD) is a phytocannabinoid with potential as a therapy for a variety of diseases. CBD may act via cannabinoid receptors but also via other G-protein-coupled receptors (GPCRs), including the adenosine A2A receptor. Homogenous binding and signaling assays in Chinese hamster ovary (CHO) cells expressing the human version of the A2A receptor were performed to address the effect of CBD on receptor functionality. CBD was not able to compete for the binding of a SCH 442416 derivative labeled with a red emitting fluorescent probe that is a selective antagonist that binds to the orthosteric site of the receptor. However, CBD reduced the effect of the selective A2A receptor agonist, CGS 21680, on Gs-coupling and on the activation of the mitogen activated kinase signaling pathway. It is suggested that CBD is a negative allosteric modulator of the A2A receptor.

Ligand entry pathways control the chemical space recognized by GPR183

The G protein-coupled receptor GPR183 is a chemotactic receptor with an important function in the immune system and association with a variety of diseases. It recognizes ligands with diverse physicochemical properties as both the endogenous oxysterol ligand 7α,25-OHC and synthetic molecules can activate the G protein pathway of the receptor. To better understand the ligand promiscuity of GPR183, we utilized both molecular dynamics simulations and cell-based validation experiments. Our work reveals that the receptor possesses two ligand entry channels: one lateral between transmembrane helices 4 and 5 facing the membrane, and one facing the extracellular environment. Using enhanced sampling, we provide a detailed structural model of 7α,25-OHC entry through the lateral membrane channel. Importantly, the first ligand recognition point at the receptor surface has been captured in diverse experimentally solved structures of different GPCRs. The proposed ligand binding pathway is supported by in vitro data employing GPR183 mutants with a sterically blocked lateral entrance, which display diminished binding and signaling. In addition, computer simulations and experimental validation confirm the existence of a polar water channel which might serve as an alternative entrance gate for less lipophilic ligands from the extracellular milieu. Our study reveals knowledge to understand GPR183 functionality and ligand recognition with implications for the development of drugs for this receptor. Beyond, our work provides insights into a general mechanism GPCRs may use to respond to chemically diverse ligands.

The Leu/Val6.51 Side Chain of Cannabinoid Receptors Regulates the Binding Mode of the Alkyl Chain of Δ9-Tetrahydrocannabinol

(-)-Δ9-trans-tetrahydrocannabinol (THC), which is the principal psychoactive constituent of Cannabis, mediates its action by binding to two members of the G-protein-coupled receptor (GPCR) family: the cannabinoid CB1 (CB1R) and CB2 (CB2R) receptors. Molecular dynamics simulations showed that the pentyl chain of THC could adopts an I-shape conformation, filling an intracellular cavity between Phe3.36 and Trp6.48 for initial agonist-induced receptor activation, in CB1R but not in CB2R. This cavity opens to the five-carbon chain of THC by the conformational change of the γ-branched, flexible, Leu6.51 side chain of CB1R, which is not feasible by the β-branched, mode rigid, Val6.51 side chain of CB2R. In agreement with our computational results, THC could not decrease the forskolin-induced cAMP levels in cells expressing mutant CB1RL6.51V receptor but could activate the mutant CB2RV6.51L receptor as efficiently as wild-type CB1R. Additionally, JWH-133, a full CB2R agonist, contains a branched dimethyl moiety in the ligand chain that bridges Phe3.36 and Val6.51 for receptor activation. In this case, the substitution of Val6.51 to Leu in CB2R makes JWH-133 unable to activate CB2RV6.51L. In conclusion, our combined computational and experimental results have shown that the amino acid at position 6.51 is a key additional player in the initial mechanism of activation of GPCRs that recognize signaling molecules derived from lipid species.

Structural basis of selective cannabinoid CB2 receptor activation

Cannabinoid CB₂ receptor (CB₂R) agonists are investigated as therapeutic agents in the clinic. However, their molecular mode-of-action is not fully understood. Here, we report the discovery of LEI-102, a CB₂R agonist, used in conjunction with three other CBR ligands (APD371, HU308, and CP55,940) to investigate the selective CB₂R activation by binding kinetics, site-directed mutagenesis, and cryo-EM studies. We identify key residues for CB₂R activation. Highly lipophilic HU308 and the endocannabinoids, but not the more polar LEI-102, APD371, and CP55,940, reach the binding pocket through a membrane channel in TM1-TM7. Favorable physico-chemical properties of LEI-102 enable oral efficacy in a chemotherapy-induced nephropathy model. This study delineates the molecular mechanism of CB₂R activation by selective agonists and highlights the role of lipophilicity in CB₂R engagement. This may have implications for GPCR drug design and sheds light on their activation by endogenous ligands.

A Single Point Mutation Blocks the Entrance of Ligands to the Cannabinoid CB2 Receptor via the Lipid Bilayer

Molecular dynamic (MD) simulations have become a common tool to study the pathway of ligand entry to the orthosteric binding site of G protein-coupled receptors. Here, we have combined MD simulations and site-directed mutagenesis to study the binding process of the potent JWH-133 agonist to the cannabinoid CB₂ receptor (CB₂R). In CB₂R, the N-terminus and extracellular loop 2 fold over the ligand binding pocket, blocking access to the binding cavity from the extracellular environment. We, thus, hypothesized that the binding pathway is a multistage process consisting of the hydrophobic ligand diffusing in the lipid bilayer to contact a lipid-facing vestibule, from which the ligand enters an allosteric site inside the transmembrane bundle through a tunnel formed between TMs 1 and 7 and finally moving from the allosteric to the orthosteric binding cavity. This pathway was experimentally validated by the Ala282^7.36Phe mutation that blocks the entrance of the ligand, as JWH-133 was not able to decrease the forskolin-induced cAMP levels in cells expressing the mutant receptor. This proposed ligand entry pathway defines transient binding sites that are potential cavities for the design of synthetic modulators.

Medium-Chain Lipid Conjugation Facilitates Cell-Permeability and Bioactivity

The majority of bioactive molecules act on membrane proteins or intracellular targets and therefore needs to partition into or cross biological membranes. Natural products often exhibit lipid modifications to facilitate critical molecule-membrane interactions, and in many cases their bioactivity is markedly reduced upon removal of a lipid group. However, despite its importance in nature, lipid-conjugation of small molecules is not commonly used in chemical biology and medicinal chemistry, and the effect of such conjugation has not been systematically studied. To understand the composition of lipids found in natural products, we carried out a chemoinformatic characterization of the "natural product lipidome". According to this analysis, lipidated natural products predominantly contain saturated medium-chain lipids (MCLs), which are significantly shorter than the long-chain lipids (LCLs) found in membranes and lipidated proteins. To study the usefulness of such modifications in probe design, we systematically explored the effect of lipid conjugation on five different small molecule chemotypes and find that permeability, cellular retention, subcellular localization, and bioactivity can be significantly modulated depending on the type of lipid tail used. We demonstrate that MCL conjugation can render molecules cell-permeable and modulate their bioactivity. With all explored chemotypes, MCL-conjugates consistently exhibited superior uptake or bioactivity compared to LCL-conjugates and either comparable or superior uptake or bioactivity to short-chain lipid (SCL)-conjugates. Together, our findings suggest that conjugation of small molecules with MCLs could be a powerful strategy for the design of probes and drugs.

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.

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.

Fly casting with ligand sliding and orientational selection supporting complex formation of a GPCR and a middle sized flexible molecule

A GA-guided multidimensional virtual-system coupled molecular dynamics (GA-mD-VcMD) simulation was conducted to elucidate binding mechanisms of a middle-sized flexible molecule, bosentan, to a GPCR protein, human endothelin receptor type B (hETB). GA-mD-VcMD is a generalized ensemble method that produces a free-energy landscape of the ligand-receptor binding by searching large-scale motions accompanied with stable maintenance of the fragile cell-membrane structure. All molecular components (bosentan, hETB, membrane, and solvent) were represented with an all-atom model. Then sampling was conducted from conformations where bosentan was distant from the binding site in the hETB binding pocket. The deepest basin in the resultant free-energy landscape was assigned to native-like complex conformation. The following binding mechanism was inferred. First, bosentan fluctuating randomly in solution is captured using a tip region of the flexible N-terminal tail of hETB via nonspecific attractive interactions (fly casting). Bosentan then slides occasionally from the tip to the root of the N-terminal tail (ligand–sliding). During this sliding, bosentan passes the gate of the binding pocket from outside to inside of the pocket with an accompanying rapid reduction of the molecular orientational variety of bosentan (orientational selection). Last, in the pocket, ligand–receptor attractive native contacts are formed. Eventually, the native-like complex is completed. The bosentan-captured conformations by the tip-region and root-region of the N-terminal tail correspond to two basins in the free-energy landscape. The ligand-sliding corresponds to overcoming of a free-energy barrier between the basins.

Impact of Inert Crowders on Host-Guest Recognition Process

Biological environments typically contain high concentrations (300-400 mg/mL) of different macromolecules at volume fractions as large as 30%-40%. Biomolecular recognition processes, a ubiquitous biological phenomena, occurring in such crowded heterogeneous media would differ significantly compared to the dilute buffer solutions. Here we quantify the potential impact of inert crowders on prototypical host-guest recognition process by explicit-solvent molecular dynamics (MD) simulations in atomic resolution. We demonstrate that the crowders, when smaller in size, would facilitate the binding process of the guest molecule by decreasing the free energy barrier for binding via excluded volume effect and desolvation of the host receptor. However, the extent of crowder-induced stabilization of a host-guest complex is found to be significantly higher when the guest molecule is sterically constricted to approach the host along a centrosymmetric direction, compared to its unrestricted, freely diffusive movement. A kinetic analysis of the recognition process reveals that the origin of a relatively stronger crowder impact during constricted movement of guest molecule lies in the relatively enhanced residence time of the guest inside the host by crowders. Together, our results suggest that the extent of impact of crowding on recognition processes would be contingent upon the presence or absence of constriction on ligand movement.

Covalent cannabinoid receptor ligands - structural insight and selectivity challenges

X-ray crystallography and cryogenic electronic microscopy have provided significant advancement in the knowledge of GPCR structure and have allowed the rational design of GPCR ligands. The class A GPCRs cannabinoid receptor type 1 and type 2 are implicated in many pathophysiological processes and thus rational design of drug and tool compounds is of great interest. Recent structural insight into cannabinoid receptors has already led to a greater understanding of ligand binding sites and receptor residues that likely contribute to ligand selectivity. Herein, classes of heterocyclic covalent cannabinoid receptor ligands are reviewed in light of the recent advances in structural knowledge of cannabinoid receptors, with particular discussion regarding covalent ligand selectivity and rationale design.

Differences in ligand-induced protein dynamics extracted from an unsupervised deep learning approach correlate with protein-ligand binding affinities

Prediction of protein–ligand binding affinity is a major goal in drug discovery. Generally, free energy gap is calculated between two states (e.g., ligand binding and unbinding). The energy gap implicitly includes the effects of changes in protein dynamics induced by ligand binding. However, the relationship between protein dynamics and binding affinity remains unclear. Here, we propose a method that represents ligand-binding-induced protein behavioral change with a simple feature that can be used to predict protein–ligand affinity. From unbiased molecular simulation data, an unsupervised deep learning method measures the differences in protein dynamics at a ligand-binding site depending on the bound ligands. A dimension reduction method extracts a dynamic feature that strongly correlates to the binding affinities. Moreover, the residues that play important roles in protein–ligand interactions are specified based on their contribution to the differences. These results indicate the potential for binding dynamics-based drug discovery.

Differences in ligand-induced protein dynamics extracted as a single feature from a deep learning-based analysis of MD simulations correlate with ligand binding affinity.

Two entry tunnels in mouse TAAR9 suggest the possibility of multi-entry tunnels in olfactory receptors

Orthosteric binding sites of olfactory receptors have been well understood for ligand-receptor interactions. However, a lack of explanation for subtle differences in ligand profile of olfactory receptors even with similar orthosteric binding sites promotes more exploration into the entry tunnels of the receptors. An important question regarding entry tunnels is the number of entry tunnels, which was previously believed to be one. Here, we used TAAR9 that recognizes important biogenic amines such as cadaverine, spermine, and spermidine as a model for entry tunnel study. We identified two entry tunnels in TAAR9 and described the residues that form the tunnels. In addition, we found two vestibular binding pockets, each located in one tunnel. We further confirmed the function of two tunnels through site-directed mutagenesis. Our study challenged the existing views regarding the number of entry tunnels in the subfamily of olfactory receptors and demonstrated the possible mechanism how the entry tunnels function in odorant recognition.

Recent Advances in Structure, Function, and Pharmacology of Class A Lipid GPCRs: Opportunities and Challenges for Drug Discovery

Great progress has been made over the past decade in understanding the structural, functional, and pharmacological diversity of lipid GPCRs. From the first determination of the crystal structure of bovine rhodopsin in 2000, much progress has been made in the field of GPCR structural biology. The extraordinary progress in structural biology and pharmacology of GPCRs, coupled with rapid advances in computational approaches to study receptor dynamics and receptor-ligand interactions, has broadened our comprehension of the structural and functional facets of the receptor family members and has helped usher in a modern age of structure-based drug design and development. First, we provide a primer on lipid mediators and lipid GPCRs and their role in physiology and diseases as well as their value as drug targets. Second, we summarize the current advancements in the understanding of structural features of lipid GPCRs, such as the structural variation of their extracellular domains, diversity of their orthosteric and allosteric ligand binding sites, and molecular mechanisms of ligand binding. Third, we close by collating the emerging paradigms and opportunities in targeting lipid GPCRs, including a brief discussion on current strategies, challenges, and the future outlook.

Characterisation of bis(4-aminoquinoline)s as α1A adrenoceptor allosteric modulators

The development of sub-type selective α₁ adrenoceptor ligands has been hampered by the high sequence similarity of the amino acids forming the orthosteric binding pocket of the three α₁ adrenoceptor subtypes, along with other biogenic amine receptors. One possible approach to overcome this issue is to target allosteric sites on the α₁ adrenoceptors. Previous docking studies suggested that one of the quinoline moieties of a bis(4-aminoquinoline), comprising a 9-carbon methylene linker attached via the amine groups, could interact with residues outside of the orthosteric binding site while, simultaneously, the other quinoline moiety bound within the orthosteric site. We therefore hypothesized that this compound could act in a bitopic manner, displaying both orthosteric and allosteric binding properties. To test this proposition, we investigated the allosteric activity of a series of bis(4-aminoquinoline)s with linker lengths ranging from 2 to 12 methylene units (designated C2-C12). A linear trend of increasing [³H]prazosin dissociation rate with increasing linker length between C7 and C11 was observed, confirming their action as allosteric modulators. These data suggest that the optimal linker length for the bis(4-aminoquinoline)s to occupy the allosteric site of the α_1A adrenoceptor is between 7 and 11 methylene units. In addition, the ability of C9 bis(4-aminoquinoline) to modulate the activation of the α_1A adrenoceptor by norepinephrine was subsequently examined, showing that C9 acted as a non-competitive antagonist. Our findings indicate that the bis(4-aminoquinolines) are acting as allosteric modulators of orthosteric ligand binding, but not efficacy, in a bitopic manner.

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.

Design of Negative and Positive Allosteric Modulators of the Cannabinoid CB2 Receptor Derived from the Natural Product Cannabidiol

Cannabidiol (CBD), the second most abundant of the active compounds found in the Cannabis sativa plant, is of increasing interest because it is approved for human use and is neither euphorizing nor addictive. Here, we design and synthesize novel compounds taking into account that CBD is both a partial agonist, when it binds to the orthosteric site, and a negative allosteric modulator, when it binds to the allosteric site of the cannabinoid CB₂ receptor. Molecular dynamic simulations and site-directed mutagenesis studies have identified the allosteric site near the receptor entrance. This knowledge has permitted to perform structure-guided design of negative and positive allosteric modulators of the CB₂ receptor with potential therapeutic utility.

A boosted unbiased molecular dynamics method for predicting ligands binding mechanisms: probing the binding pathway of dasatinib to Src-kinase

SUMMARY

Small molecules such as metabolites and drugs play essential roles in biological processes and pharmaceutical industry. Knowing their interactions with biomacromolecular targets demands a deep understanding of binding mechanisms. Dozens of papers have suggested that discovering of the binding event by means of conventional unbiased molecular dynamics (MD) simulation urges considerable amount of computational resources, therefore, only one who holds a cluster or a supercomputer can afford such extensive simulations. Thus, many researchers who do not own such resources are reluctant to take the benefits of running unbiased MD simulation, in full atomistic details, when studying a ligand binding pathway. Many researchers are impelled to be content with biased MD simulations which seek its validation due to its intrinsic preconceived framework. In this work, we have presented a workable stratagem to encourage everyone to perform unbiased (unguided) MD simulations, in this case a protein-ligand binding process, by typical desktop computers and so achieve valuable results in nanosecond time scale. Here, we have described a dynamical binding's process of an anticancer drug, the dasatinib, to the c-Src kinase in full atomistic details for the first time, without applying any biasing force or potential which may lead the drug to artificial interactions with the protein. We have attained multiple independent binding events which occurred in the nanosecond time scales, surprisingly as little as ∼30 ns. Both the protonated and deprotonated forms of the dasatinib reached the crystallographic binding mode without having any major intermediate state during induction.

AVAILABILITY AND IMPLEMENTATION

The links of the tutorial and technical documents are accessible in the article.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Markov State Models to Elucidate Ligand Binding Mechanism

Molecular dynamics simulations can now routinely access the microsecond timescale, making feasible direct sampling of ligand association events. While Markov State Model (MSM) approaches offer a useful framework for analyzing such trajectory data to gain insight into binding mechanisms, accurate modeling of ligand association pathways and kinetics must be done carefully. We describe methods and good practices for constructing MSMs of ligand binding from unbiased trajectory data and discuss how to use time-lagged independent component analysis (tICA) to build informative models, using as an example recent simulation work to model the binding of phenylalanine to the regulatory ACT domain dimer of phenylalanine hydroxylase. We describe a variety of methods for estimating association rates from MSMs and discuss how to distinguish between conformational selection and induced-fit mechanisms using MSMs. In addition, we review some examples of MSMs constructed to elucidate the mechanisms by which p53 transactivation domain (TAD) and related peptides bind the oncoprotein MDM2.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Discovery of Homobivalent Bitopic Ligands of the Cannabinoid CB2 Receptor*

Single chemical entities with potential to simultaneously interact with two binding sites are emerging strategies in medicinal chemistry. We have designed, synthesized and functionally characterized the first bitopic ligands for the CB2 receptor. These compounds selectively target CB2 versus CB1 receptors. Their binding mode was studied by molecular dynamic simulations and site-directed mutagenesis.

Addressing free fatty acid receptor 1 (FFAR1) activation using supervised molecular dynamics

The free fatty acid receptor 1 (FFAR1, formerly GPR40), is a potential G protein-coupled receptor (GPCR) target for the treatment of type 2 diabetes mellitus (T2DM), as it enhances glucose-dependent insulin secretion upon activation by endogenous long-chain free fatty acids. The presence of two allosterically communicating binding sites and the lack of the conserved GPCR structural motifs challenge the general knowledge of its activation mechanism. To date, four X-ray crystal structures are available for computer-aided drug design. In this study, we employed molecular dynamics (MD) and supervised molecular dynamics (SuMD) to deliver insights into the (un)binding mechanism of the agonist MK-8666, and the allosteric communications between the two experimentally determined FFAR1 binding sites. We found that FFAR1 extracellular loop 2 (ECL2) mediates the binding of the partial agonist MK-8666. Moreover, simulations showed that the agonists MK-8666 and AP8 are reciprocally stabilized and that AP8 influences MK-8666 unbinding from FFAR1.

Rational drug repurposing for cancer by inclusion of the unbiased molecular dynamics simulation in the structure-based virtual screening approach: Challenges and breakthroughs

Managing cancer is now one of the biggest concerns of health organizations. Many strategies have been developed in drug discovery pipelines to help rectify this problem and two of the best ones are drug repurposing and computational methods. The combination of these approaches can have immense impact on the course of drug discovery. In silico drug repurposing can significantly reduce the time, the cost and the effort of drug development. Computational methods such as structure-based drug design (SBDD) and virtual screening can predict the potentials of small molecule binders, such as drugs, for having favorable effect on a particular molecular target. However, the demand for accuracy and efficiency of SBDD requires more sophisticated and complicated approaches such as unbiased molecular dynamics (UMD) simulation that has been recently introduced. As a complementary strategy, the knowledge acquired from UMD simulations can increase the chance of finding the right candidates and the pipeline of its administration is introduced and discussed in this review. An elaboration of this pipeline is also made by detailing an example, the binding and unbinding pathways of dasatinib-c-Src kinase complex, which shows that how influential this method can be in rational drug repurposing in cancer treatment.

Membrane Phospholipid Analogues as Molecular Rulers to Probe the Position of the Hydrophobic Contact Point of Lysophospholipid Ligands on the Surface of G-Protein-Coupled Receptor during Membrane Approach

When lipid mediators bind to G-protein-coupled receptors (GPCRs), the ligand first enters the lipid bilayer, then diffuses laterally in the cell membrane to make hydrophobic contact with the receptor protein, and finally enters the receptor's binding pocket. In this process, the location of the hydrophobic contact point on the surface of the receptor has been little discussed even in cases in which the crystal structure has been determined, because the ligand binding pocket is buried inside the transmembrane (TM) domains. Here, we coupled an activator ligand to a series of membrane phospholipid surrogates, which constrain the depth of entry of the ligand into the lipid bilayer. Consequently, via measurement of the receptor-activating activity as a function of the depth of entry into the membrane, these surrogates can be used as molecular rulers to estimate the location of the hydrophobic contact point on the surface of GPCR. We focused on lysophosphatidylserine (LysoPS) receptor GPR34 and prepared a series of simplified membrane-lipid-surrogate-conjugated lysophospholipid analogues by attaching alkoxy amine chains of varying lengths to the hydrophobic tail of a potent GPR34 agonist. As expected, the activity of these lipid-conjugated LysoPS analogues was dependent on chain length. The predicted contact position matches the position of the terminal benzene ring of a nonlipidic ligand that protrudes between TMs 4 and 5 of the receptor. We further found that the nature of the terminal hydrophilic functional group of the conjugated membrane lipid surrogate strongly influences the activity, suggesting that lateral hydrophilic contact of LysoPS analogues with the receptor's surface is also crucial for ligand-GPCR binding.

Monte Carlo simulations using PELE to identify a protein-protein inhibitor binding site and pose

In silico binding site location and pose prediction for a molecule targeted at a large protein surface is a challenging task. We report a blind test with two peptidomimetic molecules that bind the flu virus hemagglutinin (HA) surface antigen, JNJ7918 and JNJ4796 (recently disclosed in van Dongen et al., Science, 2019, 363). Tests with a series of conventional approaches such as rigid (receptor) docking against available X-ray crystal structures or against an ensemble of structures generated by quick methodologies (NMA, homology modeling) gave mixed results, due to the shallowness and flexibility of the binding site and the sheer size of the target. However, tests with our Monte Carlo platform PELE in two protocols involving either exploration of the whole protein surface (global exploration), or the latter followed by refinement of best solutions (local exploration) yielded remarkably good results by locating the actual binding site and generating binding modes that recovered all native contacts found in the X-ray structures. Thus, the Monte Carlo scheme of PELE seems promising as a quick methodology to overcome the challenge of identifying entirely unknown binding sites and modes for protein–protein disruptors.

PELE prospectively unveils the binding site and mode of a protein–protein disruptor.

Exploring the role of the membrane bilayer in the recognition of candesartan by its GPCR AT1 receptor

Cardiovascular diseases and hypertension in particular are major health risks worldwide and the improvement on their treatment will be beneficial for the human health. AT1R antagonists belong to the sartans family that targets the renin-angiotensin aldosterone system (RAAS) through blocking the hormone angiotensin II to exert its detrimental effects in pathological states. As a consequence, they are beneficial to treat hypertension, diabetes related kidney failure and hyperaemic episodes. Long unbiased Molecular Dynamics (MD) simulations are performed in order to explore candesartan's possible 2D and 3D diffusion mechanisms towards AT1R receptor. 3D diffusion mechanism is referred to the direct binding of the AT1 antagonist candesartan to the AT1R 3D structure (PDB ID: 4YAY). 2D diffusion mechanism involves first, the incorporation of candesartan in the bilayer core and then its localization on the AT1R binding cavity, through a diffusion mechanism. The obtained results indicate that membranes interact significantly with the neutral form of candesartan, which is indeed approaching the receptors' active site through diffusion via the lipids. On the other hand, the deprotonated form of the drug is interacting with AT1R's extracellular loop and fails to enter the membrane, pointing out the importance of the pH microenvironment around the receptor. To validate the calculated diffusion coefficients of the drug in the lipid bilayers 2D DOSY NMR experiments were recorded and they were in good agreement. Information on the impact that has the interaction of candesartan with the membrane is very important for the rationally design and development of potent ARBs. Thus, its conformational features as well as its localization in the membrane core have to be thoroughly explored.

Probing the Druggablility on the Interface of the Protein-Protein Interaction and Its Allosteric Regulation Mechanism on the Drug Screening for the CXCR4 Homodimer

Modulating protein–protein interactions (PPIs) with small drug-like molecules targeting it exhibits great promise in modern drug discovery. G protein-coupled receptors (GPCRs) are the largest family of targeted proteins and could form dimers in living biological cells through PPIs. However, compared to drug development of the orthosteric site, there has been lack of investigations on the druggability of the PPI interface for GPCRs and its functional implication on experiments. Thus, in order to address these issues, we constructed a novel computational strategy, which involved in molecular dynamics simulation, virtual screening and protein structure network (PSN), to study one representative GPCR homodimer (CXCR4). One druggable pocket was identified in the PPI interface and one small molecule targeting it was screened, which could strengthen PPI mainly through hydrophobic interaction between the benzene rings of the PPI molecule and TM4 of the receptor. The PSN results further reveals that the PPI molecule could increase the number of the allosteric regulation pathways between the druggable pocket of the dimer interface to the orthostatic site for the subunit A but only play minor role for the other subunit B, leading to the asymmetric change in the volume of the binding pockets for the two subunits (increase for the subunit A and minor change for the subunit B). Consequently, the screening performance of the subunit A to the antagonists is enhanced while the subunit B is unchanged nearly, implying that the PPI molecule may be beneficial to enhance the drug efficacies of the antagonists. In addition, one main regulation pathway with the highest frequency was identified for the subunit A, which consists of Trp195^5.34–Tyr190^ECL2–Val196^5.35–Gln200^5.39–Asp262^6.58–Cys28^N-term, revealing their importance in the allosteric regulation from the PPI molecule. The observations from the work could provide valuable information for the development of the PPI drug-like molecule for GPCRs.

Molecular Mechanism of S1P Binding and Activation of the S1P1 Receptor

Sphingosine-1-phosphate (S1P) is a lipidic mediator in mammals that functions either as a second messenger or as a ligand. In the latter case, it is transported by its HDL-associated apoM carrier and circulated in blood where it binds to specific S1P receptors on cell membranes and induces downstream reactions. Although S1P signaling pathways are essential for many biological processes, they are poorly understood at the molecular level. Here, the solved crystal structures of the S1P₁ receptor were used to evaluate molecular dynamics (MD) simulations to generate greater detailed molecular insights into the mechanism of S1P signaling. The MD simulations provided observations at the coarse-grained and atomic levels indicating that S1P may access the receptor binding pocket directly from solvents. Lifting of the bulky N-terminal cap region of the receptor precedes initial S1P binding. Glu121^3.29 guides S1P penetration, and together with Arg292^7.34 is responsible for the stabilization of S1P in the binding pocket, which is consistent with experimental predictions. The complete binding of S1P is followed by receptor activation, wherein Trp269^6.48 moves toward the transmembrane helix (TM) 7, resulting in the formation of an enhanced hydrogen bond network in the lower region of TM7. The distance between TM3 and TM6 is subsequently increased, resulting in the opening of the intracellular binding pocket that enables G protein binding. Further analysis of the force distribution network in the receptor yielded a detailed molecular understanding of the signal transmission network that is activated upon agonist binding.

The Structural Binding Mode of the Four Autotaxin Inhibitor Types that Differentially Affect Catalytic and Non-Catalytic Functions

Autotaxin (ATX) is a secreted lysophospholipase D, catalysing the conversion of lysophosphatidylcholine (LPC) to bioactive lysophosphatidic acid (LPA). LPA acts through two families of G protein-coupled receptors (GPCRs) controlling key cellular responses, and it is implicated in many physiological processes and pathologies. ATX, therefore, has been established as an important drug target in the pharmaceutical industry. Structural and biochemical studies of ATX have shown that it has a bimetallic nucleophilic catalytic site, a substrate-binding (orthosteric) hydrophobic pocket that accommodates the lipid alkyl chain, and an allosteric tunnel that can accommodate various steroids and LPA. In this review, first, we revisit what is known about ATX-mediated catalysis, crucially in light of allosteric regulation. Then, we present the known ATX catalysis-independent functions, including binding to cell surface integrins and proteoglycans. Next, we analyse all crystal structures of ATX bound to inhibitors and present them based on the four inhibitor types that are established based on the binding to the orthosteric and/or the allosteric site. Finally, in light of these data we discuss how mechanistic differences might differentially modulate the activity of the ATX-LPA signalling axis, and clinical applications including cancer.

The size matters? A computational tool to design bivalent ligands

Motivation

Bivalent ligands are increasingly important such as for targeting G protein-coupled receptor (GPCR) dimers or proteolysis targeting chimeras (PROTACs). They contain two pharmacophoric units that simultaneously bind in their corresponding binding sites, connected with a spacer chain. Here, we report a molecular modelling tool that links the pharmacophore units via the shortest pathway along the receptors van der Waals surface and then scores the solutions providing prioritization for the design of new bivalent ligands.

Results

Bivalent ligands of known dimers of GPCRs, PROTACs and a model bivalent antibody/antigen system were analysed. The tool could rapidly assess the preferred linker length for the different systems and recapitulated the best reported results. In the case of GPCR dimers the results suggest that in some cases these ligands might bind to a secondary binding site at the extracellular entrance (vestibule or allosteric site) instead of the orthosteric binding site.

Availability and implementation

Freely accessible from the Molecular Operating Environment svl exchange server (https://svl.chemcomp.com/).

Supplementary information

Supplementary data are available at Bioinformatics online.

Investigating targets for neuropharmacological intervention by molecular dynamics simulations

Medical research has identified over 500 brain disorders. Among these, there are still only very few neuropathologies whose causes are fully understood and, consequently, very few drugs whose mechanism of action is known. No FDA drug has been identified for major neurodegenerative diseases, such as Alzheimer's and Parkinson's. We still lack effective treatments and strategies for modulating progression or even early neurodegenerative disease onset diagnostic tools. A great support toward the highly needed identification of neuroactive drugs comes from computer simulation methods and, in particular, from molecular dynamics (MD). This provides insight into structure-function relationship of a target and predicts structure, dynamics and energetics of ligand/target complexes under biologically relevant conditions like temperature and physiological saline concentration. Here, we present examples of the predictive power of MD for neuroactive ligands/target complexes. This brief survey from our own research shows the usefulness of partnerships between academia and industry, and from joint efforts between experimental and theoretical groups.

Structural Insights into CB1 Receptor Biased Signaling

The endocannabinoid system has emerged as a promising target for the treatment of numerous diseases, including cancer, neurodegenerative disorders, and metabolic syndromes. Thus far, two cannabinoid receptors, CB1 and CB2, have been discovered, which are found predominantly in the central nervous system (CB1) or the immune system (CB2), among other organs and tissues. CB1 receptor ligands have been shown to induce a complex pattern of intracellular effects. The binding of a ligand induces distinct conformational changes in the receptor, which will eventually translate into distinct intracellular signaling pathways through coupling to specific intracellular effector proteins. These proteins can mediate receptor desensitization, trafficking, or signaling. Ligand specificity and selectivity, complex cellular components, and the concomitant expression of other proteins (which either regulate the CB1 receptor or are regulated by the CB1 receptor) will affect the therapeutic outcome of its targeting. With an increased interest in G protein-coupled receptors (GPCR) research, in-depth studies using mutations, biological assays, and spectroscopic techniques (such as NMR, EPR, MS, FRET, and X-ray crystallography), as well as computational modelling, have begun to reveal a set of concerted structural features in Class A GPCRs which relate to signaling pathways and the mechanisms of ligand-induced activation, deactivation, or activity modulation. This review will focus on the structural features of the CB1 receptor, mutations known to bias its signaling, and reported studies of CB1 receptor ligands to control its specific signaling.

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.

Lipid-Based Inhibitors Act Directly on GlyT2

The endogenous lipids N-arachidonylglycine and oleoyl-l-carnitine are potential therapeutic leads in the treatment of chronic pain through their inhibition of the glycine transporter GlyT2. However, their mechanism of action is unknown. It has been hypothesized that these "bioactive" lipids either inhibit GlyT2 indirectly, by significantly perturbing the biophysical properties of the membrane, or directly, by binding directly to the transporter (either from a membrane-exposed or solvent-exposed binding site). Here, we used molecular dynamics simulations to study the effects of the lipids anandamide, N-arachidonylglycine, and oleoyl-l-carnitine on (a) the biophysical properties of the bilayer and (b) direct binding interactions with GlyT2. During the simulations, the biophysical properties of the bilayer itself, for example, the area per lipid, bilayer thickness, and order parameters, were not significantly altered by the presence or type of bioactive lipid, regardless of the presence of GlyT2. Our work, together with previous computational and experimental data, suggests that these acyl-inhibitors of GlyT2 inhibit the transporter by directly binding to it. However, these bioactive lipids bound to various parts of GlyT2 and did not prefer a single binding site during 4.5 μs of simulation. We postulate that the binding site is located at the solvent-exposed regions of GlyT2. Understanding the mechanism of action of these and related bioactive lipids is essential in effectively developing high-affinity GlyT2 inhibitors for the treatment of pain.

Influence of Cholesterol and Its Stereoisomers on Members of the Serotonin Receptor Family

Despite the ubiquity of cholesterol within the cell membrane, the mechanism by which it influences embedded proteins remains elusive. Numerous G-protein coupled receptors exhibit dramatic responses to membrane cholesterol with regard to the ligand-binding affinity and functional properties, including the 5-HT receptor family. Here, we use over 25 μs of unbiased atomistic molecular dynamics simulations to identify cholesterol interaction sites in the 5-HT_1B and 5-HT_2B receptors and evaluate their impact on receptor structure. Susceptibility to membrane cholesterol is shown to be subtype dependent and determined by the quality of interactions between the extracellular loops. Charged residues are essential for maintaining the arrangement of the extracellular surface in 5-HT_2B; in the absence of such interactions, the extracellular surface of the 5-HT_1B is malleable, populating a number of distinct conformations. Elevated cholesterol density near transmembrane helix 4 is considered to be conducive to the conformation of extracellular loop 2. Occupation of this site is also shown to be stereospecific, illustrated by differential behavior of nat-cholesterol isomers, ent- and epi-cholesterol. In simulations containing the endogenous agonist, serotonin, cholesterol binding at transmembrane helix 4 biases bound serotonin molecules toward an unexpected binding mode in the extended binding pocket. The results highlight the capability of membrane cholesterol to influence the mobility of the extracellular surface in the 5-HT₁ receptor family and manipulate the architecture of the extracellular ligand-binding pocket.

Ligand Binding Mechanisms in Human Cone Visual Pigments

Vertebrate vision starts with light absorption by visual pigments in rod and cone photoreceptor cells of the retina. Rhodopsin, in rod cells, responds to dim light, whereas three types of cone opsins (red, green, and blue) function under bright light and mediate color vision. Cone opsins regenerate with retinal much faster than rhodopsin, but the molecular mechanism of regeneration is still unclear. Recent advances in the area pinpoint transient intermediate opsin conformations, and a possible secondary retinal-binding site, as determinant factors for regeneration. In this Review, we compile previous and recent findings to discuss possible mechanisms of ligand entry in cone opsins, involving a secondary binding site, which may have relevant functional and evolutionary implications.

Mammalian class I odorant receptors exhibit a conserved vestibular-binding pocket

Odorant receptors represent the largest family of mammalian G protein-coupled receptors. Phylogenetically, they are split into two classes (I and II). By analyzing the entire subclass I odorant receptors sequences, we identified two class I-specific and highly conserved motifs. These are predicted to face each other at the extra-cellular portion of the transmembrane domain, forming a vestibular site at the entrance to the orthosteric-binding cavity. Molecular dynamics simulation combined with site-directed mutagenesis and in vitro functional assays confirm the functional role of this vestibular site in ligand-driven activation. Mutations at this part of the receptor differentially affect the receptor response to four agonists. Since this vestibular site is involved in ligand recognition, it could serve ligand design that targets specifically this sub-genome of mammalian odorant receptors.

Peeking at G-protein-coupled receptors through the molecular dynamics keyhole

Molecular dynamics is a state of the art computational tool for the investigation of biophysics phenomenon at a molecular scale, as it enables the modeling of dynamic processes, such as conformational motions, molecular solvation and ligand binding. The recent advances in structural biology have led to a bloom in published G-protein-coupled receptor structures, representing a solid and valuable resource for molecular dynamics studies. During the last decade, indeed, a plethora of physiological and pharmacological facets of this membrane protein superfamily have been addressed by means of molecular dynamics simulations, including the activation mechanism, allosterism and, very recently, biased signaling. Here, we try to recapitulate some of the main contributions that molecular dynamics has recently produced in the field.

How Effectively Can Adaptive Sampling Methods Capture Spontaneous Ligand Binding?

Molecular dynamics (MD) simulations that capture the spontaneous binding of drugs and other ligands to their target proteins can reveal a great deal of useful information, but most drug-like ligands bind on time scales longer than those accessible to individual MD simulations. Adaptive sampling methods―in which one performs multiple rounds of simulation, with the initial conditions of each round based on the results of previous rounds―offer a promising potential solution to this problem. No comprehensive analysis of the performance gains from adaptive sampling is available for ligand binding, however, particularly for protein–ligand systems typical of those encountered in drug discovery. Moreover, most previous work presupposes knowledge of the ligand’s bound pose. Here we outline existing methods for adaptive sampling of the ligand-binding process and introduce several improvements, with a focus on methods that do not require prior knowledge of the binding site or bound pose. We then evaluate these methods by comparing them to traditional, long MD simulations for realistic protein–ligand systems. We find that adaptive sampling simulations typically fail to reach the bound pose more efficiently than traditional MD. However, adaptive sampling identifies multiple potential binding sites more efficiently than traditional MD and also provides better characterization of binding pathways. We explain these results by showing that protein–ligand binding is an example of an exploration–exploitation dilemma. Existing adaptive sampling methods for ligand binding in the absence of a known bound pose vastly favor the broad exploration of protein–ligand space, sometimes failing to sufficiently exploit intermediate states as they are discovered. We suggest potential avenues for future research to address this shortcoming.

'Crystal' Clear? Lysophospholipid Receptor Structure Insights and Controversies

Lysophospholipids (LPLs), particularly sphingosine 1-phosphate (S1P) and lysophosphatidic acid (LPA), are bioactive lipid modulators of cellular homeostasis and pathology. The discovery and characterization of five S1P- and six LPA-specific G protein-coupled receptors (GPCRs), S1P1-5 and LPA1-6, have expanded their known involvement in all mammalian physiological systems. Resolution of the S1P1, LPA1, and LPA6 crystal structures has fueled the growing interest in these receptors and their ligands as targets for pharmacological manipulation. In this review, we have attempted to provide an integrated overview of the three crystallized LPL GPCRs with biochemical and physiological structure-function data. Finally, we provide a novel discussion of how chaperones for LPLs may be considered when extrapolating crystallographic and computational data toward understanding actual biological interactions and phenotypes.

Role of Extracellular Loops and Membrane Lipids for Ligand Recognition in the Neuronal Adenosine Receptor Type 2A: An Enhanced Sampling Simulation Study

Human G-protein coupled receptors (GPCRs) are important targets for pharmaceutical intervention against neurological diseases. Here, we use molecular simulation to investigate the key step in ligand recognition governed by the extracellular domains in the neuronal adenosine receptor type 2A (hA_2AR), a target for neuroprotective compounds. The ligand is the high-affinity antagonist (4-(2-(7-amino-2-(furan-2-yl)-[1,2,4]triazolo[1,5-a][1,3,5]triazin-5-ylamino)ethyl)phenol), embedded in a neuronal membrane mimic environment. Free energy calculations, based on well-tempered metadynamics, reproduce the experimentally measured binding affinity. The results are consistent with the available mutagenesis studies. The calculations identify a vestibular binding site, where lipids molecules can actively participate to stabilize ligand binding. Bioinformatic analyses suggest that such vestibular binding site and, in particular, the second extracellular loop, might drive the ligand toward the orthosteric binding pocket, possibly by allosteric modulation. Taken together, these findings point to a fundamental role of the interaction between extracellular loops and membrane lipids for ligands’ molecular recognition and ligand design in hA_2AR.

Simulations of the regulatory ACT domain of human phenylalanine hydroxylase (PAH) unveil its mechanism of phenylalanine binding

Phenylalanine hydroxylase (PAH) regulates phenylalanine (Phe) levels in mammals to prevent neurotoxicity resulting from high Phe concentrations as observed in genetic disorders leading to hyperphenylalaninemia and phenylketonuria. PAH senses elevated Phe concentrations by transient allosteric Phe binding to a protein-protein interface between ACT domains of different subunits in a PAH tetramer. This interface is present in an activated PAH (A-PAH) tetramer and absent in a resting-state PAH (RS-PAH) tetramer. To investigate this allosteric sensing mechanism, here we used the GROMACS molecular dynamics simulation suite on the Folding@home computing platform to perform extensive molecular simulations and Markov state model (MSM) analysis of Phe binding to ACT domain dimers. These simulations strongly implicated a conformational selection mechanism for Phe association with ACT domain dimers and revealed protein motions that act as a gating mechanism for Phe binding. The MSMs also illuminate a highly mobile hairpin loop, consistent with experimental findings also presented here that the PAH variant L72W does not shift the PAH structural equilibrium toward the activated state. Finally, simulations of ACT domain monomers are presented, in which spontaneous transitions between resting-state and activated conformations are observed, also consistent with a mechanism of conformational selection. These mechanistic details provide detailed insight into the regulation of PAH activation and provide testable hypotheses for the development of new allosteric effectors to correct structural and functional defects in PAH.

Entry from the Lipid Bilayer: A Possible Pathway for Inhibition of a Peptide G Protein-Coupled Receptor by a Lipophilic Small Molecule

The pathways that G protein-coupled receptor (GPCR) ligands follow as they bind to or dissociate from their receptors are largely unknown. Protease-activated receptor-1 (PAR1) is a GPCR activated by intramolecular binding of a tethered agonist peptide that is exposed by thrombin cleavage. By contrast, the PAR1 antagonist vorapaxar is a lipophilic drug that binds in a pocket almost entirely occluded from the extracellular solvent. The binding and dissociation pathway of vorapaxar is unknown. Starting with the crystal structure of vorapaxar bound to PAR1, we performed temperature-accelerated molecular dynamics simulations of ligand dissociation. In the majority of simulations, vorapaxar exited the receptor laterally into the lipid bilayer through openings in the transmembrane helix (TM) bundle. Prior to full dissociation, vorapaxar paused in metastable intermediates stabilized by interactions with the receptor and lipid headgroups. Derivatives of vorapaxar with alkyl chains predicted to extend between TM6 and TM7 into the lipid bilayer inhibited PAR1 with apparent on rates similar to that of the parent compound in cell signaling assays. These data are consistent with vorapaxar binding to PAR1 via a pathway that passes between TM6 and TM7 from the lipid bilayer, in agreement with the most consistent pathway observed by molecular dynamics. While there is some evidence of entry of the ligand into rhodopsin and lipid-activated GPCRs from the cell membrane, our study provides the first such evidence for a peptide-activated GPCR and suggests that metastable intermediates along drug binding and dissociation pathways can be stabilized by specific interactions between lipids and the ligand.