Communication over the network of binary switches regulates the activation of A2A adenosine receptor

Kinetic Model for the Desensitization of G Protein-Coupled Receptor

Desensitization of G-protein-coupled receptors (GPCR) is a general regulatory mechanism adopted by biological organisms against overstimulation of G protein signaling. Although the details of the mechanism are extensively studied, it is not easy to gain an overarching understanding of the process constituted by a multitude of molecular events with vastly differing time scales. To offer a semiquantitative yet predictive understanding of the mechanism, we formulate a kinetic model for the G protein signaling and desensitization by considering essential biochemical steps from ligand binding to receptor internalization. The internalization, followed by receptor depletion from the plasma membrane, attenuates the downstream signal. Together with the kinetic model and its full numerics of the expression derived for the dose-response relation, an approximated form of the expression clarifies the role played by the individual biochemical processes and allows us to identify four distinct regimes for the downregulation that emerge from the balance between phosphorylation, dephosphorylation, and the cellular level of β-arrestin.

General Chemical Reaction Network Theory for Olfactory Sensing Based on G-Protein-Coupled Receptors: Elucidation of Odorant Mixture Effects and Agonist-Synergist Threshold

This work presents a general chemical reaction network theory for olfactory sensing processes that employ G-protein-coupled receptors as olfactory receptors (ORs). The theory can be applied to general mixtures of odorants and an arbitrary number of ORs. Reactions of ORs with G-proteins, in both the presence and absence of odorants, are explicitly considered. A unique feature of the theory is the definition of an odor activity vector consisting of strengths of odorant-induced signals from ORs relative to those due to background G-protein activity in the absence of odorants. It is demonstrated that each component of the odor activity defined this way reduces to a Michaelis-Menten form capable of accounting for cooperation or competition effects between different odorants. The main features of the theory are illustrated for a two-odorant mixture. Known and potential mixture effects, such as suppression, shadowing, inhibition, and synergy, are quantitatively described. Effects of relative values of rate constants, basal activity, and G-protein concentration are also demonstrated.

MD Simulations Revealing Special Activation Mechanism of Cannabinoid Receptor 1

Cannabinoid receptor 1 (CB1) is a G protein-coupled receptor (GPCR) that is gaining much interest for its regulating role in the central nervous system and its value as a drug target. Structures of CB1 in inactive and active states have revealed conformational change details that are not common in other GPCRs. Here, we performed molecular dynamics simulations of CB1 in different ligand binding states and with mutations to reveal its activation mechanism. The conformational change of the “twin toggle switch” residues F200^3.36 and W356^6.48 that correlates with ligand efficacy is identified as a key barrier step in CB1 activation. Similar conformational change of residues 3.36/6.48 is also observed in melanocortin receptor 4, showing this “twin toggle switch” residue pair is crucial for the activation of multiple GPCR members.

Activation-Induced Reorganization of Energy Transport Networks in the β2 Adrenergic Receptor

We compute energy exchange networks (EENs) through the β₂ adrenergic receptor (β₂AR), a G-protein coupled receptor (GPCR), in inactive and active states, based on the results of molecular dynamics simulations of this membrane bound protein. We introduce a new definition for the reorganization of EENs upon activation that depends on the relative change in rates of energy transfer across noncovalent contacts throughout the protein. On the basis of the reorganized network that we obtain for β₂AR upon activation, we identify a branched pathway between the agonist binding site and the cytoplasmic region, where a G-protein binds to the receptor when activated. The pathway includes all of the motifs containing molecular switches previously identified as contributing to the allosteric transition of β₂AR upon agonist binding. EENs and their reorganization upon activation are compared with structure-based contact networks computed for the inactive and active states of β₂AR.

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

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

Transmembrane 4 L Six Family Member 5 Senses Arginine for mTORC1 Signaling

The mechanistic target of rapamycin complex (mTORC1) is a signaling hub on the lysosome surface, responding to lysosomal amino acids. Although arginine is metabolically important, the physiological arginine sensor that activates mTOR remains unclear. Here, we show that transmembrane 4 L six family member 5 (TM4SF5) translocates from plasma membrane to lysosome upon arginine sufficiency and senses arginine, culminating in mTORC1/S6K1 activation. TM4SF5 bound active mTOR upon arginine sufficiency and constitutively bound amino acid transporter SLC38A9. TM4SF5 binding to the cytosolic arginine sensor Castor1 decreased upon arginine sufficiency, thus allowing TM4SF5-mediated sensing of metabolic amino acids. TM4SF5 directly bound free L-arginine via its extracellular loop possibly for the efflux, being supported by mutant study and homology and molecular docking modeling. Therefore, we propose that lysosomal TM4SF5 senses and enables arginine efflux for mTORC1/S6K1 activation, and arginine-auxotroph in hepatocellular carcinoma may be targeted by blocking the arginine sensing using anti-TM4SF5 reagents.

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

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

Approaches for Differentiation and Interconverting GPCR Agonists and Antagonists

Predicting the functional preferences of the ligands was always a highly demanding task, much harder that predicting whether a ligand can bind to the receptor. This is because of significant similarities of agonists, antagonists and inverse agonists which are binding usually in the same binding site of the receptor and only small structural changes can push receptor toward a particular activation state. For G protein-coupled receptors, due to a large progress in crystallization techniques and also in receptor thermal stabilization, it was possible to obtain a large number of high-quality structures of complexes of these receptors with agonists and non-agonists. Additionally, the long-time-scale molecular dynamics simulations revealed how the activation processes of GPCRs can take place. Using both theoretical and experimental knowledge it was possible to employ many clever and sophisticated methods which can help to differentiate agonists and non-agonists, so one can interconvert them in search of the optimal drug.

Implications for human odor sensing revealed from the statistics of odorant-receptor interactions

Binding of odorants to olfactory receptors (ORs) elicits downstream chemical and neural signals, which are further processed to odor perception in the brain. Recently, Mainland and colleagues have measured more than 500 pairs of odorant-OR interaction by a high-throughput screening assay method, opening a new avenue to understanding the principles of human odor coding. Here, using a recently developed minimal model for OR activation kinetics, we characterize the statistics of OR activation by odorants in terms of three empirical parameters: the half-maximum effective concentration EC₅₀, the efficacy, and the basal activity. While the data size of odorants is still limited, the statistics offer meaningful information on the breadth and optimality of the tuning of human ORs to odorants, and allow us to relate the three parameters with the microscopic rate constants and binding affinities that define the OR activation kinetics. Despite the stochastic nature of the response expected at individual OR-odorant level, we assess that the confluence of signals in a neuron released from the multitude of ORs is effectively free of noise and deterministic with respect to changes in odorant concentration. Thus, setting a threshold to the fraction of activated OR copy number for neural spiking binarizes the electrophysiological signal of olfactory sensory neuron, thereby making an information theoretic approach a viable tool in studying the principles of odor perception.

Despite the decades of research, quantitative details of human olfaction have remained largely unexplored. However, a high-throughput measurement has recently been carried out to produce dose-response data between a set of odorants and a repertoire of human olfactory receptors. We characterized each pair of odorant-receptor interaction in terms of EC50, efficacy, and basal level, a strategy often adopted in biochemical, pharmacological sciences to describe the response of receptors to cognate ligands. The distributions of EC50 values and efficacies acquired from the analysis provide glimpses into how human olfactory receptors are tuned to odorants. Specifically, the response of human ORs is optimized around ∼ 100μM of odorant. Next, the efficacies of OR responses to odorants are bi-exponentially distributed, which indicates that the strength of odorant-OR interaction is classified into strong and weak subgroups. By showing that the stochastic response of individual receptor to odorant can effectively be binarized at cellular level through olfactory processes, we also provide a theoretical basis for an information theoretical approach in studying the principles of odor perception.

Understanding G Protein-Coupled Receptor Allostery via Molecular Dynamics Simulations: Implications for Drug Discovery

Unraveling the mystery of protein allostery has been one of the greatest challenges in both structural and computational biology. However, recent advances in computational methods, particularly molecular dynamics (MD) simulations, have led to its utility as a powerful and popular tool for the study of protein allostery. By capturing the motions of a protein's constituent atoms, simulations can enable the discovery of allosteric hot spots and the determination of the mechanistic basis for allostery. These structural and dynamic studies can provide a foundation for a wide range of applications, including rational drug design and protein engineering. In our laboratory, the use of MD simulations and network analysis assisted in the elucidation of the allosteric hotspots and intracellular signal transduction of G protein-coupled receptors (GPCRs), primarily on one of the adenosine receptor subtypes, A_2A adenosine receptor (A_2AAR). In this chapter, we describe a method for calculating the map of allosteric signal flow in different GPCR conformational states and illustrate how these concepts have been utilized in understanding the mechanism of GPCR allostery. These structural studies will provide valuable insights into the allosteric and orthosteric modulations that would be of great help to design novel drugs targeting GPCRs in pathological states.

Identification of a Different Agonist-Binding Site and Activation Mechanism of the Human P2Y1 Receptor

The human P2Y1 receptor (P2Y1R) is a purinergic G-protein-coupled receptor (GPCR) that functions as a receptor for adenosine 5'-diphosphate (ADP). An antagonist of P2Y1R might potentially have antithrombotic effects, whereas agonists might serve as antidiabetic agents. On the basis of the antagonist-bound MRS2500-P2Y1R crystal structure, we constructed computational models of apo-P2Y1R and the agonist-receptor complex 2MeSADP-P2Y1R. We then performed conventional molecular dynamics (cMD) and accelerated molecular dynamics (aMD) simulations to study the conformational dynamics after binding with agonist/antagonist as well as the P2Y1R activation mechanism. We identified a new agonist-binding site of P2Y1R that is consistent with previous mutagenesis data. This new site is deeper than those of the agonist ADP in the recently simulated ADP-P2Y1R structure and the antagonist MRS2500 in the MRS2500-P2Y1R crystal structure. During P2Y1R activation, the cytoplasmic end of helix VI shifts outward 9.1 Å, the Ser1463.47-Tyr2375.58 hydrogen bond breaks, a Tyr2375.58-Val2626.37 hydrogen bond forms, and the conformation of the χ1 rotamer of Phe2696.44 changes from parallel to perpendicular to helix VI. The apo-P2Y1R system and the MRS2500-P2Y1R system remain inactive. The newly identified agonist binding site and activation mechanism revealed in this study may aid in the design of P2Y1R antagonists/agonists as antithrombotic/antidiabetic agents, respectively.

Activation-Induced Conformational Changes of Dopamine D3 Receptor Promote the Formation of the Internal Water Channel

The atomic-level dopamine activation mechanism for transmitting extracellular ligand binding events through transmembrane helices to the cytoplasmic G protein remains unclear. In the present study, the complete dopamine D3 receptor (D3R), with a homology-modeled N-terminus, was constructed to dock different ligands to simulate conformational alterations in the receptor’s active and inactive forms during microsecond-timescale molecular dynamic simulations. In agonist-bound systems, the D3R N-terminus formed a “lid-like” structure and lay flat on the binding site opening, whereas in antagonist and inverse agonist-bound systems, the N-terminus exposed the binding cavity. Receptor activation was characterized using the different molecular switch residue distances, and G protein-binding site volumes. A continuous water pathway was observed only in the dopamine-G_αi-bound system. In the inactive D3Rs, water entry was hindered by the hydrophobic layers. Finally, a complete activation mechanism of D3R was proposed. Upon agonist binding, the “lid-like” conformation of the N-terminus induces a series of molecular switches to increase the volume of the D3R cytoplasmic binding part for G protein association. Meanwhile, water enters the transmembrane region inducing molecular switches to assist in opening the hydrophobic layers to form a continuous water channel, which is crucial for maintaining a fully active conformation for signal transduction.

Signaling within Allosteric Machines: Signal Transmission Pathways Inside G Protein-Coupled Receptors

In recent years, our understanding of function of G protein-coupled receptors (GPCRs) has changed from a picture of simple signal relays, transmitting only a particular signal to a particular G protein heterotrimer, to versatile machines, capable of various responses to different stimuli and being modulated by various factors. Some recent reports provide not only the data on ligands/modulators and resultant signals induced by them, but also deeper insights into exact pathways of signal migration and mechanisms of signal transmission through receptor structure. Combination of these computational and experimental data sheds more light on underlying mechanisms of signal transmission and signaling bias in GPCRs. In this review we focus on available clues on allosteric pathways responsible for complex signal processing within GPCRs structures, with particular emphasis on linking compatible in silico- and in vitro-derived data on the most probable allosteric connections.

Ultraslow Water-Mediated Transmembrane Interactions Regulate the Activation of A2A Adenosine Receptor

Water molecules inside a G-protein coupled receptor (GPCR) have recently been spotlighted in a series of crystal structures. To decipher the dynamics and functional roles of internal water molecules in GPCR activity, we studied the A2A adenosine receptor using microsecond molecular-dynamics simulations. Our study finds that the amount of water flux across the transmembrane (TM) domain varies depending on the receptor state, and that the water molecules of the TM channel in the active state flow three times more slowly than those in the inactive state. Depending on the location in solvent-protein interface as well as the receptor state, the average residence time of water in each residue varies from ∼O(10(2)) ps to ∼O(10(2)) ns. Especially, water molecules, exhibiting ultraslow relaxation (∼O(10(2)) ns) in the active state, are found around the microswitch residues that are considered activity hotspots for GPCR function. A continuous allosteric network spanning the TM domain, arising from water-mediated contacts, is unique in the active state, underscoring the importance of slow water molecules in the activation of GPCRs.

Biased Allostery

G-protein-coupled receptors (GPCRs) constitute a large group of integral membrane proteins that transduce extracellular signals from a wide range of agonists into targeted intracellular responses. Although the responses can vary depending on the category of G-proteins activated by a particular receptor, responses were also found to be triggered by interactions of the receptor with β-arrestins. It was subsequently discovered that for the same receptor molecule (e.g., the β-adrenergic receptor), some agonists have a propensity to specifically favor responses by G-proteins, others by β-arrestins, as has now been extensively studied. This feature of the GPCR system is known as biased agonism and is subject to various interpretations, including agonist-induced conformational change versus selective stabilization of preexisting active conformations. Here, we explore a complete allosteric framework for biased agonism based on alternative preexisting conformations that bind more strongly, but nonexclusively, either G-proteins or β-arrestins. The framework incorporates reciprocal effects among all interacting molecules. As a result, G-proteins and β-arrestins are in steric competition for binding to the cytoplasmic surface of either the G-protein-favoring or β-arrestin-favoring GPCR conformation. Moreover, through linkage relations, the strength of the interactions of G-proteins or β-arrestins with the corresponding active conformation potentiates the apparent affinity for the agonist, effectively equating these two proteins to allosteric modulators. The balance between response alternatives can also be influenced by the physiological concentrations of either G-proteins or β-arrestins, as well as by phosphorylation or interactions with positive or negative allosteric modulators. The nature of the interactions in the simulations presented suggests novel experimental tests to distinguish more fully among alternative mechanisms.

Gorge Motions of Acetylcholinesterase Revealed by Microsecond Molecular Dynamics Simulations

Acetylcholinesterase, with a deep, narrow active-site gorge, attracts enormous interest due to its particularly high catalytic efficiency and its inhibitors used for treatment of Alzheimer’s disease. To facilitate the massive pass-through of the substrate and inhibitors, “breathing” motions to modulate the size of the gorge are an important prerequisite. However, the molecular mechanism that governs such motions is not well explored. Here, to systematically investigate intrinsic motions of the enzyme, we performed microsecond molecular dynamics simulations on the monomer and dimer of Torpedo californica acetylcholinesterase (TcAChE) as well as the complex of TcAChE bound with the drug E2020. It has been revealed that protein-ligand interactions and dimerization both keep the gorge in bulk, and opening events of the gorge increase dramatically compared to the monomer. Dynamics of three subdomains, S3, S4 and the Ω-loop, are tightly associated with variations of the gorge size while the dynamics can be changed by ligand binding or protein dimerization. Moreover, high correlations among these subdomains provide a basis for remote residues allosterically modulating the gorge motions. These observations are propitious to expand our understanding of protein structure and function as well as providing clues for performing structure-based drug design.

Kinetic Model for the Activation of Mammalian Olfactory Receptor

The sense of smell is triggered by binding of odorants to a set of olfactory receptors (ORs), the activation of which generates specific patterns of neuronal signals in olfactory bulbs. Despite a long history of research and speculations, very little is known about the actual mechanism of OR activation. In particular, there is virtually no theoretical framework capable of describing the kinetics of olfactory activation at a quantitative level. Based on the fact that mammalian ORs belong to a class of G-protein coupled receptors (GPCRs) and utilizing the information available from recent studies on other types of GPCRs with known structural data, we construct a minimal kinetic model for mammalian olfactory activation, obtaining a new expression for the signal strength as a function of odorant and G-protein concentrations and defining this as odor activity (OA). The parametric dependence of OA on equilibrium dissociation and rate constants provides a new comprehensive means to describe how odorant-OR binding kinetics affects the odor signal, and offers new quantitative criteria for classifying agonistic, partially agonistic, and antagonistic (or inverse agonistic) behavior. The dependence of OA on the concentration of G-proteins also suggests a new experimental method to determine key equilibrium constants for odorant-OR and G-protein-OR association/dissociation processes.

Revealing Atomic-Level Mechanisms of Protein Allostery with Molecular Dynamics Simulations

Molecular dynamics (MD) simulations have become a powerful and popular method for the study of protein allostery, the widespread phenomenon in which a stimulus at one site on a protein influences the properties of another site on the protein. By capturing the motions of a protein's constituent atoms, simulations can enable the discovery of allosteric binding sites and the determination of the mechanistic basis for allostery. These results can provide a foundation for applications including rational drug design and protein engineering. Here, we provide an introduction to the investigation of protein allostery using molecular dynamics simulation. We emphasize the importance of designing simulations that include appropriate perturbations to the molecular system, such as the addition or removal of ligands or the application of mechanical force. We also demonstrate how the bidirectional nature of allostery-the fact that the two sites involved influence one another in a symmetrical manner-can facilitate such investigations. Through a series of case studies, we illustrate how these concepts have been used to reveal the structural basis for allostery in several proteins and protein complexes of biological and pharmaceutical interest.

DIRECT-ID: An automated method to identify and quantify conformational variations--application to β2 -adrenergic GPCR

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT-ID, was used to identify relevant ligand-modulated structural variations in the β2 -adrenergic (β2 AR) G-protein coupled receptor. Micro-second MD simulations of the β2 AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI-167107 (agonist), epinephrine (agonist), salbutamol (long-acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β2 AR differently and DIRECT-ID analysis of the inverse-agonist vs. agonist-modulated β2 AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand-specific microswitches across the GPCR. This study demonstrates the utility of DIRECT-ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems.