Identification and validation of novel microtubule suppressors with an imidazopyridine scaffold through structure-based virtual screening and docking.

Targeting the colchicine binding site of α/ß tubulin microtubules can lead to suppression of microtubule dynamics, cell cycle arrest and apoptosis. Therefore, the development of microtubule (MT) inhibitors is considered a promising route to anticancer agents. Our approach to identify novel scaffolds as MT inhibitors depends on a 3D-structure-based pharmacophore approach and docking using three programs MOE, Autodock and BUDE (Bristol University Docking Engine) to screen a library of virtual compounds. From this work we identified the compound 7-(3-hydroxy-4-methoxy-phenyl)-3-(3-trifluoromethyl-phenyl)-6,7-dihydro-3H-imidazo[4,5-b]pyridin-5-ol (6) as a novel inhibitor scaffold. This compound inhibited several types of cancer cell proliferation at low micromolar concentrations with low toxicity. Compound 6 caused cell cycle arrest in the G2/M phase and blocked tubulin polymerization at low micromolar concentration (IC50 = 6.1 ±0.1 µM), inducing apoptosis via activation of caspase 9, increasing the level of the pro-apoptotic protein Bax and decreasing the level of the anti-apoptotic protein Bcl2. In summary, our approach identified a lead compound with potential antimitotic and antiproliferative activity.

A potential interaction between the SARS-CoV-2 spike protein and nicotinic acetylcholine receptors.

Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4ß2 and α7 subtypes and the muscle-like αßγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4ß2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4ß2 and αßγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.

De Novo Designed Peptide and Protein Hairpins Self-Assemble into Sheets and Nanoparticles.

The design and assembly of peptide-based materials has advanced considerably, leading to a variety of fibrous, sheet, and nanoparticle structures. A remaining challenge is to account for and control different possible supramolecular outcomes accessible to the same or similar peptide building blocks. Here a de novo peptide system is presented that forms nanoparticles or sheets depending on the strategic placement of a "disulfide pin" between two elements of secondary structure that drive self-assembly. Specifically, homodimerizing and homotrimerizing de novo coiled-coil α-helices are joined with a flexible linker to generate a series of linear peptides. The helices are pinned back-to-back, constraining them as hairpins by a disulfide bond placed either proximal or distal to the linker. Computational modeling indicates, and advanced microscopy shows, that the proximally pinned hairpins self-assemble into nanoparticles, whereas the distally pinned constructs form sheets. These peptides can be made synthetically or recombinantly to allow both chemical modifications and the introduction of whole protein cargoes as required.

A conserved arginine with non-conserved function is a key determinant of agonist selectivity in α7 nicotinic ACh receptors.

BACKGROUND AND PURPOSE: The α7 and α4ß2* ("*" denotes possibly assembly with another subunit) nicotinic acetylcholine receptors (nAChRs) are the most abundant nAChRs in the mammalian brain. These receptors are the most targeted nAChRs in drug discovery programmes for brain disorders. However, the development of subtype-specific agonists remains challenging due to the high degree of sequence homology and conservation of function in nAChRs. We have developed C(10) variants of cytisine, a partial agonist of α4ß2 nAChR that has been used for smoking cessation. The C(10) methyl analogue used in this study displays negligible affinity for α7 nAChR, while retaining high affinity for α4ß2 nAChR. EXPERIMENTAL APPROACH: The structural underpinning of the selectivity of 10-methylcytisine for α7 and α4ß2 nAChRs was investigated using molecular dynamic simulations, mutagenesis and whole-cell and single-channel current recordings. KEY RESULTS: We identified a conserved arginine in the ß3 strand that exhibits a non-conserved function in nAChRs. In α4ß2 nAChR, the arginine forms a salt bridge with an aspartate residue in loop B that is necessary for receptor expression, whereas in α7 nAChR, this residue is not stabilised by electrostatic interactions, making its side chain highly mobile. This lack of constrain produces steric clashes with agonists and affects the dynamics of residues involved in agonist binding and the coupling network. CONCLUSION AND IMPLICATIONS: We conclude that the high mobility of the ß3-strand arginine in the α7 nAChR influences agonist binding and possibly gating network and desensitisation. The findings have implications for rational design of subtype-selective nAChR agents.

Neuropilin-1 is a host factor for SARS-CoV-2 infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19), uses the viral spike (S) protein for host cell attachment and entry. The host protease furin cleaves the full-length precursor S glycoprotein into two associated polypeptides: S1 and S2. Cleavage of S generates a polybasic Arg-Arg-Ala-Arg carboxyl-terminal sequence on S1, which conforms to a C-end rule (CendR) motif that binds to cell surface neuropilin-1 (NRP1) and NRP2 receptors. We used x-ray crystallography and biochemical approaches to show that the S1 CendR motif directly bound NRP1. Blocking this interaction by RNA interference or selective inhibitors reduced SARS-CoV-2 entry and infectivity in cell culture. NRP1 thus serves as a host factor for SARS-CoV-2 infection and may potentially provide a therapeutic target for COVID-19.

Identification of the Initial Steps in Signal Transduction in the α4ß2 Nicotinic Receptor: Insights from Equilibrium and Nonequilibrium Simulations.

Nicotinic acetylcholine receptors (nAChRs) modulate synaptic transmission in the nervous system. These receptors have emerged as therapeutic targets in drug discovery for treating several conditions, including Alzheimer's disease, pain, and nicotine addiction. In this in silico study, we use a combination of equilibrium and nonequilibrium molecular dynamics simulations to map dynamic and structural changes induced by nicotine in the human α4ß2 nAChR. They reveal a striking pattern of communication between the extracellular binding pockets and the transmembrane domains (TMDs) and show the sequence of conformational changes associated with the initial steps in this process. We propose a general mechanism for signal transduction for Cys-loop receptors: the mechanistic steps for communication proceed firstly through loop C in the principal subunit, and are subsequently transmitted, gradually and cumulatively, to loop F of the complementary subunit, and then to the TMDs through the M2-M3 linker.

Emergence of a Thrombospondin Superfamily at the Origin of Metazoans.

Extracellular matrix (ECM) is considered central to the evolution of metazoan multicellularity; however, the repertoire of ECM proteins in nonbilaterians remains unclear. Thrombospondins (TSPs) are known to be well conserved from cnidarians to vertebrates, yet to date have been considered a unique family, principally studied for matricellular functions in vertebrates. Through searches utilizing the highly conserved C-terminal region of TSPs, we identify undisclosed new families of TSP-related proteins in metazoans, designated mega-TSP, sushi-TSP, and poriferan-TSP, each with a distinctive phylogenetic distribution. These proteins share the TSP C-terminal region domain architecture, as determined by domain composition and analysis of molecular models against known structures. Mega-TSPs, the only form identified in ctenophores, are typically >2,700 aa and are also characterized by N-terminal leucine-rich repeats and central cadherin/immunoglobulin domains. In cnidarians, which have a well-defined ECM, Mega-TSP was expressed throughout embryogenesis in Nematostella vectensis, with dynamic endodermal expression in larvae and primary polyps and widespread ectodermal expression in adult Nematostella vectensis and Hydra magnipapillata polyps. Hydra Mega-TSP was also expressed during regeneration and siRNA-silencing of Mega-TSP in Hydra caused specific blockade of head regeneration. Molecular phylogenetic analyses based on the conserved TSP C-terminal region identified each of the TSP-related groups to form clades distinct from the canonical TSPs. We discuss models for the evolution of the newly defined TSP superfamily by gene duplications, radiation, and gene losses from a debut in the last metazoan common ancestor. Together, the data provide new insight into the evolution of ECM and tissue organization in metazoans.

Antiproliferative and Antimigratory Effects of a Novel YAP-TEAD Interaction Inhibitor Identified Using in Silico Molecular Docking.

The Hippo pathway is an important regulator of cell growth, proliferation, and migration. TEAD transcription factors, which lie at the core of the Hippo pathway, are essential for regulation of organ growth and wound repair. Dysregulation of TEAD and its regulatory cofactor Yes-associated protein (YAP) have been implicated in numerous human cancers and hyperproliferative pathological processes. Hence, the YAP-TEAD complex is a promising therapeutic target. Here, we use in silico molecular docking using Bristol University Docking Engine to screen a library of more than 8 million druglike molecules for novel disrupters of the YAP-TEAD interaction. We report the identification of a novel compound (CPD3.1) with the ability to disrupt YAP-TEAD protein-protein interaction and inhibit TEAD activity, cell proliferation, and cell migration. The YAP-TEAD complex is a viable drug target, and CPD3.1 is a lead compound for the development of more potent TEAD inhibitors for treating cancer and other hyperproliferative pathologies.

Intraring allostery controls the function and assembly of a hetero-oligomeric class II chaperonin.

Class II chaperonins are essential multisubunit complexes that aid the folding of nonnative proteins in the cytosol of archaea and eukarya. They use energy derived from ATP to drive a series of structural rearrangements that enable polypeptides to fold within their central cavity. These events are regulated by an elaborate allosteric mechanism in need of elucidation. We employed mutagenesis and experimental analysis in concert with in silico molecular dynamics simulations and interface-binding energy calculations to investigate the class II chaperonin from Thermoplasma acidophilum. Here we describe the effects on the asymmetric allosteric mechanism and on hetero-oligomeric complex formation in a panel of mutants in the ATP-binding pocket of the α and ß subunits. Our observations reveal a potential model for a nonconcerted folding mechanism optimized for protecting and refolding a range of nonnative substrates under different environmental conditions, starting to unravel the role of subunit heterogeneity in this folding machine and establishing important links with the behavior of the most complex eukaryotic chaperonins.-Shoemark, D. K., Sessions, R. B., Brancaccio, A., Bigotti, M. G. Intraring allostery controls the function and assembly of a hetero-oligomeric class II chaperonin.

Identification of novel small molecule TGF-ß antagonists using structure-based drug design.

Aberrant transforming growth factor-ß (TGF-ß) signalling has been associated with a number of disease pathologies, such as the development of fibrosis in the heart, lung and liver, cardiovascular disease and cancer, hence the TGF-ß pathway represents a promising target for a variety of diseases. However, highly specific ways to inhibit TGF-ß signalling need to be developed to prevent cross-talk with related receptors and minimise unwanted side effects. We have used used virtual screening and molecular docking to identify small molecule inhibitors of TGF-ß binding to TßRII. The crystal structure of TGF-ß3 in complex with the extracellular domain of the type II TGF-ß receptor was taken as a starting point for molecular docking and we developed a structure-based pharmacophore model to identify compounds that competitively inhibit the binding of TGF-ß to TßRII and antogonize TGF-ß signalling. We have experimentally tested 67 molecules suggested by in silico screening and similarity searching for their ability to inhibit TGF-ß signalling in TGF-ß-dependent luciferase assays in vitro and the molecule with the strongest inhibition had an IC50 of 18 µM. These compounds were selected to bind to the SS1 subsite (composed of F30, C31, D32, I50, T51 S52, I53, C54 and E55) of TßRII and all share the general property of being aromatic and fairly flat. Molecular dynamics simulations confirmed that this was the most likely binding mode. The computational methods used and the hits identified in this study provide an excellent guide to medicinal chemistry efforts to design tighter binding molecules to disrupt the TGF-ß/TßRII interaction.

Structure of bacterial glutathione-S-transferase maleyl pyruvate isomerase and implications for mechanism of isomerisation.

Maleyl pyruvate isomerase (MPI) is a bacterial glutathione S-transferase (GST) from the pathway for degradation of naphthalene via gentisate that enables the bacterium Ralstonia to use polyaromatic hydrocarbons as a sole carbon source. Genome sequencing projects have revealed the presence of large numbers of GSTs in bacterial genomes, often located within gene clusters encoding the degradation of different aromatic compounds. This structure is therefore an example of this under-represented class of enzymes. Unlike many glutathione transferases, the reaction catalysed by MPI is an isomerisation of an aromatic ring breakdown product, and glutathione is a true cofactor rather than a substrate in the reaction. We have solved the structure of the enzyme in complex with dicarboxyethyl glutathione, an analogue of a proposed reaction intermediate, at a resolution of 1.3 A. The structure provides direct evidence that the glutathione thiolate attacks the substrate in the C2 position, with the terminal carboxylate buried at the base of the active site cleft. Our structures suggest that the C1-C2 bond remains fixed so when rotation occurs around the C2-C3 bond the atoms from C4 onwards actually move. We identified a conserved arginine that is likely to stabilize the enolate form of the substrate during the isomerisation. Arginines at either side of the active site cleft can interact with the end of the substrate/product and preferentially stabilise the product. MPI has significant sequence similarity to maleylacetoacetate isomerase (MAAI), which performs an analogous reaction in the catabolism of phenylalanine and tyrosine. The proposed mechanism therefore has relevance to the MAAIs. Significantly, whilst the overall sequence identity is 40% only one of the five residues from the Zeta motif in the active site is conserved. We re-examined the roles of the residues in the active site of both enzymes and the Zeta motif itself.