Modeling activated states of GPCRs: the rhodopsin template

Use of Molecular Dynamics Simulations in Structure-Based Drug Discovery

BACKGROUND

Traditional drug discovery is a lengthy process which involves a huge amount of resources. Modern-day drug discovers various multidisciplinary approaches amongst which, computational ligand and structure-based drug designing methods contribute significantly. Structure-based drug designing techniques require the knowledge of structural information of drug target and drug-target complexes. Proper understanding of drug-target binding requires the flexibility of both ligand and receptor to be incorporated. Molecular docking refers to the static picture of the drug-target complex(es). Molecular dynamics, on the other hand, introduces flexibility to understand the drug binding process.

OBJECTIVE

The aim of the present study is to provide a systematic review on the usage of molecular dynamics simulations to aid the process of structure-based drug design.

METHOD

This review discussed findings from various research articles and review papers on the use of molecular dynamics in drug discovery. All efforts highlight the practical grounds for which molecular dynamics simulations are used in drug designing program. In summary, various aspects of the use of molecular dynamics simulations that underline the basis of studying drug-target complexes were thoroughly explained.

RESULTS

This review is the result of reviewing more than a hundred papers. It summarizes various problems that use molecular dynamics simulations.

CONCLUSION

The findings of this review highlight how molecular dynamics simulations have been successfully implemented to study the structure-function details of specific drug-target complexes. It also identifies the key areas such as stability of drug-target complexes, ligand binding kinetics and identification of allosteric sites which have been elucidated using molecular dynamics simulations.

Characterization of gonadotropin receptors Fshr and Lhr in Japanese medaka, Oryzias latipes

Reproduction in vertebrates is controlled by the brain-pituitary-gonad axis, where the two gonadotropins follicle-stimulating hormone (Fsh) and luteinizing hormone (Lh) play vital parts by activating their cognate receptors in the gonads. The main purpose of this work was to study intra- and interspecies ligand promiscuity of teleost gonadotropin receptors, since teleost receptor specificity is unclear, in contrast to mammalian receptors. Receptor activation was investigated by transfecting COS-7 cells with either Fsh receptor (mdFshr, tiFshr) or Lh receptor (mdLhr, tiLhr), and tested for activation by recombinant homologous and heterologous ligands (mdFshβα, mdLhβα, tiFshβα, tiLhβα) from two representative fish orders, Japanese medaka (Oryzias latipes, Beloniformes) and Nile tilapia (Oreochromis niloticus, Cichliformes). Results showed that each gonadotropin preferentially activates its own cognate receptor. Cross-reactivity was detected to some extent as mdFshβα was able to activate the mdLhr, and mdLhβα the mdFshr. Medaka pituitary extract (MPE) stimulated CRE-LUC activity in COS-7 cells expressing mdlhr, but could not stimulate cells expressing mdfshr. Recombinant tiLhβα, tiFshβα and tilapia pituitary extract (TPE) could activate the mdLhr, suggesting cross-species reactivity for mdLhr. Cross-species reactivity was also detected for mdFshr due to activation by tiFshβα, tiLhβα, and TPE, as well as for tiFshr and tiLhr due to stimulation by mdFshβα, mdLhβα, and MPE. Tissue distribution analysis of gene expression revealed that medaka receptors, fshr and lhr, are highly expressed in both ovary and testis. High expression levels were found for lhr also in brain, while fshr was expressed at low levels. Both fshr and lhr mRNA levels increased significantly during testis development. Amino acid sequence alignment and three-dimensional modelling of ligands and receptors highlighted conserved beta sheet domains of both Fsh and Lh between Japanese medaka and Nile tilapia. It also showed a higher structural homology and similarity of transmembrane regions of Lhr between both species, in contrast to Fshr, possibly related to the substitution of the conserved cysteine residue in the transmembrane domain 6 in medaka Fshr with glycine. Taken together, this is the first characterization of medaka Fshr and Lhr using homologous ligands, enabling to better understand teleost hormone-receptor interactions and specificities. The data suggest partial ligand promiscuity and cross-species reactivity between gonadotropins and their receptors in medaka and tilapia.

Joint-based description of protein structure: its application to the geometric characterization of membrane proteins

A macroscopic description of a protein structure allows an understanding of the protein conformations in a more simplistic manner. Here, a new macroscopic approach that utilizes the joints of the protein secondary structures as a basic descriptor for the protein structure is proposed and applied to study the arrangement of secondary structures in helical membrane proteins. Two types of dihedral angle, Ω and λ, were defined based on the joint points of the transmembrane (TM) helices and loops, and employed to analyze 103 non-homologous membrane proteins with 3 to 14 TM helices. The Ω-λ plot, which is a distribution plot of the dihedral angles of the joint points, identified the allowed and disallowed regions of helical arrangement. Analyses of consecutive dihedral angle patterns indicated that there are preferred patterns in the helical alignment and extension of TM proteins, and helical extension pattern in TM proteins is varied as the size of TM proteins increases. Finally, we could identify some symmetric protein pairs in TM proteins under the joint-based coordinate and 3-dimensional coordinates. The joint-based approach is expected to help better understand and model the overall conformational features of complicated large-scale proteins, such as membrane proteins.

Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study

Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.

Comparing Class A GPCRs to bitter taste receptors: Structural motifs, ligand interactions and agonist-to-antagonist ratios

G protein-coupled receptors (GPCRs) are seven transmembrane (TM) proteins that play a key role in human physiology. The GPCR superfamily comprises about 800 members, classified into several classes, with rhodopsin-like Class A being the largest and most studied thus far. A huge component of the human repertoire consists of the chemosensory GPCRs, including ∼400 odorant receptors, 25 bitter taste receptors (TAS2Rs), which are thought to guard the organism from consuming poisons, and sweet and umami TAS1R heteromers, which indicate the nutritive value of food. The location of the binding site of TAS2Rs is similar to that of Class A GPCRs. However, most of the known bitter ligands are agonists, with only a few antagonists documented thus far. The agonist-to-antagonist ratios of Class A GPCRs vary, but in general are much lower than for TAS2Rs. For a set of well-studied GPCRs, a gradual change in agonists-to-antagonists ratios is observed when comparing low (10 μM)- and high (10 nM)-affinity ligand sets from ChEMBL and the DrugBank set of drugs. This shift reflects pharmaceutical bias toward the therapeutically desirable pharmacology for each of these GPCRs, while the 10 μM sets possibly represent the native tendency of the receptors toward either agonists or antagonists. Analyzing ligand-GPCR interactions in 56 X-ray structures representative of currently available structural data, we find that the N-terminus, TM1 and TM2 are more involved in binding of antagonists than of agonists. On the other hand, ECL2 tends to be more involved in binding of agonists. This is of interest, since TAS2Rs harbor variations on the typical Class A sequence motifs, including the absence of the ECL2-TM3 disulfide bridge. This suggests an alternative mode of regulation of conformational states for TAS2Rs, with potentially less stabilized inactive state. The comparison of TAS2Rs and Class A GPCRs structural features and the pharmacology of the their ligands highlights the intricacies of GPCR architecture and provides a framework for rational design of new ligands.

Milestones in the Activation of a G Protein-Coupled Receptor. Insights from Molecular-Dynamics Simulations into the Human Cholecystokinin Receptor-1

Activation of G protein-coupled receptors (GPCRs) obeys an allosteric mechanism triggered by ligand binding. To understand how the signal is transduced in the cell, identification of the milestones paving the pathway between the active and the inactive states of the receptor is necessary. A model of the human cholecystokinin receptor-1 (CCK1R) has been proposed recently. The complex formed by CCK1R and an agonist ligand will serve as a paradigm of an active conformation to capture milestones in GPCR activation. To reach this goal, assuming microreversibility, the initial step toward the inactivation of CCK1R was modeled using free energy calculations, whereby the ligand is removed from the binding pocket. However accurate the reproduction of the experimental affinity constant, this simulation only represents an embryonic stage of the inactivation process. Starting from the apo receptor, an unprecedented 0.1-μs molecular dynamics trajectory was generated, bereft of experimental biases, bringing into the light key events in the inactivation of CCK1R, chief among which the hydration of its internal cavity, concomitant with the spatial rearrangement of the transmembrane helical segments. Hydration is intimately related to the isomerization of the highly conserved residue W326 of helix VI, acting as a two-state toggle switch, and of residue M121 of helix III. In the active state, the former residue obstructs the crevice, thereby preventing water leakage, which would otherwise trigger the disruption of an ionic lock between helices II and III involving the signature E/DRY motif ubiquitous to GPCRs.

Elastic network normal mode dynamics reveal the GPCR activation mechanism

G-protein-coupled receptors (GPCR) are a family of membrane-embedded metabotropic receptors which translate extracellular ligand binding into an intracellular response. Here, we calculate the motion of several GPCR family members such as the M2 and M3 muscarinic acetylcholine receptors, the A2A adenosine receptor, the β2 -adrenergic receptor, and the CXCR4 chemokine receptor using elastic network normal modes. The normal modes reveal a dilation and a contraction of the GPCR vestibule associated with ligand passage, and activation, respectively. Contraction of the vestibule on the extracellular side is correlated with cavity formation of the G-protein binding pocket on the intracellular side, which initiates intracellular signaling. Interestingly, the normal modes of rhodopsin do not correlate well with the motion of other GPCR family members. Electrostatic potential calculation of the GPCRs reveal a negatively charged field around the ligand binding site acting as a siphon to draw-in positively charged ligands on the membrane surface. Altogether, these results expose the GPCR activation mechanism and show how conformational changes on the cell surface side of the receptor are allosterically translated into structural changes on the inside.

Structural features of the G-protein/GPCR interactions

BACKGROUND

The details of the functional interaction between G proteins and the G protein coupled receptors (GPCRs) have long been subjected to extensive investigations with structural and functional assays and a large number of computational studies.

SCOPE OF REVIEW

The nature and sites of interaction in the G-protein/GPCR complexes, and the specificities of these interactions selecting coupling partners among the large number of families of GPCRs and G protein forms, are still poorly defined.

MAJOR CONCLUSIONS

Many of the contact sites between the two proteins in specific complexes have been identified, but the three dimensional molecular architecture of a receptor-Gα interface is only known for one pair. Consequently, many fundamental questions regarding this macromolecular assembly and its mechanism remain unanswered.

GENERAL SIGNIFICANCE

In the context of current structural data we review the structural details of the interfaces and recognition sites in complexes of sub-family A GPCRs with cognate G-proteins, with special emphasis on the consequences of activation on GPCR structure, the prevalence of preassembled GPCR/G-protein complexes, the key structural determinants for selective coupling and the possible involvement of GPCR oligomerization in this process.

Genome-wide Membrane Protein Structure Prediction

Transmembrane proteins allow cells to extensively communicate with the external world in a very accurate and specific way. They form principal nodes in several signaling pathways and attract large interest in therapeutic intervention, as the majority pharmaceutical compounds target membrane proteins. Thus, according to the current genome annotation methods, a detailed structural/functional characterization at the protein level of each of the elements codified in the genome is also required. The extreme difficulty in obtaining high-resolution three-dimensional structures, calls for computational approaches. Here we review to which extent the efforts made in the last few years, combining the structural characterization of membrane proteins with protein bioinformatics techniques, could help describing membrane proteins at a genome-wide scale. In particular we analyze the use of comparative modeling techniques as a way of overcoming the lack of high-resolution three-dimensional structures in the human membrane proteome.

Homology model-assisted elucidation of binding sites in GPCRs

G protein-coupled receptors (GPCRs) are important mediators of cell signaling and a major family of drug targets. Despite recent breakthroughs, experimental elucidation of GPCR structures remains a formidable challenge. Homology modeling of 3D structures of GPCRs provides a practical tool for elucidating the structural determinants governing the interactions of these important receptors with their ligands. The working model of the binding site can then be used for virtual screening of additional ligands that may fit this site, for determining and comparing specificity profiles of related receptors, and for structure-based design of agonists and antagonists. The current review presents the protocol and enumerates the steps for modeling and validating the residues involved in ligand binding. The main stages include (a) modeling the receptor structure using an automated fragment-based approach, (b) predicting potential binding pockets, (c) docking known binders, (d) analyzing predicted interactions and comparing with positions that have been shown to bind ligands in other receptors, (e) validating the structural model by mutagenesis.

Novel computational methodologies for structural modeling of spacious ligand binding sites of G-protein-coupled receptors: development and application to human leukotriene B4 receptor

This paper describes a novel method to predict the activated structures of G-protein-coupled receptors (GPCRs) with high accuracy, while aiming for the use of the predicted 3D structures in in silico virtual screening in the future. We propose a new method for modeling GPCR thermal fluctuations, where conformation changes of the proteins are modeled by combining fluctuations on multiple time scales. The core idea of the method is that a molecular dynamics simulation is used to calculate average 3D coordinates of all atoms of a GPCR protein against heat fluctuation on the picosecond or nanosecond time scale, and then evolutionary computation including receptor-ligand docking simulations functions to determine the rotation angle of each helix of a GPCR protein as a movement on a longer time scale. The method was validated using human leukotriene B4 receptor BLT1 as a sample GPCR. Our study demonstrated that the proposed method was able to derive the appropriate 3D structure of the active-state GPCR which docks with its agonists.

Structure-based molecular modeling approaches to GPCR oligomerization

Classical structure-based drug design techniques using G-protein-coupled receptors (GPCRs) as targets focus nearly exclusively on binding at the orthosteric site of a single receptor. Dimerization and oligomerization of GPCRs, proposed almost 30 years ago, have, however, crucial relevance for drug design. Targeting these complexes selectively or designing small molecules that affect receptor-receptor interactions might provide new opportunities for novel drug discovery. In order to study the mechanisms and dynamics that rule GPCRs oligomerization, it is essential to understand the dynamic process of receptor-receptor association and to identify regions that are suitable for selective drug binding, which may be determined with experimental methods such as Förster resonance energy transfer (FRET) or Bioluminescence resonance energy transfer (BRET) and computational sequence- and structure-based approaches. The aim of this chapter is to provide a comprehensive description of the structure-based molecular modeling methods for studying GPCR dimerization, that is, protein-protein docking, molecular dynamics, normal mode analysis, and electrostatics studies.

Experimental and computational study of inter- and intra- species specificity of gonadotropins for various gonadotropin receptors

The gonadotropins follicle-stimulating hormone (FSH) and luteinizing hormone (LH) and their receptors play critical roles in vertebrate reproduction. In order to study intra- and interspecies ligand promiscuity of gonadotropins, COS-7 cells were transiently transfected with one of the gonadotropin receptor genes, FSHR or LHR, and tested for activation by gonadotropins from representative fish orders: Aquilliformes (eel; e), Salmoniformes (trout; tr), and Perciformes (tilapia; ta), and of mammalian origin: porcine (p), bovine (b) and human (h). The study reveals complex relations between the gonadotropin hormones and their receptors. Each gonadotropin activated its own cognate receptor. However, taLHR was also activated by hCG and eLHR was activated by hFSH, hCG, and trFSH. For FSHR, the only cross-reactivity detected was for hFSHR, which was activated by pFSH and bFSH. These findings are of great interest and applicability in the context of activation of various GTHRs by their ligands and by ligands from other vertebrates. Analysis of the three-dimensional models of the structures highlights the importance of residues outside of the currently established hormone-receptor interface region. In addition, the interface residues in taFSHR and the effect of exon duplication, which causes an insert in the LRR domain, are suggested to affect the interaction and binding of taFSH.

Action of molecular switches in GPCRs--theoretical and experimental studies

G protein coupled receptors (GPCRs), also called 7TM receptors, form a huge superfamily of membrane proteins that, upon activation by extracellular agonists, pass the signal to the cell interior. Ligands can bind either to extracellular N-terminus and loops (e.g. glutamate receptors) or to the binding site within transmembrane helices (Rhodopsin-like family). They are all activated by agonists although a spontaneous auto-activation of an empty receptor can also be observed. Biochemical and crystallographic methods together with molecular dynamics simulations and other theoretical techniques provided models of the receptor activation based on the action of so-called "molecular switches" buried in the receptor structure. They are changed by agonists but also by inverse agonists evoking an ensemble of activation states leading toward different activation pathways. Switches discovered so far include the ionic lock switch, the 3-7 lock switch, the tyrosine toggle switch linked with the nPxxy motif in TM7, and the transmission switch. The latter one was proposed instead of the tryptophan rotamer toggle switch because no change of the rotamer was observed in structures of activated receptors. The global toggle switch suggested earlier consisting of a vertical rigid motion of TM6, seems also to be implausible based on the recent crystal structures of GPCRs with agonists. Theoretical and experimental methods (crystallography, NMR, specific spectroscopic methods like FRET/BRET but also single-molecule-force-spectroscopy) are currently used to study the effect of ligands on the receptor structure, location of stable structural segments/domains of GPCRs, and to answer the still open question on how ligands are binding: either via ensemble of conformational receptor states or rather via induced fit mechanisms. On the other hand the structural investigations of homoand heterodimers and higher oligomers revealed the mechanism of allosteric signal transmission and receptor activation that could lead to design highly effective and selective allosteric or ago-allosteric drugs.

Modeling of mammalian olfactory receptors and docking of odorants

Olfactory receptors (ORs) belong to the superfamily of G protein-coupled receptors (GPCRs), the second largest class of genes after those related to immunity, and account for about 3 % of mammalian genomes. ORs are present in all multicellular organisms and represent more than half the GPCRs in mammalian species (e.g., the mouse OR repertoire contains >1,000 functional genes). ORs are mainly expressed in the olfactory epithelium where they detect odorant molecules, but they are also expressed in a number of other cells, such as sperm cells, although their functions in these cells remain mostly unknown. It has recently been reported that ORs are present in tumoral tissues where they are expressed at different levels than in healthy tissues. A specific OR is over-expressed in prostate cancer cells, and activation of this OR has been shown to inhibit the proliferation of these cells. Odorant stimulation of some of these receptors results in inhibition of cell proliferation. Even though their biological role has not yet been elucidated, these receptors might constitute new targets for diagnosis and therapeutics. It is important to understand the activation mechanism of these receptors at the molecular level, in particular to be able to predict which ligands are likely to activate a particular receptor ('deorphanization') or to design antagonists for a given receptor. In this review, we describe the in silico methodologies used to model the three-dimensional (3D) structure of ORs (in the more general framework of GPCR modeling) and to dock ligands into these 3D structures.

Computational and experimental analysis of the transmembrane domain 4/5 dimerization interface of the serotonin 5-HT(1A) receptor

Experimental evidence suggests that most members of class A G-protein coupled receptors (GPCRs) can form homomers and heteromers in addition to functioning as single monomers. In particular, serotonin (5-HT) receptors were shown to homodimerize and heterodimerize with other GPCRs, although the details and the physiological role of the oligomerization has not yet been fully elucidated. Here we used computational modeling of the 5-HT(1A) receptor monomer and dimer to predict residues important for dimerization. Based on these results, we carried out rationally designed site-directed mutagenesis. The ability of the mutants to dimerize was evaluated using different FRET-based approaches. The reduced levels of acceptor photobleaching-Förster resonance energy transfer (FRET) and the lower number of monomers participating in oligomers, as assessed by lux-FRET, confirmed the decreased ability of the mutants to dimerize and the involvement of the predicted contacts (Trp175(4.64), Tyr198(5.41), Arg151(4.40), and Arg152(4.41)) at the interface. This information was reintroduced as constraints for computational protein-protein docking to obtain a high-quality dimer model. Analysis of the refined model as well as molecular dynamics simulations of wild-type (WT) and mutant dimers revealed compensating interactions in dimers composed of WT and W175A mutant. This provides an explanation for the requirement of mutations of Trp175(4.64) in both homomers for disrupting dimerization. Our iterative computational-experimental study demonstrates that transmembrane domains TM4/TM5 can form an interaction interface in 5-HT(1A) receptor dimers and indicates that specific amino acid interactions maintain this interface. The mutants and the optimized model of the dimer structure may be used in functional studies of serotonin dimers.

FINDSITE(X): a structure-based, small molecule virtual screening approach with application to all identified human GPCRs

We have developed FINDSITE(X), an extension of FINDSITE, a protein threading based algorithm for the inference of protein binding sites, biochemical function and virtual ligand screening, that removes the limitation that holo protein structures (those containing bound ligands) of a sufficiently large set of distant evolutionarily related proteins to the target be solved; rather, predicted protein structures and experimental ligand binding information are employed. To provide the predicted protein structures, a fast and accurate version of our recently developed TASSER(VMT), TASSER(VMT)-lite, for template-based protein structural modeling applicable up to 1000 residues is developed and tested, with comparable performance to the top CASP9 servers. Then, a hybrid approach that combines structure alignments with an evolutionary similarity score for identifying functional relationships between target and proteins with binding data has been developed. By way of illustration, FINDSITE(X) is applied to 998 identified human G-protein coupled receptors (GPCRs). First, TASSER(VMT)-lite provides updates of all human GPCR structures previously modeled in our lab. We then use these structures and the new function similarity detection algorithm to screen all human GPCRs against the ZINC8 nonredundant (TC < 0.7) ligand set combined with ligands from the GLIDA database (a total of 88,949 compounds). Testing (excluding GPCRs whose sequence identity > 30% to the target from the binding data library) on a 168 human GPCR set with known binding data, the average enrichment factor in the top 1% of the compound library (EF(0.01)) is 22.7, whereas EF(0.01) by FINDSITE is 7.1. For virtual screening when just the target and its native ligands are excluded, the average EF(0.01) reaches 41.4. We also analyze off-target interactions for the 168 protein test set. All predicted structures, virtual screening data and off-target interactions for the 998 human GPCRs are available at http://cssb.biology.gatech.edu/skolnick/webservice/gpcr/index.html .

Showcasing modern molecular dynamics simulations of membrane proteins through G protein-coupled receptors

Despite many years of dedicated efforts, high-resolution structural determination of membrane proteins lags far behind that of soluble proteins. Computational methods in general, and molecular dynamics (MD) simulations in particular, have represented important alternative resources over the years to advance understanding of membrane protein structure and function. However, it is only recently that much progress has been achieved owing to new high-resolution membrane protein structures, specialized parallel computer architectures, and efficient simulation algorithms. This has definitely been the case for G protein-coupled receptors (GPCRs), which have assumed a leading role in the area of structural biology with several new structures appearing in the literature during the past five years. We provide here a concise overview of recent developments in computational biophysics of membrane proteins, using GPCRs as an example to showcase important information that can be derived from modern MD simulations.

In silico identification of new ligands for GPR17: a promising therapeutic target for neurodegenerative diseases

GPR17, a previously orphan receptor responding to both uracil nucleotides and cysteinyl-leukotrienes, has been proposed as a novel promising target for human neurodegenerative diseases. Here, in order to specifically identify novel potent ligands of GPR17, we first modeled in silico the receptor by using a multiple template approach, in which extracellular loops of the receptor, quite complex to treat, were modeled making reference to the most similar parts of all the class-A GPCRs crystallized so far. A high-throughput virtual screening exploration of GPR17 binding site with more than 130,000 lead-like compounds was then applied, followed by the wet functional and pharmacological validation of the top-scoring chemical structures. This approach revealed successful for the proposed aim, and allowed us to identify five agonists or partial agonists with very diverse chemical structure. None of these compounds could have been expected 'a priori' to act on a GPCR, and all of them behaved as much more potent ligands than GPR17 endogenous activators.

Modeling GPCR active state conformations: the β(2)-adrenergic receptor

The recent publication of several G protein-coupled receptor (GPCR) structures has increased the information available for homology modeling inactive class A GPCRs. Moreover, the opsin crystal structure shows some active features. We have therefore combined information from these two sources to generate an extensively validated model of the active conformation of the β(2)-adrenergic receptor. Experimental information on fully active GPCRs from zinc binding studies, site-directed spin labeling, and other spectroscopic techniques has been used in molecular dynamics simulations. The observed conformational changes reside mainly in transmembrane helix 6 (TM6), with additional small but significant changes in TM5 and TM7. The active model has been validated by manual docking and is in agreement with a large amount of experimental work, including site-directed mutagenesis information. Virtual screening experiments show that the models are selective for β-adrenergic agonists over other GPCR ligands, for (R)- over (S)-β-hydroxy agonists and for β(2)-selective agonists over β(1)-selective agonists. The virtual screens reproduce interactions similar to those generated by manual docking. The C-terminal peptide from a model of the stimulatory G protein, readily docks into the active model in a similar manner to which the C-terminal peptide from transducin, docks into opsin, as shown in a recent opsin crystal structure. This GPCR-G protein model has been used to explain site-directed mutagenesis data on activation. The agreement with experiment suggests a robust model of an active state of the β(2)-adrenergic receptor has been produced. The methodology used here should be transferable to modeling the active state of other GPCRs.

Progress in elucidating the structural and dynamic character of G Protein-Coupled Receptor oligomers for use in drug discovery

G Protein-Coupled Receptors (GPCRs) are the most targeted group of proteins for the development of therapeutic drugs. Until the last decade, structural information about this family of membrane proteins was relatively scarce, and their mechanisms of ligand binding and signal transduction were modeled on the assumption that GPCRs existed and functioned as monomeric entities. New crystal structures of native and engineered GPCRs, together with important biochemical and biophysical data that reveal structural details of the activation mechanism(s) of this receptor family, provide a valuable framework to improve dynamic molecular models of GPCRs with the ultimate goal of elucidating their allostery and functional selectivity. Since the dynamic movements of single GPCR protomers are likely to be affected by the presence of neighboring interacting subunits, oligomeric arrangements should be taken into account to improve the predictive ability of computer-assisted structural models of GPCRs for effective use in drug design.

2,3-Dihydro-1-benzofuran derivatives as a series of potent selective cannabinoid receptor 2 agonists: design, synthesis, and binding mode prediction through ligand-steered modeling

We recently discovered and reported a series of N-alkyl-isatin acylhydrazone derivatives that are potent cannabinoid receptor 2 (CB(2)) agonists. In an effort to improve the druglike properties of these compounds and to better understand and improve the treatment of neuropathic pain, we designed and synthesized a new series of 2,3-dihydro-1-benzofuran derivatives bearing an asymmetric carbon atom that behave as potent selective CB(2) agonists. We used a multidisciplinary medicinal chemistry approach with binding mode prediction through ligand-steered modeling. Enantiomer separation and configuration assignment were carried out for the racemic mixture for the most selective compound, MDA7 (compound 18). It appeared that the S enantiomer, compound MDA104 (compound 33), was the active enantiomer. Compounds MDA42 (compound 19) and MDA39 (compound 30) were the most potent at CB(2). MDA42 was tested in a model of neuropathic pain and exhibited activity in the same range as that of MDA7. Preliminary ADMET studies for MDA7 were performed and did not reveal any problems.

BlockMaster: partitioning protein kinase structures using normal-mode analysis

Protein kinases are key signaling enzymes which are dysregulated in many health disorders and therefore represent major targets of extensive drug discovery efforts. Their regulation in the cell is exerted via various mechanisms, including control of the 3D conformation of their catalytic domains. We developed a procedure, BlockMaster, for partitioning protein structures into semirigid blocks and flexible regions based on residue-residue correlations calculated from normal modes. BlockMaster provided correct partitioning into domains and subdomains of several test set proteins for which documented expert annotation of subdomains exists. When applied to representative structures of protein kinases, BlockMaster identified semirigid blocks within the traditional N-terminal and C-terminal lobes of the kinase domain. In general, the block regions had elevated helical content and reduced, but significant, coil content compared to the nonblock (flexible) regions. The specificity-determining regions, previously used to derive inhibitory peptides, were found to be more flexible in the tyrosine kinases than in serine/threonine kinases. Two blocks were identified which spanned both lobes. The first, which we termed the "pivot" block, included the alphaC-beta4 loop in the N-terminal lobe and part of the activation loop in the C-terminal lobe and appeared in both the active and inactive conformations of the kinases. The second, which we termed the "loop" block, differed between the active and inactive conformations. In the structures of active kinases, this block included part of the activation loop in the C-terminal lobe and the alphaC helix in the N-terminal lobe, representing a known interaction that stabilizes the active conformation. In the inactive structures, this block included G loop residues instead of the alphaC residues. This novel inactive "loop" block may stabilize the inactive conformation and thus downregulate kinase activity.

Activation of the ghrelin receptor is described by a privileged collective motion: a model for constitutive and agonist-induced activation of a sub-class A G-protein coupled receptor (GPCR)

Three homology models of the human ghrelin receptor (GHS-R1a) have been generated from the available X-ray structures of rhodopsin (RHO model), opsin (OPS model) and beta-2 adrenergic receptor (B2 model). The latter was used as a starting point for combined molecular dynamics simulation (MDS) and full atom normal modes analysis (NMA). A low-frequency normal mode (mode 16) perfectly reproduced the intracellular motions observed between B2 and RHO models; in the opposite direction along the same mode, the generated structures are closer to the OPS model, suggesting a direct link with GHS-R1a activation. This was in agreement with motions of the seven transmembranous segments, increase of the solvent accessibility of the 140-ERY-142 sequence, and flip of the Trp276 (C WLP) residue, some features related to GPCRs activation. According to our model, His280 was proposed to stabilize Trp276 in the active state; this was verified by site-directed mutagenesis and biochemical characterization of the resulting H280A and H280S mutants, which were fully functional but sharing an important decrease of their basal activities. Docking performed with short ghrelin derivatives Gly-Ser-Ser ([octa])-Phe-NH (2) and Gly-Ser-Ser ([octa])-Phe-Leu-NH (2) allowed the identification of a robust position of these peptides in the active site of the receptor. This model was refined by MDS and validated by docking experiments performed on a set of 55 ghrelin receptor ligands based on the 1,2,4- triazole scaffold. Finally, NMA performed on the obtained peptide-receptor complex suggested stabilization of the Trp276 residue and of the whole receptor in the active state, preventing the motion observed along mode 16 computed for the unbound receptor. Our results show that NMA offers a powerful approach to study the conformational diversity and the activation mechanism of GPCRs.

Computational analysis of membrane proteins: the largest class of drug targets

Given the key roles of integral membrane proteins as transporters and channels, it is necessary to understand their structures and, hence, mechanisms and regulation at the molecular level. Membrane proteins represent approximately 30% of all proteins of currently sequenced genomes. Paradoxically, however, only approximately 2% of crystal structures deposited in the protein data bank are of membrane proteins, and very few of these are at high resolution (better than 2A). The great disparity between our understanding of soluble proteins and our understanding of membrane proteins is because of the practical problems of working with membrane proteins - specifically, difficulties in expression, purification and crystallization. Thus, computational modeling has been utilized extensively to make crucial advances in understanding membrane protein structure and function.

Interactions of histamine H1-receptor agonists and antagonists with the human histamine H4-receptor

The human histamine H(4)-receptor (hH(4)R) possesses high constitutive activity and, like the human H(1)-receptor (hH(1)R), is involved in the pathogenesis of type-I allergic reactions. The study aims were to explore the value of dual H(1)/H(4)R antagonists as antiallergy drugs and to address the question of whether H(1)R ligands bind to hH(4)R. In an acute murine asthma model, the H(1)R antagonist mepyramine and the H(4)R antagonist 1-[(5-chloro-1H-indol-2-yl)carbonyl]-4-methyl-piperazine (JNJ 7777120) exhibited synergistic inhibitory effects on eosinophil accumulation in the bronchoalveolar lavage fluid. At the hH(4)R expressed in Sf9 insect cells, 18 H(1)R antagonists and 22 H(1)R agonists showed lower affinity to hH(4)R than to hH(1)R as assessed in competition binding experiments. For a small number of H(1)R antagonists, hH(4)R partial agonism was observed in the steady-state GTPase assay. Most compounds were neutral antagonists or inverse agonists. Twelve phenylhistamine-type hH(1)R partial agonists were also hH(4)R partial agonists. Four histaprodifen-type hH(1)R partial agonists were hH(4)R inverse agonists. Dimeric histaprodifen was a more efficacious hH(4)R inverse agonist than the reference compound thioperamide. Suprahistaprodifen was the only histaprodifen acting as hH(4)R partial agonist. Suprahistaprodifen was docked into the binding pocket of inactive and active hH(4)R models in two different orientations, predominantly stabilizing the active state of hH(4)R. Collectively, the synergistic effects of H(1)R and H(4)R antagonists in an acute asthma model and the overlapping interaction of structurally diverse H(1)R ligands with hH(1)R and hH(4)R indicate that the development of dual H(1)R/H(4)R antagonists is a worthwhile and technically feasible goal for the treatment of type-I allergic reactions.

Allosteric communication between protomers of dopamine class A GPCR dimers modulates activation

A major obstacle to understanding the functional importance of dimerization between Class A G protein-coupled receptors (GPCRs) has been the methodological limitation in achieving control of the identity of the components comprising the signaling unit. We have developed a functional complementation assay that enables such control and illustrate it for the human dopamine D2 receptor. The minimal signaling unit, two receptors and a single G protein, is maximally activated by agonist binding to a single protomer, which suggests an asymmetrical activated dimer. Inverse agonist binding to the second protomer enhances signaling, whereas agonist binding to the second protomer blunts signaling. Ligand-independent constitutive activation of the second protomer also inhibits signaling. Thus, GPCR dimer function can be modulated by the activity state of the second protomer, which for a heterodimer may be altered in pathological states. Our novel methodology also makes possible the characterization of signaling from a defined heterodimer unit.

Rational redesign of neutral endopeptidase binding to merlin and moesin proteins

Neutral endopeptidase (NEP) is a 90- to 110-kDa cell-surface peptidase that is normally expressed by numerous tissues but whose expression is lost or reduced in a variety of malignancies. The anti-tumorigenic function of NEP is mediated not only by its catalytic activity but also through direct protein-protein interactions of its cytosolic region with several binding partners, including Lyn kinase, PTEN, and ezrin/radixin/moesin (ERM) proteins. We have previously shown that mutation of the K(19)K(20)K(21) basic cluster in NEPs' cytosolic region to residues QNI disrupts binding to the ERM proteins. Here we show that the ERM-related protein merlin (NF2) does not bind NEP or its cytosolic region. Using experimental data, threading, and sequence analysis, we predicted the involvement of moesin residues E(159)Q(160) in binding to the NEP cytosolic domain. Mutation of these residues to NL (to mimic the corresponding N(159)L(160) residues in the nonbinder merlin) disrupted moesin binding to NEP. Mutation of residues N(159)L(160)Y(161)K(162)M(163) in merlin to the corresponding moesin residues resulted in NEP binding to merlin. This engineered NEP peptide-merlin interaction was diminished by the QNI mutation in NEP, supporting the role of the NEP basic cluster in binding. We thus identified the region of interaction between NEP and moesin, and engineered merlin into a NEP-binding protein. These data form the basis for further exploration of the details of NEP-ERM binding and function.

Modeling small molecule-compound binding to G-protein-coupled receptors

G-protein-coupled receptors (GPCRs) form a superfamily of membrane proteins that play a crucial role in mediating physiological processes as well as pathogenesis of many critical diseases. They are one of the most successful drug targets, accounting for more than 30% of prescription drugs on the market today. Three-dimensional structural information on GPCRs will greatly aid the drug design process, and great strides are being made in obtaining crystallographic information on GPCRs. Since this process is both tedious and risky, a combination of computational methods and biophysical experiments is a useful approach to rapidly obtain information on a wide variety of GPCRs. In this review, we describe the methods/protocols involved in these computational techniques, as well as methods for site-directed mutagenesis and ligand-binding assays that are currently being used for validating structural-model and small-molecule-ligand binding to GPCRs. We discuss the merits and pitfalls of the various methods used in obtaining structural and dynamic information for ligand binding to GPCRs. Another important factor to consider in drug design is the conformational flexibility of GPCRs since it has been shown that small-molecule ligands of varied efficacy stabilize different receptor conformations leading to functional selectivity of ligands. We discuss the computational methods used to study this specific ligand-induced state.

Potential modes of interaction of 9-aminomethyl-9,10-dihydroanthracene (AMDA) derivatives with the 5-HT2A receptor: a ligand structure-affinity relationship, receptor mutagenesis and receptor modeling investigation

The effects of 3-position substitution of 9-aminomethyl-9,10-dihydroanthracene (AMDA) on 5-HT 2A receptor affinity were determined and compared to a parallel series of DOB-like 1-(2,5-dimethoxyphenyl)-2-aminopropanes substituted at the 4-position. The results were interpreted within the context of 5-HT 2A receptor models that suggest that members of the DOB-like series can bind to the receptor in two distinct modes that correlate with the compounds' functional activity. Automated ligand docking and molecular dynamics suggest that all of the AMDA derivatives, the parent of which is a 5-HT 2A antagonist, bind in a fashion analogous to that for the sterically demanding antagonist DOB-like compounds. The failure of the F340 (6.52)L mutation to adversely affect the affinity of AMDA and the 3-bromo derivative is consistent with the proposed modes of orientation. Evaluation of ligand-receptor complex models suggest that a valine/threonine exchange between the 5-HT 2A and D 2 receptors may be the origin of selectivity for AMDA and two substituted derivatives.

Towards a quantitative representation of the cell signaling mechanisms of hallucinogens: measurement and mathematical modeling of 5-HT1A and 5-HT2A receptor-mediated ERK1/2 activation

Through a multidisciplinary approach involving experimental and computational studies, we address quantitative aspects of signaling mechanisms triggered in the cell by the receptor targets of hallucinogenic drugs, the serotonin 5-HT2A receptors. To reveal the properties of the signaling pathways, and the way in which responses elicited through these receptors alone and in combination with other serotonin receptors' subtypes (the 5-HT1AR), we developed a detailed mathematical model of receptor-mediated ERK1/2 activation in cells expressing the 5-HT1A and 5-HT2A subtypes individually, and together. In parallel, we measured experimentally the activation of ERK1/2 by the action of selective agonists on these receptors expressed in HEK293 cells. We show here that the 5-HT1AR agonist Xaliproden HCl elicited transient activation of ERK1/2 by phosphorylation, whereas 5-HT2AR activation by TCB-2 led to higher, and more sustained responses. The 5-HT2AR response dominated the MAPK signaling pathway when co-expressed with 5-HT1AR, and diminution of the response by the 5-HT2AR antagonist Ketanserin could not be rescued by the 5-HT1AR agonist. Computational simulations produced qualitative results in good agreement with these experimental data, and parameter optimization made this agreement quantitative. In silico simulation experiments suggest that the deletion of the positive regulators PKC in the 5-HT2AR pathway, or PLA2 in the combined 5-HT1A/2AR model greatly decreased the basal level of active ERK1/2. Deletion of negative regulators of MKP and PP2A in 5-HT1AR and 5-HT2AR models was found to have even stronger effects. Under various parameter sets, simulation results implied that the extent of constitutive activity in a particular tissue and the specific drug efficacy properties may determine the distinct dynamics of the 5-HT receptor-mediated ERK1/2 activation pathways. Thus, the mathematical models are useful exploratory tools in the ongoing efforts to establish a mechanistic understanding and an experimentally testable representation of hallucinogen-specific signaling in the cellular machinery, and can be refined with quantitative, function-related information.

Atomistic insights into rhodopsin activation from a dynamic model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.

Influence of oligomerization on the dynamics of G-protein coupled receptors as assessed by normal mode analysis

The recently discovered impact of oligomerization on G-protein coupled receptor (GPCR) function further complicates the already challenging goal of unraveling the molecular and dynamic mechanisms of these receptors. To help understand the effect of oligomerization on the dynamics of GPCRs, we have compared the motion of monomeric, dimeric, and tetrameric arrangements of the prototypic GPCR rhodopsin, using an approximate-yet powerful-normal mode analysis (NMA) technique termed elastic network model (ENM). Moreover, we have used ENM to discriminate between putative dynamic mechanisms likely to account for the recently observed conformational rearrangement of the TM4,5-TM4,5 dimerization interface of GPCRs that occurs upon activation. Our results indicate: (1) significant perturbation of the normal modes (NMs) of the rhodopsin monomer upon oligomerization, which is mainly manifested at interfacial regions; (2) increased positive correlation among the transmembrane domains (TMs) and between the extracellular loop (EL) and TM regions of the rhodopsin protomer; (3) highest interresidue positive correlation at the interfaces between protomers; and (4) experimentally testable hypotheses of differential motional changes within different putative oligomeric arrangements.

Modelling the structures of G protein-coupled receptors aided by three-dimensional validation

Background

G protein-coupled receptors (GPCRs) are abundant, activate complex signalling and represent the targets for up to ~60% of pharmaceuticals but there is a paucity of structural data. Bovine rhodopsin is the first GPCR for which high-resolution structures have been completed but significant variations in structure are likely to exist among the GPCRs. Because of this, considerable effort has been expended on developing in silico tools for refining structures of individual GPCRs. We have developed REPIMPS, a modification of the inverse-folding software Profiles-3D, to assess and predict the rotational orientation and vertical position of helices within the helix bundle of individual GPCRs. We highlight the value of the method by applying it to the Baldwin GPCR template but the method can, in principle, be applied to any low- or high-resolution membrane protein template or structure.

Results

3D models were built for transmembrane helical segments of 493 GPCRs based on the Baldwin template, and the models were then scored using REPIMPS and Profiles-3D. The compatibility scores increased significantly using REPIMPS because it takes into account the physicochemical properties of the (lipid) environment surrounding the helix bundle. The arrangement of helices in the helix bundle of the 493 models was then altered systematically by rotating the individual helices. For most GPCRs in the set, changes in the rotational position of one or more helices resulted in significant improvement in the compatibility scores. In particular, for most GPCRs, a rotation of helix VII by 240–300° resulted in improved scores. Bovine rhodopsin modelled using this method showed 3.31 Å RMSD to its crystal structure for 198 C^α atom pairs, suggesting the utility of the method even when starting with idealised structures such as the Baldwin template.

Conclusion

We have developed an in silico tool which can be used to test the validity of, and refine, models of GPCRs with respect to helix rotation and vertical position based on the physicochemical properties of amino acids and the surrounding environment. The method can be applied to any multi-pass membrane protein and potentially can be used in combination with other high-throughput methodologies to generate and refine models of membrane proteins.

Mechanism of signal propagation upon retinal isomerization: insights from molecular dynamics simulations of rhodopsin restrained by normal modes

As one of the best studied members of the pharmaceutically relevant family of G-protein-coupled receptors, rhodopsin serves as a prototype for understanding the mechanism of G-protein-coupled receptor activation. Here, we aim at exploring functionally relevant conformational changes and signal transmission mechanisms involved in its photoactivation brought about through a cis-trans photoisomerization of retinal. For this exploration, we propose a molecular dynamics simulation protocol that utilizes normal modes derived from the anisotropic network model for proteins. Deformations along multiple low-frequency modes of motion are used to efficiently sample collective conformational changes in the presence of explicit membrane and water environment, consistent with interresidue interactions. We identify two highly stable regions in rhodopsin, one clustered near the chromophore, the other near the cytoplasmic ends of transmembrane helices H1, H2, and H7. Due to redistribution of interactions in the neighborhood of retinal upon stabilization of the trans form, local structural rearrangements in the adjoining H3-H6 residues are efficiently propagated to the cytoplasmic end of these particular helices. In the structures obtained by our simulations, all-trans retinal interacts with Cys(167) on H4 and Phe(203) on H5, which were not accessible in the dark state, and exhibits stronger interactions with H5, while some of the contacts made (in the cis form) with H6 are lost.

Advances in the Development and Application of Computational Methodologies for Structural Modeling of G-Protein Coupled Receptors

BACKGROUND: Despite the large amount of experimental data accumulated in the past decade on G-protein coupled receptor (GPCR) structure and function, understanding of the molecular mechanisms underlying GPCR signaling is still far from being complete, thus impairing the design of effective and selective pharmaceuticals. OBJECTIVE: Understanding of GPCR function has been challenged even further by more recent experimental evidence that several of these receptors are organized in the cell membrane as homo- or hetero-oligomers, and that they may exhibit unique pharmacological properties. Given the complexity of these new signaling systems, researcher's efforts are turning increasingly to molecular modeling, bioinformatics and computational simulations for mechanistic insights of GPCR functional plasticity. METHODS: We review here current advances in the development and application of computational approaches to improve prediction of GPCR structure and dynamics, thus enhancing current understanding of GPCR signaling. RESULTS/CONCLUSIONS: Models resulting from use of these computational approaches further supported by experiments are expected to help elucidate the complex allosterism that propagates through GPCR complexes, ultimately aiming at successful structure-based rational drug design.

Terahertz spectroscopy of bacteriorhodopsin and rhodopsin: similarities and differences

We studied the low-frequency terahertz spectroscopy of two photoactive protein systems, rhodopsin and bacteriorhodopsin, as a means to characterize collective low-frequency motions in helical transmembrane proteins. From this work, we found that the nature of the vibrational motions activated by terahertz radiation is surprisingly similar between these two structurally similar proteins. Specifically, at the lowest frequencies probed, the cytoplasmic loop regions of the proteins are highly active; and at the higher terahertz frequencies studied, the extracellular loop regions of the protein systems become vibrationally activated. In the case of bacteriorhodopsin, the calculated terahertz spectra are compared with the experimental terahertz signature. This work illustrates the importance of terahertz spectroscopy to identify vibrational degrees of freedom which correlate to known conformational changes in these proteins.