Structural basis for modulation of a G-protein-coupled receptor by allosteric drugs.

The design of G-protein-coupled receptor (GPCR) allosteric modulators, an active area of modern pharmaceutical research, has proved challenging because neither the binding modes nor the molecular mechanisms of such drugs are known. Here we determine binding sites, bound conformations and specific drug-receptor interactions for several allosteric modulators of the M2 muscarinic acetylcholine receptor (M2 receptor), a prototypical family A GPCR, using atomic-level simulations in which the modulators spontaneously associate with the receptor. Despite substantial structural diversity, all modulators form cation-π interactions with clusters of aromatic residues in the receptor extracellular vestibule, approximately 15 Å from the classical, 'orthosteric' ligand-binding site. We validate the observed modulator binding modes through radioligand binding experiments on receptor mutants designed, on the basis of our simulations, either to increase or to decrease modulator affinity. Simulations also revealed mechanisms that contribute to positive and negative allosteric modulation of classical ligand binding, including coupled conformational changes of the two binding sites and electrostatic interactions between ligands in these sites. These observations enabled the design of chemical modifications that substantially alter a modulator's allosteric effects. Our findings thus provide a structural basis for the rational design of allosteric modulators targeting muscarinic and possibly other GPCRs.

Pathway and mechanism of drug binding to G-protein-coupled receptors.

How drugs bind to their receptors--from initial association, through drug entry into the binding pocket, to adoption of the final bound conformation, or "pose"--has remained unknown, even for G-protein-coupled receptor modulators, which constitute one-third of all marketed drugs. We captured this pharmaceutically critical process in atomic detail using the first unbiased molecular dynamics simulations in which drug molecules spontaneously associate with G-protein-coupled receptors to achieve final poses matching those determined crystallographically. We found that several beta blockers and a beta agonist all traverse the same well-defined, dominant pathway as they bind to the ß(1)- and ß(2)-adrenergic receptors, initially making contact with a vestibule on each receptor's extracellular surface. Surprisingly, association with this vestibule, at a distance of 15 Å from the binding pocket, often presents the largest energetic barrier to binding, despite the fact that subsequent entry into the binding pocket requires the receptor to deform and the drug to squeeze through a narrow passage. The early barrier appears to reflect the substantial dehydration that takes place as the drug associates with the vestibule. Our atomic-level description of the binding process suggests opportunities for allosteric modulation and provides a structural foundation for future optimization of drug-receptor binding and unbinding rates.

Activation mechanism of the ß2-adrenergic receptor.

A third of marketed drugs act by binding to a G-protein-coupled receptor (GPCR) and either triggering or preventing receptor activation. Although recent crystal structures have provided snapshots of both active and inactive functional states of GPCRs, these structures do not reveal the mechanism by which GPCRs transition between these states. Here we propose an activation mechanism for the ß(2)-adrenergic receptor, a prototypical GPCR, based on atomic-level simulations in which an agonist-bound receptor transitions spontaneously from the active to the inactive crystallographically observed conformation. A loosely coupled allosteric network, comprising three regions that can each switch individually between multiple distinct conformations, links small perturbations at the extracellular drug-binding site to large conformational changes at the intracellular G-protein-binding site. Our simulations also exhibit an intermediate that may represent a receptor conformation to which a G protein binds during activation, and suggest that the first structural changes during receptor activation often take place on the intracellular side of the receptor, far from the drug-binding site. By capturing this fundamental signaling process in atomic detail, our results may provide a foundation for the design of drugs that control receptor signaling more precisely by stabilizing specific receptor conformations.

Principles of conduction and hydrophobic gating in K+ channels.

We present the first atomic-resolution observations of permeation and gating in a K(+) channel, based on molecular dynamics simulations of the Kv1.2 pore domain. Analysis of hundreds of simulated permeation events revealed a detailed conduction mechanism, resembling the Hodgkin-Keynes "knock-on" model, in which translocation of two selectivity filter-bound ions is driven by a third ion; formation of this knock-on intermediate is rate determining. In addition, at reverse or zero voltages, we observed pore closure by a novel "hydrophobic gating" mechanism: A dewetting transition of the hydrophobic pore cavity-fastest when K(+) was not bound in selectivity filter sites nearest the cavity-caused the open, conducting pore to collapse into a closed, nonconducting conformation. Such pore closure corroborates the idea that voltage sensors can act to prevent pore collapse into the intrinsically more stable, closed conformation, and it further suggests that molecular-scale dewetting facilitates a specific biological function: K(+) channel gating. Existing experimental data support our hypothesis that hydrophobic gating may be a fundamental principle underlying the gating of voltage-sensitive K(+) channels. We suggest that hydrophobic gating explains, in part, why diverse ion channels conserve hydrophobic pore cavities, and we speculate that modulation of cavity hydration could enable structural determination of both open and closed channels.

Enabling structure-based drug design of Tyk2 through co-crystallization with a stabilizing aminoindazole inhibitor.

BACKGROUND: Structure-based drug design (SBDD) can accelerate inhibitor lead design and optimization, and efficient methods including protein purification, characterization, crystallization, and high-resolution diffraction are all needed for rapid, iterative structure determination. Janus kinases are important targets that are amenable to structure-based drug design. Here we present the first mouse Tyk2 crystal structures, which are complexed to 3-aminoindazole compounds. RESULTS: A comprehensive construct design effort included N- and C-terminal variations, kinase-inactive mutations, and multiple species orthologs. High-throughput cloning and expression methods were coupled with an abbreviated purification protocol to optimize protein solubility and stability. In total, 50 Tyk2 constructs were generated. Many displayed poor expression, inadequate solubility, or incomplete affinity tag processing. One kinase-inactive murine Tyk2 construct, complexed with an ATP-competitive 3-aminoindazole inhibitor, provided crystals that diffracted to 2.5-2.6 Å resolution. This structure revealed initial "hot-spot" regions for SBDD, and provided a robust platform for ligand soaking experiments. Compared to previously reported human Tyk2 inhibitor crystal structures (Chrencik et al. (2010) J Mol Biol 400:413), our structures revealed a key difference in the glycine-rich loop conformation that is induced by the inhibitor. Ligand binding also conferred resistance to proteolytic degradation by thermolysin. As crystals could not be obtained with the unliganded enzyme, this enhanced stability is likely important for successful crystallization and inhibitor soaking methods. CONCLUSIONS: Practical criteria for construct performance and prioritization, the optimization of purification protocols to enhance protein yields and stability, and use of high-throughput construct exploration enable structure determination methods early in the drug discovery process. Additionally, specific ligands stabilize Tyk2 protein and may thereby enable crystallization.

The future of molecular dynamics simulations in drug discovery.

Molecular dynamics simulations can now track rapid processes--those occurring in less than about a millisecond--at atomic resolution for many biologically relevant systems. These simulations appear poised to exert a significant impact on how new drugs are found, perhaps even transforming the very process of drug discovery. We predict here future results we can expect from, and enhancements we need to make in, molecular dynamics simulations over the coming 25 years, and in so doing set out several Grand Challenges for the field. In the context of the problems now facing the pharmaceutical industry, we ask how we can best address drug discovery needs of the next quarter century using molecular dynamics simulations, and we suggest some possible approaches.

Identification of two distinct inactive conformations of the beta2-adrenergic receptor reconciles structural and biochemical observations.

Fully understanding the mechanisms of signaling proteins such as G protein-coupled receptors (GPCRs) will require the characterization of their conformational states and the pathways connecting those states. The recent crystal structures of the beta(2)- and beta(1)-adrenergic receptors in a nominally inactive state constituted a major advance toward this goal, but also raised new questions. Although earlier biochemical observations had suggested that these receptors possessed a set of contacts between helices 3 and 6, known as the ionic lock, which was believed to form a molecular switch for receptor activation, the crystal structures lacked these contacts. The unexpectedly broken ionic lock has raised questions about the true conformation(s) of the inactive state and the role of the ionic lock in receptor activation and signaling. To address these questions, we performed microsecond-timescale molecular dynamics simulations of the beta(2)-adrenergic receptor (beta(2)AR) in multiple wild-type and mutant forms. In wild-type simulations, the ionic lock formed reproducibly, bringing the intracellular ends of helices 3 and 6 together to adopt a conformation similar to that found in inactive rhodopsin. Our results suggest that inactive beta(2)AR exists in equilibrium between conformations with the lock formed and the lock broken, whether or not the cocrystallized ligand is present. These findings, along with the formation of several secondary structural elements in the beta(2)AR loops during our simulations, may provide a more comprehensive picture of the inactive state of the beta-adrenergic receptors, reconciling the crystal structures with biochemical studies.

Dynamic control of slow water transport by aquaporin 0: implications for hydration and junction stability in the eye lens.

Aquaporin 0 (AQP0), the most abundant membrane protein in mammalian lens fiber cells, not only serves as the primary water channel in this tissue but also appears to mediate the formation of thin junctions between fiber cells. AQP0 is remarkably less water permeable than other aquaporins, but the structural basis and biological significance of this low permeability remain uncertain, as does the permeability of the protein in a reported junctional form. To address these issues, we performed molecular dynamics (MD) simulations of water transport through membrane-embedded AQP0 in both its (octameric) junctional and (tetrameric) nonjunctional forms. From our simulations, we measured an osmotic permeability for the nonjunctional form that agrees with experiment and found that the distinct dynamics of the conserved, lumen-protruding side chains of Tyr-23 and Tyr-149 modulate water passage, accounting for the slow permeation. The junctional and nonjunctional forms conducted water equivalently, in contrast to a previous suggestion based on static crystal structures that water conduction is lost on junction formation. Our analysis suggests that the low water permeability of AQP0 may help maintain the mechanical stability of the junction. We hypothesize that the structural features leading to low permeability may have evolved in part to allow AQP0 to form junctions that both conduct water and contribute to the organizational structure of the fiber cell tissue and microcirculation within it, as required to maintain transparency of the lens.

2,4-Diaminopyrimidine MK2 inhibitors. Part I: Observation of an unexpected inhibitor binding mode.

MK2 is a Ser/Thr kinase of significant interest as an anti-inflammatory drug discovery target. Here we describe the development of in vitro tools for the identification and characterization of MK2 inhibitors, including validation of inhibitor interactions with the crystallography construct and determination of the unique binding mode of 2,4-diaminopyrimidine inhibitors in the MK2 active site. Use of these tools in the optimization of a potent and selective inhibitor lead series is described in the accompanying Letter.

2,4-Diaminopyrimidine MK2 inhibitors. Part II: Structure-based inhibitor optimization.

We describe structure-based optimization of a series of novel 2,4-diaminopyrimidine MK2 inhibitors. Co-crystal structures (see accompanying Letter) demonstrated a unique inhibitor binding mode. Resulting inhibitors had IC(50) values as low as 19nM and moderate selectivity against a kinase panel. Compounds 15, 31a, and 31b inhibit TNFalpha production in peripheral human monocytes.

Rational mutagenesis to support structure-based drug design: MAPKAP kinase 2 as a case study.

BACKGROUND: Structure-based drug design (SBDD) can provide valuable guidance to drug discovery programs. Robust construct design and expression, protein purification and characterization, protein crystallization, and high-resolution diffraction are all needed for rapid, iterative inhibitor design. We describe here robust methods to support SBDD on an oral anti-cytokine drug target, human MAPKAP kinase 2 (MK2). Our goal was to obtain useful diffraction data with a large number of chemically diverse lead compounds. Although MK2 structures and structural methods have been reported previously, reproducibility was low and improved methods were needed. RESULTS: Our construct design strategy had four tactics: N- and C-terminal variations; entropy-reducing surface mutations; activation loop deletions; and pseudoactivation mutations. Generic, high-throughput methods for cloning and expression were coupled with automated liquid dispensing for the rapid testing of crystallization conditions with minimal sample requirements. Initial results led to development of a novel, customized robotic crystallization screen that yielded MK2/inhibitor complex crystals under many conditions in seven crystal forms. In all, 44 MK2 constructs were generated, ~500 crystals were tested for diffraction, and ~30 structures were determined, delivering high-impact structural data to support our MK2 drug design effort. CONCLUSION: Key lessons included setting reasonable criteria for construct performance and prioritization, a willingness to design and use customized crystallization screens, and, crucially, initiation of high-throughput construct exploration very early in the drug discovery process.

Novel boron-containing, nonclassical antifolates: synthesis and preliminary biological and structural evaluation.

Two boron-containing, ortho-icosahedral carborane lipophilic antifolates were synthesized, and the crystal structures of their ternary complexes with human dihydrofolate reductase (DHFR) and dihydronicotinamide adenine dinucleotide phosphate were determined. The compounds were screened for activity against DHFR from six sources (human, rat liver, Pneumocystis carinii, Toxoplasma gondii, Mycobacterium avium, and Lactobacillus casei) and showed good to modest activity against these enzymes. The compounds were also tested for antibacterial activity against L. casei, M. tuberculosis H37Ra, and three M. avium strains and for cytotoxic activity against seven different human tumor cell lines. Antibacterial and cytotoxic activity was modest, with one sample, the closo-carborane 4, showing about 10-fold greater activity. The less toxic nido-carborane 2 was also tested as a candidate for boron neutron capture therapy, but showed poor tumor retention and low selectivity ratios for boron distribution in tumor tissue versus normal tissue.

Structure of Mycobacterium tuberculosis FtsZ reveals unexpected, G protein-like conformational switches.

We report three crystal structures of the Mycobacterium tuberculosis cell division protein FtsZ, as the citrate, GDP, and GTPgammaS complexes, determined at 1.89, 2.60, and 2.08A resolution. MtbFtsZ crystallized as a tight, laterally oriented dimer distinct from the longitudinal polymer observed for alphabeta-tubulin. Mutational data on Escherichia coli FtsZ suggest that this dimer interface is important for proper protofilament and "Z-ring" assembly and function. An alpha-to-beta secondary structure conformational switch at the dimer interface is spatially analogous to, and has many of the hallmarks of, the Switch I conformational changes exhibited by G-proteins upon activation. The presence of a gamma-phosphate in the FtsZ active site modulates the conformation of the "tubulin" loop T3 (spatially analogous to the G-protein Switch II); T3 switching upon gamma-phosphate ligation is directly coupled to the alpha-to-beta switch by steric overlap. The dual conformational switches observed here for the first time in an FtsZ link GTP binding and hydrolysis to FtsZ (and tubulin) lateral assembly and Z-ring contraction, and they are suggestive of an underappreciated functional analogy between FtsZ, tubulin and G-proteins.

Atomic structures of human dihydrofolate reductase complexed with NADPH and two lipophilic antifolates at 1.09 a and 1.05 a resolution.

The crystal structures of two human dihydrofolate reductase (hDHFR) ternary complexes, each with bound NADPH cofactor and a lipophilic antifolate inhibitor, have been determined at atomic resolution. The potent inhibitors 6-([5-quinolylamino]methyl)-2,4-diamino-5-methylpyrido[2,3-d]pyrimidine (SRI-9439) and (Z)-6-(2-[2,5-dimethoxyphenyl]ethen-1-yl)-2,4-diamino-5-methylpyrido[2,3-d]pyrimidine (SRI-9662) were developed at Southern Research Institute against Toxoplasma gondii DHFR-thymidylate synthase. The 5-deazapteridine ring of each inhibitor adopts an unusual puckered conformation that enables the formation of identical contacts in the active site. Conversely, the quinoline and dimethoxybenzene moieties exhibit distinct binding characteristics that account for the differences in inhibitory activity. In both structures, a salt-bridge is formed between Arg70 in the active site and Glu44 from a symmetry-related molecule in the crystal lattice that mimics the binding of methotrexate to DHFR.

SIGNAL TRANSDUCTION. Structural basis for nucleotide exchange in heterotrimeric G proteins.

G protein-coupled receptors (GPCRs) relay diverse extracellular signals into cells by catalyzing nucleotide release from heterotrimeric G proteins, but the mechanism underlying this quintessential molecular signaling event has remained unclear. Here we use atomic-level simulations to elucidate the nucleotide-release mechanism. We find that the G protein α subunit Ras and helical domains-previously observed to separate widely upon receptor binding to expose the nucleotide-binding site-separate spontaneously and frequently even in the absence of a receptor. Domain separation is necessary but not sufficient for rapid nucleotide release. Rather, receptors catalyze nucleotide release by favoring an internal structural rearrangement of the Ras domain that weakens its nucleotide affinity. We use double electron-electron resonance spectroscopy and protein engineering to confirm predictions of our computationally determined mechanism.

Molecular determinants of drug-receptor binding kinetics.

It is increasingly appreciated that the rates at which drugs associate with and dissociate from receptors--the binding kinetics--directly impact drug efficacy and safety. The molecular determinants of drug-receptor binding kinetics remain poorly understood, however, especially when compared with the well-known factors that affect binding affinity. The rational modulation of kinetics during lead optimization thus remains challenging. We review some of the key factors thought to control drug-receptor binding kinetics at the molecular level--molecular size, conformational fluctuations, electrostatic interactions and hydrophobic effects--and discuss several possible approaches for the rational design of drugs with desired binding kinetics.

Mechanism of voltage gating in potassium channels.

The mechanism of ion channel voltage gating-how channels open and close in response to voltage changes-has been debated since Hodgkin and Huxley's seminal discovery that the crux of nerve conduction is ion flow across cellular membranes. Using all-atom molecular dynamics simulations, we show how a voltage-gated potassium channel (KV) switches between activated and deactivated states. On deactivation, pore hydrophobic collapse rapidly halts ion flow. Subsequent voltage-sensing domain (VSD) relaxation, including inward, 15-angstrom S4-helix motion, completes the transition. On activation, outward S4 motion tightens the VSD-pore linker, perturbing linker-S6-helix packing. Fluctuations allow water, then potassium ions, to reenter the pore; linker-S6 repacking stabilizes the open pore. We propose a mechanistic model for the sodium/potassium/calcium voltage-gated ion channel superfamily that reconciles apparently conflicting experimental data.

Comment on "A bacterium that can grow by using arsenic instead of phosphorus".

Wolfe-Simon et al. (Research Articles, 3 June 2011, p. 1163; published online 2 December 2010) reported the discovery of an unusual bacterium, strain GFAJ-1, that can grow in the presence of high concentrations of arsenate. The authors' contention, however, that this microbe can appreciably vary the elemental composition of its fundamental biomolecules by substituting arsenic for phosphorus appears premature based on the data presented.