Principles of Cell Signaling

Introduction

No cell is an island. From single-celled bacteria coordinating swarming behavior to the trillions of cells in a human body that must grow, differentiate, and die in precise harmony, cell signaling — the ability of cells to send, receive, and interpret chemical messages — is fundamental to life. A liver cell must know when blood glucose is high; an immune cell must know when a pathogen has invaded; a developing neuron must know where to send its axon. In every case, the answer depends on extracellular signals and the molecular machinery that detects and interprets them.

The principles governing cell communication are remarkably conserved. Whether the signal is a steroid hormone diffusing through the plasma membrane to activate a nuclear receptor, or a growth factor binding a cell-surface receptor to trigger an intracellular kinase cascade, the underlying logic follows a common architecture: a signal, a receptor, a series of intracellular mediators, and a cellular response. This lesson examines that architecture in depth, beginning with the general principles of cell communication (Section 15.1) and then exploring the best-understood class of cell-surface receptors — the G-protein-coupled receptors (GPCRs) — and their second messenger pathways (Section 15.2). We also introduce the bioinformatics databases and computational tools used to analyze signaling pathways, phosphoproteomics data, and GPCR pharmacology.

15.1 — General Principles of Cell Communication

Extracellular Signal Molecules Bind to Specific Receptors

Cells communicate by releasing extracellular signal molecules (also called ligands) that are detected by receptor proteins on or inside the target cell. The signal molecule and its receptor interact with high specificity and affinity, much like an enzyme binding its substrate. This molecular recognition ensures that each signal activates only its intended pathway.

Signal molecules fall into two broad categories based on their chemistry and the location of their receptors:

Hydrophobic signals can cross the plasma membrane directly and bind intracellular receptors — often members of the nuclear receptor superfamily — that function as ligand-activated transcription factors. Hydrophilic signals cannot cross the membrane and instead bind cell-surface receptors that relay the message inward through intracellular signaling cascades.

Signal Molecules Can Act Over Short or Long Distances

The distance over which a signal acts defines four major modes of cell communication:

Endocrine signaling relies on hormones that travel through the circulatory system and can reach every cell in the body. Only cells expressing the appropriate receptor respond. Paracrine signaling involves molecules that act locally; their diffusion range is limited by uptake, degradation, or binding to extracellular matrix components. Synaptic signaling is a specialized form of paracrine signaling in which neurotransmitters are released at discrete synaptic junctions and act on immediately adjacent postsynaptic cells.

Autocrine Signaling Can Coordinate Decisions by Groups of Identical Cells

In autocrine signaling, a cell secretes a signal molecule that binds to receptors on its own surface (or on neighboring cells of the same type). This creates a positive-feedback loop that can synchronize the behavior of a population of similar cells. Autocrine signaling is important in the immune system, where T cells secrete interleukins that stimulate their own proliferation, and in development, where it can drive a group of equivalent precursor cells to differentiate simultaneously. It is also co-opted in cancer: many tumor cells produce growth factors that stimulate their own proliferation, creating a self-sustaining growth signal.

Gap Junctions Allow Signaling Information to Be Shared by Neighboring Cells

Gap junctions are intercellular channels formed by the docking of connexon hemichannels (hexamers of connexin proteins) on adjacent cells. Each channel provides a direct cytoplasmic connection through which ions, metabolites, and small signaling molecules (up to ~1,000 Da) can diffuse between cells. This allows electrically and metabolically coupled cells to share second messengers such as Ca²⁺, IP₃, and cAMP without secreting them into the extracellular space. Gap-junction coupling coordinates the contraction of cardiac muscle cells and the firing of smooth muscle cells in the gut, and it helps synchronize metabolic activity in groups of cells in the liver and other tissues.

Each Cell Responds to a Limited Set of Extracellular Signals

Although a complex organism produces hundreds of distinct signal molecules, any individual cell expresses receptors for only a subset of them. The specific combination of receptors on a cell determines its signaling repertoire — the set of signals to which it can respond. Furthermore, two different cell types may express the same receptor yet respond differently because they contain different sets of intracellular signaling proteins and effectors. For example, acetylcholine causes skeletal muscle to contract but causes heart muscle to relax, because the downstream signaling machinery differs in the two cell types.

Three Major Classes of Cell-Surface Receptor Proteins

Cell-surface receptors fall into three major families, classified by their mechanism of signal transduction:

GPCRs are by far the largest family, with ~800 members in the human genome — roughly 4% of all protein-coding genes. They mediate responses to hormones, neurotransmitters, odors, tastes, and light, and are the target of approximately 35% of all approved drugs.

Cell-Surface Receptors Relay Signals via Intracellular Signaling Molecules

When a signal molecule binds a cell-surface receptor, the receptor does not carry out the cellular response itself. Instead, it activates one or more intracellular signaling molecules that relay, amplify, and distribute the signal within the cell. These intracellular mediators come in two forms:

• Intracellular signaling proteins — proteins whose activity is regulated by phosphorylation, GTP binding, or protein-protein interactions (e.g., kinases, GTPases, adaptor proteins)
• Small intracellular mediators (second messengers) — small molecules or ions generated in large numbers in response to receptor activation (e.g., cyclic AMP, Ca²⁺, IP₃, DAG)

The signaling chain from receptor to cellular response typically involves multiple steps, organized into modular signaling pathways. Each step provides an opportunity for amplification, integration, and regulation.

Intracellular Signals Must Be Specific and Precise in a Noisy Cytoplasm

The cytoplasm is a crowded environment containing thousands of different protein species at concentrations that could, in principle, generate spurious interactions. How does the cell maintain signaling specificity — ensuring that the right signal activates the right target — in this noisy environment?

Several mechanisms contribute. First, signaling proteins are organized into dedicated signaling complexes that physically bring pathway components together and exclude non-partners. Second, scaffold proteins tether sequential components of a signaling cascade to ensure that the output of one step is delivered directly to the input of the next. Third, modular interaction domains on signaling proteins mediate specific protein-protein interactions (discussed below). Finally, the cell uses spatial localization — targeting signaling components to specific subcellular compartments — to insulate parallel pathways from one another.

Intracellular Signaling Complexes Form at Activated Receptors

When a receptor is activated, it often becomes a nucleation site for a signaling complex — a transient assembly of intracellular signaling proteins recruited to the cytoplasmic face of the receptor. For example, an activated receptor tyrosine kinase autophosphorylates on multiple tyrosine residues, each of which becomes a specific docking site for a different signaling protein. The result is a multi-protein complex that simultaneously activates several downstream pathways. Activated GPCRs similarly recruit β-arrestins and other scaffold proteins that organize signaling complexes.

Modular Interaction Domains Mediate Interactions Between Intracellular Signaling Proteins

Intracellular signaling proteins are typically built from modular interaction domains — compact, independently folding protein modules that mediate specific binding interactions:

These domains allow signaling proteins to be assembled in modular fashion: a single protein may contain an SH2 domain (to bind an activated receptor), an SH3 domain (to recruit a downstream effector), and a catalytic domain (to propagate the signal). This modularity enables evolutionary flexibility — new signaling pathways can arise by shuffling existing domains into new combinations.

The Relationship Between Signal and Response Varies in Different Signaling Pathways

Signaling pathways differ in the quantitative relationship between signal strength and cellular response:

• Graded responses — the magnitude of the response increases smoothly with increasing signal concentration (e.g., gene expression proportional to hormone level)
• Switch-like (ultrasensitive) responses — the response changes abruptly over a narrow range of signal concentration, producing an effectively binary (on/off) output
• Bistable responses — the cell has two stable states and switches irreversibly from one to the other once the signal crosses a threshold

The type of input-output relationship is determined by the molecular architecture of the pathway, including the presence of cooperative binding, multisite phosphorylation, positive feedback, and double-negative feedback loops.

The Speed of a Response Depends on the Turnover of Signaling Molecules

How fast a cell can respond to a signal — and how fast it can stop responding when the signal is removed — depends on the turnover rate of the signaling molecules involved. A signaling molecule that is rapidly synthesized and rapidly degraded can change concentration quickly, enabling fast responses. A molecule with slow turnover changes concentration slowly, producing a slow, sustained response.

This principle applies to both second messengers (cAMP is rapidly destroyed by phosphodiesterases) and signaling proteins (phosphorylated proteins are rapidly dephosphorylated by protein phosphatases). The half-life of a signaling molecule determines the time resolution of the signal: fast turnover gives millisecond responses (ion channels, Ca²⁺ transients in neurons), while slow turnover gives hour-to-day responses (changes in gene expression in response to steroid hormones).

Cells Can Respond Abruptly to a Gradually Increasing Signal

Many signaling pathways exhibit ultrasensitivity — a steeper-than-linear relationship between signal strength and response. Mechanistically, ultrasensitivity can arise from cooperative binding (as in hemoglobin’s oxygen-binding curve), multisite phosphorylation (a target protein must be phosphorylated on multiple sites before it becomes active), or zero-order ultrasensitivity (when opposing kinases and phosphatases are both operating near saturation). Ultrasensitive responses filter out low-level noise and sharpen the transition between inactive and active states, enabling cells to make clear binary decisions from noisy analog inputs.

Positive Feedback Can Generate an All-or-None Response

Positive feedback occurs when the output of a pathway enhances its own production. If combined with ultrasensitivity, positive feedback can create a bistable switch — a system with two stable states (fully off and fully on) and no stable intermediate. Once the signal exceeds a threshold, the system flips to the on state and remains there even if the signal is subsequently reduced below the threshold (hysteresis). This mechanism underlies many irreversible biological decisions, including oocyte maturation (driven by MAPK positive feedback), cell-cycle entry at the restriction point, and some forms of cellular differentiation.

Negative Feedback Is a Common Motif in Signaling Systems

Negative feedback occurs when the output of a pathway inhibits an upstream step in the same pathway. It is one of the most pervasive motifs in signaling biology. Negative feedback serves multiple functions:

• Limits the duration and intensity of the response, preventing overstimulation
• Generates oscillations when feedback operates with a time delay (e.g., Ca²⁺ oscillations, NF-κB oscillations)
• Improves robustness by making the output less sensitive to fluctuations in upstream component concentrations
• Speeds up the response by operating on a faster timescale than simple deactivation

The cAMP pathway provides a classic example: cAMP activates PKA, and PKA phosphorylates and activates phosphodiesterase, which degrades cAMP — creating a self-limiting pulse of cAMP.

Cells Can Adjust Their Sensitivity to a Signal

Prolonged exposure to a signal typically causes a cell to become less responsive — a process called adaptation or desensitization. This allows cells to respond to changes in signal concentration rather than to absolute levels, enormously extending the dynamic range of the signaling system.

Desensitization occurs through several mechanisms:

• Receptor inactivation — phosphorylation of the receptor by specific kinases (e.g., GRKs for GPCRs) reduces its ability to activate downstream signaling
• Receptor sequestration — receptors are removed from the cell surface by endocytosis and sequestered in endosomes
• Receptor downregulation — endocytosed receptors are targeted to lysosomes for degradation, reducing total receptor number
• Downstream pathway inactivation — negative feedback reduces the activity of intracellular signaling components

The reverse process, sensitization, can occur when cells upregulate receptor expression or signaling components in response to chronically low signal levels.

Signaling Pathway Databases and Analysis

Pathway Databases

The complexity of cell signaling has driven the creation of curated databases that catalog known pathways, interactions, and regulatory relationships:

These databases are essential for interpreting high-throughput experimental data. For example, after a phosphoproteomics experiment identifies thousands of phosphorylation changes, pathway enrichment analysis (using tools like DAVID, g:Profiler, or Reactome’s analysis service) can determine which signaling pathways are significantly activated or suppressed.

Kinase-Substrate Databases and Phosphoproteomics

Protein phosphorylation is the most common post-translational modification in signaling, and mass spectrometry-based phosphoproteomics can now quantify tens of thousands of phosphorylation sites in a single experiment. Interpreting this data requires knowing which kinase phosphorylates which substrate:

• PhosphoSitePlus — the most comprehensive database of experimentally observed phosphorylation sites (and other PTMs), with upstream kinase annotations and information on the functional effects of each modification
• KinaseNET — kinase-substrate relationships with quantitative data

Computational tools for phosphoproteomics data analysis include kinase enrichment analysis (determining which kinases are likely responsible for the observed pattern of phosphorylation changes) and signaling network inference (reconstructing the active signaling network from phosphoproteomics data).

Signaling Network Modeling

Mathematical models of signaling networks help researchers understand emergent properties like oscillations, bistability, and ultrasensitivity that cannot be predicted from knowledge of individual components alone. Common modeling frameworks include:

• Boolean models — each node is ON or OFF; suited for large-scale qualitative analysis of pathway logic
• Ordinary differential equation (ODE) models — continuous concentrations change over time; capture quantitative dynamics (e.g., oscillation periods, dose-response curves)
• Agent-based models — individual molecules are simulated stochastically; capture spatial effects and cell-to-cell variability

Ligand-Receptor Interaction and Drug-Target Databases

Understanding which ligands bind which receptors is critical for drug discovery. Key resources include:

• CellTalkDB — a curated database of ligand-receptor pairs used in cell-cell communication analysis (especially for single-cell RNA-seq studies)
• DLRP (Database of Ligand-Receptor Partners) — literature-curated ligand-receptor interactions
• DrugBank and ChEMBL — databases linking approved and experimental drugs to their molecular targets, enabling drug-target interaction prediction

Let us explore a bioinformatics perspective on signaling proteins. The β₂-adrenergic receptor is a prototypical GPCR whose gene sequence we can analyze:

let epinephrine = "CNC[C@@H](O)c1ccc(O)c(O)c1"
let acetylcholine = "CC(=O)OCC[N+](C)(C)C"
let cortisol = "OC(=O)C1(O)CCC2C1CCC1C2CC(O)C2=CC(=O)CCC12C"
print("Epinephrine (hydrophilic - binds surface receptor):")
print(Chem.properties(epinephrine))
print("Acetylcholine (hydrophilic - binds surface receptor):")
print(Chem.properties(acetylcholine))
print("Cortisol (hydrophobic - crosses membrane to nuclear receptor):")
print(Chem.properties(cortisol))

The chemical properties of signaling molecules determine whether they bind surface or intracellular receptors. Hydrophobic molecules like cortisol can cross the membrane; hydrophilic ones like epinephrine cannot.

15.2 — Signaling Through G-Protein-Coupled Cell-Surface Receptors (GPCRs)

GPCRs are the largest and most versatile family of cell-surface receptors. All GPCRs share a common architecture: seven transmembrane α helices connected by alternating intracellular and extracellular loops, with an extracellular N-terminus and an intracellular C-terminus. Despite this conserved fold, GPCRs detect an astonishing range of signals — from photons and small-molecule odorants to large protein hormones — and they activate diverse intracellular responses through a shared transduction mechanism involving trimeric G proteins.

Trimeric G Proteins Relay Signals from GPCRs

The immediate downstream partners of GPCRs are trimeric G proteins (also called heterotrimeric G proteins), consisting of three subunits: Gα, Gβ, and Gγ. In the resting state, Gα is bound to GDP and associated with the Gβγ dimer. Signal transduction proceeds through a GTPase cycle:

1. Ligand binds GPCR → receptor undergoes conformational change
2. Activated receptor acts as a guanine-nucleotide exchange factor (GEF), promoting the exchange of GDP for GTP on Gα
3. GTP-bound Gα dissociates from Gβγ — both Gα-GTP and free Gβγ can now activate downstream effectors
4. Gα hydrolyzes GTP to GDP (intrinsic GTPase activity, accelerated by RGS proteins), triggering reassociation with Gβγ and returning the system to the resting state

The GTPase cycle acts as a built-in molecular timer: the duration of signaling is determined by how long Gα remains in the GTP-bound state. Regulators of G-protein signaling (RGS) proteins accelerate GTP hydrolysis and are thus critical negative regulators of GPCR signaling.

Different Gα subtypes activate different effector pathways:

Some G Proteins Regulate the Production of Cyclic AMP

The cAMP pathway was the first second-messenger system to be discovered (by Earl Sutherland in the 1950s, studying the action of adrenaline on glycogen breakdown). The pathway is straightforward:

1. Gαs activates adenylyl cyclase, a transmembrane enzyme that converts ATP to cyclic AMP (cAMP)
2. cAMP rapidly diffuses through the cytoplasm, amplifying the signal
3. cAMP is degraded by cyclic nucleotide phosphodiesterases (PDEs), terminating the signal

Gαi has the opposite effect: it inhibits adenylyl cyclase, reducing cAMP levels. The net cAMP concentration reflects the balance between stimulatory and inhibitory inputs, allowing the cell to integrate signals from multiple GPCRs.

Caffeine exerts its stimulatory effects in part by inhibiting phosphodiesterases, thereby prolonging cAMP signaling. The drug forskolin directly activates adenylyl cyclase and is widely used as a research tool to elevate cAMP levels experimentally.

Cyclic AMP Activates Protein Kinase A (PKA)

The principal target of cAMP in most animal cells is protein kinase A (PKA), also called cAMP-dependent protein kinase. PKA is a tetrameric enzyme consisting of two regulatory (R) subunits and two catalytic (C) subunits. In the inactive state, the R subunits bind and inhibit the C subunits. When cAMP levels rise, each R subunit binds two molecules of cAMP, causing a conformational change that releases the active C subunits.

Active PKA phosphorylates serine and threonine residues on a wide range of target proteins, mediating diverse cellular responses:

• In liver and muscle cells, PKA phosphorylates enzymes that promote glycogen breakdown and inhibit glycogen synthesis
• In adipose cells, PKA activates hormone-sensitive lipase, promoting fat mobilization
• In the nucleus, PKA phosphorylates the transcription factor CREB (cAMP response element-binding protein), which then activates genes containing CRE (cAMP response element) sequences in their promoters

The specificity of PKA action in different cell types depends on which target proteins are expressed and on A-kinase anchoring proteins (AKAPs) that tether PKA to specific subcellular locations, ensuring it phosphorylates the right substrates.

Some G Proteins Signal via Phospholipids

A second major GPCR effector pathway involves the enzyme phospholipase C-β (PLC-β), activated by Gαq. PLC-β cleaves the membrane phospholipid phosphatidylinositol 4,5-bisphosphate (PIP₂) to generate two second messengers simultaneously:

• Inositol 1,4,5-trisphosphate (IP₃) — a small, water-soluble molecule that diffuses into the cytoplasm and binds IP₃ receptors on the endoplasmic reticulum, triggering the release of stored Ca²⁺
• Diacylglycerol (DAG) — a lipid that remains in the plasma membrane and activates protein kinase C (PKC), which phosphorylates various target proteins

This bifurcating pathway simultaneously raises intracellular Ca²⁺ and activates PKC, enabling a coordinated cellular response.

Ca²⁺ Functions as a Ubiquitous Intracellular Mediator

Calcium ions are one of the most widely used intracellular signals. The cytoplasmic Ca²⁺ concentration at rest is very low (~100 nM) compared to the extracellular fluid (~1–2 mM) and the ER lumen (~500 μM). This steep gradient means that even small changes in membrane permeability to Ca²⁺ produce large and rapid increases in cytoplasmic Ca²⁺.

Ca²⁺ signals are detected by Ca²⁺-binding proteins, the most important of which is calmodulin (CaM) — a small, highly conserved protein with four Ca²⁺-binding sites (EF-hand motifs). When Ca²⁺ binds, calmodulin undergoes a conformational change that allows it to wrap around and activate target proteins, including CaM-kinase II (CaMKII), a multifunctional kinase involved in learning and memory, and myosin light-chain kinase (MLCK), which triggers smooth muscle contraction.

Other Ca²⁺-sensing proteins include troponin C in skeletal muscle, synaptotagmin (which triggers neurotransmitter vesicle fusion at synapses), and calcineurin (a Ca²⁺-activated phosphatase that activates the transcription factor NFAT in T cells — the target of the immunosuppressant drugs cyclosporin and FK506).

Feedback Generates Ca²⁺ Waves and Oscillations

Ca²⁺ signals are often not simple step increases but rather oscillations — repeated spikes of Ca²⁺ separated by intervals of low Ca²⁺. These oscillations arise from the interplay of positive and negative feedback on Ca²⁺ release channels:

• Positive feedback: Ca²⁺ released through IP₃ receptors stimulates further Ca²⁺ release from neighboring IP₃ receptors (Ca²⁺-induced Ca²⁺ release, or CICR), generating a regenerative Ca²⁺ spike
• Negative feedback: at higher concentrations, Ca²⁺ inhibits the IP₃ receptor, terminating the spike; Ca²⁺ pumps (SERCA pumps) then restore Ca²⁺ to the ER, resetting the system for another spike

The frequency of Ca²⁺ oscillations encodes information about signal strength: a stronger stimulus increases the frequency of spikes rather than their amplitude. Target proteins like CaMKII can decode this frequency-encoded signal, becoming more active at higher spike frequencies — a mechanism reminiscent of frequency modulation in telecommunications.

Ca²⁺ signals can also propagate spatially as Ca²⁺ waves that travel across the cell via sequential CICR events, coordinating responses over large distances within a single cell (e.g., in fertilized eggs, where a Ca²⁺ wave sweeping across the egg triggers the cortical reaction that prevents polyspermy).

Some G Proteins Directly Regulate Ion Channels

Not all GPCR signaling proceeds through enzymatic cascades. In some cases, the Gβγ subunit released from the trimeric G protein directly binds and opens ion channels. The best-studied example is the action of acetylcholine on the heart: acetylcholine activates muscarinic receptors coupled to Gαi, and the released Gβγ directly opens GIRK channels (G-protein-regulated inwardly rectifying K⁺ channels), hyperpolarizing the cardiac pacemaker cells and slowing the heart rate.

Smell and Vision Depend on GPCRs

Two sensory systems are built almost entirely on GPCR signaling:

Olfaction (smell): Humans express ~400 functional odorant receptor genes (from an original repertoire of ~1,000), each encoding a different GPCR. Each olfactory neuron expresses only one type of odorant receptor. When an odorant molecule binds, the receptor activates Gαolf (a variant of Gαs), which stimulates adenylyl cyclase → cAMP → opening of cyclic-nucleotide-gated (CNG) channels → depolarization and action potential.

Vision (sight): The photoreceptor protein rhodopsin in retinal rod cells is a GPCR. Its ligand is retinal, a derivative of vitamin A that is covalently bound within the receptor. Light absorption converts 11-cis-retinal to all-trans-retinal, activating rhodopsin. Activated rhodopsin stimulates transducin (Gαt), which activates cGMP phosphodiesterase, reducing cGMP levels and closing CNG channels — hyperpolarizing the rod cell. A single photon can activate ~500 transducin molecules, which in turn activate ~500 phosphodiesterase molecules, hydrolyzing ~10⁵ molecules of cGMP — a spectacular example of signal amplification.

Nitric Oxide Gas Signals by Directly Activating an Intracellular Enzyme

Nitric oxide (NO) is an unusual signal molecule: it is a gas and a free radical that diffuses freely across membranes. It is produced from arginine by the enzyme nitric oxide synthase (NOS) and has a half-life of only a few seconds, limiting its range of action.

NO acts by directly binding to and activating soluble guanylyl cyclase, an intracellular enzyme that produces cGMP from GTP. cGMP activates protein kinase G (PKG), which relaxes smooth muscle in blood vessel walls (→ vasodilation). This is the mechanism by which endothelial cells signal to underlying smooth muscle to regulate blood flow. The drug nitroglycerin, used to treat angina, works by releasing NO. The drug sildenafil (Viagra) inhibits the cGMP phosphodiesterase PDE5 in penile smooth muscle, prolonging the NO/cGMP vasodilatory signal.

Second Messengers and Enzymatic Cascades Amplify Signals

A defining feature of intracellular signaling is signal amplification. At each step in a signaling cascade, one activated molecule can activate many downstream molecules:

• One ligand-bound GPCR activates ~100 Gα molecules (catalytic amplification by the receptor)
• Each Gαs activates one adenylyl cyclase, which produces ~1,000 cAMP molecules per second
• Each cAMP activates PKA, which phosphorylates ~10 substrate molecules
• If the substrate is itself an enzyme (e.g., phosphorylase kinase), it activates ~10 molecules of the next enzyme

The net result is that a tiny extracellular signal — even a single photon or a few hormone molecules — can produce a massive intracellular response involving millions of product molecules. This amplification cascade also means that signaling systems must be tightly controlled by negative feedback, phosphatases, and desensitization mechanisms to prevent runaway activation.

We can illustrate the diversity of signaling receptors by comparing coding sequences from different GPCR families:

let camp = "C1C(COP(=O)(O)O)OC(N2C=NC3=C2N=CN=C3N)C1O"
let ip3 = "OC1C(OP(=O)(O)O)C(OP(=O)(O)O)C(O)C(OP(=O)(O)O)C1O"
print("cAMP (second messenger):")
print(Chem.properties(camp))
print("IP3 (second messenger):")
print(Chem.properties(ip3))
let similarity = Chem.tanimoto(camp, ip3)
print("Structural similarity: " + similarity)

cAMP and IP₃ are both second messengers but have very different structures — cAMP is a cyclic nucleotide while IP₃ is a phosphorylated sugar. Their structural dissimilarity reflects their activation of completely different downstream effectors.

GPCR Bioinformatics

GPCR Classification and Phylogenetics

The ~800 human GPCRs are classified into five major families using the GRAFS system (named for the five families):

GPCRdb (gpcrdb.org) is the primary bioinformatics resource for GPCR research. It provides structure-based sequence alignments using the Ballesteros-Weinstein numbering scheme (which assigns a standard residue number to each position in the seven TM helices across all GPCRs), phylogenetic trees, mutant phenotype data, ligand binding data, and structural models for all human GPCRs.

GPCR Structure Prediction and Homology Modeling

The explosion of experimentally determined GPCR structures (from fewer than 5 before 2007 to more than 500 by 2025, thanks to advances in cryo-EM and XFEL crystallography) has transformed GPCR structural biology. Homology modeling — building a three-dimensional model of a GPCR based on the experimentally determined structure of a related family member — is now feasible for most human GPCRs. AlphaFold predictions cover all human GPCRs but may not capture the conformational changes associated with activation. GPCRdb provides precomputed homology models and tools for comparative structural analysis.

GPCR-Ligand Docking and Virtual Screening

Computational docking of small molecules into GPCR binding sites is a mainstay of structure-based drug design. Programs such as AutoDock, Glide, and GOLD predict how drug candidates bind in the orthosteric (natural ligand) or allosteric (modulatory) binding pockets. Virtual screening of compound libraries against GPCR structures has successfully identified novel agonists, antagonists, and biased ligands (which selectively activate G-protein or β-arrestin pathways).

GPCR De-orphanization and Pharmacogenomics

Of the ~800 human GPCRs, roughly 100 remain orphans — their endogenous ligands are unknown. De-orphanization strategies combine bioinformatics (phylogenetic clustering with known GPCRs, co-expression analysis with candidate ligands) and experimental screening (TANGO assays, Ca²⁺ mobilization assays) to identify their natural ligands. Orphan GPCRs represent a large untapped pool of potential drug targets.

GPCR pharmacogenomics studies how genetic variation in GPCR genes affects drug response. Single-nucleotide polymorphisms (SNPs) in GPCR coding regions can alter ligand binding affinity, G-protein coupling efficiency, or desensitization kinetics. For example, variants in the β₂-adrenergic receptor gene (ADRB2) influence the response to β-agonist bronchodilators in asthma patients.

Let us compare a conserved functional motif across GPCR subfamilies. The DRY motif (Asp-Arg-Tyr) at the cytoplasmic end of TM3 is one of the most highly conserved features in class A GPCRs and plays a critical role in receptor activation and G-protein coupling:

let dose_response = '[{"x": 0.01, "y": 5, "label": "0.01 nM"}, {"x": 0.1, "y": 15, "label": "0.1 nM"}, {"x": 1, "y": 50, "label": "1 nM (EC50)"}, {"x": 10, "y": 85, "label": "10 nM"}, {"x": 100, "y": 95, "label": "100 nM"}]'
let plot = Viz.scatter(dose_response, '{"title": "Dose-Response Curve (Hill coefficient = 1)", "x_label": "Ligand concentration (nM)", "y_label": "Response (%)", "color": "#8B5CF6"}')
print(plot)

The sigmoidal shape of the dose-response curve is characteristic of receptor-ligand binding governed by the law of mass action. The EC50 (concentration producing 50% maximal response) and the Hill coefficient (steepness of the curve) are key pharmacological parameters used to characterize GPCR-ligand interactions.

Exercise: Classify Signaling Molecules

Signaling molecules are classified by whether they can cross the plasma membrane. Use chemical property analysis to determine which of these three molecules signals through a surface receptor vs. an intracellular receptor:

let insulin_fragment = "GIVEQCCTSICSLYQLENYCN"
let testosterone = "CC12CCC3C(CCC4=CC(=O)CCC34C)C1CCC2O"
let histamine = "NCCc1c[nH]cn1"
print("Insulin (peptide):")
print(Struct.protein_props(insulin_fragment))
print("Testosterone (steroid):")
print(Chem.properties(testosterone))
print("Histamine (amine):")
print(Chem.properties(histamine))
// Which molecule crosses the membrane to bind an intracellular receptor?
let answer = "testosterone"
print(answer)

Exercise: Compare Second Messenger Properties

Second messengers amplify receptor signals inside the cell. Compare the chemical properties of cAMP and diacylglycerol (DAG) to understand why cAMP diffuses through the cytosol while DAG stays in the membrane:

let camp = "C1C(COP(=O)(O)O)OC(N2C=NC3=C2N=CN=C3N)C1O"
let dag = "OCC(COC(=O)CCCCCCCCCCCCCCC)OC(=O)CCCCCCCCCCCCCCC"
print("cAMP properties:")
print(Chem.properties(camp))
print("DAG properties:")
print(Chem.properties(dag))
// Which second messenger is membrane-bound?
let answer = "DAG"
print(answer)

Exercise: Signal Amplification Cascade

A single hormone molecule can trigger the production of millions of product molecules through enzymatic amplification at each step. Visualize the amplification at each step of the epinephrine→cAMP→PKA cascade:

let cascade = '[{"label": "Epinephrine", "value": 1}, {"label": "Activated receptors", "value": 10}, {"label": "G-protein (Gs)", "value": 100}, {"label": "Adenylyl cyclase", "value": 1000}, {"label": "cAMP molecules", "value": 10000}, {"label": "PKA activated", "value": 100000}, {"label": "Phosphorylase b→a", "value": 1000000}, {"label": "Glucose released", "value": 10000000}]'
let chart = Viz.bar(cascade, '{"title": "Signal Amplification Cascade", "color": "#EF4444"}')
print(chart)
// How many glucose molecules are released per hormone molecule?
let answer = "10000000"
print(answer)

Knowledge Check

Summary

In this lesson you covered the principles of cell signaling and G-protein-coupled receptor pathways:

• Extracellular signal molecules bind specific receptors; hydrophobic signals (steroids, NO) cross the membrane to intracellular receptors, while hydrophilic signals bind cell-surface receptors
• Signaling modes include endocrine (long-range via blood), paracrine (local diffusion), synaptic (neurotransmitter at synapses), autocrine (self-stimulation), and contact-dependent (membrane-bound ligands)
• Gap junctions allow direct sharing of small signaling molecules (Ca²⁺, cAMP, IP₃) between coupled cells
• Each cell responds to a limited set of signals determined by its receptor expression; the same signal can trigger different responses in different cell types
• Three major receptor classes transduce signals at the cell surface: ion-channel-coupled, G-protein-coupled (GPCRs), and enzyme-coupled receptors
• Intracellular signaling relies on signaling proteins (kinases, GTPases, adaptors) and small second messengers (cAMP, Ca²⁺, IP₃, DAG) organized into modular pathways
• Specificity in the noisy cytoplasm is maintained by scaffold proteins, modular interaction domains (SH2, SH3, PH, PTB, PDZ), and spatial compartmentalization
• Signaling complexes assemble at activated receptors, nucleated by phosphotyrosine docking sites or β-arrestin scaffolds
• Signal-response relationships range from graded to ultrasensitive to bistable, shaped by cooperative binding, multisite phosphorylation, and feedback loops
• Response speed depends on signaling molecule turnover — fast turnover (phosphodiesterases, phosphatases) enables rapid responses; slow turnover produces sustained responses
• Positive feedback generates all-or-none switches with hysteresis; negative feedback limits response duration, generates oscillations, and improves robustness
• Desensitization (receptor phosphorylation by GRKs, endocytosis, downregulation) allows cells to respond to changes in signal strength rather than absolute levels
• GPCRs signal through trimeric G proteins (Gα, Gβ, Gγ) via a GTPase cycle; different Gα subtypes (Gαs, Gαi, Gαq, Gα12/13) activate different effectors
• cAMP signaling: Gαs → adenylyl cyclase → cAMP → PKA → phosphorylation of targets (metabolic enzymes, CREB transcription factor); terminated by phosphodiesterases
• Phospholipid signaling: Gαq → PLC-β → IP₃ (Ca²⁺ release from ER) + DAG (PKC activation)
• Ca²⁺ signaling uses the steep Ca²⁺ gradient (100 nM cytoplasmic vs. 1–2 mM extracellular) for rapid signal transduction; calmodulin and CaMKII decode Ca²⁺ signals; oscillations encode signal strength as frequency
• Gβγ can directly regulate ion channels (e.g., GIRK channels that slow the heart)
• Smell and vision are built on GPCR signaling: odorant receptors → Gαolf → cAMP → CNG channels; rhodopsin → transducin → cGMP PDE → CNG channels
• Nitric oxide is a gaseous signal that activates soluble guanylyl cyclase → cGMP → PKG → smooth muscle relaxation (vasodilation)
• Signal amplification through enzymatic cascades enables detection of minute signals (single photons, few hormone molecules)
• Pathway databases (KEGG, Reactome, SignaLink, NetPath) and kinase-substrate databases (PhosphoSitePlus) are essential for phosphoproteomics analysis and signaling network modeling
• Signaling network modeling uses Boolean, ODE, and agent-based approaches to capture emergent properties like oscillations and bistability
• Ligand-receptor databases (CellTalkDB, DLRP) and drug-target databases (DrugBank, ChEMBL) support drug-target interaction prediction
• GPCRdb provides classification, structure-based alignments (Ballesteros-Weinstein numbering), and pharmacological data for all human GPCRs
• GPCR pharmacogenomics links genetic variants in receptor genes to drug response variability; ~100 orphan GPCRs remain as potential drug targets