Sampling and Energy Evaluation Challenges in Ligand Binding Protein Design

The steroid hormone 17α-hydroxylprogesterone (17-OHP) is a biomarker for congenital adrenal hyperplasia and hence there is considerable interest in sensor development towards this compound. We used computational protein design to generate protein models with binding sites for 17-OHP containing an extended, nonpolar, shape-complementary binding pocket corresponding to the four-ring core of the compound, and hydrogen bonding residues at the base of the pocket to interact with carbonyl and hydroxyl groups at the more polar end of the ligand. Eight of sixteen designed proteins experimentally tested bind 17-OHP with micromolar affinity. A co-crystal structure of one of the designs revealed that 17-OHP is rotated 180 degrees around a pseudo-two-fold axis in the compound and displays multiple binding modes within the pocket, while still interacting with all of the designed residues in the engineered site. Subsequent rounds of mutagenesis and binding selection improved the ligand affinity to nanomolar range, while appearing to constrain the ligand to a single bound conformation that maintains the same ‘flipped’ orientation relative to the original design. We trace the discrepancy in the design calculations to two sources: first, a failure to model subtle backbone changes which alter the distribution of sidechain rotameric states and second, an underestimation of the energetic cost of desolvating the carbonyl and hydroxyl groups of the ligand. The difference between design model and crystal structure thus arises from both sampling limitations and energy function inaccuracies that are exacerbated by the near two-fold symmetry of the molecule.

Structure and function of the nicotinic acetylcholine receptor

Jörg Grandl

The subject of this thesis is the study of the structure and function of the human muscle nicotinic acetylcholine receptor (nAChR), which is a typical member of the large class of ligand-gated ion channels. Appropriate techniques, like site-directed mutagenesis, patch-clamp and fluorescence microscopy, and combinations thereof, were applied to gain insight into selected aspects of the structure-function relationship of this receptor. Special emphasis was put on the investigation of single receptors, wherever this promised scientific benefit. Experiments with the aim to study structural rearrangements of nAChR upon ligand binding by fluorescence resonance energy transfer (FRET) were undertaken. Such intramolecular motions are expected to link ligand binding to channel gating and are therefore fundamental to understand the function of nAChR. Double-cysteine mutants were produced and receptors were labelled on living cells in order to measure changes in FRET efficiency due to ligand binding. A novel class of functionally modified ligands for the nAChR was characterized by patch-clamp technique. Derivatives that carry a thiol-reactive group were shown to be full and very potent agonists. On the bases of a structural model of the binding site, residues were selected, that are likely to be in contact with the bound ligand. Equivalent single-cysteine mutants were produced to covalently bind these agonists and thereby to gain information about the mechanisms of ligand binding and receptor activation. Ligands that were coupled to fluorophores turned out to be either partial or full agonists with nanomolar efficacies for the nAChR. Their use to study the diffusion of single nAChR on living HEK293-cells for different functional states is described in detail. The dependency of the lateral diffusion on the activation state (closed/opened/desensitized) was studied. Additionally these ligands served to demonstrate the feasibility to simultaneously combine patch-clamp and fluorescent binding experiments on a single-molecule level. Acetylcholine receptors carrying a single point-mutation, related to the congenital myasthenic syndrome, were found to stabilize an intermediate conductance state and were studied in detail by single-channel electrophysiology. The dependency of the gating behavior on different ligands, ligand concentration and transmembrane voltage was investigated. Mutants, carrying equivalent mutations on different subunits were produced to reveal the origin of this effect. Quantitative characterization was obtained by dwell-time and spectral noise-analysis, in order to gain insight on sub-millisecond fluctuations of the channel lining transmembrane regions of the nAChR and the gating mechanism in general.

EPFL2005

Computational design of ligand-binding proteins with high affinity and selectivity

Kai Johnsson, Alberto Schena

The ability to design proteins with high affinity and selectivity for any given small molecule is a rigorous test of our understanding of the physiochemical principles that govern molecular recognition. Attempts to rationally design ligand-binding proteins have met with little success, however, and the computational design of protein-small-molecule interfaces remains an unsolved problem. Current approaches for designing ligand-binding proteins for medical and biotechnological uses rely on raising antibodies against a target antigen in immunized animals and/or performing laboratory-directed evolution of proteins with an existing low affinity for the desired ligand, neither of which allows complete control over the interactions involved in binding. Here we describe a general computational method for designing pre-organized and shape complementary small-molecule-binding sites, and use it to generate protein binders to the steroid digoxigenin (DIG). Of seventeen experimentally characterized designs, two bind DIG; the model of the higher affinity binder has the most energetically favourable and pre-organized interface in the design set. A comprehensive binding-fitness landscape of this design, generated by library selections and deep sequencing, was used to optimize its binding affinity to a picomolar level, and X-ray co-crystal structures of two variants show atomic-level agreement with the corresponding computational models. The optimized binder is selective for DIG over the related steroids digitoxigenin, progesterone and β-oestradiol, and this steroid binding preference can be reprogrammed by manipulation of explicitly designed hydrogen-bonding interactions. The computational design method presented here should enable the development of a new generation of biosensors, therapeutics and diagnostics.

Nature Publishing Group2013

Structure and function of the nicotinic acetylcholine receptor

Jörg Grandl

EPFL2005