Role of aggregation in rhodopsin signal transduction

Biased Activation Mechanism Induced by GPCR Heterodimerization: Observations from μOR/δOR Dimers

GPCRs regulate multiple intracellular signaling cascades. Biasedly activating one signaling pathway over the others provides additional clinical utility to optimize GPCR-based therapies. GPCR heterodimers possess different functions from their monomeric states, including their selectivity to different transducers. However, the biased signaling mechanism induced by the heterodimerization remains unclear. Motivated by the issue, we select an important GPCR heterodimer (μOR/δOR heterodimer) as a case and use microsecond Gaussian accelerated molecular dynamics simulation coupled with potential of mean force and protein structure network (PSN) to probe mechanisms regarding the heterodimerization-induced constitutive β-arrestin activity and efficacy change of the agonist DAMGO. The results show that only the lowest energy state of the μOR/δOR heterodimer, which adopts a slightly outward shift of TM6 and an ICL2 conformation close to the receptor core, can selectively accommodate β-arrestins. PSN further reveals important roles of H8, ICL1, and ICL2 in regulating the constitutive β-arrestin-biased activity for the apo μOR/δOR heterodimer. In addition, the heterodimerization can allosterically alter the binding mode of DAMGO mainly by means of W7.35. Consequently, DAMGO transmits the structural signal mainly through TM6 and TM7 in the dimer, rather than TM3 similar to the μOR monomer, thus changing the efficacy of DAMGO from a balanced agonist to the β-arrestin-biased one. On the other side, the binding of DAMGO to the heterodimer can stabilize μOR/δOR heterodimers through a stronger interaction of TM1/TM1 and H8/H8, accordingly enhancing the interaction of μOR with δOR and the binding affinity of the dimer to the β-arrestin. The agonist DAMGO does not change main compositions of the regulation network from the dimer interface to the transducer binding pocket of the μOR protomer, but induces an increase in the structural communication of the network, which should contribute to the enhanced β-arrestin coupling. Our observations, for the first time, reveal the molecular mechanism of the biased signaling induced by the heterodimerization for GPCRs, which should be beneficial to more comprehensively understand the GPCR bias signaling.

Computational and experimental approaches to probe GPCR activation and signaling

G protein-coupled receptors (GPCRs) regulate different physiological functions, e.g., sensation, growth, digestion, reproductivity, nervous and immune systems response, and many others. In eukaryotes, they are also responsible for intercellular communication in response to pathogens. The major primary messengers binding to these cell-surface receptors constitute small-molecule or peptide hormones and neurotransmitters, nucleotides, lipids as well as small proteins. The simplicity of the way how GPCR signaling can be regulated by their endogenous agonists prompted the usage of GPCRs as major drug targets in modern pharmacology. Drugs targeting GPCRs inhibit pathological processes at the very beginning. This enables to significantly reduce the occurrence of morphological changes caused by diseases. Until recently, X-ray crystallography was the method of the first choice to obtain high-resolution structural information about GPCRs. Following X-ray crystallography, cryo-EM gained attention in GPCR studies as a quick and low-cost alternative. FRET microscopy is also widely used for GPCRs in the analysis of protein-protein interactions (PPIs) in intact cells as well as for screening purposes. Regarding computational methods, molecular dynamics (MD) for many years has proven its usefulness in studying the GPCR activation. MODELLER and Rosetta were widely used to generate preliminary homology models of GPCRs for MD simulation systems. Apart from the conventional all-atom approach with explicitly defined solvent, also other techniques have been applied to GPCRs, e.g., MARTINI or hybrid methods involving the coarse-grained representation, less demanding regarding computational resources, and thus offering much larger simulation timescales.

Supramolecular organization of rhodopsin in rod photoreceptor cell membranes

Rhodopsin is the light receptor in rod photoreceptor cells that initiates scotopic vision. Studies on the light receptor span well over a century, yet questions about the organization of rhodopsin within the photoreceptor cell membrane still persist and a consensus view on the topic is still elusive. Rhodopsin has been intensely studied for quite some time, and there is a wealth of information to draw from to formulate an organizational picture of the receptor in native membranes. Early experimental evidence in apparent support for a monomeric arrangement of rhodopsin in rod photoreceptor cell membranes is contrasted and reconciled with more recent visual evidence in support of a supramolecular organization of rhodopsin. What is known so far about the determinants of forming a supramolecular structure and possible functional roles for such an organization are also discussed. Many details are still missing on the structural and functional properties of the supramolecular organization of rhodopsin in rod photoreceptor cell membranes. The emerging picture presented here can serve as a springboard towards a more in-depth understanding of the topic.

Residue-Residue Mutual Work Analysis of Retinal-Opsin Interaction in Rhodopsin: Implications for Protein-Ligand Binding

Energetic contributions at the single-residue level for retinal-opsin interactions in rhodopsin were studied by combining molecular dynamics simulations, transition path sampling, and a newly developed energy decomposition approach. The virtual work at an infinitesimal time interval was decomposed into the work components on one residue due to its interaction with another residue, which were then averaged over the transition path ensemble along a proposed reaction coordinate. Such residue-residue mutual work analysis on 62 residues within the active center of rhodopsin resulted in a very sparse interaction matrix, which is generally not symmetric but antisymmetric to some extent. Fourteen residues were identified to be major players in retinal relaxation along a plausible pathway from bathorhodopsin to the blue-shifted intermediate, which is in good agreement with an existing NMR study. Based on the matrix of mutual work, a comprehensive network was constructed to provide detailed insights into the chromophore-protein interaction from a viewpoint of energy flow.

Exploring the Activation Mechanism of a Metabotropic Glutamate Receptor Homodimer via Molecular Dynamics Simulation

Metabotropic glutamate receptors of class C GPCRs exist as constitutive dimers, which play important roles in activating excitatory synapses of the central nervous system. However, the activation mechanism induced by agonists has not been clarified in experiments. To address the problem, we used microsecond all-atom molecular dynamics (MD) simulation couple with protein structure network (PSN) to explore the glutamate-induced activation for the mGluR1 homodimer. The results indicate that glutamate binding stabilizes not only the closure of Venus flytrap domains but also the polar interaction of LB2-LB2, in turn keeping the extracelluar domain in the active state. The activation of the extracelluar domain drives transmembrane domains (TMDs) of the two protomers closer and induces asymmetric activation for the TMD domains of the two protomers. One protomer with lower binding affinity to the agonist is activated, while the other protomer with higher binding energy is still in the inactive state. The PSN analysis identifies the allosteric regulation pathway from the ligand-binding pocket in the extracellular domain to the G-protein binding site in the intracellular TMD region and further reveals that the asymmetric activation is attributed to a combination of trans-pathway and cis-pathway regulations from two glumatates, rather than a single activation pathway. These observations could provide valuable molecular information for understanding of the structure and the implications in drug efficacy for the class C GPCR dimers.

Prediction and targeting of GPCR oligomer interfaces

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.

The effect of dimerization on the activation and conformational dynamics of adenosine A1 receptor

The adenosine A1 receptor (A1R) is one of four adenosine receptors in humans, which are involved in the function of the cardiovascular, respiratory and central nervous systems. Experimental results indicate that A1R can form a homodimer and that the protomer-protomer interaction in the A1R dimer is related to certain pharmacological characteristics of A1R activation. In this work, we performed docking, metadynamics simulation, conventional molecular dynamics simulations, Gaussian-accelerated molecular dynamics simulations, potential of mean force calculations, dynamic cross-correlation motions analysis and community network analysis to study the binding mode of 5'-N-ethylcarboxamidoadenosine (NECA) to A1R and the effect of dimerization on the activation of A1R. Our results show that NECA binds to A1R in a similar mode to adenosine in the A1R crystal structure and NECA in the A2AR crystal structure. The A1R homodimer can be activated by one or two agonists with NECA occupying its orthosteric pockets in one (which we call the NECA-A1R system) or both protomers (which we call the dNECA-A1R system). In the NECA-A1R system, activation is predicated in the protomer without NECA bound. In the dNECA-A1R system, only one protomer achieves the active state. These findings suggest an asymmetrical activation mechanism of the homodimer and a negative cooperativity between the two protomers. We envision that our results may further facilitate the drug development of A1R.

Use multiscale simulation to explore the effects of the homodimerizations between different conformation states on the activation and allosteric pathway for the μ-opioid receptor

Recently, oligomers of G-protein coupled receptors (GPCRs) have been an important topic in the GPCR fields. However, knowledge about their structures and activation mechanisms is very limited due to the absence of crystal structures reported. In this work, we used multiscale simulations to study the effects of homodimerization between different conformation states on their activation, dynamic behaviors, and allosteric communication pathways for μ-OR. The results indicated that the dimerization of one inactive monomer with either one inactive monomer or one active one could enhance its constitutive activation. However, the conformation state of the other protomer (e.g., active or inactive) can influence the activated extent. The dimerization between the two inactive protomers leads to a negative cooperativity for their activation, which should contribute to the asymmetric activation of GPCR dimers observed in some experiments. On the other hand, for the active monomer, its dimerization with one inactive receptor could alleviate its deactivation, whereby negative and positive cooperativities can be observed between the two subunits of the dimer, depending on the different regions. Observations from protein structure network (PSN) analysis indicated that the dimerization of one inactive monomer with one active one would cause a significant drop in the number of main pathways from the ligand binding pocket to the G-protein coupled region for the inactive protomer, while the impact is minor for the active protomer. But, for the active monomer or the inactive one, its dimerization with one inactive monomer would significantly change the types of residues participating in the pathway with the highest frequency.

Apo-Opsin Exists in Equilibrium Between a Predominant Inactive and a Rare Highly Active State

Bleaching adaptation in rod photoreceptors is mediated by apo-opsin, which activates phototransduction with effective activity 10⁵- to 10⁶-fold lower than that of photoactivated rhodopsin (meta II). However, the mechanism that produces such low opsin activity is unknown. To address this question, we sought to record single opsin responses in mouse rods. We used mutant mice lacking efficient calcium feedback to boosts rod responses and generated a small fraction of opsin by photobleaching ∼1% of rhodopsin. The bleach produced a dramatic increase in the frequency of discrete photoresponse-like events. This activity persisted for hours, was quenched by 11-cis-retinal, and was blocked by uncoupling opsin from phototransduction, all indicating opsin as its source. Opsin-driven discrete activity was also observed in rods containing non-activatable rhodopsin, ruling out transactivation of rhodopsin by opsin. We conclude that bleaching adaptation is mediated by opsin that exists in equilibrium between a predominant inactive and a rare meta II-like state.SIGNIFICANCE STATEMENT Electrophysiological analysis is used to show that the G-protein-coupled receptor opsin exists in equilibrium between a predominant inactive and a rare highly active state that mediates bleaching adaptation in photoreceptors.

Detection of misfolded rhodopsin aggregates in cells by Förster resonance energy transfer

Rhodopsin is the light receptor in rod photoreceptor cells of the retina that plays a central role in phototransduction and rod photoreceptor cell health. Rhodopsin mutations are the leading known cause of autosomal dominant retinitis pigmentosa, a retinal degenerative disease. A majority of rhodopsin mutations cause misfolding and aggregation of the apoprotein opsin. The nature of aggregates formed by misfolded rhodopsin mutants and the associated cell toxicity is poorly understood. Misfolding rhodopsin mutants have been characterized biochemically, and categorized as either partial or complete misfolding mutants. This classification is incomplete and does not provide sufficient information to fully understand rhodopsin aggregation, disease pathogenesis, and evaluate therapeutic strategies. To better understand the aggregation of misfolded rhodopsin mutants, a Förster resonance energy transfer assay has been developed to monitor the aggregation of fluorescently tagged mutant rhodopsins expressed in live cells.

Misfolded rhodopsin mutants display variable aggregation properties

The largest class of rhodopsin mutations causing autosomal dominant retinitis pigmentosa (adRP) is mutations that lead to misfolding and aggregation of the receptor. The misfolding mutants have been characterized biochemically, and categorized as either partial or complete misfolding mutants. This classification is incomplete and does not provide sufficient information to fully understand the disease pathogenesis and evaluate therapeutic strategies. A Förster resonance energy transfer (FRET) method was utilized to directly assess the aggregation properties of misfolding rhodopsin mutants within the cell. Partial (P23H and P267L) and complete (G188R, H211P, and P267R) misfolding mutants were characterized to reveal variability in aggregation properties. The complete misfolding mutants all behaved similarly, forming aggregates when expressed alone, minimally interacting with the wild-type receptor when coexpressed, and were unresponsive to treatment with the pharmacological chaperone 9-cis retinal. In contrast, variability was observed between the partial misfolding mutants. In the opsin form, the P23H mutant behaved similarly as the complete misfolding mutants. In contrast, the opsin form of the P267L mutant existed as both aggregates and oligomers when expressed alone and formed mostly oligomers with the wild-type receptor when coexpressed. The partial misfolding mutants both reacted similarly to the pharmacological chaperone 9-cis retinal, displaying improved folding and oligomerization when expressed alone but aggregating with wild-type receptor when coexpressed. The observed differences in aggregation properties and effect of 9-cis retinal predict different outcomes in disease pathophysiology and suggest that retinoid-based chaperones will be ineffective or even detrimental.

Oligomerization and cooperativity in GPCRs from the perspective of the angiotensin AT1 and dopamine D2 receptors

G Protein-Coupled Receptors (GPCRs) can form homo- and heterodimers or constitute higher oligomeric clusters with other heptahelical GPCRs. In this article, multiscale molecular modeling approaches as well as experimental techniques which are used to study oligomerization of GPCRs are reviewed. In particular, the effect of dimerization/oligomerization to the ligand binding affinity of individual protomers and also on the efficacy of the oligomer are discussed by including diverse examples from the literature. In addition, possible allosteric effects that may emerge upon interaction of GPCRs with membrane components, like cholesterol, is also discussed. Investigation of these above-mentioned interactions may greatly contribute to the candidate molecule screening studies and development of novel therapeutics with fewer adverse effects.

Probing the cooperative mechanism of the μ–δ opioid receptor heterodimer by multiscale simulation

The activation-cooperativity of the μ–δ opioid receptor heterodimer was probed by multiscale simulation coupled with a protein structure network.

Implications of short time scale dynamics on long time processes

This review provides a comprehensive overview of the structural dynamics in topical gas- and condensed-phase systems on multiple length and time scales. Starting from vibrationally induced dissociation of small molecules in the gas phase, the question of vibrational and internal energy redistribution through conformational dynamics is further developed by considering coupled electron/proton transfer in a model peptide over many orders of magnitude. The influence of the surrounding solvent is probed for electron transfer to the solvent in hydrated I-. Next, the dynamics of a modified PDZ domain over many time scales is analyzed following activation of a photoswitch. The hydration dynamics around halogenated amino acid side chains and their structural dynamics in proteins are relevant for iodinated TyrB26 insulin. Binding of nitric oxide to myoglobin is a process for which experimental and computational analyses have converged to a common view which connects rebinding time scales and the underlying dynamics. Finally, rhodopsin is a paradigmatic system for multiple length- and time-scale processes for which experimental and computational methods provide valuable insights into the functional dynamics. The systems discussed here highlight that for a comprehensive understanding of how structure, flexibility, energetics, and dynamics contribute to functional dynamics, experimental studies in multiple wavelength regions and computational studies including quantum, classical, and more coarse grained levels are required.

An Ensemble-Based Protocol for the Computational Prediction of Helix-Helix Interactions in G Protein-Coupled Receptors using Coarse-Grained Molecular Dynamics

The accurate identification of the specific points of interaction between G protein-coupled receptor (GPCR) oligomers is essential for the design of receptor ligands targeting oligomeric receptor targets. A coarse-grained molecular dynamics computer simulation approach would provide a compelling means of identifying these specific protein–protein interactions and could be applied both for known oligomers of interest and as a high-throughput screen to identify novel oligomeric targets. However, to be effective, this in silico modeling must provide accurate, precise, and reproducible information. This has been achieved recently in numerous biological systems using an ensemble-based all-atom molecular dynamics approach. In this study, we describe an equivalent methodology for ensemble-based coarse-grained simulations. We report the performance of this method when applied to four different GPCRs known to oligomerize using error analysis to determine the ensemble size and individual replica simulation time required. Our measurements of distance between residues shown to be involved in oligomerization of the fifth transmembrane domain from the adenosine A_2A receptor are in very good agreement with the existing biophysical data and provide information about the nature of the contact interface that cannot be determined experimentally. Calculations of distance between rhodopsin, CXCR4, and β₁AR transmembrane domains reported to form contact points in homodimers correlate well with the corresponding measurements obtained from experimental structural data, providing an ability to predict contact interfaces computationally. Interestingly, error analysis enables identification of noninteracting regions. Our results confirm that GPCR interactions can be reliably predicted using this novel methodology.

Photoactivation Intermediates of a G-Protein Coupled Receptor Rhodopsin Investigated by a Hybrid Molecular Simulation

Rhodopsin is a G-protein coupled receptor functioning as a photoreceptor for vision through photoactivation of a covalently bound ligand of a retinal protonated Schiff base chromophore. Despite the availability of structural information on the inactivated and activated forms of the receptor, the transition processes initiated by the photoabsorption have not been well understood. Here we theoretically examined the photoactivation processes by means of molecular dynamics (MD) simulations and ab initio quantum mechanical/molecular mechanical (QM/MM) free energy geometry optimizations which enabled accurate geometry determination of the ligand molecule in ample statistical conformational samples of the protein. Structures of the intermediate states of the activation process, blue-shifted intermediate and Lumi, as well as the dark state first generated by MD simulations and then refined by the QM/MM free energy geometry optimizations were characterized by large displacement of the β-ionone ring of retinal along with change in the hydrogen bond of the protonated Schiff base. The ab initio calculations of vibrational and electronic spectroscopic properties of those states well reproduced the experimental observations and successfully identified the molecular origins underlying the spectroscopic features. The structural evolution in the formation of the intermediates provides a molecular insight into the efficient activation processes of the receptor.

Methods used to study the oligomeric structure of G-protein-coupled receptors

G-protein-coupled receptors (GPCRs), which constitute the largest family of cell surface receptors, were originally thought to function as monomers, but are now recognized as being able to act in a wide range of oligomeric states and indeed, it is known that the oligomerization state of a GPCR can modulate its pharmacology and function. A number of experimental techniques have been devised to study GPCR oligomerization including those based upon traditional biochemistry such as blue-native PAGE (BN-PAGE), co-immunoprecipitation (Co-IP) and protein-fragment complementation assays (PCAs), those based upon resonance energy transfer, FRET, time-resolved FRET (TR-FRET), FRET spectrometry and bioluminescence resonance energy transfer (BRET). Those based upon microscopy such as FRAP, total internal reflection fluorescence microscopy (TIRFM), spatial intensity distribution analysis (SpIDA) and various single molecule imaging techniques. Finally with the solution of a growing number of crystal structures, X-ray crystallography must be acknowledged as an important source of discovery in this field. A different, but in many ways complementary approach to the use of more traditional experimental techniques, are those involving computational methods that possess obvious merit in the study of the dynamics of oligomer formation and function. Here, we summarize the latest developments that have been made in the methods used to study GPCR oligomerization and give an overview of their application.

Impaired thromboxane receptor dimerization reduces signaling efficiency: A potential mechanism for reduced platelet function in vivo

Thromboxane A₂ is a potent mediator of inflammation and platelet aggregation exerting its effects through the activation of a G protein-coupled receptor (GPCR), termed TP. Although the existence of dimers/oligomers in Class A GPCRs is widely accepted, their functional significance still remains controversial. Recently, we have shown that TPα and TPβ homo-/hetero-dimers interact through an interface of residues in transmembrane domain 1 (TM1) whose disruption impairs dimer formation. Here, biochemical and pharmacological characterization of this dimer deficient mutant (DDM) in living cells indicates a significant impairment in its response to agonists. Interestingly, two single loss-of-function TPα variants, namely W29C and N42S recently identified in two heterozygous patients affected by bleeding disorders, match some of the residues mutated in our DDM. These two naturally occurring variants display a reduced potency to TP agonists and are characterized by impaired dimer formation in transfected HEK-293T cells. These findings provide proofs that lack of homo-dimer formation is a crucial process for reduced TPα function in vivo, and might represent one molecular mechanism through which platelet TPα receptor dysfunction affects the patient(s) carrying these mutations.

A second visual rhodopsin gene, rh1-2, is expressed in zebrafish photoreceptors and found in other ray-finned fishes

Rhodopsin (rh1) is the visual pigment expressed in rod photoreceptors of vertebrates that is responsible for initiating the critical first step of dim-light vision. Rhodopsin is usually a single copy gene; however, we previously discovered a novel rhodopsin-like gene expressed in the zebrafish retina, rh1-2, which we identified as a functional photosensitive pigment that binds 11-cis retinal and activates in response to light. Here, we localized expression of rh1-2 in the zebrafish retina to a subset of peripheral photoreceptor cells, which indicates a partially overlapping expression pattern with rh1 We also expressed, purified and characterized Rh1-2, including investigation of the stability of the biologically active intermediate. Using fluorescence spectroscopy, we found the half-life of the rate of retinal release of Rh1-2 following photoactivation to be more similar to that of the visual pigment rhodopsin than to the non-visual pigment exo-rhodopsin (exorh), which releases retinal around 5 times faster. Phylogenetic and molecular evolutionary analyses show that rh1-2 has ancient origins within teleost fishes, is under similar selective pressure to rh1, and likely experienced a burst of positive selection following its duplication and divergence from rh1 These findings indicate that rh1-2 is another functional visual rhodopsin gene, which contradicts the prevailing notion that visual rhodopsin is primarily found as a single copy gene within ray-finned fishes. The reasons for retention of this duplicate gene, as well as possible functional consequences for the visual system, are discussed.

Use of network model to explore dynamic and allosteric properties of three GPCR homodimers

We used an elastic network model and protein structure network to study three class A GPCR homodimers.

Computational approaches for modeling GPCR dimerization

Growing experimental evidences suggest that dimerization and oligomerization are important for G Protein- Coupled Receptors (GPCRs) function. The detailed structural information of dimeric/oligomeric GPCRs would be very important to understand their function. Although it is encouraging that recently several experimental GPCR structures in oligomeric forms have appeared, experimental determination of GPCR structures in oligomeric forms is still a big challenge, especially in mimicking the membrane environment. Therefore, development of computational approaches to predict dimerization of GPCRs will be highly valuable. In this review, we summarize computational approaches that have been developed and used for modeling of GPCR dimerization. In addition, we introduce a novel two-dimensional Brownian Dynamics based protein docking approach, which we have recently adapted, for GPCR dimer prediction.

Human infrared vision is triggered by two-photon chromophore isomerization

Vision relies on photoactivation of visual pigments in rod and cone photoreceptor cells of the retina. The human eye structure and the absorption spectra of pigments limit our visual perception of light. Our visual perception is most responsive to stimulating light in the 400- to 720-nm (visible) range. First, we demonstrate by psychophysical experiments that humans can perceive infrared laser emission as visible light. Moreover, we show that mammalian photoreceptors can be directly activated by near infrared light with a sensitivity that paradoxically increases at wavelengths above 900 nm, and display quadratic dependence on laser power, indicating a nonlinear optical process. Biochemical experiments with rhodopsin, cone visual pigments, and a chromophore model compound 11-cis-retinyl-propylamine Schiff base demonstrate the direct isomerization of visual chromophore by a two-photon chromophore isomerization. Indeed, quantum mechanics modeling indicates the feasibility of this mechanism. Together, these findings clearly show that human visual perception of near infrared light occurs by two-photon isomerization of visual pigments.

The CHARMM-TURBOMOLE interface for efficient and accurate QM/MM molecular dynamics, free energies, and excited state properties

The quantum mechanical (QM)/molecular mechanical (MM) interface between Chemistry at HARvard Molecular Mechanics (CHARMM) and TURBOMOLE is described. CHARMM provides an extensive set of simulation algorithms, like molecular dynamics (MD) and free energy perturbation, and support for mature nonpolarizable and Drude polarizable force fields. TURBOMOLE provides fast QM calculations using density functional theory or wave function methods and excited state properties. CHARMM-TURBOMOLE is well-suited for extended QM/MM MD simulations using first principles methods with large (triple-ζ) basis sets. We demonstrate these capabilities with a QM/MM simulation of Mg(2+) (aq), where the MM outer sphere water molecules are represented using the SWM4-NDP Drude polarizable force field and the ion and inner coordination sphere are represented using QM PBE, PBE0, and MP2 methods. The relative solvation free energies of Mg(2+) and Zn(2+) were calculated using thermodynamic integration. We also demonstrate the features for excited state properties. We calculate the time-averaged solution absorption spectrum of indole, the emission spectrum of the indole 1La excited state, and the electronic circular dichroism spectrum of an oxacepham.

Generalized QM/MM Force Matching Approach Applied to the 11-cis Protonated Schiff Base Chromophore of Rhodopsin

We extended a previously developed force matching approach to systems with covalent QM/MM boundaries and describe its user-friendly implementation in the publicly available software package CPMD. We applied this approach to the challenging case of the retinal protonated Schiff base in dark state bovine rhodopsin. We were able to develop a highly accurate force field that is able to capture subtle structural changes within the chromophore that have a pronounced influence on the optical properties. The optical absorption spectrum calculated from configurations extracted from a MD trajectory using the new force field is in excellent agreement with QM/MM and experimental references.

GPCR: G protein complexes--the fundamental signaling assembly

G protein coupled receptors (GPCR) constitute the largest group of cell surface receptors that transmit various signals across biological membranes through the binding and activation of heterotrimeric G proteins, which amplify the signal and activate downstream effectors leading to the biological responses. Thus, the first critical step in this signaling cascade is the interaction between receptor and its cognate G protein. Understanding this critical event at the molecular level is of high importance because abnormal function of GPCRs is associated with many diseases. Thus, these receptors are targets for drug development.

Retinal conformation governs pKa of protonated Schiff base in rhodopsin activation

We have explored the relationship between conformational energetics and the protonation state of the Schiff base in retinal, the covalently bound ligand responsible for activating the G protein-coupled receptor rhodopsin, using quantum chemical calculations. Guided by experimental structural determinations and large-scale molecular simulations on this system, we examined rotation about each bond in the retinal polyene chain, for both the protonated and deprotonated states that represent the dark and photoactivated states, respectively. Particular attention was paid to the torsional degrees of freedom that determine the shape of the molecule, and hence its interactions with the protein binding pocket. While most torsional degrees of freedom in retinal are characterized by large energetic barriers that minimize structural fluctuations under physiological temperatures, the C6-C7 dihedral defining the relative orientation of the β-ionone ring to the polyene chain has both modest barrier heights and a torsional energy surface that changes dramatically with protonation of the Schiff base. This surprising coupling between conformational degrees of freedom and protonation state is further quantified by calculations of the pKa as a function of the C6-C7 dihedral angle. Notably, pKa shifts of greater than two units arise from torsional fluctuations observed in molecular dynamics simulations of the full ligand-protein-membrane system. It follows that fluctuations in the protonation state of the Schiff base occur prior to forming the activated MII state. These new results shed light on important mechanistic aspects of retinal conformational changes that are involved in the activation of rhodopsin in the visual process.

Cis-retinoids and the chemistry of vision

We discuss here principal biochemical transformations of retinoid molecules in the visual cycle. We focus our analysis on the accumulating evidence of alternate pathways and functional redundancies in the cycle. The efficiency of the visual cycle depends, on one hand, on fast regeneration of the photo-bleached chromophores. On the other hand, it is crucial that the cyclic process should be highly selective to avoid accumulation of byproducts. The state-of-the-art knowledge indicates that single enzymatically active components of the cycle are not strictly selective and may require chaperones to enhance their rates. It appears that protein-protein interactions significantly improve the biological stability of the visual cycle. In particular, synthesis of thermodynamically less stable 11-cis-retinoid conformers is favored by physical interactions of the isomerases present in the retina with cellular retinaldehyde binding protein.

From atomic structures to neuronal functions of g protein-coupled receptors

G protein-coupled receptors (GPCRs) are essential mediators of signal transduction, neurotransmission, ion channel regulation, and other cellular events. GPCRs are activated by diverse stimuli, including light, enzymatic processing of their N-termini, and binding of proteins, peptides, or small molecules such as neurotransmitters. GPCR dysfunction caused by receptor mutations and environmental challenges contributes to many neurological diseases. Moreover, modern genetic technology has helped identify a rich array of mono- and multigenic defects in humans and animal models that connect such receptor dysfunction with disease affecting neuronal function. The visual system is especially suited to investigate GPCR structure and function because advanced imaging techniques permit structural studies of photoreceptor neurons at both macro and molecular levels that, together with biochemical and physiological assessment in animal models, provide a more complete understanding of GPCR signaling.

Rhodopsin Absorption from First Principles: Bypassing Common Pitfalls

Bovine rhodopsin is the most extensively studied retinal protein and is considered the prototype of this important class of photosensitive biosystems involved in the process of vision. Many theoretical investigations have attempted to elucidate the role of the protein matrix in modulating the absorption of retinal chromophore in rhodopsin, but, while generally agreeing in predicting the correct location of the absorption maximum, they often reached contradicting conclusions on how the environment tunes the spectrum. To address this controversial issue, we combine here a thorough structural and dynamical characterization of rhodopsin with a careful validation of its excited-state properties via the use of a wide range of state-of-the-art quantum chemical approaches including various flavors of time-dependent density functional theory (TDDFT), different multireference perturbative schemes (CASPT2 and NEVPT2), and quantum Monte Carlo (QMC) methods. Through extensive quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations, we obtain a comprehensive structural description of the chromophore-protein system and sample a wide range of thermally accessible configurations. We show that, in order to obtain reliable excitation properties, it is crucial to employ a sufficient number of representative configurations of the system. In fact, the common use of a single, ad hoc structure can easily lead to an incorrect model and an agreement with experimental absorption spectra due to cancelation of errors. Finally, we show that, to properly account for polarization effects on the chromophore and to quench the large blue-shift induced by the counterion on the excitation energies, it is necessary to adopt an enhanced description of the protein environment as given by a large quantum region including as many as 250 atoms.

Wagging the Tail: Essential Role of Substrate Flexibility in FAAH Catalysis

The serine hydrolase, fatty acid amide hydrolase (FAAH), is responsible for the intracellular degradation of anandamide and other bioactive fatty acid ethanolamides involved in the regulation of pain, inflammation, and other pathophysiological processes. The catalytic site of FAAH is composed of multiple cavities with mixed hydrophobic and hydrophilic properties, the role of which remains incompletely understood. Anandamide is thought to enter the active site through a "membrane-access" (MA) channel and position its flexible fatty acyl chain in a highly hydrophobic "acyl chain-binding" (AB) cavity to allow for hydrolysis to occur. Using microsecond molecular dynamics (MD) simulations of FAAH embedded in a realistic membrane/water environment, we show now that anandamide may not lock itself into the AB cavity but may rather assume catalytically significant conformations required for hydrolysis by moving its flexible arachidonoyl tail between the MA and AB cavities. This process is regulated by a phenylalanine residue (Phe432) located at the boundary between the two cavities, which may act as a "dynamic paddle." The results identify structural flexibility as a key determinant by which FAAH recognizes its primary lipid substrate.

Asymmetry of the rhodopsin dimer in complex with transducin

A large body of evidence for G-protein-coupled receptor (GPCR) oligomerization has accumulated over the past 2 decades. The smallest of these oligomers in vivo most likely is a dimer that buries 1000-Å(2) intramolecular surfaces and on stimulation forms a complex with heterotrimeric G protein in 2:1 stoichiometry. However, it is unclear whether each of the monomers adopts the same or a different conformation and function after activation of this dimer. With bovine rhodopsin (Rho) and its cognate bovine G-protein transducin (Gt) as a model system, we used the retinoid chromophores 11-cis-retinal and 9-cis-retinal to monitor each monomer of the dimeric GPCR within a stable complex with nucleotide-free Gt. We found that only 50% of Rho* in the Rho*-Gt complex is trapped in a Meta II conformation, while 50% evolves toward an opsin conformation and can be regenerated with 9-cis-retinal. We also found that all-trans-retinal can regenerate chromophore-depleted Rho*e complexed with Gt and FAK*TSA peptide containing Lys(296) with the attached all-trans retinoid (m/z of 934.5[MH](+)) was identified by mass spectrometry. Thus, our study shows that each of the monomers contributes unequally to the pentameric (2:1:1:1) complex of Rho dimer and Gt heterotrimer, validating the oligomeric structure of the complex and the asymmetry of the GPCR dimer, and revealing its structural/functional signature. This study provides a clear functional distinction between monomers of family A GPCRs in their oligomeric form.

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.

Molecular simulations and solid-state NMR investigate dynamical structure in rhodopsin activation

Rhodopsin has served as the primary model for studying G protein-coupled receptors (GPCRs)-the largest group in the human genome, and consequently a primary target for pharmaceutical development. Understanding the functions and activation mechanisms of GPCRs has proven to be extraordinarily difficult, as they are part of a complex signaling cascade and reside within the cell membrane. Although X-ray crystallography has recently solved several GPCR structures that may resemble the activated conformation, the dynamics and mechanism of rhodopsin activation continue to remain elusive. Notably solid-state ((2))H NMR spectroscopy provides key information pertinent to how local dynamics of the retinal ligand change during rhodopsin activation. When combined with molecular mechanics simulations of proteolipid membranes, a new paradigm for the rhodopsin activation process emerges. Experiment and simulation both suggest that retinal isomerization initiates the rhodopsin photocascade to yield not a single activated structure, but rather an ensemble of activated conformational states. This article is part of a Special Issue entitled: Membrane protein structure and function.

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.

Protein dynamics and enzyme catalysis: insights from simulations

The role of protein dynamics in enzyme catalysis is one of the most active and controversial areas in enzymology today. Some researchers claim that protein dynamics are at the heart of enzyme catalytic efficiency, while others state that dynamics make no significant contribution to catalysis. What is the biochemist - or student - to make of the ferocious arguments in this area? Protein dynamics are complex and fascinating, as molecular dynamics simulations and experiments have shown. The essential question is: do these complex motions have functional significance? In particular, how do they affect or relate to chemical reactions within enzymes, and how are chemical and conformational changes coupled together? Biomolecular simulations can analyse enzyme reactions and dynamics in atomic detail, beyond that achievable in experiments: accurate atomistic modelling has an essential part to play in clarifying these issues. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Trans-activation between 7TM domains: implication in heterodimeric GABAB receptor activation

Seven-transmembrane domain (7TM) receptors have important functions in cell-cell communication and can assemble into dimers or oligomers. Such complexes may allow specific functional cross-talk through trans-activation of interacting 7TMs, but this hypothesis requires further validation. Herein, we used the GABAB receptor, which is composed of two distinct subunits, GABAB1, which binds the agonist, and GABAB2, which activates G proteins, as a model system. By using a novel orthogonal-labelling approach compatible with time-resolved FRET and based on ACP- and SNAP-tag technologies to verify the heterodimerization of wild-type and mutated GABAB subunits, we demonstrate the existence of a direct allosteric coupling between the 7TMs of GABAB heterodimers. Indeed, a GABAB receptor, in which the GABAB2 extracellular domain was deleted, was still capable of activating G proteins. Furthermore, synthetic ligands for the GABAB2 7TM could increase agonist affinity at the GABAB1 subunit in this mutated receptor. In addition to bringing new information on GABAB receptor activation, these data clearly demonstrate the existence of direct trans-activation between the 7TM of two interacting proteins.