Phenotypic classification of mutants: a tool for understanding ligand binding and activation of muscarinic acetylcholine receptors

Exploration of the orthosteric/allosteric interface in human M1 muscarinic receptors by bitopic fluorescent ligands

Bitopic binding properties apply to a variety of muscarinic compounds that span and simultaneously bind to both the orthosteric and allosteric receptor sites. We provide evidence that fluorescent pirenzepine derivatives, with the M1 antagonist fused to the boron-dipyrromethene [Bodipy (558/568)] fluorophore via spacers of varying lengths, exhibit orthosteric/allosteric binding properties at muscarinic M1 receptors. This behavior was inferred from a combination of functional, radioligand, and fluorescence resonance energy transfer binding experiments performed under equilibrium and kinetic conditions on enhanced green fluorescent protein-fused M1 receptors. Although displaying a common orthosteric component, the fluorescent compounds inherit bitopic properties from a linker-guided positioning of their Bodipy moiety within the M1 allosteric vestibule. Depending on linker length, the fluorophore is allowed to reach neighboring allosteric domains, overlapping or not with the classic gallamine site, but distinct from the allosteric indolocarbazole "WIN" site. Site-directed mutagenesis, as well as molecular modeling and ligand docking studies based on recently solved muscarinic receptor structures, further support the definition of two groups of Bodipy-pirenzepine derivatives exhibiting distinct allosteric binding poses. Thus, the linker may dictate pharmacological outcomes for bitopic molecules that are hardly predictable from the properties of individual orthosteric and allosteric building blocks. Our findings also demonstrate that the fusion of a fluorophore to an orthosteric ligand is not neutral, as it may confer, unless carefully controlled, unexpected properties to the resultant fluorescent tracer. Altogether, this study illustrates the importance of a "multifacet" experimental approach to unravel and validate bitopic ligand binding mechanisms.

Role of conserved glycine in zinc-dependent medium chain dehydrogenase/reductase superfamily

The medium-chain dehydrogenase/reductase (MDR) superfamily consists of a large group of enzymes with a broad range of activities. Members of this superfamily are currently the subject of intensive investigation, but many aspects, including the zinc dependence of MDR superfamily proteins, have not yet have been adequately investigated. Using a density functional theory-based screening strategy, we have identified a strictly conserved glycine residue (Gly) in the zinc-dependent MDR superfamily. To elucidate the role of this conserved Gly in MDR, we carried out a comprehensive structural, functional, and computational analysis of four MDR enzymes through a series of studies including site-directed mutagenesis, isothermal titration calorimetry, electron paramagnetic resonance (EPR), quantum mechanics, and molecular mechanics analysis. Gly substitution by other amino acids posed a significant threat to the metal binding affinity and activity of MDR superfamily enzymes. Mutagenesis at the conserved Gly resulted in alterations in the coordination of the catalytic zinc ion, with concomitant changes in metal-ligand bond length, bond angle, and the affinity (K(d)) toward the zinc ion. The Gly mutants also showed different spectroscopic properties in EPR compared with those of the wild type, indicating that the binding geometries of the zinc to the zinc binding ligands were changed by the mutation. The present results demonstrate that the conserved Gly in the GHE motif plays a role in maintaining the metal binding affinity and the electronic state of the catalytic zinc ion during catalysis of the MDR superfamily enzymes.

Molecular modeling of the M3 acetylcholine muscarinic receptor and its binding site

The present study reports the results of a combined computational and site mutagenesis study designed to provide new insights into the orthosteric binding site of the human M3 muscarinic acetylcholine receptor. For this purpose a three-dimensional structure of the receptor at atomic resolution was built by homology modeling, using the crystallographic structure of bovine rhodopsin as a template. Then, the antagonist N-methylscopolamine was docked in the model and subsequently embedded in a lipid bilayer for its refinement using molecular dynamics simulations. Two different lipid bilayer compositions were studied: one component palmitoyl-oleyl phosphatidylcholine (POPC) and two-component palmitoyl-oleyl phosphatidylcholine/palmitoyl-oleyl phosphatidylserine (POPC-POPS). Analysis of the results suggested that residues F222 and T235 may contribute to the ligand-receptor recognition. Accordingly, alanine mutants at positions 222 and 235 were constructed, expressed, and their binding properties determined. The results confirmed the role of these residues in modulating the binding affinity of the ligand.

Molecular scanning: combining random mutagenesis, ribosome display, and bioinformatic analysis for protein engineering

Protein engineering techniques can facilitate the direct de-convolution of specific domains, regions, and particular amino acids that contribute to protein function. Many tools are available to aid this enterprise and herein we describe one such tool, a technique we term "Molecular Scanning" (MS). MS is analogous to previously described alanine scanning in that it samples potentially functional sequence space, but differs in that it uses Error-Prone polymerase chain reaction to randomly introduce all amino acids across the sequence space, as opposed to simply introducing alanine at each desired position. We commonly use MS in conjunction with ribosome-display, selecting for specific character traits (e.g., improved affinity) which allows us to sample functionally relevant diversity on a reasonably large scale. This approach is amenable to a variety of different mutational techniques and display technologies as dictated by user requirements or needs. In this chapter we present a general outline of the process as we have previously successfully applied it.

Functional and structural roles of residues in the third extramembrane segment of adrenal cytochrome b561

Several residues in the third extramembrane segment (EM3) of adrenal cytochrome b(561) have been proposed to be involved in this cytochrome's interaction with ascorbate, but there has been no systematic evaluation of residues in the segment. We used alanine scanning mutagenesis to assess the functional and structural roles of the EM3 residues and several adjacent residues (residues 70-85) in the bovine cytochrome. Each alanine mutant was expressed in a bacterial system, solubilized with detergent, and affinity-purified. The recombinant proteins contained approximately two hemes per monomer and, except for R74A, retained basic functionality (≥ 94% reduced by 20 mM ascorbate). Equilibrium spectrophotometric titrations with ascorbate were used to analyze the α-band line shape and amplitude during reduction of the high- and low-potential heme centers (b(H) and b(L), respectively) and the midpoint ascorbate concentrations for the b(H) and b(L) transitions (C(H) and C(L), respectively). Y73A and K85A markedly narrowed the b(H) α-band peak; other mutants had weaker effects or no effect on b(H) or b(L) spectra. Relative changes in C(H) for the mutants were larger than changes in C(L), with 1.5-2.9-fold increases in C(H) for L70A, L71A, Y73A, R74A, N78A, and K85A. The amounts of functional b(H) and b(L) centers in additional Arg74 mutants, assessed by ascorbate titration and EPR spectroscopy, declined in concert in the following order: wild type > R74K > R74Q > R74T and R74Y > R74E. The results of this first comprehensive experimental test of the proposed roles of EM3 residues have identified residues with a direct or indirect impact on ascorbate interactions, on the environment of the b(H) heme center, and on formation of the native b(H)-b(L) unit. Surprisingly, no individual EM3 residue was by itself indispensable for the interaction with ascorbate, and the role of the segment appears to be more subtle than previously thought. These results also support our topological model of the adrenal cytochrome, which positions b(H) near the cytoplasmic side of the membrane.

Helix 8 of the M1 muscarinic acetylcholine receptor: scanning mutagenesis delineates a G protein recognition site

We have used alanine-scanning mutagenesis followed by functional expression and molecular modeling to analyze the roles of the 14 residues, Asn422 to Cys435, C-terminal to transmembrane (TM) helix 7 of the M(1) muscarinic acetylcholine receptor. The results suggest that they form an eighth (H8) helix, associated with the cytoplasmic surface of the cell membrane in the active state of the receptor. We suggest that the amide side chain of Asn422 may act as a cap to the C terminus of TM7, stabilizing its junction with H8, whereas the side chain of Phe429 may restrict the relative movements of H8 and the C terminus of TM7 in the inactive ground state of the receptor. We have identified four residues, Phe425, Arg426, Thr428, and Leu432, which are important for G protein binding and signaling. These may form a docking site for the C-terminal helix of the G protein α subunit, and collaborate with G protein recognition residues elsewhere in the cytoplasmic domain of the receptor to form a coherent surface for G protein binding in the activated state of the receptor.

In silico saturation mutagenesis and docking screening for the analysis of protein-ligand interaction: the Endothelial Protein C Receptor case study

Background

The design of mutants in protein functional regions, such as the ligand binding sites, is a powerful approach to recognize the determinants of specific protein activities in cellular pathways. For an exhaustive analysis of selected positions of protein structure large scale mutagenesis techniques are often employed, with laborious and time consuming experimental set-up. 'In silico' mutagenesis and screening simulation represents a valid alternative to laboratory methods to drive the 'in vivo' testing toward more focused objectives.

Results

We present here a high performance computational procedure for large-scale mutant modelling and subsequent evaluation of the effect on ligand binding affinity. The mutagenesis was performed with a 'saturation' approach, where all 20 natural amino acids were tested in positions involved in ligand binding sites. Each modelled mutant was subjected to molecular docking simulation and stability evaluation. The simulated protein-ligand complexes were screened for their impairment of binding ability based on change of calculated Ki compared to the wild-type.

An example of application to the Endothelial Protein C Receptor residues involved in lipid binding is reported.

Conclusion

The computational pipeline presented in this work is a useful tool for the design of structurally stable mutants with altered affinity for ligand binding, considerably reducing the number of mutants to be experimentally tested. The saturation mutagenesis procedure does not require previous knowledge of functional role of the residues involved and allows extensive exploration of all possible substitutions and their pairwise combinations. Mutants are screened by docking simulation and stability evaluation followed by a rationally driven selection of those presenting the required characteristics. The method can be employed in molecular recognition studies and as a preliminary approach to select models for experimental testing.

Mutagenic mapping suggests a novel binding mode for selective agonists of M1 muscarinic acetylcholine receptors

Point mutations and molecular modeling have been used to study the activation of the M(1) muscarinic acetylcholine receptor (mAChR) by the functionally selective agonists 4-n-butyl-1-[4-(2-methylphenyl)-4-oxo-1-butyl]-piperidine (AC-42), and 1-[3-(4-butyl-1-piperidinyl)propyl]-3,4-dihydro-2(1H)-quinolinone (77-LH-28-1), comparing them with N-desmethylclozapine (NDMC) and acetylcholine (ACh). Unlike NDMC and ACh, the activities of AC-42 and 77-LH-28-1 were undiminished by mutations of Tyr404 and Cys407 (transmembrane helix 7), although they were reduced by mutations of Tyr408. Signaling by AC-42, 77-LH-28-1, and NDMC was reduced by L102A and abolished by D105E, suggesting that all three may interact with transmembrane helix 3 at or near the binding site Asp105 to activate the M(1) mAChR. In striking contrast to NDMC and ACh, the affinities of AC-42 and 77-LH-28-1 were increased 100-fold by W101A, and their signaling activities were abolished by Y82A. Tyr82 and Leu102 contact the indole ring of Trp101 in a structural model of the M(1) mAChR. We suggest the hypothesis that the side chain of Trp101 undergoes conformational isomerization, opening a novel binding site for the aromatic side chain of the AC-42 analogs. This may allow the positively charged piperidine nitrogen of the ligands to access the neighboring Asp105 carboxylate to activate signaling following a vector within the binding site that is distinct from that of acetylcholine. NDMC does not seem to use this mechanism. Subtype-specific differences in the free energy of rotation of the side chain and indole ring of Trp101 might underlie the M(1) selectivity of the AC-42 analogs. Tryptophan conformational isomerization may open up new avenues in selective muscarinic receptor drug design.

Family Resemblances? Ligand Binding and Activation of Family A and B G-Protein-Coupled Receptors

In April 2007, the Biochemical Society held a meeting to compare and contrast ligand binding and activation of Family A and B GPCRs (G-protein-coupled receptors). Being the largest class, Family A GPCRs usually receive the most attention, although a previous Biochemical Society meeting has focused on Family B GPCRs. The aim of the present meeting was to bring researchers of both families together in order to identify commonalities between the two. The present article introduces the proceedings of the meeting, briefly commenting on the focus of each of the following articles.