Activation and dynamic network of the M2 muscarinic receptor.

G-protein-coupled receptors (GPCRs) mediate cellular responses to various hormones and neurotransmitters and are important targets for treating a wide spectrum of diseases. Although significant advances have been made in structural studies of GPCRs, details of their activation mechanism remain unclear. The X-ray crystal structure of the M2 muscarinic receptor, a key GPCR that regulates human heart rate and contractile forces of cardiomyocytes, was determined recently in an inactive antagonist-bound state. Here, activation of the M2 receptor is directly observed via accelerated molecular dynamics simulation, in contrast to previous microsecond-timescale conventional molecular dynamics simulations in which the receptor remained inactive. Receptor activation is characterized by formation of a Tyr206(5.58)-Tyr440(7.53) hydrogen bond and ∼6-Å outward tilting of the cytoplasmic end of transmembrane α-helix 6, preceded by relocation of Trp400(6.48) toward Phe195(5.47) and Val199(5.51) and flipping of Tyr430(7.43) away from the ligand-binding cavity. Network analysis reveals that communication in the intracellular domains is greatly weakened during activation of the receptor. Together with the finding that residue motions in the ligand-binding and G-protein-coupling sites of the apo receptor are correlated, this result highlights a dynamic network for allosteric regulation of the M2 receptor activation.

Free energy landscape of G-protein coupled receptors, explored by accelerated molecular dynamics.

G-protein coupled receptors (GPCRs) mediate cellular responses to various hormones and neurotransmitters and are important targets for treating a wide spectrum of diseases. They are known to adopt multiple conformational states (e.g., inactive, intermediate and active) during their modulation of various cell signaling pathways. Here, the free energy landscape of GPCRs is explored using accelerated molecular dynamics (aMD) simulations as demonstrated on the M2 muscarinic receptor, a key GPCR that regulates human heart rate and contractile forces of cardiomyocytes. Free energy profiles of important structural motifs that undergo conformational transitions upon GPCR activation and allosteric signaling are analyzed in detail, including the Arg(3.50)-Glu(6.30) ionic lock, the Trp(6.48) toggle switch and the hydrogen interactions between Tyr(5.58)-Tyr(7.53).

Effects of histidine protonation and rotameric states on virtual screening of M. tuberculosis RmlC.

While it is well established that protonation and tautomeric states of ligands can significantly affect the results of virtual screening, such effects of ionizable residues of protein receptors are less well understood. In this study, we focus on histidine protonation and rotameric states and their impact on virtual screening of Mycobacterium tuberculosis enzyme RmlC. Depending on the net charge and the location of proton(s), a histidine can adopt three states: HIP (+1 charged, both δ- and ε-nitrogens protonated), HID (neutral, δ-nitrogen protonated), and HIE (neutral, ε-nitrogen protonated). Due to common ambiguities in X-ray crystal structures, a histidine may also be resolved as three additional states with its imidazole ring flipped. Here, we systematically investigate the predictive power of 36 receptor models with different protonation and rotameric states of two histidines in the RmlC active site by using results from a previous high-throughput screening. By measuring enrichment factors and area under the receiver operating characteristic curves, we show that virtual screening results vary depending on hydrogen bonding networks provided by the histidines, even in the cases where the ligand does not obviously interact with the side chain. Our results also suggest that, even with the help of widely used pKa prediction software, assigning histidine protonation and rotameric states for virtual screening can still be challenging and requires further examination and systematic characterization of the receptor-ligand complex.

Thermodynamic integration to predict host-guest binding affinities.

An alchemical free energy method with explicit solvent molecular dynamics simulations was applied as part of the blind prediction contest SAMPL3 to calculate binding free energies for seven guests to an acyclic cucurbit-[n]uril host. The predictions included determination of protonation states for both host and guests, docking pose generation, and binding free energy calculations using thermodynamic integration. We found a root mean square error (RMSE) of 3.6 kcal mol(-1) from the reference experimental results, with an R(2) correlation of 0.51. The agreement with experiment for the largest contributor to this error, guest 6, is improved by 1.7 kcal mol(-1) when a periodicity-induced free energy correction is applied. The corrections for the other ligands were significantly smaller, and altogether the RMSE was reduced by 0.4 kcal mol(-1). We link properties of the host-guest systems during simulation to errors in the computed free energies. Overall, we show that charged host-guest systems studied here, initialized in unconfirmed docking poses, present a challenge to accurate alchemical simulation methods.

Predictive power of molecular dynamics receptor structures in virtual screening.

Molecular dynamics (MD) simulation is a well-established method for understanding protein dynamics. Conformations from unrestrained MD simulations have yet to be assessed for blind virtual screening (VS) by docking. This study presents a critical analysis of the predictive power of MD snapshots to this regard, evaluating two well-characterized systems of varying flexibility in ligand-bound and unbound configurations. Results from such VS predictions are discussed with respect to experimentally determined structures. In all cases, MD simulations provide snapshots that improve VS predictive power over known crystal structures, possibly due to sampling more relevant receptor conformations. Additionally, MD can move conformations previously not amenable to docking into the predictive range.

Mapping of allosteric druggable sites in activation-associated conformers of the M2 muscarinic receptor.

G-protein-coupled receptors (GPCRs) are key cellular signaling proteins and have been targeted by approximately 30-40% of marketed drugs for treating many human diseases including cancer and heart failure. Recently, we directly observed activation of the M2 muscarinic receptor through long-timescale accelerated molecular dynamics (aMD) simulation, which revealed distinct inactive, intermediate and active conformers of the receptor. Here, FTMAP is applied to search for 'hot spots' in these activation-associated conformers using a library of 16 organic probe molecules that represent fragments of potential drugs. Seven allosteric (non-orthosteric) binding sites are identified in the M2 receptor through the FTMAP analysis. These sites are distributed in the solvent-exposed extracellular and intracellular mouth regions, as well as the lipid-exposed pockets formed by the transmembrane α helices TM3-TM4, TM5-TM6 and TM7-TM1/TM2. They serve as promising target sites for designing novel allosteric modulators as receptor-selective drugs.

Structure-based network analysis of an evolved G protein-coupled receptor homodimer interface.

Crystallographic structures and experimental assays of human CXC chemokine receptor type 4 (CXCR4) provide strong evidence for the capacity to homodimerize, potentially as a means of allosteric regulation. Even so, how this homodimer forms and its biological significance has yet to be fully characterized. By applying principles from network analysis, sequence-based approaches such as statistical coupling analysis to determine coevolutionary residues, can be used in conjunction with molecular dynamics simulations to identify residues relevant to dimerization. Here, the predominant coevolution sector lies along the observed dimer interface, suggesting functional relevance. Furthermore, coevolution scoring provides a basis for determining significant nodes, termed hubs, in the network formed by residues found along the interface of the homodimer. These node residues coincide with hotspots indicating potential druggability. Drug design efforts targeting such key residues could potentially result in modulation of binding and therapeutic benefits for disease states, such as lung cancers, lymphomas and latent HIV-1 infection. Furthermore, this method may be applied to any protein-protein interaction.

Rational prediction with molecular dynamics for hit identification.

Although the motions of proteins are fundamental for their function, for pragmatic reasons, the consideration of protein elasticity has traditionally been neglected in drug discovery and design. This review details protein motion, its relevance to biomolecular interactions and how it can be sampled using molecular dynamics simulations. Within this context, two major areas of research in structure-based prediction that can benefit from considering protein flexibility, binding site detection and molecular docking, are discussed. Basic classification metrics and statistical analysis techniques, which can facilitate performance analysis, are also reviewed. With hardware and software advances, molecular dynamics in combination with traditional structure-based prediction methods can potentially reduce the time and costs involved in the hit identification pipeline.

On the use of molecular dynamics receptor conformations for virtual screening.

Receptors are inherently dynamic and this flexibility is important to consider when constructing a model of molecular association. Conformations from molecular dynamics simulations, a well-established method for examining protein dynamics, can be used in virtual screening to account for flexibility in structure-based drug discovery. Different receptor configurations influence docking results. Molecular dynamics simulations can provide snapshots that improve virtual screening predictive power over known crystal structures, most likely as a result of sampling more relevant receptor conformations. Here we highlight some details and nuances of using such snapshots and evaluating them for predictive performance.

Identification of potential small molecule binding pockets on Rho family GTPases.

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (>100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies. Although the inhibition of Rho signaling with drug-like compounds is an active area of current research, very little attention has been devoted to directly inhibiting Rho by targeting potential allosteric non-nucleotide binding sites. By avoiding the nucleotide binding site, compounds may minimize the potential for undesirable off-target interactions with other ubiquitous GTP and ATP binding proteins. Here we describe the application of molecular dynamics simulations, principal component analysis, sequence conservation analysis, and ensemble small-molecule fragment mapping to provide an extensive mapping of potential small-molecule binding pockets on Rho family members. Characterized sites include novel pockets in the vicinity of the conformationaly responsive switch regions as well as distal sites that appear to be related to the conformations of the nucleotide binding region. Furthermore the use of accelerated molecular dynamics simulation, an advanced sampling method that extends the accessible time-scale of conventional simulations, is found to enhance the characterization of novel binding sites when conformational changes are important for the protein mechanism.

A novel steroidal antiandrogen targeting wild type and mutant androgen receptors.

Prostate cancer (PCa) progression is enhanced by androgen and treatment with antiandrogens represents an alternative to castration. While patients initially respond favorably to androgen ablation therapy, most experience a relapse of the disease within 1-2 years by expressing androgen receptor (AR) mutants. Such mutations, indeed, promote unfavorable agonistic behavior from classical antagonists. Here, we have synthesized and screened 37 novel compounds derived from dihydrotestosterone (DHT), cyanolutamide and hydroxyflutamide. These derivatives were tested for their potential antagonistic activity using a luciferase reporter gene assay and binding properties were determined for wild type (WT) and mutant ARs (T877A, W741C, W741L, H874Y). In the absence and presence of antiandrogens, androgen dependent cellular proliferation and prostate specific antigen (PSA) expression were assayed in the prostate cancer cell line LNCaP by crystal violet, real time PCR and by Western blots. Also, cellular proliferation and PSA expression were assayed in 22Rv1. A novel compound RB346, derived from DHT, was found to be an antagonist for all tested AR forms, preventing DHT induced proliferation and PSA expression in LNCaP and 22Rv1 cells. RB346 displayed no agonistic activity, in contrast to the non-steroidal antiandrogen bicalutamide (Casodex) with unfavorable agonistic activity for W741L-AR. Additionally, RB346 has a slightly higher binding affinity for WT-AR, T877A-AR and H874Y-AR than bicalutamide. Thus, RB346 is the first potent steroidal antiandrogen with efficacy for WT and various AR mutants.

Discovery of wild-type and Y181C mutant non-nucleoside HIV-1 reverse transcriptase inhibitors using virtual screening with multiple protein structures.

To discover non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs) that are effective against both wild-type (WT) virus and variants that encode the clinically troublesome Tyr181Cys (Y181C) RT mutation, virtual screening by docking was carried out using three RT structures and more than 2 million commercially available compounds. Two of the structures are for WT-virus with different conformations of Tyr181, while the third structure incorporates the Y181C modification. Eventually nine compounds were purchased and assayed. Three of the compounds show low-micromolar antiviral activity toward either or both the wild-type and Y181C HIV-1 strains. The study illustrates a viable protocol to seek anti-HIV agents with enhanced resistance profiles.

Exploiting structural analysis, in silico screening, and serendipity to identify novel inhibitors of drug-resistant falciparum malaria.

Plasmodium falciparum thymidylate synthase-dihydrofolate reductase (TS-DHFR) is an essential enzyme in folate biosynthesis and a major malarial drug target. This bifunctional enzyme thus presents different design approaches for developing novel inhibitors against drug-resistant mutants. We performed a high-throughput in silico screen of a database of diverse, drug-like molecules against a non-active-site pocket of TS-DHFR. The top compounds from this virtual screen were evaluated by in vitro enzymatic and cellular culture studies. Three compounds active to 20 microM IC(50)'s in both wildtype and antifolate-resistant P. falciparum parasites were identified; moreover, no inhibition of human DHFR enzyme was observed, indicating that the inhibitory effects appeared to be parasite-specific. Notably, all three compounds had a biguanide scaffold. However, relative free energy of binding calculations suggested that the compounds might preferentially interact with the active site over the screened non-active-site region. To resolve the two possible modes of binding, co-crystallization studies of the compounds complexed with TS-DHFR enzyme were performed. Surprisingly, the structural analysis revealed that these novel, biguanide compounds do indeed bind at the active site of DHFR and additionally revealed the molecular basis by which they overcome drug resistance. To our knowledge, these are the first co-crystal structures of novel, biguanide, non-WR99210 compounds that are active against folate-resistant malaria parasites in cell culture.