Ligand-dependent conformations and dynamics of the serotonin 5-HT(2A) receptor determine its activation and membrane-driven oligomerization properties

G-protein coupled receptors: advances in simulation and drug discovery

G-protein coupled receptors (GPCRs), the largest family of human membrane proteins, mediate cellular signaling and represent primary targets of about one third of currently marketed drugs. GPCRs undergo highly dynamic structural transitions during signal transduction, from binding of extracellular ligands to coupling with intracellular effector proteins. Molecular dynamics (MD) simulations have been utilized to investigate GPCR signaling mechanisms (such as pathways of ligand binding and receptor activation/deactivation) and to design novel small-molecule drug candidates. Future research directions point towards modeling cooperative binding of multiple orthosteric and allosteric ligands to GPCRs, GPCR oligomerization and interactions of GPCRs with different intracellular signaling proteins. Through methodological and supercomputing advances, MD simulations will continue to provide important insights into GPCR signaling mechanisms and further facilitate structure-based drug design.

Lipid interaction sites on channels, transporters and receptors: Recent insights from molecular dynamics simulations

Lipid molecules are able to selectively interact with specific sites on integral membrane proteins, and modulate their structure and function. Identification and characterization of these sites are of importance for our understanding of the molecular basis of membrane protein function and stability, and may facilitate the design of lipid-like drug molecules. Molecular dynamics simulations provide a powerful tool for the identification of these sites, complementing advances in membrane protein structural biology and biophysics. We describe recent notable biomolecular simulation studies which have identified lipid interaction sites on a range of different membrane proteins. The sites identified in these simulation studies agree well with those identified by complementary experimental techniques. This demonstrates the power of the molecular dynamics approach in the prediction and characterization of lipid interaction sites on integral membrane proteins. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

The membrane as the gatekeeper of infection: Cholesterol in host-pathogen interaction

The cellular plasma membrane serves as a portal for the entry of intracellular pathogens. An essential step for an intracellular pathogen to gain entry into a host cell therefore is to be able to cross the cell membrane. In this review, we highlight the role of host membrane cholesterol in regulating the entry of intracellular pathogens using insights obtained from work on the interaction of Leishmania and Mycobacterium with host cells. The entry of these pathogens is known to be dependent on host membrane cholesterol. Importantly, pathogen entry is inhibited either upon depletion (or complexation), or enrichment of membrane cholesterol. In other words, an optimum level of host membrane cholesterol is necessary for efficient infection by pathogens. In this overall context, we propose a general mechanism, based on cholesterol-induced conformational changes, involving cholesterol binding sites in host cell surface receptors that are implicated in this process. A therapeutic strategy targeting modulation of membrane cholesterol would have the advantage of avoiding the commonly encountered problem of drug resistance in tackling infection by intracellular pathogens. Insights into the role of host membrane cholesterol in pathogen entry would be instrumental in the development of novel therapeutic strategies to effectively tackle intracellular pathogenesis.

A GIPC1-Palmitate Switch Modulates Dopamine Drd3 Receptor Trafficking and Signaling

Palmitoylation is involved in several neuropsychiatric and movement disorders for which a dysfunctional signaling of the dopamine D3 receptor (Drd3) is hypothesized. Computational modeling of Drd3's homologue, Drd2, has shed some light on the putative role of palmitoylation as a reversible switch for dopaminergic receptor signaling. Drd3 is presumed to be palmitoylated, based on sequence homology with Drd2, but the functional attributes afforded by Drd3 palmitoylation have not been studied. Since these receptors are major targets of antipsychotic and anti-Parkinsonian drugs, a better characterization of Drd3 signaling and posttranslational modifications, like palmitoylation, may improve the prospects for drug development. Using molecular dynamics simulations, we evaluated in silico how Drd3 palmitoylation could elicit significant remodeling of the C-terminal cytoplasmic domain to expose docking sites for signaling proteins. We tested this model in cellulo by using the interaction of Drd3 with the G-alpha interacting protein (GAIP) C terminus 1 (GIPC1) as a template. From a series of biochemical studies, live imaging, and analyses of mutant proteins, we propose that Drd3 palmitoylation acts as a molecular switch for Drd3-biased signaling via a GIPC1-dependent route, which is likely to affect the mode of action of antipsychotic drugs.

Comparative MD Simulations Indicate a Dual Role for Arg1323.50 in Dopamine-Dependent D2R Activation

Residue Arg^3.50 belongs to the highly conserved DRY-motif of class A GPCRs, which is located at the bottom of TM3. On the one hand, Arg^3.50 has been reported to help stabilize the inactive state of GPCRs, but on the other hand has also been shown to be crucial for stabilizing active receptor conformations and mediating receptor-G protein coupling. The combined results of these studies suggest that the exact function of Arg^3.50 is likely to be receptor-dependent and must be characterized independently for every GPCR. Consequently, we now present comparative molecular-dynamics simulations that use our recently described inactive-state and Gα-bound active-state homology models of the dopamine D2 receptor (D2R), which are either bound to dopamine or ligand-free, performed to identify the function of Arg132^3.50 in D2R. Our results are consistent with a dynamic model of D2R activation in which Arg132^3.50 adopts a dual role, both by stabilizing the inactive-state receptor conformation and enhancing dopamine-dependent D2R-G protein coupling.

Progress toward rationally designed small-molecule peptide and peptidomimetic CXCR4 antagonists

Over the last 5 years, X-ray structures of CXCR4 in complex with three different ligands (the small-molecule antagonist IT1t, the polypeptide antagonist CVX15 and the viral chemokine antagonist vMIP-II) have been released. In addition to the inherent scientific value of these specific X-ray structures, they provide a reliable structural foundation for studies of the molecular interactions between CXCR4 and its key peptide ligands (CXCL12 and HIV-1 gp120), and serve as valuable templates for further development of small-molecule CXCR4 antagonists with therapeutic potential. We here review recent computational studies of the molecular interactions between CXCR4 and its peptide ligands – based on the X-ray structures of CXCR4 – and the current status of small-molecule peptide and peptidomimetic CXCR4 antagonists.

AIM for Allostery: Using the Ising Model to Understand Information Processing and Transmission in Allosteric Biomolecular Systems

In performing their biological functions, molecular machines must process and transmit information with high fidelity. Information transmission requires dynamic coupling between the conformations of discrete structural components within the protein positioned far from one another on the molecular scale. This type of biomolecular “action at a distance” is termed allostery. Although allostery is ubiquitous in biological regulation and signal transduction, its treatment in theoretical models has mostly eschewed quantitative descriptions involving the system's underlying structural components and their interactions. Here, we show how Ising models can be used to formulate an approach to allostery in a structural context of interactions between the constitutive components by building simple allosteric constructs we termed Allosteric Ising Models (AIMs). We introduce the use of AIMs in analytical and numerical calculations that relate thermodynamic descriptions of allostery to the structural context, and then show that many fundamental properties of allostery, such as the multiplicative property of parallel allosteric channels, are revealed from the analysis of such models. The power of exploring mechanistic structural models of allosteric function in more complex systems by using AIMs is demonstrated by building a model of allosteric signaling for an experimentally well-characterized asymmetric homodimer of the dopamine D2 receptor.

Structural conservation among the rhodopsin-like and other G protein-coupled receptors

Intramolecular remote coupling within the polypeptide backbones of membrane proteins is difficult to analyze owing to the limited structural information available at the atomic level. Nonetheless, recent progress in the crystallographic study of G protein-coupled receptors (GPCRs) has provided an unprecedented opportunity for understanding the sophisticated architecture of heptahelical transmembrane (7TM) bundles. These 7TM bundles can respond to a wide range of extracellular stimuli while retaining the common function of binding trimeric G proteins. Here we have systematically analyzed select sets of inactive-like 7TM bundles to highlight the structural conservation of the receptors, in terms of intramolecular C_α-C_α distances. Distances with the highest scores were found to be dominated by the intrahelical distances of helix III, regardless of the choice of bundles in the set, indicating that the intracellular half of this helix is highly conserved. Unexpectedly, the distances between the cytoplasmic side of helix I and the extracellular region of helix VI provided the largest contribution to the high score populations among the interhelical pairs in most of the selected sets, including class B, C and frizzled receptors. These findings are expected to be valuable in further studies of GPCRs with unknown structure and of other protein families.

A mechanistic role of Helix 8 in GPCRs: Computational modeling of the dopamine D2 receptor interaction with the GIPC1-PDZ-domain

Helix-8 (Hx8) is a structurally conserved amphipathic helical motif in class-A GPCRs, adjacent to the C-terminal sequence that is responsible for PDZ-domain-recognition. The Hx8 segment in the dopamine D2 receptor (D2R) constitutes the C-terminal segment and we investigate its role in the function of D2R by studying the interaction with the PDZ-containing GIPC1 using homology models based on the X-ray structures of very closely related analogs: the D3R for the D2R model, and the PDZ domain of GIPC2 for GIPC1-PDZ. The mechanism of this interaction was investigated with all-atom unbiased molecular dynamics (MD) simulations that reveal the role of the membrane in maintaining the helical fold of Hx8, and with biased MD simulations to elucidate the energy drive for the interaction with the GIPC1-PDZ. We found that it becomes more favorable energetically for Hx8 to adopt the extended conformation observed in all PDZ-ligand complexes when it moves away from the membrane, and that C-terminus palmitoylation of D2R enhanced membrane penetration by the Hx8 backbone. De-palmitoylation enables Hx8 to move out into the aqueous environment for interaction with the PDZ domain. All-atom unbiased MD simulations of the full D2R-GIPC1-PDZ complex in sphingolipid/cholesterol membranes show that the D2R carboxyl C-terminus samples the region of the conserved GFGL motif located on the carboxylate-binding loop of the GIPC1-PDZ, and the entire complex distances itself from the membrane interface. Together, these results outline a likely mechanism of Hx8 involvement in the interaction of the GPCR with PDZ-domains in the course of signaling.

Update on the Mechanism of Action of Aripiprazole: Translational Insights into Antipsychotic Strategies Beyond Dopamine Receptor Antagonism

Dopamine partial agonism and functional selectivity have been innovative strategies in the pharmacological treatment of schizophrenia and mood disorders and have shifted the concept of dopamine modulation beyond the established approach of dopamine D2 receptor (D2R) antagonism. Despite the fact that aripiprazole was introduced in therapy more than 12 years ago, many questions are still unresolved regarding the complexity of the effects of this agent on signal transduction and intracellular pathways, in part linked to its pleiotropic receptor profile. The complexity of the mechanism of action has progressively shifted the conceptualization of this agent from partial agonism to functional selectivity. From the induction of early genes to modulation of scaffolding proteins and activation of transcription factors, aripiprazole has been shown to affect multiple cellular pathways and several cortical and subcortical neurotransmitter circuitries. Growing evidence shows that, beyond the consequences of D2R occupancy, aripiprazole has a unique neurobiology among available antipsychotics. The effect of chronic administration of aripiprazole on D2R affinity state and number has been especially highlighted, with relevant translational implications for long-term treatment of psychosis. The hypothesized effects of aripiprazole on cell-protective mechanisms and neurite growth, as well as the differential effects on intracellular pathways [i.e. extracellular signal-regulated kinase (ERK)] compared with full D2R antagonists, suggest further exploration of these targets by novel and future biased ligand compounds. This review aims to recapitulate the main neurobiological effects of aripiprazole and discuss the potential implications for upcoming improvements in schizophrenia therapy based on dopamine modulation beyond D2R antagonism.

New insights into the meaning and usefulness of principal component analysis of concatenated trajectories

A comparison between different conformations of a given protein, relating both structure and dynamics, can be performed in terms of combined principal component analysis (combined-PCA). To that end, a trajectory is obtained by concatenating molecular dynamics trajectories of the individual conformations under comparison. Then, the principal components are calculated by diagonalizing the correlation matrix of the concatenated trajectory. Since the introduction of this approach in 1995 it has had a large number of applications. However, the interpretation of the eigenvectors and eigenvalues so obtained is based on intuitive foundations, because analytical expressions relating the concatenated correlation matrix with those of the individual trajectories under consideration have not been provided yet. In this article, we present such expressions for the cases of two, three, and an arbitrary number of concatenated trajectories. The formulas are simple and show what is to be expected and what is not to be expected from a combined-PCA. Their correctness and usefulness is demonstrated by discussing some representative examples. The results can be summarized in a simple sentence: the correlation matrix of a concatenated trajectory is given by the average of the individual correlation matrices plus the correlation matrix of the individual averages. From this it follows that the combined-PCA of trajectories belonging to different free energy basins provides information that could also be obtained by alternative and more straightforward means.

A functional selectivity mechanism at the serotonin-2A GPCR involves ligand-dependent conformations of intracellular loop 2

With recent progress in determination of G protein-coupled receptor (GPCR) structure with crystallography, a variety of other experimental approaches (e.g., NMR spectroscopy, fluorescent-based assays, mass spectrometry techniques) are also being used to characterize state-specific and ligand-specific conformational states. MD simulations offer a powerful complementary approach to elucidate the dynamic features associated with ligand-specific GPCR conformations. To shed light on the conformational elements and dynamics of the important aspect of GPCR functional selectivity, we carried out unbiased microsecond-length MD simulations of the human serotonin 2A receptor (5-HT(2A)R) in the absence of ligand and bound to four distinct serotonergic agonists. The 5-HT(2A)R is a suitable system to study the structural features involved in the ligand-dependent conformational heterogeneity of GPCRs because it is well-characterized experimentally and exhibits a strong agonist-specific phenotype in that some 5-HT(2A)R agonists induce LSD-like hallucinations, while others lack this psychoactive property entirely. Here we report evidence for structural and dynamic differences in 5-HT(2A)R interacting with such pharmacologically distinct ligands, hallucinogens, and nonhallucinogens obtained from all-atom MD simulations. Differential ligand binding contacts were identified for structurally similar hallucinogens and nonhallucinogens and found to correspond to different conformations in the intracellular loop 2 (ICL2). From the different ICL2 conformations, functional selective phenotypes are suggested through effects on dimerization and/or distinct direct interaction with effector proteins. The findings are presented in the context of currently proposed hallucinogenesis mechanisms, and ICL2 is proposed as a fine-tuning selective switch that can differentiates modes of 5-HT(2A)R activation.

Active-state model of a dopamine D2 receptor-Gαi complex stabilized by aripiprazole-type partial agonists

Partial agonists exhibit a submaximal capacity to enhance the coupling of one receptor to an intracellular binding partner. Although a multitude of studies have reported different ligand-specific conformations for a given receptor, little is known about the mechanism by which different receptor conformations are connected to the capacity to activate the coupling to G-proteins. We have now performed molecular-dynamics simulations employing our recently described active-state homology model of the dopamine D2 receptor-Gαi protein-complex coupled to the partial agonists aripiprazole and FAUC350, in order to understand the structural determinants of partial agonism better. We have compared our findings with our model of the D2R-Gαi-complex in the presence of the full agonist dopamine. The two partial agonists are capable of inducing different conformations of important structural motifs, including the extracellular loop regions, the binding pocket and, in particular, intracellular G-protein-binding domains. As G-protein-coupling to certain intracellular epitopes of the receptor is considered the key step of allosterically triggered nucleotide-exchange, it is tempting to assume that impaired coupling between the receptor and the G-protein caused by distinct ligand-specific conformations is a major determinant of partial agonist efficacy.

Structure-based simulations reveal concerted dynamics of GPCR activation

G protein-coupled receptors (GPCRs) are a vital class of proteins that transduce biological signals across the cell membrane. However, their allosteric activation mechanism is not fully understood; crystal structures of active and inactive receptors have been reported, but the functional pathway between these two states remains elusive. Here, we use structure-based (Gō-like) models to simulate activation of two GPCRs, rhodopsin and the β₂ adrenergic receptor (β₂AR). We used data-derived reaction coordinates that capture the activation mechanism for both proteins, showing that activation proceeds through quantitatively different paths in the two systems. Both reaction coordinates are determined from the dominant concerted motions in the simulations so the technique is broadly applicable. There were two surprising results. First, the main structural changes in the simulations were distributed throughout the transmembrane bundle, and not localized to the obvious areas of interest, such as the intracellular portion of Helix 6. Second, the activation (and deactivation) paths were distinctly nonmonotonic, populating states that were not simply interpolations between the inactive and active structures. These transitions also suggest a functional explanation for β₂AR's basal activity: it can proceed through a more broadly defined path during the observed transitions.

Simulating G protein-coupled receptors in native-like membranes: from monomers to oligomers

G protein-coupled receptors (GPCRs) are one of the most relevant superfamilies of transmembrane proteins as they participate in an important variety of biological events. Recently, the scientific community is witnessing an advent of a GPCR crystallization age along with impressive improvements achieved in the field of computer simulations during the last two decades. Computer simulation techniques such as molecular dynamics (MD) simulations are now frequent tools to study the dynamic behavior of GPCRs and, more importantly, to model the complex membrane environment where these proteins spend their lifetime. Thanks to these tools, GPCRs can be simulated not only longer but also in a more "physiological" fashion. In this scenario, scientists are taking advantage of such advances to approach certain phenomena such as GPCR oligomerization occurring only at timescales not reachable until now. Thus, despite current MD simulations having important limitations today, they have become an essential tool to study key biophysical properties of GPCRs and GPCR oligomers.

Not just an oil slick: how the energetics of protein-membrane interactions impacts the function and organization of transmembrane proteins

The membrane environment, its composition, dynamics, and remodeling, have been shown to participate in the function and organization of a wide variety of transmembrane (TM) proteins, making it necessary to study the molecular mechanisms of such proteins in the context of their membrane settings. We review some recent conceptual advances enabling such studies, and corresponding computational models and tools designed to facilitate the concerted experimental and computational investigation of protein-membrane interactions. To connect productively with the high resolution achieved by cognate experimental approaches, the computational methods must offer quantitative data at an atomistically detailed level. We show how such a quantitative method illuminated the mechanistic importance of a structural characteristic of multihelical TM proteins, that is, the likely presence of adjacent polar and hydrophobic residues at the protein-membrane interface. Such adjacency can preclude the complete alleviation of the well-known hydrophobic mismatch between TM proteins and the surrounding membrane, giving rise to an energy cost of residual hydrophobic mismatch. The energy cost and biophysical formulation of hydrophobic mismatch and residual hydrophobic mismatch are reviewed in the context of their mechanistic role in the function of prototypical members of multihelical TM protein families: 1), LeuT, a bacterial homolog of mammalian neurotransmitter sodium symporters; and 2), rhodopsin and the β1- and β2-adrenergic receptors from the G-protein coupled receptor family. The type of computational analysis provided by these examples is poised to translate the rapidly growing structural data for the many TM protein families that are of great importance to cell function into ever more incisive insights into mechanisms driven by protein-ligand and protein-protein interactions in the membrane environment.

Free energy landscape of G-protein coupled receptors, explored by accelerated molecular dynamics

G-protein coupled receptors (GPCRs) mediate cellular responses to various hormones and neurotransmitters and are important targets for treating a wide spectrum of diseases. They are known to adopt multiple conformational states (e.g., inactive, intermediate and active) during their modulation of various cell signaling pathways. Here, the free energy landscape of GPCRs is explored using accelerated molecular dynamics (aMD) simulations as demonstrated on the M2 muscarinic receptor, a key GPCR that regulates human heart rate and contractile forces of cardiomyocytes. Free energy profiles of important structural motifs that undergo conformational transitions upon GPCR activation and allosteric signaling are analyzed in detail, including the Arg(3.50)-Glu(6.30) ionic lock, the Trp(6.48) toggle switch and the hydrogen interactions between Tyr(5.58)-Tyr(7.53).

Molecular dynamics simulations of the adenosine A2a receptor in POPC and POPE lipid bilayers: effects of membrane on protein behavior

Analysis of 300 ns (ns) molecular dynamics (MD) simulations of an adenosine A2a receptor (A2a AR) model, conducted in triplicate, in 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and 1-palmitoyl-2-oleoylphosphatidylethanolamine (POPE) bilayers reveals significantly different protein dynamical behavior. Principal component analysis (PCA) shows that the dissimilarities stem from interhelical rather than intrahelical motions. The difference in the hydrophobic thicknesses of these simulated lipid bilayers is potentially a significant reason for the observed difference in results. The distinct lipid headgroups might also lead to different molecular interactions and hence different protein loop motions. Overall, the A2a AR shows higher mobility and flexibility in POPC as compared to POPE.

How the dynamic properties and functional mechanisms of GPCRs are modulated by their coupling to the membrane environment

Experimental observations of the dependence of function and organization of G protein-coupled receptors (GPCRs) on their lipid environment have stimulated new quantitative studies of the coupling between the proteins and the membrane. It is important to develop such a quantitative understanding at the molecular level because the effects of the coupling are seen to be physiologically and clinically significant. Here we review findings that offer insight into how membrane-GPCR coupling is connected to the structural characteristics of the GPCR, from sequence to 3D structural detail, and how this coupling is involved in the actions of ligands on the receptor. The application of a recently developed computational approach designed for quantitative evaluation of membrane remodeling and the energetics of membrane-protein interactions brings to light the importance of the radial asymmetry of the membrane-facing surface of GPCRs in their interaction with the surrounding membrane. As the radial asymmetry creates adjacencies of hydrophobic and polar residues at specific sites of the GPCR, the ability of membrane remodeling to achieve complete hydrophobic matching is limited, and the residual mismatch carries a significant energy cost. The adjacencies are shown to be affected by ligand-induced conformational changes. Thus, functionally important organization of GPCRs in the cell membrane can depend both on ligand-determined properties and on the lipid composition of various membrane regions with different remodeling capacities. That this functionally important reorganization can be driven by oligomerization patterns that reduce the energy cost of the residual mismatch, suggests a new perspective on GPCR dimerization and ligand-GPCR interactions. The relation between the modulatory effects on GPCRs from the binding of specific cell-membrane components, e.g., cholesterol, and those produced by the non-local energetics of hydrophobic mismatch are discussed in this context.

Structural features of the regulatory ACT domain of phenylalanine hydroxylase

Phenylalanine hydroxylase (PAH) catalyzes the conversion of L-Phe to L-Tyr. Defects in PAH activity, caused by mutations in the human gene, result in the autosomal recessively inherited disease hyperphenylalaninemia. PAH activity is regulated by multiple factors, including phosphorylation and ligand binding. In particular, PAH displays positive cooperativity for L-Phe, which is proposed to bind the enzyme on an allosteric site in the N-terminal regulatory domain (RD), also classified as an ACT domain. This domain is found in several proteins and is able to bind amino acids. We used molecular dynamics simulations to obtain dynamical and structural insights into the isolated RD of PAH. Here we show that the principal motions involve conformational changes leading from an initial open to a final closed domain structure. The global intrinsic motions of the RD are correlated with exposure to solvent of a hydrophobic surface, which corresponds to the ligand binding-site of the ACT domain. Our results strongly suggest a relationship between the Phe-binding function and the overall dynamic behaviour of the enzyme. This relationship may be affected by structure-disturbing mutations. To elucidate the functional implications of the mutations, we investigated the structural effects on the dynamics of the human RD PAH induced by six missense hyperphenylalaninemia-causing mutations, namely p.G46S, p.F39C, p.F39L, p.I65S, p.I65T and p.I65V. These studies showed that the alterations in RD hydrophobic interactions induced by missense mutations could affect the functionality of the whole enzyme.

Membrane driven spatial organization of GPCRs

Spatial organization of G-protein coupled receptors (GPCRs) into dimers and higher order oligomers has been demonstrated in vitro and in vivo. The pharmacological readout was shown to depend on the specific interfaces, but why particular regions of the GPCR structure are involved, and how ligand-determined states change them remains unknown. Here we show why protein-membrane hydrophobic matching is attained upon oligomerization at specific interfaces from an analysis of coarse-grained molecular dynamics simulations of the spontaneous diffusion-interaction of the prototypical beta2-adrenergic (β₂AR) receptors in a POPC lipid bilayer. The energy penalty from mismatch is significantly reduced in the spontaneously emerging oligomeric arrays, making the spatial organization of the GPCRs dependent on the pattern of mismatch in the monomer. This mismatch pattern is very different for β₂AR compared to the highly homologous and structurally similar β₁AR, consonant with experimentally observed oligomerization patterns of β₂AR and β₁AR. The results provide a mechanistic understanding of the structural context of oligomerization.

Molecular pharmacology and ligand docking studies reveal a single amino acid difference between mouse and human serotonin 5-HT2A receptors that impacts behavioral translation of novel 4-phenyl-2-dimethylaminotetralin ligands

During translational studies to develop 4-phenyl-2-dimethylaminotetralin (PAT) compounds for neuropsychiatric disorders, the (2R,4S)-trans-(+)- and (2S,4R)-trans-(-)-enantiomers of the analog 6-hydroxy-7-chloro-PAT (6-OH-7-Cl-PAT) demonstrated unusual pharmacology at serotonin (5-HT) 5-HT2 G protein-coupled receptors (GPCRs). The enantiomers had similar affinities (Ki) at human (h) 5-HT2A receptors (≈ 70 nM). In an in vivo mouse model of 5-HT2A receptor activation [(±)-(2,5)-dimethoxy-4-iodoamphetamine (DOI)-elicited head twitch], however, (-)-6-OH-7-Cl-PAT was about 5-fold more potent than the (+)-enantiomer at attenuating the DOI-elicited response. It was discovered that (+)-6-OH-7-Cl-PAT (only) had ≈ 40-fold-lower affinity at mouse (m) compared with h5-HT2A receptors. Molecular modeling and computational ligand docking studies indicated that the 6-OH moiety of (+)- but not (-)-6-OH-7-Cl-PAT could form a hydrogen bond with serine residue 5.46 of the h5-HT2A receptor. The m5-HT2A as well as m5-HT2B, h5-HT2B, m5-HT2C, and h5-HT2C receptors have alanine at position 5.46, obviating this interaction; (+)-6-OH-7-Cl-PAT also showed ≈ 50-fold lower affinity than (-)-6-OH-7-Cl-PAT at m5-HT2C and h5-HT2C receptors. Mutagenesis studies confirmed that 5-HT2A S5.46 is critical for (+)- but not (-)-6-OH-7-Cl-PAT binding, as well as function. The (+)-6-OH-7-Cl-PAT enantiomer showed partial agonist effects at h5-HT2A wild-type (WT) and m5-HT2A A5.46S point-mutated receptors but did not activate m5-HT2A WT and h5-HT2A S5.46A point-mutated receptors, or h5-HT2B, h5-HT2C, and m5-HT2C receptors; (-)-6-OH-7-Cl-PAT did not activate any of the 5-HT2 receptors. Experiments also included the (2R,4S)-trans-(+)- and (2S,4R)-trans-(-)-enantiomers of 6-methoxy-7-chloro-PAT to validate hydrogen bonding interactions proposed for the corresponding 6-OH analogs. Results indicate that PAT ligand three-dimensional structure impacts target receptor binding and translational outcomes, supporting the hypothesis that GPCR ligand structure governs orthosteric binding pocket molecular determinants and resulting pharmacology.

The membrane protein LeuT in micellar systems: aggregation dynamics and detergent binding to the S2 site

Structural and functional properties of integral membrane proteins are often studied in detergent micellar environments (proteomicelles), but how such proteomicelles form and organize is not well understood. This makes it difficult to evaluate the relationship between the properties of the proteins measured in such a detergent-solubilized form and under native conditions. To obtain mechanistic information about this relationship for the leucine transporter (LeuT), a prokaryotic homologue of the mammalian neurotransmitter/sodium symporters (NSSs), we studied the properties of proteomicelles formed by n-dodecyl-β,D-maltopyranoside (DDM) detergent. Extensive atomistic molecular dynamics simulations of different protein/detergent/water number ratios revealed the formation of a proteomicelle characterized by a constant-sized shell of detergents surrounding LeuT protecting its transmembrane segments from unfavorable hydrophobic/hydrophilic exposure. Regardless of the DDM content in the simulated system, this shell consisted of a constant number of DDM molecules (∼120 measured at a 4 Å cutoff distance from LeuT). In contrast, the overall number of DDMs in the proteomicelle (aggregation number) was found to depend on the detergent concentration, reaching a saturation value of 226±17 DDMs in the highest concentration regime simulated. Remarkably, we found that at high detergent-to-protein ratios we observed two independent ways of DDM penetration into LeuT, both leading to a positioning of the DDM molecule in the second substrate (S2) binding site of LeuT. Consonant with several recent experimental studies demonstrating changes in functional properties of membrane proteins due to detergent, our findings highlight how the environment in which the membrane proteins are examined may affect the outcome and interpretation of their mechanistic features.

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.

Conformational restriction of G-proteins Coupled Receptors (GPCRs) upon complexation to G-proteins: a putative activation mode of GPCRs?

GPCRs undergo large conformational changes during their activation. Starting from existing X-ray structures, we used Normal Modes Analyses to study the collective motions of the agonist-bound β2-adrenergic receptor both in its isolated "uncoupled" and G-protein "coupled" conformations. We interestingly observed that the receptor was able to adopt only one major motion in the protein:protein complex. This motion corresponded to an anti-symmetric rotation of both its extra- and intra-cellular parts, with a key role of previously identified highly conserved proline residues. Because this motion was also retrieved when performing NMA on 7 other GPCRs which structures were available, it is strongly suspected to possess a significant biological role, possibly being the "activation mode" of a GPCR when coupled to G-proteins.

Ligand-dependent activation and deactivation of the human adenosine A(2A) receptor

G-protein-coupled receptors (GPCRs) are membrane proteins with critical functions in cellular signal transduction, representing a primary class of drug targets. Acting by direct binding, many drugs modulate GPCR activity and influence the signaling pathways associated with numerous diseases. However, complete details of ligand-dependent GPCR activation/deactivation are difficult to obtain from experiments. Therefore, it remains unclear how ligands modulate a GPCR's activity. To elucidate the ligand-dependent activation/deactivation mechanism of the human adenosine A2A receptor (AA2AR), a member of the class A GPCRs, we performed large-scale unbiased molecular dynamics and metadynamics simulations of the receptor embedded in a membrane. At the atomic level, we have observed distinct structural states that resemble the active and inactive states. In particular, we noted key structural elements changing in a highly concerted fashion during the conformational transitions, including six conformational states of a tryptophan (Trp246(6.48)). Our findings agree with a previously proposed view that, during activation, this tryptophan residue undergoes a rotameric transition that may be coupled to a series of coherent conformational changes, resulting in the opening of the G-protein binding site. Further, metadynamics simulations provide quantitative evidence for this mechanism, suggesting how ligand binding shifts the equilibrium between the active and inactive states. Our analysis also proposes that a few specific residues are associated with agonism/antagonism, affinity, and selectivity, and suggests that the ligand-binding pocket can be thought of as having three distinct regions, providing dynamic features for structure-based design. Additional simulations with AA2AR bound to a novel ligand are consistent with our proposed mechanism. Generally, our study provides insights into the ligand-dependent AA2AR activation/deactivation in addition to what has been found in crystal structures. These results should aid in the discovery of more effective and selective GPCR ligands.

Computational Insights into the Mechanism of Inhibition of OASS-A by a Small Molecule Inhibitor

O-Acetylserine sulfhydrylase (isoform A, OASS-A) is a PLP-dependent enzyme involved in the last step of cysteine biosynthesis in many pathogens. Many microorganisms use cysteine as the main building block for sulfur-containing antioxidants, and cysteine depletion in several pathogens resulted in a reduced antibiotic resistance, thus leading to the identification of OASS as novel suitable molecular targets to overcome antimicrobial resistances. The precise molecular mechanism of OASS-A inhibition by small peptides or by small molecule inhibitors is still unclear. To shed more lights on the structural basis underlying the inhibition mechanism for OASS, we engaged ourselves in studying the dynamic properties of this enzyme. In this paper, we describe a computational study involving unbiased MD simulations of OASS-A from Haemophilus influenzae (HiOASS) in its inhibitor free, PLP-bound form, and in complex with a pentapeptide inhibitor and with UPAR40, a small molecule which we have recently reported as a potent OASS-A inhibitors. We proposed that UPAR40 inhibits HiOASS-A through the stabilization of a closed conformation. Moreover, preliminary docking studies and sequence analysis allow us to speculate about the non-specificity of UPAR40 toward a particular OASS enzyme species or isoforms.

The cost of living in the membrane: a case study of hydrophobic mismatch for the multi-segment protein LeuT

Many observations of the role of the membrane in the function and organization of transmembrane (TM) proteins have been explained in terms of hydrophobic mismatch between the membrane and the inserted protein. For a quantitative investigation of this mechanism in the lipid-protein interactions of functionally relevant conformations adopted by a multi-TM segment protein, the bacterial leucine transporter (LeuT), we employed a novel method, Continuum-Molecular Dynamics (CTMD), that quantifies the energetics of hydrophobic mismatch by combining the elastic continuum theory of membrane deformations with an atomistic level description of the radially asymmetric membrane-protein interface from MD simulations. LeuT has been serving as a model for structure-function studies of the mammalian neurotransmitter:sodium symporters (NSSs), such as the dopamine and serotonin transporters, which are the subject of intense research in the field of neurotransmission. The membrane models in which LeuT was embedded for these studies were composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid, or 3:1 mixture of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) lipids. The results show that deformation of the host membrane alone is not sufficient to alleviate the hydrophobic mismatch at specific residues of LeuT. The calculations reveal significant membrane thinning and water penetration due to the specific local polar environment produced by the charged K288 of TM7 in LeuT, that is membrane-facing deep inside the hydrophobic milieu of the membrane. This significant perturbation is shown to result in unfavorable polar-hydrophobic interactions at neighboring hydrophobic residues in TM1a and TM7. We show that all the effects attributed to the K288 residue (membrane thinning, water penetration, and the unfavorable polar-hydrophobic interactions at TM1a and TM7), are abolished in calculations with the K288A mutant. The involvement of hydrophobic mismatch is somewhat different in the functionally distinct conformations (outward-open, occluded, inward-open) of LeuT, and the differences are shown to connect to structural elements (e.g., TM1a) known to play key roles in transport. This finding suggests a mechanistic hypothesis for the enhanced transport activity observed for the K288A mutant, suggesting that the unfavorable hydrophobic-hydrophilic interactions hinder the motion of TM1a in the functionally relevant conformational transition to the inward-open state. Various extents of such unfavorable interactions, involving exposure to the lipid environment of adjacent hydrophobic and polar residues, are common in multi-segment transmembrane proteins, and must be considered to affect functionally relevant conformational transitions.

Behavioral and serotonergic response changes in the Dhcr7-HET mouse model of Smith-Lemli-Opitz syndrome

Smith-Lemli-Opitz syndrome (SLOS) is a developmental disorder resulting from mutations to the Dhcr7 gene, which is required for cholesterol synthesis. Patients with SLOS typically exhibit a number of severe behavioral deficits and many are diagnosed with autistic spectrum disorder. Although the molecular pathophysiology underlying behavioral changes in SLOS and autism spectrum disorders is poorly understood, there is evidence for the involvement of the serotonergic system in SLOS and autism in general. Behavioral testing was undertaken to ascertain the basal behavioral differences between Dhcr7-heterozygous (HET) and wild-type control mice and explore the utility of a Dhcr7-HET mouse line in the development of new treatments for this disorder. Dhcr7-HET mice did not differ from wild-type control mice on basic measures of locomotor activity, anxiety and neuromuscular ability. However, female Dhcr7-HET mice at 6 months of age or older were significantly more likely to win on the social dominance tube test against an unfamiliar mouse. Pharmacological testing, using the 5-HT2A agonist 1-(2,5-dimethoxy-4-iodophenyl)-2-aminopropane (DOI), showed increased head-twitch response in Dhcr7-HET mice, which was apparent from 6 months of age. No differences were found between the genotypes in testing for 5-HT1A agonist 8-OH-DPAT-induced hypothermia. These data indicate an underlying dysfunction of the 5-HT2A receptors in Dhcr7-HET mice that warrants further investigation to establish how this may relate to behavioral disturbances in human patients carrying Dhcr7 mutations.

Sphingosine-1-phosphate and its receptors: structure, signaling, and influence

The sphingosine-1-phosphate (S1P) receptor signaling system has biological and medical importance and is the first lipid G protein-coupled receptor (GPCR) structure to be solved to 2.8-Å resolution. S1P binds to five high-affinity GPCRs generating multiple downstream signals that play essential roles in vascular development and endothelial integrity, control of cardiac rhythm, and routine oral treatment of multiple sclerosis. Genetics, chemistry, and now structural biology have advanced this integrated biochemical system. The S1P receptors have a novel N-terminal fold that occludes access to the binding pocket from the extracellular environment as well as orthosteric and bitopic ligands with very different physicochemical properties. S1P receptors and metabolizing enzymes have been deleted, inducibly deleted, and knocked in as tagged or altered receptors in mice. An array of genetic models allows analysis of integrated receptor function in vivo. We can now directly understand causal relationships among protein expression, signal, and control points in physiology and pathology.

Why GPCRs behave differently in cubic and lamellar lipidic mesophases

Recent successes in the crystallographic determination of structures of transmembrane proteins in the G protein-coupled receptor (GPCR) family have established the lipidic cubic phase (LCP) environment as the medium of choice for growing structure-grade crystals by the method termed "in meso". The understanding of in meso crystallogenesis is currently at a descriptive level. To enable an eventual quantitative, energy-based description of the nucleation and crystallization mechanism, we have examined the properties of the lipidic cubic phase system and the dynamics of the GPCR rhodopsin reconstituted into the LCP with coarse-grained molecular dynamics simulations with the Martini force-field. Quantifying the differences in the hydrophobic/hydrophilic exposure of the GPCR to lipids in the cubic and lamellar phases, we found that the highly curved geometry of the cubic phase provides more efficient shielding of the protein from unfavorable hydrophobic exposure, which leads to a lesser hydrophobic mismatch and less unfavorable hydrophobic-hydrophilic interactions between the protein and lipid-water interface in the LCP, compared to the lamellar phase. Since hydrophobic mismatch is considered a driving force for oligomerization, the differences in exposure mismatch energies between the LCP and the lamellar structures suggest that the latter provide a more favorable setting in which GPCRs can oligomerize as a prelude to nucleation and crystal growth. These new findings lay the foundation for future investigations of in meso crystallization mechanisms related to the transition from the LCP to the lamellar phase and studies aimed at an improved rational approach for generating structure-quality crystals of membrane proteins.