Mapping the druggable allosteric space of G-protein coupled receptors: a fragment-based molecular dynamics approach.

To address the problem of specificity in G-protein coupled receptor (GPCR) drug discovery, there has been tremendous recent interest in allosteric drugs that bind at sites topographically distinct from the orthosteric site. Unfortunately, structure-based drug design of allosteric GPCR ligands has been frustrated by the paucity of structural data for allosteric binding sites, making a strong case for predictive computational methods. In this work, we map the surfaces of the beta1 (beta1AR) and beta2 (beta2AR) adrenergic receptor structures to detect a series of five potentially druggable allosteric sites. We employ the FTMAP algorithm to identify 'hot spots' with affinity for a variety of organic probe molecules corresponding to drug fragments. Our work is distinguished by an ensemble-based approach, whereby we map diverse receptor conformations taken from molecular dynamics (MD) simulations totaling approximately 0.5 micros. Our results reveal distinct pockets formed at both solvent-exposed and lipid-exposed cavities, which we interpret in light of experimental data and which may constitute novel targets for GPCR drug discovery. This mapping data can now serve to drive a combination of fragment-based and virtual screening approaches for the discovery of small molecules that bind at these sites and which may offer highly selective therapies

Cystic Fibrosis Transmembrane Conductance Regulator: Using Differential Reactivity toward Channel-Permeant and Channel-Impermeant Thiol-Reactive Probes To Test a Molecular Model for the Pore†

ABSTRACT: The sixth transmembrane segment (TM6) of the CFTR chloride channel has been intensively investigated. The effects of amino acid substitutions and chemical modification of engineered cysteines (cysteine scanning) on channel properties strongly suggest that TM6 is a key component of the anion-conducting pore, but previous cysteine-scanning studies of TM6 have produced conflicting results. Our aim was to resolve these conflicts by combining a screening strategy based on multiple, thiol-directed probes with molecular modeling of the pore. CFTR constructs were screened for reactivity toward both channel-permeant and channel-impermeant thiol-directed reagents, and patterns of reactivity in TM6 were mapped onto two new, molecular models of the CFTR pore: one based on homology modeling using Sav1866 as the template and a second derived from the first by molecular dynamics simulation. Comparison of the pattern of cysteine reactivity with model predictions suggests that nonreactive sites are those where the TM6 side chains are occluded by other TMs. Reactive sites, in contrast, are generally situated such that the respective amino acid side chains either project into the predicted pore or lie within a predicted extracellular loop. Sites wher

Cell cycle alterations in the mussel Mytilus galloprovincialis hemocytes caused by environmental contamination

Environmental contamination includes a mixture of organic substances that can have detrimental effects on marine organisms and should be evaluated in the quality and risk assessment of investigated marine areas. Marine areas selected for this study are a protected area, a mariculture area, a shipyard and an industrial area. Based on the toxicity of the organic seawater extracts these areas were classified as an undisturbed reference area (S1), an area with the low anthropogenic impact (S2), a potentially endangered area (S3) and an area with high anthropogenic impact (S4) respectively. The organic mixtures present in seawater samples collected at the above defined areas were tested for the induction of DNA damage and cell cycle alterations in the mussel Mytilus galloprovincialis hemocytes. Flow cytometric analyses were performed to detect changes in hemocytes DNA content distribution throughout the cell cycle. Organic seawater extracts from sampling sites S2, S3 and S4 induced an increase in the coefficient of variation of the G0/G1 peak and an increase in the number of cells in the G2/M phase reflecting the extent of DNA damage and G2/M arrest, respectively. The G2/M arrest in mussel hemocytes was concentration-dependent upon injection with organic seawater extracts from the S3 site and time dependant for S2, S3 and S4 sampling sites. The time dependence of the induction of the G2/M arrest showed a characteristic pattern for each site due to the different quantitative and qualitative composition of the organic seawater extracts. The G2/M arrest was reversible 24 or 72 hours after treatment with organic seawater extracts from S2 or S3, and S4 sites, respectively. This reversibility was time- and site-specific indicating that such DNA damage is repairable to a certain degree according to the organic seawater extract composition. Thus, the hemocytes cell cycle alterations in the mussel Mytilus galloprovincialis caused by organic seawater extracts reliably reflect the extent of organic contamination effects for selected marine areas

In silico assessment of potential druggable pockets on the surface of α1-Antitrypsin conformers

The search for druggable pockets on the surface of a protein is often performed on a single conformer, treated as a rigid body. Transient druggable pockets may be missed in this approach. Here, we describe a methodology for systematic in silico analysis of surface clefts across multiple conformers of the metastable protein α1-antitrypsin (A1AT). Pathological mutations disturb the conformational landscape of A1AT, triggering polymerisation that leads to emphysema and hepatic cirrhosis. Computational screens for small molecule inhibitors of polymerisation have generally focused on one major druggable site visible in all crystal structures of native A1AT. In an alternative approach, we scan all surface clefts observed in crystal structures of A1AT and in 100 computationally produced conformers, mimicking the native solution ensemble. We assess the persistence, variability and druggability of these pockets. Finally, we employ molecular docking using publicly available libraries of small molecules to explore scaffold preferences for each site. Our approach identifies a number of novel target sites for drug design. In particular one transient site shows favourable characteristics for druggability due to high enclosure and hydrophobicity. Hits against this and other druggable sites achieve docking scores corresponding to a Kd in the µM–nM range, comparing favourably with a recently identified promising lead. Preliminary ThermoFluor studies support the docking predictions. In conclusion, our strategy shows considerable promise compared with the conventional single pocket/single conformer approach to in silico screening. Our best-scoring ligands warrant further experimental investigation

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

A combinatorial TIR1/AFB–Aux/IAA co-receptor system for differential sensing of auxin

The plant hormone auxin regulates virtually every aspect of plant growth and development. Auxin acts by binding the F-box protein transport inhibitor response 1 (TIR1) and promotes the degradation of the AUXIN/INDOLE-3-ACETIC ACID (Aux/IAA) transcriptional repressors. Here we show that efficient auxin binding requires assembly of an auxin co-receptor complex consisting of TIR1 and an Aux/IAA protein. Heterologous experiments in yeast and quantitative IAA binding assays using purified proteins showed that different combinations of TIR1 and Aux/IAA proteins form co-receptor complexes with a wide range of auxin-binding affinities. Auxin affinity seems to be largely determined by the Aux/IAA. As there are 6 TIR1/AUXIN SIGNALING F-BOX proteins (AFBs) and 29 Aux/IAA proteins in Arabidopsis thaliana, combinatorial interactions may result in many co-receptors with distinct auxin-sensing properties. We also demonstrate that the AFB5–Aux/IAA co-receptor selectively binds the auxinic herbicide picloram. This co-receptor system broadens the effective concentration range of the hormone and may contribute to the complexity of auxin response

Introducing Catastrophe-QSAR. Application on Modeling Molecular Mechanisms of Pyridinone Derivative-Type HIV Non-Nucleoside Reverse Transcriptase Inhibitors

The classical method of quantitative structure-activity relationships (QSAR) is enriched using non-linear models, as Thom’s polynomials allow either uni- or bi-variate structural parameters. In this context, catastrophe QSAR algorithms are applied to the anti-HIV-1 activity of pyridinone derivatives. This requires calculation of the so-called relative statistical power and of its minimum principle in various QSAR models. A new index, known as a statistical relative power, is constructed as an Euclidian measure for the combined ratio of the Pearson correlation to algebraic correlation, with normalized t-Student and the Fisher tests. First and second order inter-model paths are considered for mono-variate catastrophes, whereas for bi-variate catastrophes the direct minimum path is provided, allowing the QSAR models to be tested for predictive purposes. At this stage, the max-to-min hierarchies of the tested models allow the interaction mechanism to be identified using structural parameter succession and the typical catastrophes involved. Minimized differences between these catastrophe models in the common structurally influential domains that span both the trial and tested compounds identify the “optimal molecular structural domains” and the molecules with the best output with respect to the modeled activity, which in this case is human immunodeficiency virus type 1 HIV-1 inhibition. The best molecules are characterized by hydrophobic interactions with the HIV-1 p66 subunit protein, and they concur with those identified in other 3D-QSAR analyses. Moreover, the importance of aromatic ring stacking interactions for increasing the binding affinity of the inhibitor-reverse transcriptase ligand-substrate complex is highlighted

Control of actin polymerization via the coincidence of phosphoinositides and high membrane curvature

The conditional use of actin during clathrin-mediated endocytosis in mammalian cells suggests that the cell controls whether and how actin is used. Using a combination of biochemical reconstitution and mammalian cell culture, we elucidate a mechanism by which the coincidence of PI(4,5)P2 and PI(3)P in a curved vesicle triggers actin polymerization. At clathrin-coated pits, PI(3)P is produced by the INPP4A hydrolysis of PI(3,4)P2, and this is necessary for actin-driven endocytosis. Both Cdc42⋅guanosine triphosphate and SNX9 activate N-WASP–WIP- and Arp2/3-mediated actin nucleation. Membrane curvature, PI(4,5)P2, and PI(3)P signals are needed for SNX9 assembly via its PX–BAR domain, whereas signaling through Cdc42 is activated by PI(4,5)P2 alone. INPP4A activity is stimulated by high membrane curvature and synergizes with SNX9 BAR domain binding in a process we call curvature cascade amplification. We show that the SNX9-driven actin comets that arise on human disease–associated oculocerebrorenal syndrome of Lowe (OCRL) deficiencies are reduced by inhibiting PI(3)P production, suggesting PI(3)P kinase inhibitors as a therapeutic strategy in Lowe syndrome.J.L. Gallop is supported by a Wellcome Trust Research Career Development Fellowship (grant WT095829AIA). F.  Daste, A.  Walrant, J.R. Gadsby, and J. Mason are supported by an H2020 European Research Council Starting Grant (281971) awarded to J.L. Gallop. Gurdon Institute funding is provided by the Wellcome Trust (grant 092096) and Cancer Research UK (grant C6946/A14492). The Swedish Medical Research Council and the Swedish Foundation for Strategic Research supported the work of M.R. Holst and R. Lundmark. S.F. Lee is funded by a Royal Society University Research Fellowship (grant UF120277). M. Mettlen is funded by grant MH73125 to Sandra L. Schmid (University of Texas Southwestern Medical Center)

Inter-domain Communication Mechanisms in an ABC Importer: A Molecular Dynamics Study of the MalFGK2E Complex

ATP-Binding Cassette transporters are ubiquitous membrane proteins that convert the energy from ATP-binding and hydrolysis into conformational changes of the transmembrane region to allow the translocation of substrates against their concentration gradient. Despite the large amount of structural and biochemical data available for this family, it is still not clear how the energy obtained from ATP hydrolysis in the ATPase domains is “transmitted” to the transmembrane domains. In this work, we focus our attention on the consequences of hydrolysis and inorganic phosphate exit in the maltose uptake system (MalFGK2E) from Escherichia coli. The prime goal is to identify and map the structural changes occurring during an ATP-hydrolytic cycle. For that, we use extensive molecular dynamics simulations to study three potential intermediate states (with 10 replicates each): an ATP-bound, an ADP plus inorganic phosphate-bound and an ADP-bound state. Our results show that the residues presenting major rearrangements are located in the A-loop, in the helical sub-domain, and in the “EAA motif” (especially in the “coupling helices” region). Additionally, in one of the simulations with ADP we were able to observe the opening of the NBD dimer accompanied by the dissociation of ADP from the ABC signature motif, but not from its corresponding P-loop motif. This work, together with several other MD studies, suggests a common communication mechanism both for importers and exporters, in which ATP-hydrolysis induces conformational changes in the helical sub-domain region, in turn transferred to the transmembrane domains via the “coupling helices”