Observation of "ionic lock" formation in molecular dynamics simulations of wild-type beta 1 and beta 2 adrenergic receptors

Effects of an external static EF on the conformational transition of 5-HT1A receptor: A molecular dynamics simulation study

The serotonin receptor subtype 1A (5-HT1AR), one of the G-protein-coupled receptor (GPCR) family, has been implicated in several neurological conditions. Understanding the activation and inactivation mechanism of 5-HT1AR at the molecular level is critical for discovering novel therapeutics in many diseases. Recently there has been a growing appreciation for the role of external electric fields (EFs) in influencing the structure and activity of biomolecules. In this study, we used molecular dynamics (MD) simulations to examine conformational features of active states of 5-HT1AR and investigate the effect of an external static EF with 0.02 V/nm applied on the active state of 5-HT1AR. Our results showed that the active state of 5-HT1AR maintained the native structure, while the EF led to structural modifications in 5-HT1AR, particularly inducing the inward movement of transmembrane helix 6 (TM6). Furthermore, it disturbed the conformational switches associated with activation in the CWxP, DRY, PIF, and NPxxY motifs, consequently predisposing an inclination towards the inactive-like conformation. We also found that the EF led to an overall increase in the dipole moment of 5-HT1AR, encompassing TM6 and pivotal amino acids. The analyses of conformational properties of TM6 showed that the changed secondary structure and decreased solvent exposure occurred upon the EF condition. The interaction of 5-HT1AR with the membrane lipid bilayer was also altered under the EF. Our findings reveal the molecular mechanism underlying the transition of 5-HT1AR conformation induced by external EFs, which offer potential novel insights into the prospect of employing structure-based EF applications for GPCRs.

Molecular mechanism of fatty acid activation of FFAR1

FFAR1 is a G-protein-coupled receptor (GPCR) that responds to circulating free fatty acids to enhance glucose-stimulated insulin secretion and release of incretin hormones. Due to the glucose-lowering effect of FFAR1 activation, potent agonists for this receptor have been developed for the treatment of diabetes. Previous structural and biochemical studies of FFAR1 showed multiple sites of ligand binding to the inactive state but left the mechanism of fatty acid interaction and receptor activation unknown. We used cryo-electron microscopy to elucidate structures of activated FFAR1 bound to a Gq mimetic, which were induced either by the endogenous FFA ligand docosahexaenoic acid or γ-linolenic acid and the agonist drug TAK-875. Our data identify the orthosteric pocket for fatty acids and show how both endogenous hormones and synthetic agonists induce changes in helical packing along the outside of the receptor that propagate to exposure of the G-protein-coupling site. These structures show how FFAR1 functions without the highly conserved "DRY" and "NPXXY" motifs of class A GPCRs and also illustrate how the orthosteric site of a receptor can be bypassed by membrane-embedded drugs to confer full activation of G protein signaling.

Exploring the kinetic selectivity of drugs targeting the β1 -adrenoceptor

In this study, we report the β1 -adrenoceptor binding kinetics of several clinically relevant β1/2 -adrenoceptor (β1/2 AR) agonists and antagonists. [3 H]-DHA was used to label CHO-β1 AR for binding studies. The kinetics of ligand binding was assessed using a competition association binding method. Ligand physicochemical properties, including logD7.4 and the immobilized artificial membrane partition coefficient (KIAM ), were assessed using column-based methods. Protein Data Bank (PDB) structures and hydrophobic and electrostatic surface maps were constructed in PyMOL. We demonstrate that the hydrophobic properties of a molecule directly affect its kinetic association rate (kon ) and affinity for the β1 AR. In contrast to our findings at the β2 -adrenoceptor, KIAM , reflecting both hydrophobic and electrostatic interactions of the drug with the charged surface of biological membranes, was no better predictor than simple hydrophobicity measurements such as clogP or logD7.4 , at predicting association rate. Bisoprolol proved kinetically selective for the β1 AR subtype, dissociating 50 times slower and partly explaining its higher measured affinity for the β1 AR. We speculate that the association of positively charged ligands at the β1 AR is curtailed somewhat by its predominantly neutral/positive charged extracellular surface. Consequently, hydrophobic interactions in the ligand-binding pocket dominate the kinetics of ligand binding. In comparison at the β2 AR, a combination of hydrophobicity and negative charge attracts basic, positively charged ligands to the receptor's surface promoting the kinetics of ligand binding. Additionally, we reveal the potential role kinetics plays in the on-target and off-target pharmacology of clinically used β-blockers.

Kinetics of Ligand-Protein Dissociation from All-Atom Simulations: Are We There Yet?

Large parallel gains in the development of both computational resources and sampling methods have now made it possible to simulate dissociation events in ligand-protein complexes with all-atom resolution. Such encouraging progress, together with the inherent spatiotemporal resolution associated with molecular simulations, has left their use for investigating dissociation processes brimming with potential, both in rational drug design, where it can be an invaluable tool for determining the mechanistic driving forces behind dissociation rate constants, and in force-field development, where it can provide a catalog of transient molecular structures with which to refine force fields. Although much progress has been made in making force fields more accurate, reducing their error for transient structures along a transition path could yet prove to be a critical development helping to make kinetic predictions much more accurate. In what follows, we will provide a state-of-the-art compilation of the enhanced sampling methods based on molecular dynamics (MD) simulations used to investigate the kinetics and mechanisms of ligand-protein dissociation processes. Due to the time scales of such processes being slower than what is accessible using straightforward MD simulations, several ingenious schemes are being devised at a rapid rate to overcome this obstacle. Here we provide an up-to-date compendium of such methods and their achievements and shortcomings in extracting mechanistic insight into ligand-protein dissociation. We conclude with a critical and provocative appraisal attempting to answer the title of this Perspective.

Ligand modulation of sidechain dynamics in a wild-type human GPCR

GPCRs regulate all aspects of human physiology, and biophysical studies have deepened our understanding of GPCR conformational regulation by different ligands. Yet there is no experimental evidence for how sidechain dynamics control allosteric transitions between GPCR conformations. To address this deficit, we generated samples of a wild-type GPCR (A_2AR) that are deuterated apart from ¹H/¹³C NMR probes at isoleucine δ1 methyl groups, which facilitated ¹H/¹³C methyl TROSY NMR measurements with opposing ligands. Our data indicate that low [Na⁺] is required to allow large agonist-induced structural changes in A_2AR, and that patterns of sidechain dynamics substantially differ between agonist (NECA) and inverse agonist (ZM241385) bound receptors, with the inverse agonist suppressing fast ps-ns timescale motions at the G protein binding site. Our approach to GPCR NMR creates a framework for exploring how different regions of a receptor respond to different ligands or signaling proteins through modulation of fast ps-ns sidechain dynamics.

Designing Safer Analgesics via μ-Opioid Receptor Pathways

Pain is both a major clinical and economic problem, affecting more people than diabetes, heart disease, and cancer combined. While a variety of prescribed or over-the-counter (OTC) medications are available for pain management, opioid medications, especially those acting on the μ-opioid receptor (μOR) and related pathways, have proven to be the most effective, despite some serious side effects including respiration depression, pruritus, dependence, and constipation. It is therefore imperative that both academia and industry develop novel μOR analgesics which retain their opioid analgesic properties but with fewer or no adverse effects. In this review we outline recent progress towards the discovery of safer opioid analgesics.

Understanding the molecular basis of agonist/antagonist mechanism of human mu opioid receptor through gaussian accelerated molecular dynamics method

The most powerful analgesic and addictive properties of opiate alkaloids are mediated by the μ opioid receptor (MOR). The MOR has been extensively investigated as a drug target in the twentieth century, with numerous compounds of varying efficacy being identified. We employed molecular dynamics and Gaussian accelerated molecular dynamics techniques to identify the binding mechanisms of MORs to BU72 (agonist) and β-funaltrexamine (antagonist). Our approach theoretically suggests that the 34 residues (Lys209–Phe221 and Ile301–Cys321) of the MORs were the key regions enabling the two compounds to bind to the active site of the MORs. When the MORs were in the holo form, the key region was in the open conformation. When the MORs were in the apo form, the key region was in the closed conformation. The key region might be responsible for the selectivity of new MOR agonists and antagonists.

An Ensemble-Based Protocol for the Computational Prediction of Helix-Helix Interactions in G Protein-Coupled Receptors using Coarse-Grained Molecular Dynamics

The accurate identification of the specific points of interaction between G protein-coupled receptor (GPCR) oligomers is essential for the design of receptor ligands targeting oligomeric receptor targets. A coarse-grained molecular dynamics computer simulation approach would provide a compelling means of identifying these specific protein–protein interactions and could be applied both for known oligomers of interest and as a high-throughput screen to identify novel oligomeric targets. However, to be effective, this in silico modeling must provide accurate, precise, and reproducible information. This has been achieved recently in numerous biological systems using an ensemble-based all-atom molecular dynamics approach. In this study, we describe an equivalent methodology for ensemble-based coarse-grained simulations. We report the performance of this method when applied to four different GPCRs known to oligomerize using error analysis to determine the ensemble size and individual replica simulation time required. Our measurements of distance between residues shown to be involved in oligomerization of the fifth transmembrane domain from the adenosine A_2A receptor are in very good agreement with the existing biophysical data and provide information about the nature of the contact interface that cannot be determined experimentally. Calculations of distance between rhodopsin, CXCR4, and β₁AR transmembrane domains reported to form contact points in homodimers correlate well with the corresponding measurements obtained from experimental structural data, providing an ability to predict contact interfaces computationally. Interestingly, error analysis enables identification of noninteracting regions. Our results confirm that GPCR interactions can be reliably predicted using this novel methodology.

The DRY motif and the four corners of the cubic ternary complex model

Recent structural data on GPCRs using a variety of spectroscopic approaches suggest that GPCRs adopt a dynamic conformational landscape, with ligands stabilizing subsets of these states to activate one or more downstream signaling effectors. A key outstanding question posed by this emerging dynamic structural model of GPCRs is what states, active, inactive, or intermediate are captured by the numerous crystal structures of GPCRs complexed with a variety of agonists, partial agonists, and antagonists. In the early nineties the discovery of inverse agonists and constitutive activity led to the idea that the active receptor state (R^⁎) is an intrinsic property of the receptor itself rather than of the RG complex, eventually leading to the formulation of the cubic ternary complex model (CTC). Here, by a careful analysis of a series of data obtained with a number of mutants of the highly conserved E/DRY motif, we show evidences for the existence of all the receptor states theorized by the CTC, four 'uncoupled (R, R^⁎ and HR and HR^⁎), and, consequently four 'coupled' (RG, R^⁎G, HRG and HR^⁎G). The E/DRY motif located at the cytosolic end of transmembrane helix III of Class A GPCRs has been widely studied and analyzed because it forms a network of interactions believed to lock receptors in the inactive conformation (R), and, thus, to play a key role in receptor activation. Our conclusions are supported by recent crystal and NMR spectra, as well as by results obtained with two prototypical GPCRs using a new FRET technology that de-couples G protein binding to the receptor from signal transduction. Thus, despite its complexity and limitations, we propose that the CTC is a useful framework to reconcile pharmacological, biochemical and structural data.

Can Specific Protein-Lipid Interactions Stabilize an Active State of the Beta 2 Adrenergic Receptor?

G-protein-coupled receptors are eukaryotic membrane proteins with broad biological and pharmacological relevance. Like all membrane-embedded proteins, their location and orientation are influenced by lipids, which can also impact protein function via specific interactions. Extensive simulations totaling 0.25 ms reveal a process in which phospholipids from the membrane's cytosolic leaflet enter the empty G-protein binding site of an activated β2 adrenergic receptor and form salt-bridge interactions that inhibit ionic lock formation and prolong active-state residency. Simulations of the receptor embedded in an anionic membrane show increased lipid binding, providing a molecular mechanism for the experimental observation that anionic lipids can enhance receptor activity. Conservation of the arginine component of the ionic lock among Rhodopsin-like G-protein-coupled receptors suggests that intracellular lipid ingression between receptor helices H6 and H7 may be a general mechanism for active-state stabilization.

Intramolecular allosteric communication in dopamine D2 receptor revealed by evolutionary amino acid covariation

The structural basis of allosteric signaling in G protein-coupled receptors (GPCRs) is important in guiding design of therapeutics and understanding phenotypic consequences of genetic variation. The Evolutionary Trace (ET) algorithm previously proved effective in redesigning receptors to mimic the ligand specificities of functionally distinct homologs. We now expand ET to consider mutual information, with validation in GPCR structure and dopamine D2 receptor (D2R) function. The new algorithm, called ET-MIp, identifies evolutionarily relevant patterns of amino acid covariations. The improved predictions of structural proximity and D2R mutagenesis demonstrate that ET-MIp predicts functional interactions between residue pairs, particularly potency and efficacy of activation by dopamine. Remarkably, although most of the residue pairs chosen for mutagenesis are neither in the binding pocket nor in contact with each other, many exhibited functional interactions, implying at-a-distance coupling. The functional interaction between the coupled pairs correlated best with the evolutionary coupling potential derived from dopamine receptor sequences rather than with broader sets of GPCR sequences. These data suggest that the allosteric communication responsible for dopamine responses is resolved by ET-MIp and best discerned within a short evolutionary distance. Most double mutants restored dopamine response to wild-type levels, also suggesting that tight regulation of the response to dopamine drove the coevolution and intramolecular communications between coupled residues. Our approach provides a general tool to identify evolutionary covariation patterns in small sets of close sequence homologs and to translate them into functional linkages between residues.

Comparative MD Simulations Indicate a Dual Role for Arg1323.50 in Dopamine-Dependent D2R Activation

Residue Arg^3.50 belongs to the highly conserved DRY-motif of class A GPCRs, which is located at the bottom of TM3. On the one hand, Arg^3.50 has been reported to help stabilize the inactive state of GPCRs, but on the other hand has also been shown to be crucial for stabilizing active receptor conformations and mediating receptor-G protein coupling. The combined results of these studies suggest that the exact function of Arg^3.50 is likely to be receptor-dependent and must be characterized independently for every GPCR. Consequently, we now present comparative molecular-dynamics simulations that use our recently described inactive-state and Gα-bound active-state homology models of the dopamine D2 receptor (D2R), which are either bound to dopamine or ligand-free, performed to identify the function of Arg132^3.50 in D2R. Our results are consistent with a dynamic model of D2R activation in which Arg132^3.50 adopts a dual role, both by stabilizing the inactive-state receptor conformation and enhancing dopamine-dependent D2R-G protein coupling.

Comparing Class A GPCRs to bitter taste receptors: Structural motifs, ligand interactions and agonist-to-antagonist ratios

G protein-coupled receptors (GPCRs) are seven transmembrane (TM) proteins that play a key role in human physiology. The GPCR superfamily comprises about 800 members, classified into several classes, with rhodopsin-like Class A being the largest and most studied thus far. A huge component of the human repertoire consists of the chemosensory GPCRs, including ∼400 odorant receptors, 25 bitter taste receptors (TAS2Rs), which are thought to guard the organism from consuming poisons, and sweet and umami TAS1R heteromers, which indicate the nutritive value of food. The location of the binding site of TAS2Rs is similar to that of Class A GPCRs. However, most of the known bitter ligands are agonists, with only a few antagonists documented thus far. The agonist-to-antagonist ratios of Class A GPCRs vary, but in general are much lower than for TAS2Rs. For a set of well-studied GPCRs, a gradual change in agonists-to-antagonists ratios is observed when comparing low (10 μM)- and high (10 nM)-affinity ligand sets from ChEMBL and the DrugBank set of drugs. This shift reflects pharmaceutical bias toward the therapeutically desirable pharmacology for each of these GPCRs, while the 10 μM sets possibly represent the native tendency of the receptors toward either agonists or antagonists. Analyzing ligand-GPCR interactions in 56 X-ray structures representative of currently available structural data, we find that the N-terminus, TM1 and TM2 are more involved in binding of antagonists than of agonists. On the other hand, ECL2 tends to be more involved in binding of agonists. This is of interest, since TAS2Rs harbor variations on the typical Class A sequence motifs, including the absence of the ECL2-TM3 disulfide bridge. This suggests an alternative mode of regulation of conformational states for TAS2Rs, with potentially less stabilized inactive state. The comparison of TAS2Rs and Class A GPCRs structural features and the pharmacology of the their ligands highlights the intricacies of GPCR architecture and provides a framework for rational design of new ligands.

Multi-scaled normal mode analysis method for dynamics simulation of protein-membrane complexes: A case study of potassium channel gating motion correlations

Membrane proteins play critically important roles in many cellular activities such as ions and small molecule transportation, signal recognition, and transduction. In order to fulfill their functions, these proteins must be placed in different membrane environments and a variety of protein-lipid interactions may affect the behavior of these proteins. One of the key effects of protein-lipid interactions is their ability to change the dynamics status of membrane proteins, thus adjusting their functions. Here, we present a multi-scaled normal mode analysis (mNMA) method to study the dynamics perturbation to the membrane proteins imposed by lipid bi-layer membrane fluctuations. In mNMA, channel proteins are simulated at all-atom level while the membrane is described with a coarse-grained model. mNMA calculations clearly show that channel gating motion can tightly couple with a variety of membrane deformations, including bending and twisting. We then examined bi-channel systems where two channels were separated with different distances. From mNMA calculations, we observed both positive and negative gating correlations between two neighboring channels, and the correlation has a maximum as the channel center-to-center distance is close to 2.5 times of their diameter. This distance is larger than recently found maximum attraction distance between two proteins embedded in membrane which is 1.5 times of the protein size, indicating that membrane fluctuation might impose collective motions among proteins within a larger area. The hybrid resolution feature in mNMA provides atomic dynamics information for key components in the system without costing much computer resource. We expect it to be a conventional simulation tool for ordinary laboratories to study the dynamics of very complicated biological assemblies. The source code is available upon request to the authors.

Computational study on the different ligands induced conformation change of β2 adrenergic receptor-Gs protein complex

β₂ adrenergic receptor (β₂AR) regulated many key physiological processes by activation of a heterotrimeric GTP binding protein (Gs protein). This process could be modulated by different types of ligands. But the details about this modulation process were still not depicted. Here, we performed molecular dynamics (MD) simulations on the structures of β₂AR-Gs protein in complex with different types of ligands. The simulation results demonstrated that the agonist BI-167107 could form hydrogen bonds with Ser203^5.42, Ser207^5.46 and Asn293^6.55 more than the inverse agonist ICI 118,551. The different binding modes of ligands further affected the conformation of β₂AR. The energy landscape profiled the energy contour map of the stable and dissociated conformation of Gαs and Gβγ when different types of ligands bound to β₂AR. It also showed the minimum energy pathway about the conformational change of Gαs and Gβγ along the reaction coordinates. By using interactive essential dynamics analysis, we found that Gαs and Gβγ domain of Gs protein had the tendency to separate when the inverse agonist ICI 118,551 bound to β₂AR. The α5-helix had a relatively quick movement with respect to transmembrane segments of β₂AR when the inverse agonist ICI 118,551 bound to β₂AR. Besides, the analysis of the centroid distance of Gαs and Gβγ showed that the Gαs was separated from Gβγ during the MD simulations. Our results not only could provide details about the different types of ligands that induced conformational change of β₂AR and Gs protein, but also supplied more information for different efficacies of drug design of β₂AR.

Study of structural dynamics of ligand-activated membrane receptors by means of principal component analysis

The structural dynamics of three different ligand-activated G-protein coupled receptors (GPCRs) and the photoreactive receptor rhodopsin from mammals were comparatively studied. As a result, diagrams demonstrating the main structural differences between the studied membrane receptors were obtained. These diagrams represent the projection of the crystal structures of rhodopsin photointermediates and ligand-activated receptors onto the plane defined by the principal components. Thus, we were able to associate the activation process of the receptors with large-scale movements of their individual transmembrane (TM) domains. In addition, the dynamics of extracellular loops of ligand-activated receptors responsible for recognition and initial binding of ligands was studied. Based on these results, two parameters of functionally significant structural dynamics of membrane receptors can be thoroughly analyzed simultaneously - movements of individual TM helices and of extracellular loops.

Biomolecular simulation: a computational microscope for molecular biology

Molecular dynamics simulations capture the behavior of biological macromolecules in full atomic detail, but their computational demands, combined with the challenge of appropriately modeling the relevant physics, have historically restricted their length and accuracy. Dramatic recent improvements in achievable simulation speed and the underlying physical models have enabled atomic-level simulations on timescales as long as milliseconds that capture key biochemical processes such as protein folding, drug binding, membrane transport, and the conformational changes critical to protein function. Such simulation may serve as a computational microscope, revealing biomolecular mechanisms at spatial and temporal scales that are difficult to observe experimentally. We describe the rapidly evolving state of the art for atomic-level biomolecular simulation, illustrate the types of biological discoveries that can now be made through simulation, and discuss challenges motivating continued innovation in this field.

Full and partial agonists of thromboxane prostanoid receptor unveil fine tuning of receptor superactive conformation and G protein activation

The intrahelical salt bridge between E/D^3.49 and R^3.50 within the E/DRY motif on helix 3 (H3) and the interhelical hydrogen bonding between the E/DRY and residues on H6 are thought to be critical in stabilizing the class A G protein-coupled receptors in their inactive state. Removal of these interactions is expected to generate constitutively active receptors. This study examines how neutralization of E^3.49/6.30 in the thromboxane prostanoid (TP) receptor alters ligand binding, basal, and agonist-induced activity and investigates the molecular mechanisms of G protein activation. We demonstrate here that a panel of full and partial agonists showed an increase in affinity and potency for E129V and E240V mutants. Yet, even augmenting the sensitivity to detect constitutive activity (CA) with overexpression of the receptor or the G protein revealed resistance to an increase in basal activity, while retaining fully the ability to cause agonist-induced signaling. However, direct G protein activation measured through bioluminescence resonance energy transfer (BRET) indicates that these mutants more efficiently communicate and/or activate their cognate G proteins. These results suggest the existence of additional constrains governing the shift of TP receptor to its active state, together with an increase propensity of these mutants to agonist-induced signaling, corroborating their definition as superactive mutants. The particular nature of the TP receptor as somehow “resistant” to CA should be examined in the context of its pathophysiological role in the cardiovascular system. Evolutionary forces may have favored regulation mechanisms leading to low basal activity and selected against more highly active phenotypes.

Characterization of the dynamic events of GPCRs by automated computational simulations

The recent advances in membrane protein crystallography have provided extremely valuable structural information of the superfamily of GPCRs (G-protein-coupled receptors). This has been particularly true for a few receptors whose structure was solved several times under different biochemical conditions. It follows that the mechanisms of receptor conformational equilibrium and related dynamic events can be explored by computational simulations. In the present article, we summarize our recent understanding of several dynamic features of GPCRs, accomplished through the use of MD (molecular dynamics) simulations. Our pipeline for the MD simulations of GPCRs, implemented in the web service http://gpcr.usc.es, is updated in the present paper and illustrated by recent applications. Special emphasis is put on the A2A adenosine receptor, one of the selected cases where crystal structures in several conformations and conditions exist, and on the dimerization process of the CXCR4 (CXC chemokine receptor 4).

Wagging the Tail: Essential Role of Substrate Flexibility in FAAH Catalysis

The serine hydrolase, fatty acid amide hydrolase (FAAH), is responsible for the intracellular degradation of anandamide and other bioactive fatty acid ethanolamides involved in the regulation of pain, inflammation, and other pathophysiological processes. The catalytic site of FAAH is composed of multiple cavities with mixed hydrophobic and hydrophilic properties, the role of which remains incompletely understood. Anandamide is thought to enter the active site through a "membrane-access" (MA) channel and position its flexible fatty acyl chain in a highly hydrophobic "acyl chain-binding" (AB) cavity to allow for hydrolysis to occur. Using microsecond molecular dynamics (MD) simulations of FAAH embedded in a realistic membrane/water environment, we show now that anandamide may not lock itself into the AB cavity but may rather assume catalytically significant conformations required for hydrolysis by moving its flexible arachidonoyl tail between the MA and AB cavities. This process is regulated by a phenylalanine residue (Phe432) located at the boundary between the two cavities, which may act as a "dynamic paddle." The results identify structural flexibility as a key determinant by which FAAH recognizes its primary lipid substrate.

Structure-function of the G protein-coupled receptor superfamily

During the past few years, crystallography of G protein-coupled receptors (GPCRs) has experienced exponential growth, resulting in the determination of the structures of 16 distinct receptors-9 of them in 2012 alone. Including closely related subtype homology models, this coverage amounts to approximately 12% of the human GPCR superfamily. The adrenergic, rhodopsin, and adenosine receptor systems are also described by agonist-bound active-state structures, including a structure of the receptor-G protein complex for the β(2)-adrenergic receptor. Biochemical and biophysical techniques, such as nuclear magnetic resonance and hydrogen-deuterium exchange coupled with mass spectrometry, are providing complementary insights into ligand-dependent dynamic equilibrium between different functional states. Additional details revealed by high-resolution structures illustrate the receptors as allosteric machines that are controlled not only by ligands but also by ions, lipids, cholesterol, and water. This wealth of data is helping redefine our knowledge of how GPCRs recognize such a diverse array of ligands and how they transmit signals 30 angstroms across the cell membrane; it also is shedding light on a structural basis of GPCR allosteric modulation and biased signaling.

Action of molecular switches in GPCRs--theoretical and experimental studies

G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called "molecular switches" buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homoand heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

The significance of G protein-coupled receptor crystallography for drug discovery

Crucial as molecular sensors for many vital physiological processes, seven-transmembrane domain G protein-coupled receptors (GPCRs) comprise the largest family of proteins targeted by drug discovery. Together with structures of the prototypical GPCR rhodopsin, solved structures of other liganded GPCRs promise to provide insights into the structural basis of the superfamily's biochemical functions and assist in the development of new therapeutic modalities and drugs. One of the greatest technical and theoretical challenges to elucidating and exploiting structure-function relationships in these systems is the emerging concept of GPCR conformational flexibility and its cause-effect relationship for receptor-receptor and receptor-effector interactions. Such conformational changes can be subtle and triggered by relatively small binding energy effects, leading to full or partial efficacy in the activation or inactivation of the receptor system at large. Pharmacological dogma generally dictates that these changes manifest themselves through kinetic modulation of the receptor's G protein partners. Atomic resolution information derived from increasingly available receptor structures provides an entrée to the understanding of these events and practically applying it to drug design. Supported by structure-activity relationship information arising from empirical screening, a unified structural model of GPCR activation/inactivation promises to both accelerate drug discovery in this field and improve our fundamental understanding of structure-based drug design in general. This review discusses fundamental problems that persist in drug design and GPCR structural determination.

Comparison of dynamics of extracellular accesses to the β(1) and β(2) adrenoceptors binding sites uncovers the potential of kinetic basis of antagonist selectivity

From the molecular mechanism of antagonist unbinding in the β(1) and β(2) adrenoceptors investigated by steered molecular dynamics, we attempt to provide further possibilities of ligand subtype and subspecies selectivity. We have simulated unbinding of β(1)-selective Esmolol and β(2)-selective ICI-118551 from both receptors to the extracellular environment and found distinct molecular features of unbinding. By calculating work profiles, we show different preference in antagonist unbinding pathways between the receptors, in particular, perpendicular to the membrane pathway is favourable in the β(1) adrenoceptor, whereas the lateral pathway involving helices 5, 6 and 7 is preferable in the β(2) adrenoceptor. The estimated free energy change of unbinding based on the preferable pathway correlates with the experimental ligand selectivity. We then show that the non-conserved K347 (6.58) appears to facilitate in guiding Esmolol to the extracellular surface via hydrogen bonds in the β(1) adrenoceptor. In contrast, hydrophobic and aromatic interactions dominate in driving ICI-118551 through the easiest pathway in the β(2) adrenoceptor. We show how our study can stimulate design of selective antagonists and discuss other possible molecular reasons of ligand selectivity, involving sequential binding of agonists and glycosylation of the receptor extracellular surface.

Structural modelling and dynamics of proteins for insights into drug interactions

Proteins are the workhorses of biomolecules and their function is affected by their structure and their structural rearrangements during ligand entry, ligand binding and protein-protein interactions. Hence, the knowledge of protein structure and, importantly, the dynamic behaviour of the structure are critical for understanding how the protein performs its function. The predictions of the structure and the dynamic behaviour can be performed by combinations of structure modelling and molecular dynamics simulations. The simulations also need to be sensitive to the constraints of the environment in which the protein resides. Standard computational methods now exist in this field to support the experimental effort of solving protein structures. This review presents a comprehensive overview of the basis of the calculations and the well-established computational methods used to generate and understand protein structure and function and the study of their dynamic behaviour with the reference to lung-related targets.

Membrane Protein Simulations Using AMBER Force Field and Berger Lipid Parameters

AMBER force fields are among the most commonly used in molecular dynamics (MD) simulations of proteins. Unfortunately, they lack a specific set of lipid parameters, thus limiting its use in membrane protein simulations. In order to overcome this limitation we assessed whether the widely used united-atom lipid parameters described by Berger and co-workers could be used in conjunction with AMBER force fields in simulations of membrane proteins. Thus, free energies of solvation in water and in cyclohexane, and free energies of water to cyclohexane transfer, were computed by thermodynamic integration procedures for neutral amino acid side-chains employing AMBER99, AMBER03, and OPLS-AA amino acid force fields. In addition, MD simulations of three membrane proteins in a POPC lipid bilayer, the β2 adrenergic G protein-coupled receptor, Aquaporin-1, and the outer membrane protein Omp32, were performed with the aim of comparing the AMBER99SB/Berger combination of force fields with the OPLS-AA/Berger combination. We have shown that AMBER99SB and Berger force fields are compatible, they provide reliable free energy estimations relative to experimental values, and their combination properly describes both membrane and protein structural properties. We then suggest that the AMBER99SB/Berger combination is a reliable choice for the simulation of membrane proteins, which links the easiness of ligand parametrization and the ability to reproduce secondary structure of AMBER99SB force field with the largely validated Berger lipid parameters.

Heterologous expression of functional G-protein-coupled receptors in Caenorhabditis elegans

New strategies for expression, purification, functional characterization, and structural determination of membrane-spanning G-protein-coupled receptors (GPCRs) are constantly being developed because of their importance to human health. Here, we report a Caenorhabditis elegans heterologous expression system able to produce milligram amounts of functional native and engineered GPCRs. Both bovine opsin [(b)opsin] and human adenosine A(2A) subtype receptor [(h)A(2A)R] expressed in neurons or muscles of C. elegans were localized to cell membranes. Worms expressing these GPCRs manifested changes in motor behavior in response to light and ligands, respectively. With a newly devised protocol, 0.6-1 mg of purified homogenous 9-cis-retinal-bound bovine isorhodopsin [(b)isoRho] and ligand-bound (h)A(2A)R were obtained from C. elegans from one 10-L fermentation at low cost. Purified recombinant (b)isoRho exhibited its signature absorbance spectrum and activated its cognate G-protein transducin in vitro at a rate similar to native rhodopsin (Rho) obtained from bovine retina. Generally high expression levels of 11 native and mutant GPCRs demonstrated the potential of this C. elegans system to produce milligram quantities of high-quality GPCRs and possibly other membrane proteins suitable for detailed characterization.

Predicting novel binding modes of agonists to β adrenergic receptors using all-atom molecular dynamics simulations

Understanding the binding mode of agonists to adrenergic receptors is crucial to enabling improved rational design of new therapeutic agents. However, so far the high conformational flexibility of G protein-coupled receptors has been an obstacle to obtaining structural information on agonist binding at atomic resolution. In this study, we report microsecond classical molecular dynamics simulations of β(1) and β(2) adrenergic receptors bound to the full agonist isoprenaline and in their unliganded form. These simulations show a novel agonist binding mode that differs from the one found for antagonists in the crystal structures and from the docking poses reported by in silico docking studies performed on rigid receptors. Internal water molecules contribute to the stabilization of novel interactions between ligand and receptor, both at the interface of helices V and VI with the catechol group of isoprenaline as well as at the interface of helices III and VII with the ethanolamine moiety of the ligand. Despite the fact that the characteristic N-C-C-OH motif is identical in the co-crystallized ligands and in the full agonist isoprenaline, the interaction network between this group and the anchor site formed by Asp(3.32) and Asn(7.39) is substantially different between agonists and inverse agonists/antagonists due to two water molecules that enter the cavity and contribute to the stabilization of a novel network of interactions. These new binding poses, together with observed conformational changes in the extracellular loops, suggest possible determinants of receptor specificity.

Protein dynamics and enzyme catalysis: insights from simulations

The role of protein dynamics in enzyme catalysis is one of the most active and controversial areas in enzymology today. Some researchers claim that protein dynamics are at the heart of enzyme catalytic efficiency, while others state that dynamics make no significant contribution to catalysis. What is the biochemist - or student - to make of the ferocious arguments in this area? Protein dynamics are complex and fascinating, as molecular dynamics simulations and experiments have shown. The essential question is: do these complex motions have functional significance? In particular, how do they affect or relate to chemical reactions within enzymes, and how are chemical and conformational changes coupled together? Biomolecular simulations can analyse enzyme reactions and dynamics in atomic detail, beyond that achievable in experiments: accurate atomistic modelling has an essential part to play in clarifying these issues. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Putative active states of a prototypic g-protein-coupled receptor from biased molecular dynamics

A major current focus of structural work on G-protein-coupled receptors (GPCRs) pertains to the investigation of their active states. However, for virtually all GPCRs, active agonist-bound intermediate states have been difficult to characterize experimentally owing to their higher conformational flexibility, and thus intrinsic instability, as compared to inactive inverse agonist-bound states. In this work, we explored possible activation pathways of the prototypic GPCR bovine rhodopsin by means of biased molecular dynamics simulations. Specifically, we used an explicit atomistic representation of the receptor and its environment, and sampled the conformational transition from the crystal structure of a photoactivated deprotonated state of rhodopsin to the low pH crystal structure of opsin in the presence of 11-trans-retinal, using adiabatic biased molecular dynamics simulations. We then reconstructed the system free-energy landscape along the predetermined transition trajectories using a path collective variable approach based on metadynamics. Our results suggest that the two experimental endpoints of rhodopsin/opsin are connected by at least two different pathways, and that the conformational transition is populated by at least four metastable states of the receptor, characterized by a different amplitude of the outward movement of transmembrane helix 6.

Role of aggregation in rhodopsin signal transduction

Many G protein-coupled receptors (GPCRs) are known to form dimers or even oligomers, and these aggregated states have been proposed as functional units responsible for signal transduction and G protein activation. However, the nature of their involvement has remained elusive. Here, we have investigated the role of aggregation in the signal transduction for dimeric forms of the prototypical GPCR rhodopsin using molecular dynamics simulations. The early steps after photoexcitation are characterized by a tandem mechanism in which one monomer is responsible for light detection while the other serves as the G protein activation site. Dimerization ensures efficient cross-talk between the two units within a few tens of nanoseconds following photoexcitation. This interface-mediated pathway suggests oligomerization-aided signal transduction as a crucial biological strategy to enhance activation efficiency across the entire family of GPCRs.

A lipid pathway for ligand binding is necessary for a cannabinoid G protein-coupled receptor

Recent isothiocyanate covalent labeling studies have suggested that a classical cannabinoid, (-)-7'-isothiocyanato-11-hydroxy-1',1'dimethylheptyl-hexahydrocannabinol (AM841), enters the cannabinoid CB2 receptor via the lipid bilayer (Pei, Y., Mercier, R. W., Anday, J. K., Thakur, G. A., Zvonok, A. M., Hurst, D., Reggio, P. H., Janero, D. R., and Makriyannis, A. (2008) Chem. Biol. 15, 1207-1219). However, the sequence of steps involved in such a lipid pathway entry has not yet been elucidated. Here, we test the hypothesis that the endogenous cannabinoid sn-2-arachidonoylglycerol (2-AG) attains access to the CB2 receptor via the lipid bilayer. To this end, we have employed microsecond time scale all-atom molecular dynamics (MD) simulations of the interaction of 2-AG with CB2 via a palmitoyl-oleoyl-phosphatidylcholine lipid bilayer. Results suggest the following: 1) 2-AG first partitions out of bulk lipid at the transmembrane alpha-helix (TMH) 6/7 interface; 2) 2-AG then enters the CB2 receptor binding pocket by passing between TMH6 and TMH7; 3) the entrance of the 2-AG headgroup into the CB2 binding pocket is sufficient to trigger breaking of the intracellular TMH3/6 ionic lock and the movement of the TMH6 intracellular end away from TMH3; and 4) subsequent to protonation at D3.49/D6.30, further 2-AG entry into the ligand binding pocket results in both a W6.48 toggle switch change and a large influx of water. To our knowledge, this is the first demonstration via unbiased molecular dynamics that a ligand can access the binding pocket of a class A G protein-coupled receptor via the lipid bilayer and the first demonstration via molecular dynamics of G protein-coupled receptor activation triggered by a ligand binding event.