Computational modeling of human coreceptor CCR5 antagonist as a HIV-1 entry inhibitor: using an integrated homology modeling,docking, and membrane molecular dynamics simulation analysis approach

Chemokine receptor 5 (CCR5) is an integral membrane protein that is utilized during human immunodeficiency virus type-1 entry into host cells. CCR5 is a G-protein coupled receptor that contains seven transmembrane (TM) helices. However, the crystal structure of CCR5 has not been reported. A homology model of CCR5 was developed based on the recently reported CXCR4 structure as template. Automated docking of the most potent (14), medium potent (37), and least potent (25) CCR5 antagonists was performed using the CCR5 model. To characterize the mechanism responsible for the interactions between ligands (14, 25, and 37) and CCR5, membrane molecular dynamic (MD) simulations were performed. The position and orientation of ligands (14, 25, and 37) were found to be changed after MD simulations, which demonstrated the ability of this technique to identify binding modes. Furthermore, at the end of simulation, it was found that residues identified by docking were changed and some new residues were introduced in the proximity of ligands. Our results are in line with the majority of previous mutational reports. These results show that hydrophobicity is the determining factor of CCR5 antagonism. In addition, salt bridging and hydrogen bond contacts between ligands (14, 25, and 37) and CCR5 are also crucial for inhibitory activity. The residues newly identified by MD simulation are Ser160, Phe166, Ser180, His181, and Trp190, and so far no site-directed mutagenesis studies have been reported. To determine the contributions made by these residues, additional mutational studies are suggested. We propose a general binding mode for these derivatives based on the MD simulation results of higher (14), medium (37), and lower (25) potent inhibitors. Interestingly, we found some trend for these inhibitors such as, salt bridge interaction between basic nitrogen of ligand and acidic Glu283 seemed necessary for inhibitory activity. Also, two aromatic pockets (pocket I – TM1-3 and pocket II – TM3-6) were linked by the central polar region in TM7, and the simulated inhibitors show important interactions with the Trp86, Tyr89, Tyr108, Phe112, Ile198, Tyr251, Leu255, and Gln280 and Glu283 residues. These results shed light on the usage of MD simulation to identify more stable, optimal binding modes of the inhibitors.     

Molecular Mechanisms of Bitopic Ligand Engagement with the M1 Muscarinic Acetylcholine Receptor

TBPB and 77-LH-28-1 are selective agonists of the M₁ muscarinic acetylcholine receptor (mAChR) that may gain their selectivity through a bitopic mechanism, interacting concomitantly with the orthosteric site and part of an allosteric site. The current study combined site-directed mutagenesis, analytical pharmacology,and molecular modeling to gain further insights into the structural basis underlying binding and signaling by these agonists. Mutations within the orthosteric binding site caused similar reductions in affinity and signaling efficacy for both selective and prototypical orthosteric ligands. In contrast, the mutation of residues within transmembrane helix (TM) 2 and the second extracellular loop (ECL2) discriminated between the different classes of ligand. In particular, ECL2 appears to be involved in the selective binding of bitopic ligands and in coordinating biased agonism between intracellular calcium mobilization and ERK1/2 phosphorylation. Molecular modeling of the interaction between TBPB and the M₁ mAChR revealed a binding pose predicted to extend from the orthosteric site up toward a putative allosteric site bordered by TM2, TM3, and TM7, thus consistent with a bitopic mode of binding. Overall, these findings provide valuable structural and mechanistic insights into bitopic ligand actions and receptor activation and support a role for ECL2 in dictating the active states that can be adopted by a G protein-coupled receptor. This may enable greater selective ligand design and development for mAChRs and facilitate improved identification of bitopic ligands.     

Molecular Characterization and Expression of Cloned Human Galanin Receptors GALR2 and GALR3

Abstract: Galanin is a 29- or 30-amino acid peptide with wide-ranging effects on hormone release, feeding behavior, smooth muscle contractility, and somatosensory neuronal function. Three distinct galanin receptor (GALR) subtypes, designated GALR1, 2, and 3, have been cloned from the rat. We report here the cloning of the human GALR2 and GALR3 genes, an initial characterization of their pharmacology with respect to radioligand binding and signal transduction pathways, and a profile of their expression in brain and peripheral tissues. Human GALR2 and GALR3 show, respectively, 92 and 89% amino acid sequence identity with their rat homologues. Radioligand binding studies with ¹²⁵I-galanin show that recombinant human GALR2 binds with high affinity to human galanin (K_D = 0.3 nM). Human GALR3 binds galanin with less affinity (IC₅₀ of 12 nM for porcine galanin and 75 nM for human galanin). Human GALR2 was shown to couple to phospholipase C and elevation of intracellular calcium levels as assessed by aequorin luminescence in HEK-293 cells and by Xenopus melanophore pigment aggregation and dispersion assays, in contrast to human GALR1 and human GALR3, which signal predominantly through inhibition of adenylate cyclase. GALR2 mRNA shows a wide distribution in the brain (mammillary nuclei, dentate gyrus, cingulate gyrus, and posterior hypothalamic, supraoptic, and arcuate nuclei), and restricted peripheral tissue distribution with highest mRNA levels detected in human small intestine. In comparison, whereas GALR3 mRNA was expressed in many areas of the rat brain, there was abundant expression in the primary olfactory cortex, olfactory tubercle, the islands of Calleja, the hippocampal CA regions of Ammon's horn, and the dentate gyrus. GALR3 mRNA was highly expressed in human testis and was detectable in adrenal gland and pancreas. The genes for human GALR2 and 3 were localized to chromosomes 17q25 and 22q12.2–13.1, respectively.     

Mechanism of μ-Opioid Receptor-Magnesium Interaction and Positive Allosteric Modulation

In the era of opioid abuse epidemics, there is an increased demand for understanding how opioid receptors can be allosterically modulated to guide the development of more effective and safer opioid therapies. Among the modulators of the μ-opioid (MOP) receptor, which is the pharmacological target for the majority of clinically used opioid drugs, are monovalent and divalent cations. Specifically, the monovalent sodium cation (Na⁺) has been known for decades to affect MOP receptor signaling by reducing agonist binding, whereas the divalent magnesium cation (Mg²⁺) has been shown to have the opposite effect, notwithstanding the presence of sodium chloride. Although ultra-high-resolution opioid receptor crystal structures have revealed a specific Na⁺ binding site and molecular dynamics (MD) simulation studies have supported the idea that this monovalent ion reduces agonist binding by stabilizing the receptor inactive state, the putative binding site of Mg²⁺ on the MOP receptor, as well as the molecular determinants responsible for its positive allosteric modulation of the receptor, are unknown. In this work, we carried out tens of microseconds of all-atom MD simulations to investigate the simultaneous binding of Mg²⁺ and Na⁺ cations to inactive and active crystal structures of the MOP receptor embedded in an explicit lipid-water environment and confirmed adequate sampling of Mg²⁺ ion binding with a grand canonical Monte Carlo MD method. Analyses of these simulations shed light on 1) the preferred binding sites of Mg²⁺ on the MOP receptor, 2) details of the competition between Mg²⁺ and Na⁺ cations for specific sites, 3) estimates of binding affinities, and 4) testable hypotheses of the molecular mechanism underlying the positive allosteric modulation of the MOP receptor by the Mg²⁺ cation.     

Structure of the human galanin receptor 2 bound to galanin and Gq reveals the basis of ligand specificity and how binding affects the G-protein interface

Galanin is a neuropeptide expressed in the central and peripheral nervous systems, where it regulates various processes including neuroendocrine release, cognition, and nerve regeneration. Three G-protein coupled receptors (GPCRs) for galanin have been discovered, which is the focus of efforts to treat diseases including Alzheimer’s disease, anxiety, and addiction. To understand the basis of the ligand preferences of the receptors and to assist structure-based drug design, we used cryo-electron microscopy (cryo-EM) to solve the molecular structure of GALR2 bound to galanin and a cognate heterotrimeric G-protein, providing a molecular view of the neuropeptide binding site. Mutant proteins were assayed to help reveal the basis of ligand specificity, and structural comparison between the activated GALR2 and inactive hβ₂AR was used to relate galanin binding to the movements of transmembrane (TM) helices and the G-protein interface.

Galanin is a neuropeptide expressed in the central and peripheral nervous systems, where it regulates various processes including neuroendocrine release, cognition, and nerve regeneration. This cryo-electron microscopy study shows how galanin interacts with one of its three human receptor proteins, GALR2, and reveals the basis of the selectivity of this GPCR for Gq.     

Use of Cysteine Trapping to Map Spatial Approximations between Residues Contributing to the Helix N-capping Motif of Secretin and Distinct Residues within Each of the Extracellular Loops of Its Receptor

Amino-terminal regions of secretin-family peptides contain key determinants for biological activity and binding specificity, although the nature of interactions with receptors is unclear. A helix N-capping motif within this region has been postulated to directly contribute to agonist activity while also stabilizing formation of a helix extending toward the peptide carboxyl terminus and docking within the receptor amino terminus. We used cysteine trapping to systematically explore spatial approximations between cysteines replacing each residue in this motif of secretin (sec), Phe⁶, Thr⁷, and Leu¹⁰, and cysteines incorporated into the extracellular face of the receptor. Each peptide was a full agonist for cAMP, but had a lower binding affinity than natural hormone. These bound to COS cells expressing 61 receptor constructs incorporating cysteines in every position along each extracellular loop (ECL) and adjacent parts of transmembrane (TM) segments. Patterns of covalent labeling were distinct for each probe, with Cys⁶-sec labeling multiple residues in the carboxyl-terminal half of ECL2 and throughout ECL3, Cys⁷-sec predominantly labeling only single residues in the carboxyl-terminal end of ECL2 and the amino-terminal end of ECL3, and Cys¹⁰-sec not efficiently labeling any of these residues. These spatial constraints were used to refine our model of secretin bound to its receptor, now bringing ECL3 above the amino terminus of the ligand and revealing possible charge-charge interactions between this part of secretin and receptor residues in TM5, TM6, ECL2, and ECL3, which can orient and stabilize the peptide-receptor complex. This was validated by testing predicted approximations by mutagenesis and residue-residue complementation studies.     

Mechanistic Insights into Allosteric Structure-Function Relationships at the M1 Muscarinic Acetylcholine Receptor

Benzylquinolone carboxylic acid (BQCA) is the first highly selective positive allosteric modulator (PAM) for the M₁ muscarinic acetylcholine receptor (mAChR), but it possesses low affinity for the allosteric site on the receptor. More recent drug discovery efforts identified 3-((1S,2S)-2-hydroxycyclohexyl)-6-((6-(1-methyl-1H-pyrazol-4-yl)pyridin-3-yl)methyl)benzo[h]quinazolin-4(3H)-one (referred to herein as benzoquinazolinone 12) as a more potent M₁ mAChR PAM with a structural ancestry originating from BQCA and related compounds. In the current study, we optimized the synthesis of and fully characterized the pharmacology of benzoquinazolinone 12, finding that its improved potency derived from a 50-fold increase in allosteric site affinity as compared with BQCA, while retaining a similar level of positive cooperativity with acetylcholine. We then utilized site-directed mutagenesis and molecular modeling to validate the allosteric binding pocket we previously described for BQCA as a shared site for benzoquinazolinone 12 and provide a molecular basis for its improved activity at the M₁ mAChR. This includes a key role for hydrophobic and polar interactions with residues Tyr-179, in the second extracellular loop (ECL2) and Trp-400^7.35 in transmembrane domain (TM) 7. Collectively, this study highlights how the properties of affinity and cooperativity can be differentially modified on a common structural scaffold and identifies molecular features that can be exploited to tailor the development of M₁ mAChR-targeting PAMs.     

Negative allosteric modulators of cannabinoid receptor 2: protein modeling,binding site identification and molecular dynamics simulations in the presence of an orthosteric agonist

Abstract

Selective activation of the cannabinoid receptor subtype 2 (CB2) shows promise for treating pain, inflammation, multiple sclerosis, cancer, ischemic/reperfusion injury and osteoporosis. Target selectivity and off-target side effects are two major limiting factors for orthosteric ligands, and therefore, the search for allosteric modulators (AMs) is a widely used drug discovery approach. To date, only a limited number of negative CB2 AMs have been identified, possessing only micromolar activity at best, and the CB2 receptor’s allosteric site(s) are not well characterized. Herein, we used computational approaches including receptor modeling, site mapping, docking, molecular dynamics (MD) simulations and binding free energy calculations to predict, characterize and validate allosteric sites within the complex of the CB2 receptor with bound orthosteric agonist CP55,940. After docking of known negative CB2 allosteric modulators (NAMs), dihydro-gambogic acid (DHGA) and trans-β-caryophyllene (TBC) (note that TBC also shows agonist activity), at the predicted allosteric sites, the best total complex with CB2, CP55,940 and NAM was embedded into a hydrated lipid bilayer and subjected to a 200 ns MD simulation. The presence of an AM affected the CB2–CP55,940 complex, altering the relative positioning of the toggle switch residues and promoting a strong π–π interaction between Phe117^3.36 and Trp258^6.48. Binding of either TBC or DHGA to a putative allosteric pocket directly adjacent to the orthosteric ligand reduced the binding free energy of CP55,940, which is consistent with the expected effect of a negative AM. The identified allosteric sites present immense scope for the discovery of novel classes of CB2 AMs.     

Identification of a 3rd Na+ Binding Site of the Glycine Transporter,GlyT2

The Na⁺/Cl^- dependent glycine transporters GlyT1 and GlyT2 regulate synaptic glycine concentrations. Glycine transport by GlyT2 is coupled to the co-transport of three Na⁺ ions, whereas transport by GlyT1 is coupled to the co-transport of only two Na⁺ ions. These differences in ion-flux coupling determine their respective concentrating capacities and have a direct bearing on their functional roles in synaptic transmission. The crystal structures of the closely related bacterial Na⁺-dependent leucine transporter, LeuT_Aa, and the Drosophila dopamine transporter, dDAT, have allowed prediction of two Na⁺ binding sites in GlyT2, but the physical location of the third Na⁺ site in GlyT2 is unknown. A bacterial betaine transporter, BetP, has also been crystallized and shows structural similarity to LeuT_Aa. Although betaine transport by BetP is coupled to the co-transport of two Na⁺ ions, the first Na⁺ site is not conserved between BetP and LeuT_Aa, the so called Na1'' site. We hypothesized that the third Na⁺ binding site (Na3 site) of GlyT2 corresponds to the BetP Na1'' binding site. To identify the Na3 binding site of GlyT2, we performed molecular dynamics (MD) simulations. Surprisingly, a Na⁺ placed at the location consistent with the Na1'' site of BetP spontaneously dissociated from its initial location and bound instead to a novel Na3 site. Using a combination of MD simulations of a comparative model of GlyT2 together with an analysis of the functional properties of wild type and mutant GlyTs we have identified an electrostatically favorable novel third Na⁺ binding site in GlyT2 formed by Trp263 and Met276 in TM3, Ala481 in TM6 and Glu648 in TM10.     

Long Range Effect of Mutations on Specific Conformational Changes in the Extracellular Loop 2 of Angiotensin II Type 1 Receptor

The topology of the second extracellular loop (ECL2) and its interaction with ligands is unique in each G protein-coupled receptor. When the orthosteric ligand pocket located in the transmembrane (TM) domain is occupied, ligand-specific conformational changes occur in the ECL2. In more than 90% of G protein-coupled receptors, ECL2 is tethered to the third TM helix via a disulfide bond. Therefore, understanding the extent to which the TM domain and ECL2 conformations are coupled is useful. To investigate this, we examined conformational changes in ECL2 of the angiotensin II type 1 receptor (AT1R) by introducing mutations in distant sites that alter the activation state equilibrium of the AT1R. Differential accessibility of reporter cysteines introduced at four conformation-sensitive sites in ECL2 of these mutants was measured. Binding of the agonist angiotensin II (AngII) and inverse agonist losartan in wild-type AT1R changed the accessibility of reporter cysteines, and the pattern was consistent with ligand-specific “lid” conformations of ECL2. Without agonist stimulation, the ECL2 in the gain of function mutant N111G assumed a lid conformation similar to AngII-bound wild-type AT1R. In the presence of inverse agonists, the conformation of ECL2 in the N111G mutant was similar to the inactive state of wild-type AT1R. In contrast, AngII did not induce a lid conformation in ECL2 in the loss of function D281A mutant, which is consistent with the reduced AngII binding affinity in this mutant. However, a lid conformation was induced by [Sar¹,Gln²,Ile⁸] AngII, a specific analog that binds to the D281A mutant with better affinity than AngII. These results provide evidence for the emerging paradigm of domain coupling facilitated by long range interactions at distant sites on the same receptor.     

Phenotypic analysis of boldit Galr1 knockout mice reveals a role for GALR1 galanin receptor in modulating seizure activity but not nerve regeneration

The GALR1 galanin receptor is expressed at high levels withinthe central nervous system and is hypothesised to play asignificant role in many of the central actions of galanin. Todetermine which specific actions of galanin are mediated byGALR1, we have developed mice that carry an insertionalinactivating mutation within the first coding exon of the geneencoding GALR1 (Galr1). HomozygousGalr1 ^-/-mice are viable. Both male and female mice exhibit reducedcirculating levels of insulin-like growth factor-I (IGF-I) butno significant difference in growth rate relative to Galr1 ^+/+ controls. Female homozygousGalr1 ^-/-mice are capable of breeding and nursing offspring. Functionalrecovery after sciatic nerve crush is not significantlydifferent in Galr1 ^-/- mice relative to Galr1 ^+/+ controls, indicating that GALR1 does not mediate the nerve regenerative effects of galanin. However, homozygous Galr1 ^-/- mice exhibit spontaneous seizures, identifying a critical role for GALR1 in mediating the anti-seizure activity of galanin.     

Phenotypic analysis of boldit Galr1 knockout mice reveals a role for GALR1 galanin receptor in modulating seizure activity but not nerve regeneration

Summary The GALR1 galanin receptor is expressed at high levels within the central nervous system and is hypothesised to play a significant role in many of the central actions of galanin. To determine which specific actions of galanin are mediated by GALR1, we have developed mice that carry an insertional inactivating mutation within the first coding exon of the gene encoding GALR1 (Galr1). HomozygousGalr1 ^−/− mice are viable. Both male and female mice exhibit reduced circulating levels of insulin-like growth factor-I (IGF-I) but no significant difference in growth rate relative toGalr1 ^+/+ controls. Female homozygousGalr1 ^−/− mice are capable of breeding and nursing offspring. Functional recovery after sciatic nerve crush is not significantly different inGalr1 ^−/− mice relative toGalr1 ^+/+ controls, indicating that GALR1 does not mediate the nerve regenerative effects of galanin. However, homozygousGalr1 ^−/− mice exhibit spontaneous seizures, identifying a critical role for GALR1 in mediating the anti-seizure activity of galanin.     

Molecular modelling and site-directed mutagenesis of human GALR1 galanin receptor defines determinants of receptor subtype specificity

Human galanin is a 30 amino acid neuropeptide that elicits a range of biological activities by interaction with G protein-coupled receptors. We have generated a model of the human GALR1 galanin receptor subtype (hGALR1) based on the alpha carbon maps of frog rhodopsin and investigated the significance of potential contact residues suggested by the model using site-directed mutagenesis. Mutation of Phe186 within the second extracellular loop to Ala resulted in a 6-fold decrease in affinity for galanin, representing a change in free energy consistent with hydrophobic interaction. Our model suggests interaction between Phe186 of hGALR1 and Ala7 or Leu11 of galanin. Receptor subtype specificity was investigated by replacement of residues in hGALR1 with the corresponding residues in hGALR2 and use of the hGALR2-specific ligands hGalanin(2-30) and [D-Trp2]hGalanin(1-30). The His267Ile mutant receptor exhibited a pharmacological profile corresponding to that of hGALR1, suggesting that His267 is not involved in a receptor-ligand interaction. The mutation Phe115Ala resulted in a decreased binding affinity for hGalanin and for hGALR2-specific analogues, indicating Phe115 to be of structural importance to the ligand binding pocket of hGALR1 but not involved in direct ligand interaction. Analysis of Glu271Trp suggested that Glu271 of hGALR1 interacts with the N-terminus of galanin and that the Trp residue in the corresponding position in hGALR2 is involved in receptor subtype specificity of binding. Our model supports previous reports of Phe282 of hGALR1 interacting with Trp2 of galanin and His264 of hGALR1 interacting with Tyr9 of galanin.     

Isolation and cDNA cloning of a novel galanin-like peptide (GALP) from porcine hypothalamus

Galanin is a widely distributed neuropeptide with a variety of physiological functions. Three galanin receptor subtypes, GALR1, GALR2, and GALR3, have been reported. We isolated a novel galanin-like peptide (GALP) from porcine hypothalamus by observing its activity for increasing [(35)S]GTPgammaS binding to a membrane preparation of GALR2-transfected cells. The peptide had 60 amino acid residues and a non-amidated C terminus. The amino acid sequence of GALP-(9-21) was completely identical to that of galanin-(1-13). A cloned porcine GALP cDNA indicated that GALP was processed from a 120-amino acid GALP precursor protein. The structures of rat and human GALP-(1-60) were deduced from cloned cDNA, which indicated that the amino acid sequences 1-24 and 41-53 were highly conserved between humans, rats, and pigs. Receptor binding studies revealed that porcine GALP-(1-60) had a high affinity for the GALR2 receptor (IC(50) = 0.24 nM) and a lower affinity for the GALR1 receptor (IC(50) = 4.3 nM). In contrast, galanin showed high affinity for the GALR1 (IC(50) = 0.097 nM) and GALR2 receptors (IC(50) = 0.48 nM). GALP is therefore an endogenous ligand that preferentially binds the GALR2 receptor, whereas galanin is relatively non-selective.     

Functional Analysis of the Na+,K+/H+ Antiporter PeNHX3 from the Tree Halophyte Populus euphratica in Yeast by Model-Guided Mutagenesis

Na⁺,K⁺/H⁺ antiporters are H⁺-coupled cotransporters that are crucial for cellular homeostasis. Populus euphratica, a well-known tree halophyte, contains six Na⁺/H⁺ antiporter genes (PeNHX1-6) that have been shown to function in salt tolerance. However, the catalytic mechanisms governing their ion transport remain largely unknown. Using the crystal structure of the Na⁺/H⁺ antiporter from the Escherichia coli (EcNhaA) as a template, we built the three-dimensional structure of PeNHX3 from P. euphratica. The PeNHX3 model displays the typical TM4-TM11 assembly that is critical for ion binding and translocation. The PeNHX3 structure follows the ‘positive-inside’ rule and exhibits a typical physicochemical property of the transporter proteins. Four conserved residues, including Tyr149, Asn187, Asp188, and Arg356, are indentified in the TM4-TM11 assembly region of PeNHX3. Mutagenesis analysis showed that these reserved residues were essential for the function of PeNHX3: Asn187 and Asp188 (forming a ND motif) controlled ion binding and translocation, and Tyr149 and Arg356 compensated helix dipoles in the TM4-TM11 assembly. PeNHX3 mediated Na⁺, K⁺ and Li⁺ transport in a yeast growth assay. Domain-switch analysis shows that TM11 is crucial to Li⁺ transport. The novel features of PeNHX3 in ion binding and translocation are discussed.     

Binding of histamine to the H1 receptor—a molecular dynamics study

Binding of histamine to the G-protein coupled histamine H₁ receptor plays an important role in the context of allergic reactions; however, no crystal structure of the resulting complex is available yet. To deduce the histamine binding site, we performed unbiased molecular dynamics (MD) simulations on a microsecond time scale, which allowed to monitor one binding event, in which particularly the residues of the extracellular loop 2 were involved in the initial recognition process. The final histamine binding pose in the orthosteric pocket is characterized by interactions with Asp107^3.32, Tyr108^3.33, Thr194^5.43, Asn198^5.46, Trp428^6.48, Tyr431^6.51, Phe432^6.52, and Phe435^6.55, which is in agreement with existing mutational data. The conformational stability of the obtained complex structure was subsequently confirmed in 2 μs equilibrium MD simulations, and a metadynamics simulation proved that the detected binding site represents an energy minimum. A complementary investigation of a D107A mutant, which has experimentally been shown to abolish ligand binding, revealed that this exchange results in a significantly weaker interaction and enhanced ligand dynamics. This finding underlines the importance of the electrostatic interaction between the histamine ammonium group and the side chain of Asp107^3.32 for histamine binding.     

Ligand Entry and Exit Pathways in the β2-Adrenergic Receptor

The recently determined crystal structure of the human β₂-adrenergic (β₂AR) G-protein-coupled receptor provides an excellent structural basis for exploring β₂AR-ligand binding and dissociation process. Based on this crystal structure, we simulated ligand exit from the β₂AR receptor by applying the random acceleration molecular dynamics (RAMD) simulation method. The simulation results showed that the extracellular opening on the receptor surface was the most frequently observed egress point (referred to as pathway A), and a few other pathways through interhelical clefts were also observed with significantly lower frequencies. In the egress trajectories along pathway A, the D192-K305 salt bridge between the extracellular loop 2 (ECL2) and the apex of the transmembrane helix 7 (TM7) was exclusively broken. The spatial occupancy maps of the ligand computed from the 100 RAMD simulation trajectories indicated that the receptor-ligand interactions that restrained the ligand in the binding pocket were the major resistance encountered by the ligand during exit and no second barrier was notable. We next performed RAMD simulations by using a putative ligand-free conformation of the receptor as input structure. This conformation was obtained in a standard molecular dynamics simulation in the absence of the ligand and it differed from the ligand-bound conformation in a hydrophobic patch bridging ECL2 and TM7 due to the rotation of F193 of ECL2. Results from the RAMD simulations with this putative ligand-free conformation suggest that the cleft formed by the hydrophobic bridge, TM2, TM3, and TM7 on the extracellular surface likely serves as a more specific ligand-entry site and the ECL2-TM7 hydrophobic junction can be partially interrupted upon the entry of ligand that pushes F193 to rotate, resulting in a conformation as observed in the ligand-bound crystal structure. These results may help in the design of β₂AR-targeting drugs with improved efficacy, as well as in understanding the receptor subtype selectivity of ligand binding in the β family of the adrenergic receptors that share almost identical ligand-binding pockets, but show notable amino acid sequence divergence in the putative ligand-entry site, including ECL2 and the extracellular end of TM7.     

The Role of Water in Activation Mechanism of Human N-Formyl Peptide Receptor 1 (FPR1) Based on Molecular Dynamics Simulations

The Formyl Peptide Receptor 1 (FPR1) is an important chemotaxis receptor involved in various aspects of host defense and inflammatory processes. We constructed a model of FPR1 using as a novel template the chemokine receptor CXCR4 from the same branch of the phylogenetic tree of G-protein-coupled receptors. The previously employed template of rhodopsin contained a bulge at the extracellular part of TM2 which directly influenced binding of ligands. We also conducted molecular dynamics (MD) simulations of FPR1 in the apo form as well as in a form complexed with the agonist fMLF and the antagonist tBocMLF in the model membrane. During all MD simulation of the fMLF-FPR1 complex a water molecule transiently bridged the hydrogen bond between W254^6.48 and N108^3.35 in the middle of the receptor. We also observed a change in the cytoplasmic part of FPR1 of a rotamer of the Y301^7.53 residue (tyrosine rotamer switch). This effect facilitated movement of more water molecules toward the receptor center. Such rotamer of Y301^7.53 was not observed in any crystal structures of GPCRs which can suggest that this state is temporarily formed to pass the water molecules during the activation process. The presence of a distance between agonist and residues R201^5.38 and R205^5.42 on helix TM5 may suggest that the activation of FPR1 is similar to the activation of β-adrenergic receptors since their agonists are separated from serine residues on helix TM5. The removal of water molecules bridging these interactions in FPR1 can result in shrinking of the binding site during activation similarly to the shrinking observed in β-ARs. The number of GPCR crystal structures with agonists is still scarce so the designing of new ligands with agonistic properties is hampered, therefore homology modeling and docking can provide suitable models. Additionally, the MD simulations can be beneficial to outline the mechanisms of receptor activation and the agonist/antagonist sensing.     

Extracellular Allosteric Na Binding to the Na,K-ATPase in Cardiac Myocytes

Whole-cell patch-clamp measurements of the current, I_p, produced by the Na⁺,K⁺-ATPase across the plasma membrane of rabbit cardiac myocytes show an increase in I_p over the extracellular Na⁺ concentration range 0–50 mM. This is not predicted by the classical Albers-Post scheme of the Na⁺,K⁺-ATPase mechanism, where extracellular Na⁺ should act as a competitive inhibitor of extracellular K⁺ binding, which is necessary for the stimulation of enzyme dephosphorylation and the pumping of K⁺ ions into the cytoplasm. The increase in I_p is consistent with Na⁺ binding to an extracellular allosteric site, independent of the ion transport sites, and an increase in turnover via an acceleration of the rate-determining release of K⁺ to the cytoplasm, E2(K⁺)₂ → E1 + 2K⁺. At normal physiological concentrations of extracellular Na⁺ of 140 mM, it is to be expected that binding of Na⁺ to the allosteric site would be nearly saturated. Its purpose would seem to be simply to optimize the enzyme’s ion pumping rate under its normal physiological conditions. Based on published crystal structures, a possible location of the allosteric site is within a cleft between the α- and β-subunits of the enzyme.