Evolutionary conserved Tyr169 stabilizes the ß2-α2 loop of the prion protein.

Experimental evidence indicates that the primary structure of the ß2-α2 loop region (residues 165-175) in mammalian prion proteins (PrP) influences the conversion from the cellular species (PrP(C)) to the ß-sheet-rich aggregate. Here, we captured the transition of the ß2-α2 loop from 310-helical turn to ß turn by unbiased molecular dynamics simulations of the single-point mutant Y169G. Multiple conformations along the spontaneous transition of the mutant were then used as starting point for sampling of the free-energy surface of the wild type and other single-point mutants. Using two different methods for the determination of free energy profiles, we found that the barrier for the 310-helical turn to ß turn transition of the wild type is higher by about 2.5 kcal/mol than for the Y169G mutant, which is due to favorable stacking of the aromatic rings of Y169 and F175, and a stable hydrogen bond between the side chains of Y169 and D178. The transition of the ß2-α2 loop to ß turn increases the solvent-exposure of the hydrophobic stretch 169-YSNQNNF-175. The simulations indicate that the strictly conserved Y169 in mammalian prion proteins stabilizes the 310-helical turn in the ß2-α2 loop, thus hindering the conversion to an aggregation-prone conformation.

The catalytic mechanism of the RNA methyltransferase METTL3.

The complex of methyltransferase-like proteins 3 and 14 (METTL3-14) is the major enzyme that deposits N6-methyladenosine (m6A) modifications on messenger RNA (mRNA) in humans. METTL3-14 plays key roles in various biological processes through its methyltransferase (MTase) activity. However, little is known about its substrate recognition and methyl transfer mechanism from its cofactor and methyl donor S-adenosylmethionine (SAM). Here, we study the MTase mechanism of METTL3-14 by a combined experimental and multiscale simulation approach using bisubstrate analogues (BAs), conjugates of a SAM-like moiety connected to the N6-atom of adenosine. Molecular dynamics simulations based on crystal structures of METTL3-14 with BAs suggest that the Y406 side chain of METTL3 is involved in the recruitment of adenosine and release of m6A. A crystal structure with a BA representing the transition state of methyl transfer shows a direct involvement of the METTL3 side chains E481 and K513 in adenosine binding which is supported by mutational analysis. Quantum mechanics/molecular mechanics (QM/MM) free energy calculations indicate that methyl transfer occurs without prior deprotonation of adenosine-N6. Furthermore, the QM/MM calculations provide further support for the role of electrostatic contributions of E481 and K513 to catalysis. The multidisciplinary approach used here sheds light on the (co)substrate binding mechanism, catalytic step, and (co)product release, and suggests that the latter step is rate-limiting for METTL3. The atomistic information on the substrate binding and methyl transfer reaction of METTL3 can be useful for understanding the mechanisms of other RNA MTases and for the design of transition state analogues as their inhibitors.

Structure-Based Design of Inhibitors of the m6A-RNA Writer Enzyme METTL3.

Methyltransferase-like 3 (METTL3) and METTL14 form a heterodimeric complex that catalyzes the most abundant internal mRNA modification, N6-methyladenosine (m6A). METTL3 is the catalytic subunit that binds the co-substrate S-adenosyl methionine (SAM), while METTL14 is involved in mRNA binding. The m6A modification provides post-transcriptional level control over gene expression as it affects almost all stages of the mRNA life cycle, including splicing, nuclear export, translation, and decay. There is increasing evidence for an oncogenic role of METTL3 in acute myeloid leukemia. Here, we use structural and dynamic details of the catalytic subunit METTL3 for developing small-molecule inhibitors that compete with SAM. Starting from a hit identified by high-throughput docking, protein crystallography and molecular dynamics simulations were employed to guide the optimization of inhibitory activity. The potency was successfully improved by 8000-fold as measured by a homogeneous time-resolved fluorescence assay. The optimized compound is selective against the off-targets RNA methyltransferases METTL1 and METTL16.

Rapid ex vivo reverse genetics identifies the essential determinants of prion protein toxicity.

The cellular prion protein PrPC mediates the neurotoxicity of prions and other protein aggregates through poorly understood mechanisms. Antibody-derived ligands against the globular domain of PrPC (GDL) can also initiate neurotoxicity by inducing an intramolecular R208 -H140 hydrogen bond ("H-latch") between the α2-α3 and ß2-α2 loops of PrPC . Importantly, GDL that suppresses the H-latch prolong the life of prion-infected mice, suggesting that GDL toxicity and prion infections exploit convergent pathways. To define the structural underpinnings of these phenomena, we transduced 19 individual PrPC variants to PrPC -deficient cerebellar organotypic cultured slices using adenovirus-associated viral vectors (AAV). We report that GDL toxicity requires a single N-proximal cationic residue (K27 or R27 ) within PrPC . Alanine substitution of K27 also prevented the toxicity of PrPC mutants that induce Shmerling syndrome, a neurodegenerative disease that is suppressed by co-expression of wild-type PrPC . K27 may represent an actionable target for compounds aimed at preventing prion-related neurodegeneration.

The catalytic mechanism of the RNA methyltransferase METTL3.

The complex of methyltransferase-like proteins 3 and 14 (METTL3-14) is the major enzyme that deposits N6-methyladenosine (m6A) modifications on mRNA in humans. METTL3-14 plays key roles in various biological processes through its methyltransferase (MTase) activity. However, little is known about its substrate recognition and methyl transfer mechanism from its cofactor and methyl donor S-adenosylmethionine (SAM). Here, we study the MTase mechanism of METTL3-14 by a combined experimental and multiscale simulation approach using bisubstrate analogues (BAs), conjugates of a SAM-like moiety connected to the N6-atom of adenosine. Molecular dynamics simulations based on crystal structures of METTL3-14 with BAs suggest that the Y406 side chain of METTL3 is involved in the recruitment of adenosine and release of m6A. A crystal structure with a bisubstrate analogue representing the transition state of methyl transfer shows a direct involvement of the METTL3 side chains E481 and K513 in adenosine binding which is supported by mutational analysis. Quantum mechanics/molecular mechanics (QM/MM) free energy calculations indicate that methyl transfer occurs without prior deprotonation of adenosine-N6. Furthermore, the QM/MM calculations provide further support for the role of electrostatic contributions of E481 and K513 to catalysis. The multidisciplinary approach used here sheds light on the (co)substrate binding mechanism, catalytic step, and (co)product release catalysed by METTL3, and suggests that the latter step is rate-limiting. The atomistic information on the substrate binding and methyl transfer reaction of METTL3 can be useful for understanding the mechanisms of other RNA MTases and for the design of transition state analogues as their inhibitors.

The free energy landscape of small molecule unbinding.

The spontaneous dissociation of six small ligands from the active site of FKBP (the FK506 binding protein) is investigated by explicit water molecular dynamics simulations and network analysis. The ligands have between four (dimethylsulphoxide) and eleven (5-diethylamino-2-pentanone) non-hydrogen atoms, and an affinity for FKBP ranging from 20 to 0.2 mM. The conformations of the FKBP/ligand complex saved along multiple trajectories (50 runs at 310 K for each ligand) are grouped according to a set of intermolecular distances into nodes of a network, and the direct transitions between them are the links. The network analysis reveals that the bound state consists of several subbasins, i.e., binding modes characterized by distinct intermolecular hydrogen bonds and hydrophobic contacts. The dissociation kinetics show a simple (i.e., single-exponential) time dependence because the unbinding barrier is much higher than the barriers between subbasins in the bound state. The unbinding transition state is made up of heterogeneous positions and orientations of the ligand in the FKBP active site, which correspond to multiple pathways of dissociation. For the six small ligands of FKBP, the weaker the binding affinity the closer to the bound state (along the intermolecular distance) are the transition state structures, which is a new manifestation of Hammond behavior. Experimental approaches to the study of fragment binding to proteins have limitations in temporal and spatial resolution. Our network analysis of the unbinding simulations of small inhibitors from an enzyme paints a clear picture of the free energy landscape (both thermodynamics and kinetics) of ligand unbinding.

PARP1 ADP-ribosylates lysine residues of the core histone tails.

The chromatin-associated enzyme PARP1 has previously been suggested to ADP-ribosylate histones, but the specific ADP-ribose acceptor sites have remained enigmatic. Here, we show that PARP1 covalently ADP-ribosylates the amino-terminal histone tails of all core histones. Using biochemical tools and novel electron transfer dissociation mass spectrometric protocols, we identify for the first time K13 of H2A, K30 of H2B, K27 and K37 of H3, as well as K16 of H4 as ADP-ribose acceptor sites. Multiple explicit water molecular dynamics simulations of the H4 tail peptide into the catalytic cleft of PARP1 indicate that two stable intermolecular salt bridges hold the peptide in an orientation that allows K16 ADP-ribosylation. Consistent with a functional cross-talk between ADP-ribosylation and other histone tail modifications, acetylation of H4K16 inhibits ADP-ribosylation by PARP1. Taken together, our computational and experimental results provide strong evidence that PARP1 modifies important regulatory lysines of the core histone tails.

Kinase selectivity potential for inhibitors targeting the ATP binding site: a network analysis.

UNLABELLED: MOTIVATION AND METHOD: Small-molecule inhibitors targeting the adenosine triphosphate (ATP) binding pocket of the catalytic domain of protein kinases have potential to become drugs devoid of (major) side effects, particularly if they bind selectively. Here, the sequences of the 518 human kinases are first mapped onto the structural alignment of 116 kinases of known three-dimensional structure. The multiple structure alignment is then used to encode the known strategies for developing selective inhibitors into a fingerprint. Finally, a network analysis is used to partition the kinases into clusters according to similarity of their fingerprints, i.e. physico-chemical characteristics of the residues responsible for selective binding. RESULTS: For each kinase the network analysis reveals the likelihood to find selective inhibitors targeting the ATP binding site. Systematic guidelines are proposed to develop selective inhibitors. Importantly, the network analysis suggests that the tyrosine kinase EphB4 has high selectivity potential, which is consistent with the selectivity profile of two novel EphB4 inhibitors. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Atomistic and Thermodynamic Analysis of N6-Methyladenosine (m6A) Recognition by the Reader Domain of YTHDC1.

N6-Methyladenosine (m6A) is the most frequent modification in eukaryotic messenger RNA (mRNA) and its cellular processing and functions are regulated by the reader proteins YTHDCs and YTHDFs. However, the mechanism of m6A recognition by the reader proteins is still elusive. Here, we investigate this recognition process by combining atomistic simulations, site-directed mutagenesis, and biophysical experiments using YTHDC1 as a model. We find that the N6 methyl group of m6A contributes to the binding through its specific interactions with an aromatic cage (formed by Trp377 and Trp428) and also by favoring the association-prone conformation of m6A-containing RNA in solution. The m6A binding site dynamically equilibrates between multiple metastable conformations with four residues being involved in the regulation of m6A binding (Trp428, Met438, Ser378, and Thr379). Trp428 switches between two conformational states to build and dismantle the aromatic cage. Interestingly, mutating Met438 and Ser378 to alanine does not alter m6A binding to the protein but significantly redistributes the binding enthalpy and entropy terms, i.e., enthalpy-entropy compensation. Such compensation is reasoned by different entropy-enthalpy transduction associated with both conformational changes of the wild-type and mutant proteins and the redistribution of water molecules. In contrast, the point mutant Thr379Val significantly changes the thermal stability and binding capability of YTHDC1 to its natural ligand. Additionally, thermodynamic analysis and free energy calculations shed light on the role of a structural water molecule that synergistically binds to YTHDC1 with m6A and acts as the hub of a hydrogen-bond network. Taken together, the experimental data and simulation results may accelerate the discovery of chemical probes, m6A-editing tools, and drug candidates against reader proteins.

1,4,9-Triazaspiro[5.5]undecan-2-one Derivatives as Potent and Selective METTL3 Inhibitors.

N6-methyladenosine (m6A) is the most frequent of the 160 RNA modifications reported so far. Accumulating evidence suggests that the METTL3/METTL14 protein complex, part of the m6A regulation machinery, is a key player in a variety of diseases including several types of cancer, type 2 diabetes, and viral infections. Here we report on a protein crystallography-based medicinal chemistry optimization of a METTL3 hit compound that has resulted in a 1400-fold potency improvement (IC50 of 5 nM for the lead compound 22 (UZH2) in a time-resolved Förster resonance energy transfer (TR-FRET) assay). The series has favorable ADME properties as physicochemical characteristics were taken into account during hit optimization. UZH2 shows target engagement in cells and is able to reduce the m6A/A level of polyadenylated RNA in MOLM-13 (acute myeloid leukemia) and PC-3 (prostate cancer) cell lines.

METTL3 Inhibitors for Epitranscriptomic Modulation of Cellular Processes.

The methylase METTL3 is the writer enzyme of the N6 -methyladenosine (m6 A) modification of RNA. Using a structure-based drug discovery approach, we identified a METTL3 inhibitor with potency in a biochemical assay of 280 nM, while its enantiomer is 100 times less active. We observed a dose-dependent reduction in the m6 A methylation level of mRNA in several cell lines treated with the inhibitor already after 16 h of treatment, which lasted for at least 6 days. Importantly, the prolonged incubation (up to 6 days) with the METTL3 inhibitor did not alter levels of other RNA modifications (i. e., m1 A, m6 Am , m7 G), suggesting selectivity of the developed compound towards other RNA methyltransferases.

Library screening by fragment-based docking.

We review our computational tools for high-throughput screening by fragment-based docking of large collections of small molecules. Applications to six different enzymes, four proteases, and two protein kinases, are presented. Remarkably, several low-micromolar inhibitors were discovered in each of the high-throughput docking campaigns. Probable reasons for the lack of submicromolar inhibitors are the tiny fraction of chemical space covered by the libraries of available compounds, as well as the approximations in the methods employed for scoring, and the use of a rigid conformation of the target protein.

Small-Molecule Inhibitors of METTL3, the Major Human Epitranscriptomic Writer.

The RNA methylase METTL3 catalyzes the transfer of a methyl group from the cofactor S-adenosyl-L-methionine (SAM) to the N6 atom of adenine. We have screened a library of 4000 analogues and derivatives of the adenosine moiety of SAM by high-throughput docking into METTL3. Two series of adenine derivatives were identified in silico, and the binding mode of six of the predicted inhibitors was validated by protein crystallography. Two compounds, one for each series, show good ligand efficiency. We propose a route for their further development into potent and selective inhibitors of METTL3.

Selectively Disrupting m6A-Dependent Protein-RNA Interactions with Fragments.

We report a crystallographic analysis of small-molecule ligands of the human YTHDC1 domain that recognizes N6-methylated adenine (m6A) in RNA. The 30 binders are fragments (molecular weight < 300 g mol-1) that represent 10 different chemotypes identified by virtual screening. Despite the structural disorder of the binding site loop (residues 429-439), most of the 30 fragments emulate the two main interactions of the -NHCH3 group of m6A. These interactions are the hydrogen bond to the backbone carbonyl of Ser378 and the van der Waals contacts with the tryptophan cage. Different chemical groups are involved in the conserved binding motifs. Some of the fragments show favorable ligand efficiency for YTHDC1 and selectivity against other m6A reader domains. The structural information is useful for the design of modulators of m6A recognition by YTHDC1.

Hitting a Moving Target: Simulation and Crystallography Study of ATAD2 Bromodomain Blockers.

Small molecule ligand binding to the ATAD2 bromodomain is investigated here through the synergistic combination of molecular dynamics and protein crystallography. A previously unexplored conformation of the binding pocket upon rearrangement of the gatekeeper residue Ile1074 has been found. Further, our investigations reveal how minor structural differences in the ligands result in binding with different plasticity of the ZA loop for this difficult-to-drug bromodomain.

A fluorescence quenching assay to discriminate between specific and nonspecific inhibitors of dengue virus protease.

In drug discovery, the occurrence of false positives is a major hurdle in the search for lead compounds that can be developed into drugs. A small-molecular-weight compound that inhibits dengue virus protease at low micromolar levels was identified in a screening campaign. Binding to the enzyme was confirmed by isothermal titration calorimetry (ITC) and nuclear magnetic resonance (NMR). However, a structure-activity relationship study that ensued did not yield more potent leads. To further characterize the parental compound and its analogues, we developed a high-speed, low-cost, quantitative fluorescence quenching assay. We observed that specific analogues quenched dengue protease fluorescence and showed variation in IC(50) values. In contrast, nonspecifically binding compounds did not quench its fluorescence and showed similar IC(50) values with steep dose-response curves. We validated the assay using single Trp-to-Ala protease mutants and the competitive protease inhibitor aprotinin. Specific compounds detected in the binding assay were further analyzed by competitive ITC, NMR, and surface plasmon resonance, and the assay's utility in comparison with these biophysical methods is discussed. The sensitivity of this assay makes it highly useful for hit finding and validation in drug discovery. Furthermore, the technique can be readily adapted for studying other protein-ligand interactions.

Flexible Binding of m6A Reader Protein YTHDC1 to Its Preferred RNA Motif.

N6-Methyladenosine (m6A) is the most prevalent chemical modification in human mRNAs. Its recognition by reader proteins enables many cellular functions, including splicing and translation of mRNAs. However, the binding mechanisms of m6A-containing RNAs to their readers are still elusive due to the unclear roles of m6A-flanking ribonucleotides. Here, we use a model system, YTHDC1 with its RNA motif 5'-G-2G-1(m6A)C+1U+2-3', to investigate the binding mechanisms by atomistic simulations, X-ray crystallography, and isothermal titration calorimetry. The experimental data and simulation results show that m6A is captured by an aromatic cage of YTHDC1 and the 3' terminus nucleotides are stabilized by cation-π-π interactions, while the 5' terminus remains flexible. Notably, simulations of unbound RNA motifs reveal that the methyl group of m6A and the 5' terminus shift the conformational preferences of the oligoribonucleotide to the bound-like conformation, thereby facilitating the association process. The binding mechanisms may help in the discovery of chemical probes against m6A reader proteins.

Structure-based tailoring of compound libraries for high-throughput screening: discovery of novel EphB4 kinase inhibitors.

High-throughput docking is a computational tool frequently used to discover small-molecule inhibitors of enzymes or receptors of known three-dimensional structure. Because of the large number of molecules in chemical libraries, automatic procedures to prune multimillion compound collections are useful for high-throughput docking and necessary for in vitro screening. Here, we propose an anchor-based library tailoring approach (termed ALTA) to focus a chemical library by docking and prioritizing molecular fragments according to their binding energy which includes continuum electrostatics solvation. In principle, ALTA does not require prior knowledge of known inhibitors, but receptor-based pharmacophore information (hydrogen bonds with the hinge region) is additionally used here to identify molecules with optimal anchor fragments for the ATP-binding site of the EphB4 receptor tyrosine kinase. The 21,418 molecules of the focused library (from an initial collection of about 730,000) are docked into EphB4 and ranked by force-field-based energy including electrostatic solvation. Among the 43 compounds tested in vitro, eight molecules originating from two different anchors show low-micromolar activity in a fluorescence-based enzymatic assay. Four of them are active in a cell-based assay and are potential anti-angiogenic compounds.

Discovery of kinase inhibitors by high-throughput docking and scoring based on a transferable linear interaction energy model.

The linear interaction energy method with continuum electrostatics (LIECE) is evaluated in depth on five kinases. The two multiplicative coefficients for the van der Waals energy and electrostatic free energy are shown to be transferable among different kinases. Moreover, good enrichment factors are obtained for a library of 40375 diverse compounds seeded with 73 known inhibitors of CDK2. Therefore, a general two-parameter LIECE model for kinases is derived by combining large data sets of inhibitors of CDK2, Lck, and p38. This two-parameter model is cross-validated on two kinases not used for fitting; it shows an average error of about 1.5 kcal/mol for the prediction of absolute binding affinity of 37 and 128 known inhibitors of EphB4 and EGFR, respectively. High-throughput docking and ranking by two-parameter LIECE models are shown to be able to identify novel low-micromolar EphB4 and CDK2 inhibitors of low-molecular weight (< or =355 g/mol).

Discovery of cell-permeable non-peptide inhibitors of beta-secretase by high-throughput docking and continuum electrostatics calculations.

A fragment-based docking procedure followed by substructure search were used to identify active-site beta-secretase inhibitors from a composite set of about 300 000 and a library of nearly 180 000 small molecules, respectively. EC(50) values less than 10 microM were measured in at least one of two different mammalian cell-based assays for 12 of the 72 purchased compounds. In particular, the phenylureathiadiazole 2 and the diphenylurea derivative 3 are promising lead compounds for beta-secretase inhibition.