State-dependent dynamics of extramembrane domains in the ß2 -adrenergic receptor.

G protein-coupled receptors (GPCRs) are membrane-bound signaling proteins that play an essential role in cellular signaling processes. Due to their intrinsic function of transmitting internal signals in response to external cues, these receptors are adapted to be highly dynamic in nature. The ß2 -adrenergic receptor (ß2 AR) is a representative member of the family that has been extensively analyzed in terms of its structure and activation. Although the structure of the transmembrane domain has been characterized in the different functional states of the receptor, the conformational dynamics of the extramembrane domains, especially the intrinsically disordered regions are still emerging. In this study, we analyze the state-dependent dynamics of extramembrane domains of ß2 AR using atomistic molecular dynamics simulations. We introduce a parameter, the residue excess dynamics that allows us to better quantify receptor dynamics. Using this measure, we show that the dynamics of the extramembrane domains are sensitive to the receptor state. Interestingly, the ligand-bound intermediate R ' state shows the maximal dynamics compared to either the active R*G or inactive R states. Ligand binding appears to be correlated with high residue excess dynamics that are dampened upon G protein coupling. The intracellular loop-3 (ICL3) domain has a tendency to flip towards the membrane upon ligand binding, which could contribute to receptor "priming." We highlight an important ICL1-helix-8 interplay that is broken in the ligand-bound state but is retained in the active state. Overall, our study highlights the importance of characterizing the functional dynamics of the GPCR loop domains.

Molecular determinants of GPCR pharmacogenetics: Deconstructing the population variants in ß2-adrenergic receptor.

G protein-coupled receptors (GPCRs) are membrane proteins that play a central role in cell signaling and constitute one of the largest classes of drug targets. The molecular mechanisms underlying GPCR function have been characterized by several experimental and computational methods and provide an understanding of their role in physiology and disease. Population variants arising from nsSNPs affect the native function of GPCRs and have been implicated in differential drug response. In this chapter, we provide an overview on GPCR structure and activation, with a special focus on the ß2-adrenergic receptor (ß2-AR). First, we discuss the current understanding of the structural and dynamic features of the wildtype receptor. Subsequently, the population variants identified in this receptor from clinical and large-scale genomic studies are described. We show how computational approaches such as bioinformatics tools and molecular dynamics simulations can be used to characterize the variant receptors in comparison to the wildtype receptor. In particular, we discuss three examples of clinically important variants and discuss how the structure and function of these variants differ from the wildtype receptor at a molecular level. Overall, the chapter provides an overview of structure and function of GPCR variants and is a step towards the study of inter-individual differences and personalized medicine.

Molecular Characterization of the ß2-like Octopamine Receptor of Helicoverpa armigera.

Helicoverpa armigera is a devastating polyphagous and cosmopolitan crop pest. There are reports of this insect being resistant to a variety of pesticides raising concern worldwide. The Octopamine (OA) binding ß2-like receptor (OAR), a GPCR, is widely distributed in the nervous system of the insect and plays essential roles in the physiology and development and thus is an important target for insecticides. Yet, the molecular characterization of the H. armigera OAR (HarmOAR) and rational design of compounds based on this receptor is lacking. As a first step, we performed multiple sequence alignment of all insect OARs, which revealed that the sequences contained all conserved class A GPCR motifs. Phylogenetic studies showed clade-specific variations in the protein sequences primarily arising owing to differences in the ICL3 loop region. Further, a structural model of HarmOAR was built using the inactive human ß2AR as a template. 0.9 µs atomistic simulations revealed conserved inter helical contacts and water molecules of HarmOAR. The detailed binding of octopamine was studied using molecular docking and 0.3 µs atomistic simulations. Twenty-two insecticides active against octopamine receptors of other insects were compiled and docked to HarmOAR followed by rescoring with binding free energies to prioritize them for H. armigera. Our study suggests α-terpineol to be a good candidate as an insecticide or insect repellent for Helicoverpa armigera.

Loss of a water-mediated network results in reduced agonist affinity in a ß2-adrenergic receptor clinical variant.

The ß2-adrenergic receptor (ß2AR) is a member of the G protein-coupled receptor (GPCR) family that is an important drug target for asthma and COPD. Clinical studies coupled with biochemical data have identified a critical receptor variant, Thr164Ile, to have a reduced response to agonist-based therapy, although the molecular mechanism underlying this seemingly "non-deleterious" substitution is not clear. Here, we couple molecular dynamics simulations with network analysis and free-energy calculations to identify the molecular determinants underlying the differential drug response. We are able to identify hydration sites in the transmembrane domain that are essential to maintain the integrity of the binding site but are absent in the variant. The loss of these hydration sites in the variant correlates with perturbations in the intra-protein interaction network and rearrangements in the orthosteric ligand binding site. In conjunction, we observe an altered binding and reduced free energy of a series of agonists, in line with experimental trends. Our work identifies a functional allosteric pathway connected by specific hydration sites in ß2AR that has not been reported before and provides insight into water-mediated networks in GPCRs in general. Overall, the work is one of the first step towards developing variant-specific potent and selective agonists.

In vitro and in silico studies on membrane interactions of diverse Capsicum annuum flower γ-thionin peptides.

Thionins are small, cysteine-rich peptides that play an important role in plant defense, primarily through their interactions with membranes. Eight novel γ-thionin peptides (CanThio1-8) were isolated from the flower of Capsicum annuum. Sequence analysis revealed that the peptides cluster into three groups. A representative peptide from each group (CanThio1, 2, and 3) was used for experimental characterization. Interestingly, peptides were found to possess some cytotoxic activity against normal human embryonic kidney cell line but higher cytotoxicity against cancer cell line MCF-7. CanThio3 peptide was chosen as a representative peptide to study the molecular mechanism of action on membranes. Microsecond timescale atomistic simulations of CanThio3 were performed in the presence of a POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) lipid bilayer. Simulations revealed that CanThio3 interacts with the bilayer and causes lipid thinning in the vicinity. Nonpolar amino acids specific to the α-core region of CanThio3 along with nonpolar residues in the γ-core region are seen to interact with the lipid tails. The differences in the amino acid sequence of CanThio peptides in these regions explain the variability in cytotoxic activities. In summary, our results demonstrate the membrane-mediated activity of a novel series of γ-thionin peptides from C. annuum.

Differential Dynamics Underlying the Gln27Glu Population Variant of the ß2-Adrenergic Receptor.

The ß2-adrenergic receptor (ß2AR) is a membrane-bound G-protein-coupled receptor and an important drug target for asthma. Clinical studies report that the population variant Gln27Glu is associated with a differential response to common asthma drugs, such as albuterol, isoproterenol and terbutaline. Interestingly, the 27th amino acid is positioned on the N-terminal region that is the most flexible and consequently the least studied part of the receptor. In this study, we probe the molecular origin of the differential drug binding by performing structural modeling and simulations of the wild-type (Gln) and variant (Glu) receptors followed by ensemble docking with the ligands, albuterol, isoproterenol and terbutaline. In line with clinical studies, the ligands were observed to interact preferentially with the Glu variant. Our results indicate that the Glu residue at the 27th position perturbs the network of electrostatic interactions that connects the N-terminal region to the binding site in the wild-type receptor. As a result, the Glu variant is observed to bind better to the three ligands tested in this study. Our study provides a structural basis to explain the variable drug response associated with the 27th position polymorphism in the ß2AR and is a starting step to identify genotype-specific therapeutics.