Oscillatory Activity of P-Type Membrane Adenosine Triphosphatases: a Kinetic Model

A kinetic model for membrane P-type adenosine triphosphatases is considered, the main application being to the erythrocyte Ca²⁺-ATPase. It is shown that a simple modification of the known catalytic mechanism of the ATPase by addition of a self-inhibition step and the steady calcium influx leads to damped oscillations in the system discussed. In this way, the model can explain the kinetic experimental results obtained for the purified enzyme in solution as well as for the enzyme incorporated into liposome membranes. The estimated kinetic parameters are close to the experimental ones. Alternative changes in time, demonstrated by the kinetic model for the conformational enzyme states, E₁ and E₂, confirm the model of two alternatively functioning gates in the ion pumping Ca²⁺-ATPase.__________Translated from Biokhimiya, Vol. 70, No. 4, 2005, pp. 533–538.Original Russian Text Copyright © 2005 by Goldstein, Mayevsky, Zakrjevskaya.     

Conformational dynamics of partially denatured myoglobin studied by time-resolved electrospray mass spectrometry with online hydrogen-deuterium exchange

This study demonstrates the use of electrospray mass spectrometry in conjunction with rapid online mixing ("time-resolved" ESI-MS) for monitoring protein conformational dynamics under equilibrium conditions. The hydrogen/deuterium exchange (HDX) kinetics of mildly denatured myoglobin (Mb) at pD 9.3, in the presence of 27% acetonitrile, were studied with millisecond time resolution. Analytical ultracentrifugation indicates that the average protein compactness under these solvent conditions is similar to that of native holomyoglobin (hMb). The mass spectrum shows protein ions in a wide array of charge and heme binding states, indicating the presence of multiple coexisting conformations. The experimental approach used allows the HDX kinetics of all of these species to be monitored separately. A combination of EX1 and EX2 behavior was observed for hMb ions in charge states 7+ to 9+, which predominantly represent nativelike hMb in solution. The EX1 kinetics are biphasic, indicating the presence of two protein populations that undergo conformational opening events with different rate constants. The EX2 kinetics observed for nativelike hMb are biphasic as well. All other charge and heme binding states represent non-native protein conformations that are involved in rapid interconversion processes, thus leading to monoexponential EX2 kinetics with a common rate constant. Burst phase labeling for these non-native proteins occurs at 125 sites. In contrast, the nativelike protein conformation shows burst phase labeling only for 88 sites. A kinetic model is developed which is based on the assumption of three distinct (un)folding units in Mb. The model implies that the free energy landscape of the protein exhibits a major barrier. The crossing of this barrier is most likely associated with slow, cooperative opening/closing events of the heme binding pocket. Rapid conformational fluctuations on either side of the barrier give rise to the observed EX2 kinetics. Simulated HDX kinetics based on this model are in excellent agreement with the experimental data.     

Characterizing Single-Channel Behavior of GluA3 Receptors

AMPA receptors play a major role in excitatory neurotransmission in the CNS and are involved in numerous neurological disorders. Agonists bind to each of four bilobed LBDs of this tetrameric receptor, and upon binding, the lobes close to envelope the agonist, leading to channel activation. However, AMPA receptors exhibit complex activation kinetics, the mechanism of which has not yet been determined. We report here single-channel studies of a homomeric AMPA receptor (GluA3) activated by the full agonist, glutamate, and a partial agonist, fluorowillardiine. Both agonists activate the channel to the same three open conductance levels but with different open probabilities in each level. The closed probability (P_c) varied within records, particularly at low agonist concentrations. By sorting discrete segments of the record according to P_c using the X-means algorithm, we defined five modes of activity. The kinetic behavior could then be analyzed for both agonists over a range of agonist concentrations with a relatively simple model (three closed states and two open states for each open conductance level). The structural mechanism underlying the modal behavior is not clear; however, it occurs on a timescale consistent with hydrogen bonding across the lobe interface in the LBD.     

Kinetics of anesthetic-induced conformational transitions in a four-alpha-helix bundle protein

Inhaled anesthetics are thought to alter the conformational states of Cys-loop ligand-gated ion channels (LGICs) by binding within discrete cavities that are lined by portions of four alpha-helical transmembrane domains. Because Cys-loop LGICs are complex molecules that are notoriously difficult to express and purify, scaled-down models have been used to better understand the basic molecular mechanisms of anesthetic action. In this study, stopped-flow fluorescence spectroscopy was used to define the kinetics with which inhaled anesthetics interact with (Aalpha(2)-L1M/L38M)(2), a four-alpha-helix bundle protein that was designed to model anesthetic binding sites on Cys-loop LGICs. Stopped-flow fluorescence traces obtained upon mixing (Aalpha(2)-L1M/L38M)(2) with halothane revealed immediate, fast, and slow components of quenching. The immediate component, which occurred within the mixing time of the spectrofluorimeter, was attributed to direct quenching of tryptophan fluorescence upon halothane binding to (Aalpha(2)-L1M/L38M)(2). This was followed by a biexponential fluorescence decay containing fast and slow components, reflecting anesthetic-induced conformational transitions. Fluorescence traces obtained in studies using sevoflurane, isoflurane, and desflurane, which poorly quench tryptophan fluorescence, did not contain the immediate component. However, these anesthetics did produce the fast and slow components, indicating that they also alter the conformation of (Aalpha(2)-L1M/L38M)(2). Cyclopropane, an anesthetic that acts with unusually low potency on Cys-loop LGICs, acted with low apparent potency on (Aalpha(2)-L1M/L38M)(2). These results suggest that four-alpha-helix bundle proteins may be useful models of in vivo sites of action that allow the use of a wide range of techniques to better understand how anesthetic binding leads to changes in protein structure and function.     

Effect of protein kinase A-induced phosphorylation on the gating mechanism of the brain Na+ channel: model fitting to whole-cell current traces.

The activity of the voltage-gated Na+ channel is subjected to modulation through covalent modifications. It has been previously shown that brain Na+ currents are reduced following the activation of the protein kinase A (PKA) pathway, but the effect of the phosphorylation on the gating mechanism of the channel has not been demonstrated so far. In this study, we analyze the whole-cell Na+ current recorded in the absence or presence of forskolin, which stimulates the PKA pathway. A minimal molecular model of the gating mechanism of the Na+ channel is defined to fit the experimental data: it consists of three closed states, one open state, and two inactivated states. We experimentally demonstrate that the kinetics of inactivation from the closed states are not affected by phosphorylation. The results obtained by computer fitting indicate that, among all the kinetic parameters describing the transitions between states, only one parameter is significantly modified in the presence of forskolin, and corresponds to the acceleration of the inactivation from the open state. This conclusion is supported by the analysis of current traces obtained from cells in the presence of a phosphatase inhibitor or loaded with the PKA catalytic unit, and is in agreement with previously reported single channel records.     

Protein homeostasis disorders of key enzymes of amino acids metabolism: mutation-induced protein kinetic destabilization and new therapeutic strategies

Many inborn errors of amino acids metabolism are caused by single point mutations affecting the ability of proteins to fold properly (i.e., protein homeostasis), thus leading to enzyme loss-of-function. Mutations may affect protein homeostasis by altering intrinsic physical properties of the polypeptide (folding thermodynamics, and rates of folding/unfolding/misfolding) as well as the interaction of partially folded states with elements of the protein homeostasis network (such as molecular chaperones and proteolytic machineries). Understanding these mutational effects on protein homeostasis is required to develop new therapeutic strategies aimed to target specific features of the mutant polypeptide. Here, I review recent work in three different diseases of protein homeostasis associated to inborn errors of amino acids metabolism: phenylketonuria, inherited homocystinuria and primary hyperoxaluria type I. These three different genetic disorders involve proteins operating in different cell organelles and displaying different structural complexities. Mutations often decrease protein kinetic stability of the native state (i.e., its half-life for irreversible denaturation), which can be studied using simple kinetic models amenable to biophysical and biochemical characterization. Natural ligands and pharmacological chaperones are shown to stabilize mutant enzymes, thus supporting their therapeutic application to overcome protein kinetic destabilization. The role of molecular chaperones in protein folding and misfolding is also discussed as well as their potential pharmacological modulation as promising new therapeutic approaches. Since current available treatments for these diseases are either burdening or only successful in a fraction of patients, alternative treatments must be considered covering studies from protein structure and biophysics to studies in animal models and patients.     

Single-Channel Kinetic Analysis for Activation and Desensitization of Homomeric 5-HT3A Receptors

The 5-HT₃A receptor is a member of the Cys-loop family of ligand-gated ion channels. To perform kinetic analysis, we mutated the 5-HT3A subunit to obtain a high-conductance form so that single-channel currents can be detected. At all 5-HT concentrations (>0.1 μM), channel activity appears as openings in quick succession that form bursts, which coalesce into clusters. By combining single-channel and macroscopic data, we generated a kinetic model that perfectly describes activation, deactivation, and desensitization. The model shows that full activation arises from receptors with three molecules of agonist bound. It reveals an earlier conformational change of the fully liganded receptor that occurs while the channel is still closed. From this pre-open closed state, the receptor enters into an open-closed cycle involving three open states, which form the cluster whose duration parallels the time constant of desensitization. A similar model lacking the pre-open closed state can describe the data only if the opening rates are fixed to account for the slow activation rate. The application of the model to M4 mutant receptors shows that position 10′ contributes to channel opening and closing rates. Thus, our kinetic model provides a foundation for understanding structural bases of activation and drug action.     

Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding?

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplified protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits protein-like calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Our analysis indicates that a coherent, self-consistent physical account of both the thermodynamic and kinetic properties of the model leads naturally to the concept of a native state ensemble that encompasses considerable conformational fluctuations. Such a multiple-conformation native state is seen to involve conformational states similar to those revealed by native-state hydrogen exchange. Many of these conformational states are predicted to lie below native baselines commonly used in interpreting calorimetric data. Folding and unfolding kinetics are studied under a range of intrachain interaction strengths as in experimental chevron plots. Kinetically determined transition midpoints match well with their thermodynamic counterparts. Kinetic relaxations are found to be essentially single-exponential over an extended range of model interaction strengths. This includes the entire unfolding regime and a significant part of a folding regime with a chevron rollover, as has been observed for real proteins that fold with non-two-state kinetics. The transition state picture of protein folding and unfolding is evaluated by comparing thermodynamic free energy profiles with actual kinetic rates. These analyses suggest that some chevron rollovers may arise from an internal frictional effect that increasingly impedes chain motions with more native conditions, rather than being caused by discrete deadtime folding intermediates or shifts of the transition state peak as previously posited.     

A kinetic mechanism for nicotinic acetylcholine receptors based on multiple allosteric transitions

 Nicotinic acetylcholine receptors are transmembrane oligomeric proteins that mediate interconversions between open and closed channel states under the control of neurotransmitters. Fast in vitro chemical kinetics and in vivo electrophysiological recordings are consistent with the following multi-step scheme. Upon binding of agonists, receptor molecules in the closed but activatable resting state (the Basal state, B) undergo rapid transitions to states of higher affinities with either open channels (the Active state, A) or closed channels (the initial Inactivatable and fully Desensitized states, I and D). In order to represent the functional properties of such receptors, we have developed a kinetic model that links conformational interconversion rates to agonist binding and extends the general principles of the Monod-Wyman-Changeux model of allosteric transitions. The crucial assumption is that the linkage is controlled by the position of the interconversion transition states on a hypothetical linear reaction coordinate. Application of the model to the peripheral nicotinic acetylcholine receptor (nAChR) accounts for the main properties of ligand-gating, including single-channel events, and several new relationships are predicted. Kinetic simulations reveal errors inherent in using the dose-response analysis, but justify its application under defined conditions. The model predicts that (in order to overcome the intrinsic stability of the B state and to produce the appropriate cooperativity) channel activation is driven by an A state with a K_d in the 50 nM range, hence some 140-fold stronger than the apparent affinity of the open state deduced previously. According to the model, recovery from the desensitized states may occur via rapid transit through the A state with minimal channel opening, thus without necessarily undergoing a distinct recovery pathway, as assumed in the standard ‘cyclic’ model. Transitions to the desensitized states by low concentration ‘pre-pulses’ are predicted to occur without significant channel opening, but equilibrium values of IC₅₀ can be obtained only with long pre-pulse times. Predictions are also made concerning allosteric effectors and their possible role in coincidence detection. In terms of future developments, the analysis presented here provides a physical basis for constructing more biologically realistic models of synaptic modulation that may be applied to artificial neural networks. Received: 22 November 1995/Accepted in revised form: 24 July 1996     

On the mechanism of protein fold‐switching by a molecular sensor

Alternate frame folding (AFF) is a mechanism by which conformational change can be engineered into a protein. The protein structure switches from the wild‐type fold (N) to a circularly‐permuted fold (N′), or vice versa, in response to a signaling event such as ligand binding. Despite the fact that the two native states have similar structures, their interconversion involves folding and unfolding of large parts of the molecule. This rearrangement is reported by fluorescent groups whose relative proximities change as a result of the order–disorder transition. The nature of the conformational change is expected to be similar from protein to protein; thus, it may be possible to employ AFF as a general method to create optical biosensors. Toward that goal, we test basic aspects of the AFF mechanism using the AFF variant of calbindin D_9k. A simple three‐state model for fold switching holds that N and N′ interconvert through the unfolded state. This model predicts that the fundamental properties of the switch—calcium binding affinity, signal response (i.e., fluorescence change upon binding), and switching rate—can be controlled by altering the relative stabilities of N and N′. We find that selectively destabilizing N or N′ changes the equilibrium properties of the switch (binding affinity and signal response) in accordance with the model. However, kinetic data indicate that the switching pathway does not require whole‐molecule unfolding. The rate is instead limited by unfolding of a portion of the protein, possibly in concert with folding of a corresponding region. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Communication over the Network of Binary Switches Regulates the Activation of A2A Adenosine Receptor

Dynamics and functions of G-protein coupled receptors (GPCRs) are accurately regulated by the type of ligands that bind to the orthosteric or allosteric binding sites. To glean the structural and dynamical origin of ligand-dependent modulation of GPCR activity, we performed total ~ 5 μsec molecular dynamics simulations of A_2A adenosine receptor (A_2AAR) in its apo, antagonist-bound, and agonist-bound forms in an explicit water and membrane environment, and examined the corresponding dynamics and correlation between the 10 key structural motifs that serve as the allosteric hotspots in intramolecular signaling network. We dubbed these 10 structural motifs “binary switches” as they display molecular interactions that switch between two distinct states. By projecting the receptor dynamics on these binary switches that yield 2¹⁰ microstates, we show that (i) the receptors in apo, antagonist-bound, and agonist-bound states explore vastly different conformational space; (ii) among the three receptor states the apo state explores the broadest range of microstates; (iii) in the presence of the agonist, the active conformation is maintained through coherent couplings among the binary switches; and (iv) to be most specific, our analysis shows that W246, located deep inside the binding cleft, can serve as both an agonist sensor and actuator of ensuing intramolecular signaling for the receptor activation. Finally, our analysis of multiple trajectories generated by inserting an agonist to the apo state underscores that the transition of the receptor from inactive to active form requires the disruption of ionic-lock in the DRY motif.     

A nonequilibrium binary elements-based kinetic model for benzodiazepine regulation of GABAA receptors

Ion channels, like many other proteins, are composed of multiple structural domains. A stimulus that impinges on one domain, such as binding of a ligand to its recognition site, can influence the activity of another domain, such as a transmembrane channel gate, through interdomain interactions. Kinetic schemes that describe the function of interacting domains typically incorporate a minimal number of states and transitions, and do not explicitly model interactions between domains. Here, we develop a kinetic model of the GABA_A receptor, a ligand-gated ion channel modulated by numerous compounds including benzodiazepines, a class of drugs used clinically as sedatives and anxiolytics. Our model explicitly treats both the kinetics of distinct functional domains within the receptor and the interactions between these domains. The model describes not only how benzodiazepines that potentiate GABA_A receptor activity, such as diazepam, affect peak current dose–response relationships in the presence of desensitization, but also their effect on the detailed kinetics of current activation, desensitization, and deactivation in response to various stimulation protocols. Finally, our model explains positive modulation by benzodiazepines of receptor currents elicited by either full or partial agonists, and can resolve conflicting observations arguing for benzodiazepine modulation of agonist binding versus channel gating.     

Experimentally Approaching the Solvent-Accessible Surface Area of a Protein: Insights into the Acid Molten Globule of Bovine α-Lactalbumin

Each conformational state of a protein is inextricably related to a defined extent of solvent exposure that plays a key role in protein folding and protein interactions. However, accurate measurement of the solvent-accessible surface area (ASA) is difficult for any state other than the native (N) state. We address this fundamental physicochemical parameter through a new experimental approach based on the reaction of the photochemical reagent diazirine (DZN) with the polypeptide chain. By virtue of its size, DZN is a reasonable molecular mimic of aqueous solvent. Here, we structurally characterize nonnative states of the paradigmatic protein α-lactalbumin. Covalent tagging resulting from unspecific methylene (:CH₂) reaction allows one to obtain a global estimate of ASA and to map out solvent accessibility along the amino acid sequence. By its mild apolar nature, DZN also reveals a hydrophobic phase in the acid-stabilized state of α-lactalbumin, in which there is clustering of core residues accessible to the solvent. In a fashion reminiscent of the N state, this acid-stabilized state also exhibits local regions where increased :CH₂ labeling indicates its nonhomogenous nature, likely pointing to the existence of packing defects. By contrast, the virtual absence of a defined long-range organization brings about a featureless labeling pattern for the unfolded state. Overall, :CH₂ labeling emerges as a fruitful technique that is able to quantify the ASA of the polypeptide chain, thus probing conformational features such as the outer exposed surface and inner cavities, as well as revealing the existence of noncompact apolar phases in nonnative states.     

Modeling Facilitative Sugar Transporters: Transitions Between Single and Double Ligand Occupancy of Multiconformational Channel Models Explain Anomalous Kinetics

The four-state simple carrier model (SCM) is employed to describe ligand translocation by diverse passive membrane transporters. However, its application to systems like facilitative sugar transporters (GLUTs) is controversial: unidirectional fluxes under zero-trans and equilibrium-exchange experimental conditions fit a SCM, but flux data from infinite-cis and infinite-trans experiments appear not to fit the same SCM. More complex kinetic models have been proposed to explain this ``anomalous' behavior of GLUTs, but none of them accounts for all the experimental findings. We propose an alternative model in which GLUTs are channels subject to conformational transitions, and further assume that the results from zero-trans and equilibrium-exchange experiments as well as trans-effects corresponds to a single-occupancy channel regime, whereas the results from the infinite-cis and infinite-trans experiments correspond to a regime including higher channel occupancies. We test the plausibility of this hypothesis by studying a kinetic model of a two-site channel with two conformational states. In each state, the channel can bind the ligand from only one of the compartments. Under single-occupancy, for conditions corresponding to zero-trans and equilibrium-exchange experiments, the model behaves as a SCM capable of exhibiting trans-stimulations. For a regime including higher degrees of occupancy and infinite-cis and infinite-trans conditions, the same channel model can exhibit a behavior qualitatively similar to a SCM, albeit with kinetic parameters different from those for the single-occupancy regime. Numerical results obtained with our model are consistent with available experimental data on facilitative glucose transport across erythrocyte membranes. Hence, if GLUTs are multiconformational channels, their particular kinetic properties can result from transitions between single and double channel occupancies. Received: 12 April 1995/Revised: 28 August 1995     

Structural basis and kinetics of force-induced conformational changes of an αA domain-containing integrin

Background

Integrin α_Lβ₂ (lymphocyte function-associated antigen, LFA-1) bears force upon binding to its ligand intercellular adhesion molecule 1 (ICAM-1) when a leukocyte adheres to vascular endothelium or an antigen presenting cell (APC) during immune responses. The ligand binding propensity of LFA-1 is related to its conformations, which can be regulated by force. Three conformations of the LFA-1 αA domain, determined by the position of its α₇-helix, have been suggested to correspond to three different affinity states for ligand binding.

Methodology/Principal Findings

The kinetics of the force-driven transitions between these conformations has not been defined and dynamically coupled to the force-dependent dissociation from ligand. Here we show, by steered molecular dynamics (SMD) simulations, that the αA domain was successively transitioned through three distinct conformations upon pulling the C-terminus of its α₇-helix. Based on these sequential transitions, we have constructed a mathematical model to describe the coupling between the αA domain conformational changes of LFA-1 and its dissociation from ICAM-1 under force. Using this model to analyze the published data on the force-induced dissociation of single LFA-1/ICAM-1 bonds, we estimated the force-dependent kinetic rates of interstate transition from the short-lived to intermediate-lived and from intermediate-lived to long-lived states. Interestingly, force increased these transition rates; hence activation of LFA-1 was accelerated by pulling it via an engaged ICAM-1.

Conclusions/Significance

Our study defines the structural basis for mechanical regulation of the kinetics of LFA-1 αA domain conformational changes and relates these simulation results to experimental data of force-induced dissociation of single LFA-1/ICAM-1 bonds by a new mathematical model, thus provided detailed structural and kinetic characterizations for force-stabilization of LFA-1/ICAM-1 interaction.     

Computer programs for single channel kinetic analysis

Programs are represented that permit us to solve the following three problems of single channel kinetic activity: 1) to fulfill a preliminary analysis of single channel activity, i.e., to estimate the minimum number of agonist molecules that bind with a channel structure, the minimum number of open and closed kinetic states of the channel, and possible connections between these states, to check the condition of thermodynamic equilibrium and adequacy of a Markovian assumption for modeling of the channel kinetics; 2) to select the model of channel gating: for a given number of kinetic states, the potential models that permit us to fit a list of dwell-times observed in an experiment are examined; the model is selected by a criterium of the maximum likelihood; 3) to simulate the response to a rapid change of the agonist concentration of single channels that demonstrate processes of adaptation or inactivation.     

Capturing Conformational States in Proteins Using Sparse Paramagnetic NMR Data

Capturing conformational changes in proteins or protein-protein complexes is a challenge for both experimentalists and computational biologists. Solution nuclear magnetic resonance (NMR) is unique in that it permits structural studies of proteins under greatly varying conditions, and thus allows us to monitor induced structural changes. Paramagnetic effects are increasingly used to study protein structures as they give ready access to rich structural information of orientation and long-range distance restraints from the NMR signals of backbone amides, and reliable methods have become available to tag proteins with paramagnetic metal ions site-specifically and at multiple sites. In this study, we show how sparse pseudocontact shift (PCS) data can be used to computationally model conformational states in a protein system, by first identifying core structural elements that are not affected by the environmental change, and then computationally completing the remaining structure based on experimental restraints from PCS. The approach is demonstrated on a 27 kDa two-domain NS2B-NS3 protease system of the dengue virus serotype 2, for which distinct closed and open conformational states have been observed in crystal structures. By changing the input PCS data, the observed conformational states in the dengue virus protease are reproduced without modifying the computational procedure. This data driven Rosetta protocol enables identification of conformational states of a protein system, which are otherwise difficult to obtain either experimentally or computationally.     

A model for the behaviour of phosphorylase b. The generation of different binding sites via intermediate enzymatic states.

The model given in this paper can be applied to enzymatic systems which have more than two conformational states in equilibrium and which clearly exhibit heterogeneity in the binding of one ligand. The model we propose makes possible quantitative interpretation of our experimental results and of those of many other workers as well. In some cases calorimetric, dialysis and kinetic magnitudes, when plotted against ligand concentration, give multiregional or "stepwise" curves. We suggest that such a behaviour arises because total occupation of one class of binding sites completely moves the enzyme towards a different conformational state in which the affinity for the ligand is greatly increased by the formation of a new class of binding sites. Our calorimetric results for the interaction between some nucleotides and phosphorylase b closely conform to our model.     

Fractal dimension as a measure of surface roughness of G protein-coupled receptors: implications for structure and function

Protein surface roughness is a structural property associated with ligand-protein and protein-protein binding interfaces. In this work we apply for the first time the concept of surface roughness, expressed as the fractal dimension, to address structure and function of G protein-coupled receptors (GPCRs) which are an important group of drug targets. We calculate the exposure ratio and the fractal dimension for helix-forming residues of the β(2) adrenergic receptor (β(2)AR), a model system in GPCR studies, in different conformational states: in complex with agonist, antagonist and partial inverse agonists. We show that both exposure ratio and roughness exhibit periodicity which results from the helical structure of GPCRs. The pattern of roughness and exposure ratio of a protein patch depends on its environment: the residues most exposed to membrane are in general most rough whereas parts of receptors mediating interhelical contacts in a monomer or protein complex are much smoother. We also find that intracellular ends (TM3, TM5, TM6 and TM7) which are relevant for G protein binding and thus receptor signaling, are exposed but smooth. Mapping the values of residual fractal dimension onto receptor 3D structures makes it possible to conclude that the binding sites of orthosteric ligands as well as of cholesterol are characterized with significantly higher roughness than the average for the whole protein. In summary, our study suggests that identification of specific patterns of roughness could be a novel approach to spot possible binding sites which could serve as original drug targets for GPCRs modulation.