Correction: Ligand Binding at the ?4-?4 Agonist-Binding Site of the ?4?2 nAChR Triggers Receptor Activation through a Pre-Activated Conformational State.
Ontology highlight
ABSTRACT: [This corrects the article DOI: 10.1371/journal.pone.0161154.].
Correction: Ligand Binding at the α4-α4 Agonist-Binding Site of the α4β2 nAChR Triggers Receptor Activation through a Pre-Activated Conformational State.
PloS one 20161019 10
[This corrects the article DOI: 10.1371/journal.pone.0161154.]. ...[more]
Project description:The ?4?2 nicotinic acetylcholine receptor (nAChR) is the most abundant subtype in the brain and exists in two functional stoichiometries: (?4)3(?2)2 and (?4)2(?2)3. A distinct feature of the (?4)3(?2)2 receptor is the biphasic activation response to the endogenous agonist acetylcholine, where it is activated with high potency and low efficacy when two ?4-?2 binding sites are occupied and with low potency/high efficacy when a third ?4-?4 binding site is occupied. Further, exogenous ligands can bind to the third ?4-?4 binding site and potentiate the activation of the receptor by ACh that is bound at the two ?4-?2 sites. We propose that perturbations of the recently described pre-activation step when a third binding site is occupied are a key driver of these distinct activation properties. To investigate this, we used a combination of simple linear kinetic models and voltage clamp electrophysiology to determine whether transitions into the pre-activated state were increased when three binding sites were occupied. We separated the binding at the two different sites with ligands selective for the ?4-?2 site (Sazetidine-A and TC-2559) and the ?4-?4 site (NS9283) and identified that when a third binding site was occupied, changes in the concentration-response curves were best explained by an increase in transitions into a pre-activated state. We propose that perturbations of transitions into a pre-activated state are essential to explain the activation properties of the (?4)3(?2)2 receptor by acetylcholine and other ligands. Considering the widespread clinical use of benzodiazepines, this discovery of a conserved mechanism that benzodiazepines and ACh potentiate receptor activation via a third binding site can be exploited to develop therapeutics with similar properties at other cys-loop receptors.
Project description:Neurological glutamate receptors bind a variety of artificial ligands, both agonistic and antagonistic, in addition to glutamate. Studying their small molecule binding properties increases our understanding of the central nervous system and a variety of associated pathologies. The large, oligomeric multidomain membrane protein contains a large and flexible ligand binding domains which undergoes large conformational changes upon binding different ligands. A recent application of glutamate receptors is their activation or inhibition via photo-switchable ligands, making them key systems in the emerging field of optochemical genetics. In this work, we present a theoretical study on the binding mode and complex stability of a novel photo-switchable ligand, ATA-3, which reversibly binds to glutamate receptors ligand binding domains (LBDs). We propose two possible binding modes for this ligand based on flexible ligand docking calculations and show one of them to be analogues to the binding mode of a similar ligand, 2-BnTetAMPA. In long MD simulations, it was observed that transitions between both binding poses involve breaking and reforming the T686-E402 protein hydrogen bond. Simulating the ligand photo-isomerization process shows that the two possible configurations of the ligand azo-group have markedly different complex stabilities and equilibrium binding modes. A strong but slow protein response is observed after ligand configuration changes. This provides a microscopic foundation for the observed difference in ligand activity upon light-switching.
Project description:IL-17 is the hallmark cytokine of the newly described "Th17" lymphocyte population. The composition, subunit dynamics, and ligand contacts of the IL-17 receptor are poorly defined. We previously demonstrated that the IL-17RA subunit oligomerizes in the membrane without a ligand. In this study, computational modeling identified two fibronectin-III-like (FN) domains in IL-17RA connected by a nonstructured linker, which we predicted to mediate homotypic interactions. In yeast two-hybrid, the membrane-proximal FN domain (FN2), but not the membrane-distal domain (FN1), formed homomeric interactions. The ability of FN2 to drive ligand-independent multimerization was verified by coimmunoprecipitation and fluorescence resonance energy transfer microscopy. Thus, FN2 constitutes a "pre-ligand assembly domain" (PLAD). Further studies indicated that the FN2 linker domain contains the IL-17 binding site, which was never mapped. However, the FN1 domain is also required for high affinity interactions with IL-17. Therefore, although the PLAD is located entirely within FN2, effective ligand binding also involves contributions from the linker and FN1.
Project description:BackgroundModelling the ligand binding site of a protein is an important component of understanding protein-ligand interactions and is being actively studied. Even if the side chains are restricted to rotamers, a set of commonly-observed low-energy conformations, the exhaustive combinatorial search of ligand binding site conformers is known as NP-hard. Here we propose a new method, ROTAIMAGE, for modelling the plausible conformers for the ligand binding site given a fixed backbone structure.ResultsROTAIMAGE includes a procedure of selecting ligand binding site residues, exhaustively searching rotameric conformers, clustering them by dissimilarities in pocket shape, and suggesting a representative conformer per cluster. Prior to the clustering, the list of conformers generated by exhaustive search can be reduced by pruning the conformers that have near identical pocket shapes, which is done using simple bit operations. We tested our approach by modelling the active-site inhibitor binding pockets of matrix metalloproteinase-1 and -13. For both cases, analyzing the conformers based on their pocket shapes substantially reduced the 'computational complexity' (10 to 190 fold). The subsequent clustering revealed that the pocket shapes of both proteins could be grouped into approximately 10 distinct clusters. At this level of clustering, the conformational space spanned by the known crystal structures was well covered. Heatmap analysis identified a few bit blocks that combinatorially dictated the clustering pattern. Using this analytical approach, we demonstrated that each of the bit blocks was associated with a specific pocket residue. Identification of residues that influenced the shape of the pocket is an interesting feature unique to the ROTAIMAGE algorithm.ConclusionsROTAIMAGE is a novel algorithm that was efficient in exploring the conformational space of the ligand binding site. Its ability to identify 'key' pocket residues also provides further insight into conformational flexibility with specific implications for protein-ligand interactions.
Project description:Low pH-induced ligand release and receptor recycling are important steps for endocytosis. The transmembrane protein sortilin, a β-propeller containing endocytosis receptor, internalizes a diverse set of ligands with roles in cell differentiation and homeostasis. The molecular mechanisms of pH-mediated ligand release and sortilin recycling are unresolved. Here we present crystal structures that show the sortilin luminal segment (s-sortilin) undergoes a conformational change and dimerizes at low pH. The conformational change, within all three sortilin luminal domains, provides an altered surface and the dimers sterically shield a large interface while bringing the two s-sortilin C-termini into close proximity. Biophysical and cell-based assays show that members of two different ligand families, (pro)neurotrophins and neurotensin, preferentially bind the sortilin monomer. This indicates that sortilin dimerization and conformational change discharges ligands and triggers recycling. More generally, this work may reveal a double mechanism for low pH-induced ligand release by endocytosis receptors.
Project description:The anesthetic propofol inhibits the currents of the homopentameric ligand-gated ion channel GLIC, yet the crystal structure of GLIC with five propofol molecules bound symmetrically shows an open-channel conformation. To address this dilemma and determine if the symmetry of propofol binding sites affects the channel conformational transition, we performed a total of 1.5 ?s of molecular dynamics simulations for different GLIC systems with propofol occupancies of 0, 1, 2, 3, and 5. GLIC without propofol binding or with five propofol molecules bound symmetrically, showed similar channel conformation and hydration status over multiple replicates of 100-ns simulations. In contrast, asymmetric binding to one, two or three equivalent sites in different subunits accelerated the channel dehydration, increased the conformational heterogeneity of the pore-lining TM2 helices, and shifted the lateral and radial tilting angles of TM2 toward a closed-channel conformation. The results differentiate two groups of systems based on the propofol binding symmetry. The difference between symmetric and asymmetric groups is correlated with the variance in the propofol-binding cavity adjacent to the hydrophobic gate and the force imposed by the bound propofol. Asymmetrically bound propofol produced greater variance in the cavity size that could further elevate the conformation heterogeneity. The force trajectory generated by propofol in each subunit over the course of a simulation exhibits an ellipsoidal shape, which has the larger component tangential to the pore. Asymmetric propofol binding creates an unbalanced force that expedites the channel conformation transitions. The findings from this study not only suggest that asymmetric binding underlies the propofol functional inhibition of GLIC, but also advocate for the role of symmetry breaking in facilitating channel conformational transitions.
Project description:The conformation adopted by a ligand on binding to a receptor may differ from its lowest-energy conformation in solution. In addition, the bound ligand is more conformationally restricted, which is associated with a configurational entropy loss. The free energy change due to these effects is often neglected or treated crudely in current models for predicting binding affinity. We present a method for estimating this contribution, based on perturbation theory using the quasi-harmonic model of Karplus and Kushick as a reference system. The consistency of the method is checked for small model systems. Subsequently we use the method, along with an estimate for the enthalpic contribution due to ligand-receptor interactions, to calculate relative binding affinities. The AMBER force field and generalized Born implicit solvent model is used. Binding affinities were estimated for a test set of 233 protein-ligand complexes for which crystal structures and measured binding affinities are available. In most cases, the ligand conformation in the bound state was significantly different from the most favorable conformation in solution. In general, the correlation between measured and calculated ligand binding affinities including the free energy change due to ligand conformational change is comparable to or slightly better than that obtained by using an empirically-trained docking score. Both entropic and enthalpic contributions to this free energy change are significant.
Project description:Understanding the activation mechanism of Cys loop ion channel receptors is key to understanding their physiological and pharmacological properties under normal and pathological conditions. The ligand-binding domains of these receptors comprise inner and outer beta-sheets and structural studies indicate that channel opening is accompanied by conformational rearrangements in both beta-sheets. In an attempt to resolve ligand-dependent movements in the ligand-binding domain, we employed voltage-clamp fluorometry on alpha1 glycine receptors to compare changes mediated by the agonist, glycine, and by the antagonist, strychnine. Voltage-clamp fluorometry involves labeling introduced cysteines with environmentally sensitive fluorophores and inferring structural rearrangements from ligand-induced fluorescence changes. In the inner beta-sheet, we labeled residues in loop 2 and in binding domain loops D and E. At each position, strychnine and glycine induced distinct maximal fluorescence responses. The pre-M1 domain responded similarly; at each of four labeled positions glycine produced a strong fluorescence signal, whereas strychnine did not. This suggests that glycine induces conformational changes in the inner beta-sheet and pre-M1 domain that may be important for activation, desensitization, or both. In contrast, most labeled residues in loops C and F yielded fluorescence changes identical in magnitude for glycine and strychnine. A notable exception was H201C in loop C. This labeled residue responded differently to glycine and strychnine, thus underlining the importance of loop C in ligand discrimination. These results provide an important step toward mapping the domains crucial for ligand discrimination in the ligand-binding domain of glycine receptors and possibly other Cys loop receptors.
Project description:A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.
Project description:Thrombin participates in coagulation, anticoagulation, and initiation of platelet activation. To fulfill its diverse roles and maintain hemostasis, this serine protease is regulated via the extended active site region and anion-binding exosites (ABEs) I and II. For the current project, amide proton hydrogen-deuterium exchange coupled with MALDI-TOF mass spectrometry was used to characterize ligand binding to individual exosites and to investigate the presence of exosite-active site and exosite-exosite interactions. PAR3(44-56) and PAR1(49-62) were observed to bind to thrombin ABE I and then to exhibit long range effects over to ABE II. By contrast, Hirudin(54-65) focused more on ABE I and did not transmit influences over to ABE II. Although these three ligands were each directed to ABE I, they did not promote the same conformational consequences. D-Phe-Pro-Arg-chloromethyl ketone inhibition at the thrombin active site led to further local and long range consequences to thrombin-ABE I ligand complexes with the autolysis loop often most affected. When Hirudin(54-65) was bound to ABE I, it was still possible to bind GpIbα(269-286) or fibrinogen γ'(410-427) to ABE II. Each ligand exerted its predominant influences on thrombin and also allowed interexosite communication. The results obtained support the proposal that thrombin is a highly dynamic protein. The transmission of ligand-specific local and long range conformational events is proposed to help regulate this multifunctional enzyme.