Dataset Information

Homology modeling of a Class A GPCR in the inactive conformation: A quantitative analysis of the correlation between model/template sequence identity and model accuracy.

ABSTRACT: With the present work we quantitatively studied the modellability of the inactive state of Class A G protein-coupled receptors (GPCRs). Specifically, we constructed models of one of the Class A GPCRs for which structures solved in the inactive state are available, namely the ?2 AR, using as templates each of the other class members for which structures solved in the inactive state are also available. Our results showed a detectable linear correlation between model accuracy and model/template sequence identity. This suggests that the likely accuracy of the homology models that can be built for a given receptor can be generally forecasted on the basis of the available templates. We also probed whether sequence alignments that allow for the presence of gaps within the transmembrane domains to account for structural irregularities afford better models than the classical alignment procedures that do not allow for the presence of gaps within such domains. As our results indicated, although the overall differences are very subtle, the inclusion of internal gaps within the transmembrane domains has a noticeable a beneficial effect on the local structural accuracy of the domain in question.

SUBMITTER: Costanzi S

PROVIDER: S-EPMC5138091 | biostudies-literature | 2016 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Homology modeling of a Class A GPCR in the inactive conformation: A quantitative analysis of the correlation between model/template sequence identity and model accuracy.

Costanzi Stefano S Skorski Matthew M Deplano Alessandro A Habermehl Brett B Mendoza Mary M Wang Keyun K Biederman Michelle M Dawson Jessica J Gao Jia J

Journal of molecular graphics & modelling 20161004

With the present work we quantitatively studied the modellability of the inactive state of Class A G protein-coupled receptors (GPCRs). Specifically, we constructed models of one of the Class A GPCRs for which structures solved in the inactive state are available, namely the β<sub>2</sub> AR, using as templates each of the other class members for which structures solved in the inactive state are also available. Our results showed a detectable linear correlation between model accuracy and model/t ...[more]

PMID: 27723562

Similar Datasets

Project description:G protein-coupled receptors (GPCR) are integral membrane proteins of considerable interest as targets for drug development. GPCR ligand interaction studies often have a starting point with either crystal structures or comparative models. The majority of GPCR do not have experimentally-characterized 3-dimensional structures, so comparative modeling, also called homology modeling, is a good structure-based starting point. Comparative modeling is a widely used method for generating models of proteins with unknown structures by analogy to crystallized proteins that are expected to exhibit structural conservation. Traditionally, comparative modeling template selection is based on global sequence identity and shared function. However high sequence identity localized to the ligand binding pocket may produce better models to examine protein-ligand interactions. This in silico benchmark study examined the performance of a global versus local similarity measure applied to comparative modeling template selection for 6 previously crystallized, class A GCPR (CXCR4, FFAR1, NOP, P2Y12, OPRK, and M1) with the long-term goal of optimizing GPCR ligand identification efforts. Comparative models were generated from templates selected using both global and local similarity measures. Similarity to reference crystal structures was reflected in RMSD values between atom positions throughout the structure or localized to the ligand binding pocket. Overall, models deviated from the reference crystal structure to a similar degree regardless of whether the template was selected using a global or local similarity measure. Ligand docking simulations were performed to assess relative performance in predicting protein-ligand complex structures and interaction networks. Calculated RMSD values between ligand poses from docking simulations and crystal structures indicate that models based on locally selected templates give docked poses that better mimic crystallographic ligand positions than those based on globally-selected templates in five of the six benchmark cases. However, protein model refinement strategies in advance of ligand docking applications are clearly essential as the average RMSD between crystallographic poses and poses docked into local template models was 9.7 Å and typically less than half of the ligand interaction sites are shared between the docked and crystallographic poses. These data support the utilization of local similarity measures to guide template selection in protocols using comparative models to investigate ligand-receptor interactions.

Project description:G protein-coupled receptors (GPCR) comprise the largest family of membrane proteins and are of considerable interest as targets for drug development. However, many GPCR structures remain unsolved. To address the structural ambiguity of these receptors, computational tools such as homology modeling and loop modeling are often employed to generate predictive receptor structures. Here we combined both methods to benchmark a protocol incorporating homology modeling based on a locally selected template and extracellular loop modeling that additionally evaluates the presence of template ligands during these modeling steps. Ligands were also docked using three docking methods and two pose selection methods to elucidate an optimal ligand pose selection method. Results suggest that local template-based homology models followed by loop modeling produce more accurate and predictive receptor models than models produced without loop modeling, with decreases in average receptor and ligand RMSD of 0.54 Å and 2.91 Å, respectively. Ligand docking results showcased the ability of MOE induced fit docking to produce ligand poses with atom root-mean-square deviation (RMSD) values at least 0.20 Å lower (on average) than the other two methods benchmarked in this study. In addition, pose selection methods (software-based scoring, ligand complementation) selected lower RMSD poses with MOE induced fit docking than either of the other methods (averaging at least 1.57 Å lower), indicating that MOE induced fit docking is most suited for docking into GPCR homology models in our hands. In addition, target receptor models produced with a template ligand present throughout the modeling process most often produced target ligand poses with RMSD values ≤ 4.5 Å and Tanimoto coefficients > 0.6 after selection based on ligand complementation than target receptor models produced in the absence of template ligands. Overall, the findings produced by this study support the use of local template homology modeling in combination with de novo ECL2 modeling in the presence of a ligand from the template crystal structure to generate GPCR models intended to study ligand binding interactions.

Project description:BACKGROUND: G protein-coupled receptors (GPCRs) transduce a wide variety of extracellular signals to within the cell and therefore have a key role in regulating cell activity and physiological function. GPCR malfunction is responsible for a wide range of diseases including cancer, diabetes and hyperthyroidism and a large proportion of drugs on the market target these receptors. The three dimensional structure of GPCRs is important for elucidating the molecular mechanisms underlying these diseases and for performing structure-based drug design. Although structural data are restricted to only a handful of GPCRs, homology models can be used as a proxy for those receptors not having crystal structures. However, many researchers working on GPCRs are not experienced homology modellers and are therefore unable to benefit from the information that can be gleaned from such three-dimensional models. Here, we present a comprehensive database called the GPCR-SSFE, which provides initial homology models of the transmembrane helices for a large variety of family A GPCRs. DESCRIPTION: Extending on our previous theoretical work, we have developed an automated pipeline for GPCR homology modelling and applied it to a large set of family A GPCR sequences. Our pipeline is a fragment-based approach that exploits available family A crystal structures. The GPCR-SSFE database stores the template predictions, sequence alignments, identified sequence and structure motifs and homology models for 5025 family A GPCRs. Users are able to browse the GPCR dataset according to their pharmacological classification or search for results using a UniProt entry name. It is also possible for a user to submit a GPCR sequence that is not contained in the database for analysis and homology model building. The models can be viewed using a Jmol applet and are also available for download along with the alignments. CONCLUSIONS: The data provided by GPCR-SSFE are useful for investigating general and detailed sequence-structure-function relationships of GPCRs, performing structure-based drug design and for better understanding the molecular mechanisms underlying disease-associated mutations in GPCRs. The effectiveness of our multiple template and fragment approach is demonstrated by the accuracy of our predicted homology models compared to recently published crystal structures.

Project description:The human genome encodes 850 G protein-coupled receptors (GPCRs), half of which are considered potential drug targets. GPCRs transduce extracellular stimuli into a plethora of vital physiological processes. Consequently, GPCRs are an attractive drug target class. This is underlined by the fact that approximately 40% of marketed drugs modulate GPCRs. Intriguingly 60% of non-olfactory GPCRs have no drugs or candidates in clinical development, highlighting the continued potential of GPCRs as drug targets. The discovery of small molecules targeting these GPCRs by conventional high throughput screening (HTS) campaigns is challenging. Although the definition of success varies per company, the success rate of HTS for GPCRs is low compared to other target families (Fujioka and Omori, 2012; Dragovich et al., 2022). Beyond this, GPCR structure determination can be difficult, which often precludes the application of structure-based drug design approaches to arising HTS hits. GPCR structural studies entail the resource-demanding purification of native receptors, which can be challenging as they are inherently unstable when extracted from the lipid matrix. Moreover, GPCRs are flexible molecules that adopt distinct conformations, some of which need to be stabilized if they are to be structurally resolved. The complexity of targeting distinct therapeutically relevant GPCR conformations during the early discovery stages contributes to the high attrition rates for GPCR drug discovery programs. Multiple strategies have been explored in an attempt to stabilize GPCRs in distinct conformations to better understand their pharmacology. This review will focus on the use of camelid-derived immunoglobulin single variable domains (VHHs) that stabilize disease-relevant pharmacological states (termed ConfoBodies by the authors) of GPCRs, as well as GPCR:signal transducer complexes, to accelerate drug discovery. These VHHs are powerful tools for supporting in vitro screening, deconvolution of complex GPCR pharmacology, and structural biology purposes. In order to demonstrate the potential impact of ConfoBodies on translational research, examples are presented of their role in active state screening campaigns and structure-informed rational design to identify de novo chemical space and, subsequently, how such matter can be elaborated into more potent and selective drug candidates with intended pharmacology.

Dataset Information

Homology modeling of a Class A GPCR in the inactive conformation: A quantitative analysis of the correlation between model/template sequence identity and model accuracy.

Publications

Homology modeling of a Class A GPCR in the inactive conformation: A quantitative analysis of the correlation between model/template sequence identity and model accuracy.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets