Effect of intracellular loop 3 on intrinsic dynamics of human β2-adrenergic receptor

Elucidation of a dynamic interplay between a beta-2 adrenergic receptor, its agonist, and stimulatory G protein

G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (β2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to β2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from β2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-β2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and β2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of β2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.

Stepwise Frontal Analysis Coupled with Affinity Chromatography: A Fast and Reliable Method for Potential Ligand Isolation and Evaluation from Mahuang-Fuzi-Xixin Decoction

Mahuang-Fuzi-Xixin Decoction (MFXD) is widely used in the treatment of asthma, however, the functional components in the decoction targeting beta2-adrenoceptor (β₂ -AR) remain unclear. Herein, we immobilized the haloalkane dehalogenase (Halo)-tagged β₂ -AR on the 6-chlorocaproic acid-modified microspheres. Using the affinity stationary phase, the interactions of four ligands with the receptor were analyzed by stepwise frontal analysis. The association constants were (4.75±0.28)×10⁴  M^-1 for salbutamol, (2.93±0.15)×10⁴  M^-1 for terbutaline, (1.23±0.03)×10⁴  M^-1 for methoxyphenamine, (5.67±0.38)×10⁴  M^-1 for clorprenaline at high-affinity binding site, and (2.73±0.05)×10³  M^-1 at low-affinity binding site. These association constants showed the same rank order as the radioligand binding assay, demonstrating that immobilized β₂ -AR had capacity to screen bioactive compounds binding to the receptor while stepwise frontal analysis could predict their binding affinities. Application of the immobilized receptor in analysis of MFXD by chromatographic method revealed that ephedrine, aconifine, karakoline, and chasmanine were the bioactive compounds targeting β₂ -AR. Among them, ephedrine and chasmanine exhibited association constants of (2.94±0.02)×10⁴ M^-1 and (4.60±0.15)×10⁴  M^-1 to the receptor by stepwise frontal analysis. Molecular docking analysis demonstrated that ephedrine, chasmanine, and the other two compounds interact with β₂ -AR through the same pocket involving the key amino acids such as Asn312, Asp113, Phe289, Trp286, Tyr316, and Val114. As such, we reasoned that the four compounds dominate the therapeutic effect of MFXD against asthma through β₂ -AR mediating pathway. This work shed light on the potential of immobilized β₂ -AR for drug discovery and provided a valuable methodology for rapid screening.

Autoregulation of GPCR signalling through the third intracellular loop

The third intracellular loop (ICL3) of the G protein-coupled receptor (GPCR) fold is important for the signal transduction process downstream of receptor activation1-3. Despite this, the lack of a defined structure of ICL3, combined with its high sequence divergence among GPCRs, complicates characterization of its involvement in receptor signalling4. Previous studies focusing on the β2 adrenergic receptor (β2AR) suggest that ICL3 is involved in the structural process of receptor activation and signalling5-7. Here we derive mechanistic insights into the role of ICL3 in β2AR signalling, observing that ICL3 autoregulates receptor activity through a dynamic conformational equilibrium between states that block or expose the receptor's G protein-binding site. We demonstrate the importance of this equilibrium for receptor pharmacology, showing that G protein-mimetic effectors bias the exposed states of ICL3 to allosterically activate the receptor. Our findings additionally reveal that ICL3 tunes signalling specificity by inhibiting receptor coupling to G protein subtypes that weakly couple to the receptor. Despite the sequence diversity of ICL3, we demonstrate that this negative G protein-selection mechanism through ICL3 extends to GPCRs across the superfamily, expanding the range of known mechanisms by which receptors mediate G protein subtype selective signalling. Furthermore, our collective findings suggest ICL3 as an allosteric site for receptor- and signalling pathway-specific ligands.

Activation and Speciation Mechanisms in Class A GPCRs

Accurate development of allosteric modulators of GPCRs require a thorough assessment of their sequence, structure, and dynamics, toward gaining insights into their mechanisms of actions shared by family members, as well as dynamic features that distinguish subfamilies. Building on recent progress in the characterization of the signature dynamics of proteins, we analyzed here a dataset of 160 Class A GPCRs to determine their sequence similarities, structural landscape, and dynamic features across different species (human, bovine, mouse, squid, and rat), different activation states (active/inactive), and different subfamilies. The two dominant directions of variability across experimentally resolved structures, identified by principal component analysis of the dataset, shed light to cooperative mechanisms of activation, subfamily differentiation, and speciation of Class A GPCRs. The analysis reveals the functional significance of the conformational flexibilities of specific structural elements, including: the dominant role of the intracellular loop 3 (ICL3) together with the cytoplasmic ends of the adjoining helices TM5 and TM6 in enabling allosteric activation; the role of particular structural motifs at the extracellular loop 2 (ECL2) connecting TM4 and TM5 in binding ligands specific to different subfamilies; or even the differentiation of the N-terminal conformation across different species. Detailed analyses of the modes of motions accessible to the members of the dataset and their variations across members demonstrate how the active and inactive states of GPCRs obey distinct conformational dynamics. The collective fluctuations of the GPCRs are robustly defined in the active state, while the inactive conformers exhibit broad variance among members.

Lipids and Phosphorylation Conjointly Modulate Complex Formation of β2-Adrenergic Receptor and β-arrestin2

G protein-coupled receptors (GPCRs) are the largest class of human membrane proteins that bind extracellular ligands at their orthosteric binding pocket to transmit signals to the cell interior. Ligand binding evokes conformational changes in GPCRs that trigger the binding of intracellular interaction partners (G proteins, G protein kinases, and arrestins), which initiate diverse cellular responses. It has become increasingly evident that the preference of a GPCR for a certain intracellular interaction partner is modulated by a diverse range of factors, e.g., ligands or lipids embedding the transmembrane receptor. Here, by means of molecular dynamics simulations of the β₂-adrenergic receptor and β-arrestin2, we study how membrane lipids and receptor phosphorylation regulate GPCR-arrestin complex conformation and dynamics. We find that phosphorylation drives the receptor’s intracellular loop 3 (ICL3) away from a native negatively charged membrane surface to interact with arrestin. If the receptor is embedded in a neutral membrane, the phosphorylated ICL3 attaches to the membrane surface, which widely opens the receptor core. This opening, which is similar to the opening in the G protein-bound state, weakens the binding of arrestin. The loss of binding specificity is manifested by shallower arrestin insertion into the receptor core and higher dynamics of the receptor-arrestin complex. Our results show that receptor phosphorylation and the local membrane composition cooperatively fine-tune GPCR-mediated signal transduction. Moreover, the results suggest that deeper understanding of complex GPCR regulation mechanisms is necessary to discover novel pathways of pharmacological intervention.

DEER Analysis of GPCR Conformational Heterogeneity

G protein-coupled receptors (GPCRs) represent a large class of transmembrane helical proteins which are involved in numerous physiological signaling pathways and therefore represent crucial pharmacological targets. GPCR function and the action of therapeutic molecules are defined by only a few parameters, including receptor basal activity, ligand affinity, intrinsic efficacy and signal bias. These parameters are encoded in characteristic receptor conformations existing in equilibrium and their populations, which are thus of paramount interest for the understanding of receptor (mal-)functions and rational design of improved therapeutics. To this end, the combination of site-directed spin labeling and EPR spectroscopy, in particular double electron-electron resonance (DEER), is exceedingly valuable as it has access to sub-Angstrom spatial resolution and provides a detailed picture of the number and populations of conformations in equilibrium. This review gives an overview of existing DEER studies on GPCRs with a focus on the delineation of structure/function frameworks, highlighting recent developments in data analysis and visualization. We introduce "conformational efficacy" as a parameter to describe ligand-specific shifts in the conformational equilibrium, taking into account the loose coupling between receptor segments observed for different GPCRs using DEER.

Structural and functional characterization of G protein-coupled receptors with deep mutational scanning

The >800 human G protein–coupled receptors (GPCRs) are responsible for transducing diverse chemical stimuli to alter cell state- and are the largest class of drug targets. Their myriad structural conformations and various modes of signaling make it challenging to understand their structure and function. Here, we developed a platform to characterize large libraries of GPCR variants in human cell lines with a barcoded transcriptional reporter of G protein signal transduction. We tested 7800 of 7828 possible single amino acid substitutions to the beta-2 adrenergic receptor (β₂AR) at four concentrations of the agonist isoproterenol. We identified residues specifically important for β₂AR signaling, mutations in the human population that are potentially loss of function, and residues that modulate basal activity. Using unsupervised learning, we identify residues critical for signaling, including all major structural motifs and molecular interfaces. We also find a previously uncharacterized structural latch spanning the first two extracellular loops that is highly conserved across Class A GPCRs and is conformationally rigid in both the inactive and active states of the receptor. More broadly, by linking deep mutational scanning with engineered transcriptional reporters, we establish a generalizable method for exploring pharmacogenomics, structure and function across broad classes of drug receptors.

How Do Molecular Dynamics Data Complement Static Structural Data of GPCRs

G protein-coupled receptors (GPCRs) are implicated in nearly every physiological process in the human body and therefore represent an important drug targeting class. Advances in X-ray crystallography and cryo-electron microscopy (cryo-EM) have provided multiple static structures of GPCRs in complex with various signaling partners. However, GPCR functionality is largely determined by their flexibility and ability to transition between distinct structural conformations. Due to this dynamic nature, a static snapshot does not fully explain the complexity of GPCR signal transduction. Molecular dynamics (MD) simulations offer the opportunity to simulate the structural motions of biological processes at atomic resolution. Thus, this technique can incorporate the missing information on protein flexibility into experimentally solved structures. Here, we review the contribution of MD simulations to complement static structural data and to improve our understanding of GPCR physiology and pharmacology, as well as the challenges that still need to be overcome to reach the full potential of this technique.

Distinctive communication networks in inactive states of β2 -adrenergic receptor: Mutual information and entropy transfer analysis

Mutual information and entropy transfer analysis employed on two inactive states of human beta-2 adrenergic receptor (β₂ -AR) unraveled distinct communication pathways. Previously, a so-called "highly" inactive state of the receptor was observed during 1.5 microsecond long molecular dynamics simulation where the largest intracellular loop (ICL3) was swiftly packed onto the G-protein binding cavity, becoming entirely inaccessible. Mutual information quantifying the degree of correspondence between backbone-C_α fluctuations was mostly shared between intra- and extra-cellular loop regions in the original inactive state, but shifted to entirely different regions in this latest inactive state. Interestingly, the largest amount of mutual information was always shared among the mobile regions. Irrespective of the conformational state, polar residues always contributed more to mutual information than hydrophobic residues, and also the number of polar-polar residue pairs shared the highest degree of mutual information compared to those incorporating hydrophobic residues. Entropy transfer, quantifying the correspondence between backbone-C_α fluctuations at different timesteps, revealed a distinctive pathway directed from the extracellular site toward intracellular portions in this recently exposed inactive state for which the direction of information flow was the reverse of that observed in the original inactive state where the mobile ICL3 and its intracellular surroundings drove the future fluctuations of extracellular regions.

Intrinsic Dynamics and Causality in Correlated Motions Unraveled in Two Distinct Inactive States of Human β2-Adrenergic Receptor

The alternative inactive state of the human β₂-adrenergic receptor originally exposed in molecular dynamics simulations was investigated using various analysis tools to evaluate causality between correlated residue-pair fluctuations and suggest allosteric communication pathways. A major conformational shift observed in the third intracellular loop (ICL3) displayed a novel inactive state, featuring an inaccessible G protein binding site blocked by ICL3 and an expanded orthosteric ligand binding site. Residue-based mean-square fluctuation and stiffness calculations revealed a significant mobility decrease in ICL3, which induced a mobility increase in the remaining loop regions. This indicates conformational entropy loss in one mobile region being compensated by residual intermolecular motions in other mobile regions. Moreover, the extent of significantly correlated motions decreased, and correlations that once existed between transmembrane helices shifted toward regions with increased mobility. Conditional time-delayed cross-correlation analysis identified distinct driver-follower relationship profiles. Prior to its packing, freely moving ICL3 was markedly driven by transmembrane helix-8 whereas once packed, ICL3 controlled future fluctuations of nearby helices. Moreover, two transmembrane helices, (H5 and H6), started to control future fluctuations of a remote site, the extracellular loop, ECL2. This clearly suggests that allosteric coupling between extra- and intracellular parts intensified, in agreement with the receptor's well recognized feature, which is the inverse proportionality between activity and the degree of coupling.

Ligand-binding affinity of alternative conformers of human β2 -adrenergic receptor in the presence of intracellular loop 3 (ICL3) and their potential use in virtual screening studies

This study investigates the structural distinctiveness of orthosteric ligand-binding sites of several human β₂ adrenergic receptor (β₂ -AR) conformations that have been obtained from a set of independent molecular dynamics (MD) simulations in the presence of intracellular loop 3 (ICL3). A docking protocol was established in order to classify each receptor conformation via its binding affinity to selected ligands with known efficacy. This work's main goal was to reveal many subtle features of the ligand-binding site, presenting alternative conformations, which might be considered as either active- or inactive-like but mostly specific for that ligand. Agonists, inverse agonists, and antagonists were docked to each MD conformer with distinct binding pockets, using different docking tools and scoring functions. Mostly favored receptor conformation persistently observed in all docking/scoring evaluations was classified as active or inactive based on the type of ligand's biological effect. Classified MD conformers were further tested for their ability to discriminate agonists from inverse agonists/antagonists, and several conformers were proposed as important targets to be used in virtual screening experiments that were often limited to a single X-ray structure.

Exploring the Ligand Efficacy of Cannabinoid Receptor 1 (CB1) using Molecular Dynamics Simulations

Cannabinoid receptor 1 (CB1) is a promising therapeutic target for a variety of disorders. Distinct efficacy profiles showed different therapeutic effects on CB1 dependent on three classes of ligands: agonists, antagonists, and inverse agonists. To discriminate the distinct efficacy profiles of the ligands, we carried out molecular dynamics (MD) simulations to identify the dynamic behaviors of inactive and active conformations of CB1 structures with the ligands. In addition, the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method was applied to analyze the binding free energy decompositions of the CB1-ligand complexes. With these two methods, we found the possibility that the three classes of ligands can be discriminated. Our findings shed light on the understanding of different efficacy profiles of ligands by analyzing the structural behaviors of intact CB1 structures and the binding energies of ligands, thereby yielding insights that are useful for the design of new potent CB1 drugs.

Constrained dynamics of the sole tryptophan in the third intracellular loop of the serotonin1A receptor

G protein-coupled receptors (GPCRs) are major signaling proteins in eukaryotic cells and are important drug targets. In spite of their role in GPCR function, the extramembranous regions of GPCRs are relatively less appreciated. The third intracellular loop (ICL3), which connects transmembrane helices V and VI, is important in this context since its crucial role in signaling has been documented for a number of GPCRs. Unfortunately, the structure of this loop is generally not visualized in x-ray crystallographic studies since this flexible loop is either stabilized using a monoclonal antibody or replaced with lysozyme. In this work, we expressed and purified the ICL3 region of the serotonin_1A receptor and monitored its motional restriction and organization utilizing red edge excitation shift (REES) of its sole tryptophan and circular dichroism (CD) spectroscopy. Our results show that the tryptophan in ICL3 exhibits REES of 4 nm, implying that it is localized in a restricted microenvironment. These results are further supported by wavelength-selective changes in fluorescence anisotropy and lifetime. This constrained dynamics was relaxed upon denaturation of the peptide, thereby suggesting the involvement of the peptide secondary structure in the observed motional restriction, as evident from CD spectroscopy and apparent rotational correlation time. To the best of our knowledge, these results constitute one of the first measurements of motional constraint in the ICL3 region of GPCRs. Our results are relevant in the context of the reported intrinsically disordered nature of ICL3 and its role in providing functional diversity to GPCRs due to conformational plasticity.

Structural insights into positive and negative allosteric regulation of a G protein-coupled receptor through protein-lipid interactions

Lipids are becoming known as essential allosteric modulators of G protein-coupled receptor (GPCRs). However, how they exert their effects on GPCR conformation at the atomic level is still unclear. In light of recent experimental data, we have performed several long-timescale molecular dynamics (MD) simulations, totalling 24 μs, to rigorously map allosteric modulation and conformational changes in the β₂ adrenergic receptor (β2AR) that occur as a result of interactions with three different phospholipids. In particular, we identify different sequential mechanisms behind receptor activation and deactivation, respectively, mediated by specific lipid interactions with key receptor regions. We show that net negatively charged lipids stabilize an active-like state of β2AR that is able to dock G_sα protein. Clustering of anionic lipids around the receptor with local distortion of membrane thickness is also apparent. On the other hand, net-neutral zwitterionic lipids inactivate the receptor, generating either fully inactive or intermediate states, with kinetics depending on lipid headgroup charge distribution and hydrophobicity. These chemical differences alter membrane thickness and density, which differentially destabilize the β2AR active state through lateral compression effects.

In Silico Studies Targeting G-protein Coupled Receptors for Drug Research Against Parkinson's Disease

Parkinson's Disease (PD) is a long-term neurodegenerative brain disorder that mainly affects the motor system. The causes are still unknown, and even though currently there is no cure, several therapeutic options are available to manage its symptoms. The development of novel antiparkinsonian agents and an understanding of their proper and optimal use are, indeed, highly demanding. For the last decades, L-3,4-DihydrOxyPhenylAlanine or levodopa (L-DOPA) has been the gold-standard therapy for the symptomatic treatment of motor dysfunctions associated to PD. However, the development of dyskinesias and motor fluctuations (wearing-off and on-off phenomena) associated with long-term L-DOPA replacement therapy have limited its antiparkinsonian efficacy. The investigation for non-dopaminergic therapies has been largely explored as an attempt to counteract the motor side effects associated with dopamine replacement therapy. Being one of the largest cell membrane protein families, G-Protein-Coupled Receptors (GPCRs) have become a relevant target for drug discovery focused on a wide range of therapeutic areas, including Central Nervous System (CNS) diseases. The modulation of specific GPCRs potentially implicated in PD, excluding dopamine receptors, may provide promising non-dopaminergic therapeutic alternatives for symptomatic treatment of PD. In this review, we focused on the impact of specific GPCR subclasses, including dopamine receptors, adenosine receptors, muscarinic acetylcholine receptors, metabotropic glutamate receptors, and 5-hydroxytryptamine receptors, on the pathophysiology of PD and the importance of structure- and ligand-based in silico approaches for the development of small molecules to target these receptors.

Mechanism of allosteric regulation of β2-adrenergic receptor by cholesterol

There is evidence that lipids can be allosteric regulators of membrane protein structure and activation. However, there are no data showing how exactly the regulation emerges from specific lipid-protein interactions. Here we show in atomistic detail how the human β₂-adrenergic receptor (β₂AR) - a prototypical G protein-coupled receptor - is modulated by cholesterol in an allosteric fashion. Extensive atomistic simulations show that cholesterol regulates β₂AR by limiting its conformational variability. The mechanism of action is based on the binding of cholesterol at specific high-affinity sites located near the transmembrane helices 5-7 of the receptor. The alternative mechanism, where the β₂AR conformation would be modulated by membrane-mediated interactions, plays only a minor role. Cholesterol analogues also bind to cholesterol binding sites and impede the structural flexibility of β₂AR, however cholesterol generates the strongest effect. The results highlight the capacity of lipids to regulate the conformation of membrane receptors through specific interactions.

Can Specific Protein-Lipid Interactions Stabilize an Active State of the Beta 2 Adrenergic Receptor?

G-protein-coupled receptors are eukaryotic membrane proteins with broad biological and pharmacological relevance. Like all membrane-embedded proteins, their location and orientation are influenced by lipids, which can also impact protein function via specific interactions. Extensive simulations totaling 0.25 ms reveal a process in which phospholipids from the membrane's cytosolic leaflet enter the empty G-protein binding site of an activated β2 adrenergic receptor and form salt-bridge interactions that inhibit ionic lock formation and prolong active-state residency. Simulations of the receptor embedded in an anionic membrane show increased lipid binding, providing a molecular mechanism for the experimental observation that anionic lipids can enhance receptor activity. Conservation of the arginine component of the ionic lock among Rhodopsin-like G-protein-coupled receptors suggests that intracellular lipid ingression between receptor helices H6 and H7 may be a general mechanism for active-state stabilization.

Investigation of allosteric coupling in human β2-adrenergic receptor in the presence of intracellular loop 3

Background

This study investigates the allosteric coupling that exists between the intra- and extracellular parts of human β₂-adrenergic receptor (β₂-AR), in the presence of the intracellular loop 3 (ICL3), which is missing in all crystallographic experiments and most of the simulation studies reported so far. Our recent 1 μs long MD run has revealed a transition to the so-called very inactive state of the receptor, in which ICL3 packed under the G protein’s binding cavity and completely blocked its accessibility to G protein. Simultaneously, an outward tilt of transmembrane helix 5 (TM5) caused an expansion of the extracellular ligand-binding site. In the current study, we performed independent runs with a total duration of 4 μs to further investigate the very inactive state with packed ICL3 and the allosteric coupling event (three unrestrained runs and five runs with bond restraints at the ligand-binding site).

Results

In all three independent unrestrained runs (each 500 ns long), ICL3 preserved its initially packed/closed conformation within the studied time frame, suggesting an inhibition of the receptor’s activity. Specific bond restraints were later imposed between some key residues at the ligand-binding site, which have been experimentally determined to interact with the ligand. Restraining the binding site region to an open state facilitated ICL3 closure, whereas a relatively constrained/closed binding site hindered ICL3 packing. However, the reverse operation, i.e. opening of the packed ICL3, could not be realized by restraining the binding site region to a closed state. Thus, any attempt failed to free the ICL3 from its locked state due to the presence of persistent hydrogen bonds.

Conclusions

Overall, our simulations indicated that starting with very inactive states, the receptor stayed almost irreversibly inhibited, which in turn decreased the overall mobility of the receptor. Bond restraints which represented the geometric restrictions caused by ligands of various sizes when bound at the ligand-binding site, induced the expected conformational changes in TM5, TM6 and consequently, ICL3. Still, once ICL3 was packed, the allosteric coupling became ineffective due to strong hydrogen bonds connecting ICL3 to the core of the receptor.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-016-0061-9) contains supplementary material, which is available to authorized users.

How To Minimize Artifacts in Atomistic Simulations of Membrane Proteins, Whose Crystal Structure Is Heavily Engineered: β₂-Adrenergic Receptor in the Spotlight

Atomistic molecular dynamics (MD) simulations are used extensively to elucidate membrane protein properties. These simulations are based on three-dimensional protein structures that in turn are often based on crystallography. The protein structures resolved in crystallographic studies typically do not correspond to pristine proteins, however. Instead the crystallized proteins are commonly engineered, including structural modifications (mutations, replacement of protein sequences by antibodies, bound ligands, etc.) whose impact on protein structure and dynamics is largely unknown. Here we explore this issue through atomistic MD simulations (∼5 μs in total), focusing on the β2-adrenergic receptor (β2AR) that is one of the most studied members of the G-protein coupled receptor superfamily. Starting from an inactive-state crystal structure of β2AR, we remove the many modifications in β2AR systematically one at a time, in six consecutive steps. After each step, we equilibrate the system and simulate it quite extensively. The results of this step-by-step approach highlight that the structural modifications used in crystallization can affect ligand and G-protein binding sites, packing at the transmembrane-helix interface region, and the dynamics of connecting loops in β2AR. When the results of the systematic step-by-step approach are compared to an all-at-once technique where all modifications done on β2AR are removed instantaneously at the same time, it turns out that the step-by-step method provides results that are superior in terms of maintaining protein structural stability. The results provide compelling evidence that for membrane proteins whose 3D structure is based on structural engineering, the preparation of protein structure for atomistic MD simulations is a delicate and sensitive process. The results show that most valid results are found when the structural modifications are reverted slowly, one at a time.

Structural insights and functional implications of inter-individual variability in β2-adrenergic receptor

The human β₂-adrenergic receptor (β₂AR) belongs to the G protein-coupled receptor (GPCR) family and due to its central role in bronchodilation, is an important drug target. The inter-individual variability in β₂AR has been implicated in disease susceptibility and differential drug response. In this work, we identified nine potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) using a consensus approach. The deleterious nsSNPs were found to cluster near the ligand binding site and towards the G-protein binding site. To assess their molecular level effects, we built structural models of these receptors and performed atomistic molecular dynamics simulations. Most notably, in the Phe290Ser variant we observed the rotameric flip of Trp286^6.48, a putative activation switch that has not been reported in β₂AR thus far. In contrast, the variant Met82Lys was found to be the most detrimental to epinephrine binding. Additionally, a few of the nsSNPs were seen to cause perturbations to the lipid bilayer, while a few lead to differences at the G-protein coupling site. We are thus able to classify the variants as ranging from activating to damaging, prioritising them for experimental studies.

Ligand Binding Pathways of Clozapine and Haloperidol in the Dopamine D2 and D3 Receptors

The binding of a small molecule ligand to its protein target is most often characterized by binding affinity and is typically viewed as an on/off switch. The more complex reality is that binding involves the ligand passing through a series of intermediate states between the solution phase and the fully bound pose. We have performed a set of 29 unbiased molecular dynamics simulations to model the binding pathways of the dopamine receptor antagonists clozapine and haloperidol binding to the D2 and D3 dopamine receptors. Through these simulations we have captured the binding pathways of clozapine and haloperidol from the extracellular vestibule to the orthosteric binding site and thereby, we also predict the bound pose of each ligand. These are the first long time scale simulations of haloperidol or clozapine binding to dopamine receptors. From these simulations, we have identified several important stages in the binding pathway, including the involvement of Tyr7.35 in a "handover" mechanism that transfers the ligand between the extracellular vestibule and Asp3.32. We have also performed interaction and cluster analyses to determine differences in binding pathways between the D2 and D3 receptors and identified metastable states that may be of use in drug design.