Electrostatic interactions determine entrance/release order of substrates in the catalytic cycle of adenylate kinase.

Adenylate kinase is a monomeric phosphotransferase with important biological function in regulating concentration of adenosine triphosphate (ATP) in cells, by transferring the terminal phosphate group from ATP to adenosine monophosphate (AMP) and forming two adenosine diphosphate (ADP) molecules. During this reaction, the kinase may undergo a large conformational transition, forming different states with its substrates. Although many structures of the protein are available, atomic details of the whole process remain unclear. In this article, we use both conventional molecular dynamics (MD) simulation and an enhanced sampling technique called parallel cascade selection MD simulation to explore different conformational states of the Escherichia coli adenylate kinase. Based on the simulation results, we propose a possible entrance/release order of substrates during the catalytic cycle. The substrate-free protein prefers an open conformation, but changes to a closed state once ATP·Mg enters into its binding pocket first and then AMP does. After the reaction of ATP transferring the terminal phosphate group to AMP, ADP·Mg and ADP are released sequentially, and finally the whole catalyze cycle is completed. Detailed contact and distance analysis reveals that the entrance/release order of substrates may be largely controlled by electrostatic interactions between the protein and the substrates.

NMR Fragment-Based Screening against Tandem RNA Recognition Motifs of TDP-43.

The TDP-43 is originally a nuclear protein but translocates to the cytoplasm in the pathological condition. TDP-43, as an RNA-binding protein, consists of two RNA Recognition Motifs (RRM1 and RRM2). RRMs are known to involve both protein-nucleotide and protein-protein interactions and mediate the formation of stress granules. Thus, they assist the entire TDP-43 protein with participating in neurodegenerative and cancer diseases. Consequently, they are potential therapeutic targets. Protein-observed and ligand-observed nuclear magnetic resonance (NMR) spectroscopy were used to uncover the small molecule inhibitors against the tandem RRM of TDP-43. We identified three hits weakly binding the tandem RRMs using the ligand-observed NMR fragment-based screening. The binding topology of these hits is then depicted by chemical shift perturbations (CSP) of the 15N-labeled tandem RRM and RRM2, respectively, and modeled by the CSP-guided High Ambiguity Driven biomolecular DOCKing (HADDOCK). These hits mainly bind to the RRM2 domain, which suggests the druggability of the RRM2 domain of TDP-43. These hits also facilitate further studies regarding the hit-to-lead evolution against the TDP-43 RRM domain.

NMR characterization of weak interactions between RhoGDI2 and fragment screening hits.

BACKGROUND: The delineation of intrinsically weak interactions between novel targets and fragment screening hits has long limited the pace of hit-to-lead evolution. Rho guanine-nucleotide dissociation inhibitor 2 (RhoGDI2) is a novel target that lacks any chemical probes for the treatment of tumor metastasis. METHODS: Protein-observed and ligand-observed NMR spectroscopy was used to characterize the weak interactions between RhoGDI2 and fragment screening hits. RESULTS: We identified three hits of RhoGDI2 using streamlined NMR fragment-based screening. The binding site residues were assigned using non-uniformly sampled Cα- and Hα-based three dimensional NMR spectra. The molecular docking to the proposed geranylgeranyl binding pocket of RhoGDI2 was guided by NMR restraints of chemical shift perturbations and ligand-observed transferred paramagnetic relaxation enhancement. We further validated the weak RhoGDI2-hit interactions using mutagenesis and structure-affinity analysis. CONCLUSIONS: Weak interactions between RhoGDI2 and fragment screening hits were delineated using an integrated NMR approach. GENERAL INTERESTS: Binders to RhoGDI2 as a potential anti-cancer target have been first reported, and their weak interactions were depicted using NMR spectroscopy. Our work highlights the powerfulness and the versatility of the integrative NMR techniques to provide valuable structural insight into the intrinsically weak interactions between RhoGDI2 and the fragment screening hits, which could hardly be conceived using other biochemical techniques.

Fluorine Pseudocontact Shifts Used for Characterizing the Protein-Ligand Interaction Mode in the Limit of NMR Intermediate Exchange.

The characterization of protein-ligand interaction modes becomes recalcitrant in the NMR intermediate exchange regime as the interface resonances are broadened beyond detection. Here, we determined the 19 F low-populated bound-state pseudocontact shifts (PCSs) of mono- and di-fluorinated inhibitors of the BRM bromodomain using a highly skewed protein/ligand ratio. The bound-state 19 F PCSs were retrieved from 19 F chemical exchange saturation transfer (CEST) in the presence of the lanthanide-labeled protein, which was termed the 19 F PCS-CEST approach. These PCSs enriched in spatial information enabled the identification of best-fitting poses, which agree well with the crystal structure of a more soluble analog in complex with the BRM bromodomain. This approach fills the gap of the NMR structural characterization of lead-like inhibitors with moderate affinities to target proteins, which are essential for structure-guided hit-to-lead evolution.

Proteins of well-defined structures can be designed without backbone readjustment by a statistical model.

We report that using mainly a statistical energy model, protein sequence design for designable backbones can be carried out with high confidence without considering backbone relaxation. A recently-developed statistical energy function for backbone-based protein sequence design has been rationally revised to improve its accuracy. As a demonstrative example, this revised model is applied to design a de novo protein for a target backbone for which the previous model had relied on after-design directed evolution to produce a well-folded protein. The actual backbone structure of the newly designed protein agrees excellently with the corresponding target. Besides presenting a new protein design protocol with experimentally verifications on different backbone types, our study implies that with an energy model of an appropriate resolution, proteins of well-defined structures instead of molten globules can be designed without the explicit consideration of backbone variations due to side chain changes, even if the side chain changes correspond to complete sequence redesigns.

Crystal structure of tRNA m1G9 methyltransferase Trm10: insight into the catalytic mechanism and recognition of tRNA substrate.

Transfer RNA (tRNA) methylation is necessary for the proper biological function of tRNA. The N(1) methylation of guanine at Position 9 (m(1)G9) of tRNA, which is widely identified in eukaryotes and archaea, was found to be catalyzed by the Trm10 family of methyltransferases (MTases). Here, we report the first crystal structures of the tRNA MTase spTrm10 from Schizosaccharomyces pombe in the presence and absence of its methyl donor product S-adenosyl-homocysteine (SAH) and its ortholog scTrm10 from Saccharomyces cerevisiae in complex with SAH. Our crystal structures indicated that the MTase domain (the catalytic domain) of the Trm10 family displays a typical SpoU-TrmD (SPOUT) fold. Furthermore, small angle X-ray scattering analysis reveals that Trm10 behaves as a monomer in solution, whereas other members of the SPOUT superfamily all function as homodimers. We also performed tRNA MTase assays and isothermal titration calorimetry experiments to investigate the catalytic mechanism of Trm10 in vitro. In combination with mutational analysis and electrophoretic mobility shift assays, our results provide insights into the substrate tRNA recognition mechanism of Trm10 family MTases.

Process of Fragment-Based Lead Discovery-A Perspective from NMR.

Fragment-based lead discovery (FBLD) has proven fruitful during the past two decades for a variety of targets, even challenging protein-protein interaction (PPI) systems. Nuclear magnetic resonance (NMR) spectroscopy plays a vital role, from initial fragment-based screening to lead generation, because of its power to probe the intrinsically weak interactions between targets and low-molecular-weight fragments. Here, we review the NMR FBLD process from initial library construction to lead generation. We describe technical aspects regarding fragment library design, ligand- and protein-observed screening, and protein-ligand structure model generation. For weak binders, the initial hit-to-lead evolution can be guided by structural information retrieved from NMR spectroscopy, including chemical shift perturbation, transferred pseudocontact shifts, and paramagnetic relaxation enhancement. This perspective examines structure-guided optimization from weak fragment screening hits to potent leads for challenging PPI targets.

Repurposing Low-Molecular-Weight Drugs against the Main Protease of Severe Acute Respiratory Syndrome Coronavirus 2.

The coronavirus disease pandemic caused by infection with the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has affected the global healthcare system. As low-molecular-weight drugs have high potential to completely match interactions with essential SARS-CoV-2 targets, we propose a strategy to identify such drugs using the fragment-based approach. Herein, using ligand- and protein-observed fragment screening approaches, we identified niacin and hit 1 binding to the catalytic pocket of the main protease (Mpro) of SARS-CoV-2, thereby modestly inhibiting the enzymatic activity of Mpro. We further searched for low-molecular-weight drugs containing niacin or hit 1 pharmacophores with enhanced inhibiting activity, e.g., carmofur, bendamustine, triclabendazole, emedastine, and omeprazole, in which omeprazole is the only one binding to the C-terminal domain of SARS-CoV-2 Mpro. Our study demonstrates that the fragment-based approach is a feasible strategy for identifying low-molecular-weight drugs against the SARS-CoV-2 and other potential targets lacking specific drugs.

Conformational Selection in Ligand Recognition by the First Tudor Domain of PHF20L1.

The first Tudor domain (Tudor1) of PHF20L1 recognizes (non)histone methylation to play versatile roles. However, the underlying ligand-recognition mechanism remains unknown as a closed state revealed in the free-form structure. NMR relaxation dispersion and molecular dynamics simulations suggest a pre-existing low-population conformation with a remarkable rearrangement of aromatic cage residues of PHF20L1 Tudor1. Such an open-form conformation is utilized to recognize lysine 142 methylated DNMT1, a cosolvent, and an NMR fragment screening hit, as revealed by the complex crystal structures. Intriguingly, the ligand binding capacity was enhanced by mutation that tunes up the open-state population only. The recognition of DNMT1 by PHF20L1 was further validated in cancer cells. This conformational selection mechanism will enable the discovery of small molecule inhibitors against the seemingly "undruggable" PHF20L1 Tudor1.

Backmapping from Multiresolution Coarse-Grained Models to Atomic Structures of Large Biomolecules by Restrained Molecular Dynamics Simulations Using Bayesian Inference.

Coarse-grained (CG) simulations have allowed access to larger length scales and longer time scales in the study of the dynamic processes of large biomolecules than all-atom (AA) molecular dynamics (MD) simulations. Backmapping from CG models to AA structures has long been studied because it enables us to gain detailed structure insights from CG simulations. Many methods first construct an AA structure from the CG model by fragments, random placement, or geometrical rules and subsequently optimize the solution via energy minimization, simulated annealing or position-restrained simulations. However, such methods may only work well on residue-level CG models and cannot consider the deviations of CG models. In this work, we describe, to the best of our knowledge, a new backmapping method based on Bayesian inference and restrained MD simulations. Restraints with log harmonic energy terms are defined according to the target CG model using the Bayesian inference in which the CG deviations can be estimated. From an initial AA structure obtained from either high-resolution experiments or homology modeling, a MD simulation with the aforementioned restraints is performed to obtain a final AA structure that is a backmapping of the target CG model. The method was validated using multiresolution CG models of the soluble extracellular region of the human epidermal growth factor receptor and was further applied to construct AA structures from CG simulations of the nucleosome core particle. The results demonstrate that our method can generate accurate AA structures of different types of biomolecules from multiple CG models with either residue-level resolution or much lower resolution than one-site-per-residue.

Structural basis for the complete resistance of the human prion protein mutant G127V to prion disease.

Prion diseases are caused by the propagation of misfolded cellular prion proteins (PrPs). A completely prion disease-resistant genotype, V127M129, has been identified in Papua New Guinea and verified in transgenic mice. To disclose the structural basis of the disease-resistant effect of the G127V mutant, we determined and compared the structural and dynamic features of the G127V-mutated human PrP (residues 91-231) and the wild-type PrP in solution. HuPrP(G127V) contains α1, α2 and α3 helices and a stretch-strand (SS) pattern comprising residues Tyr128-Gly131 (SS1) and Val161-Arg164 (SS2), with extending atomic distances between the SS1 and SS2 strands, and a structural rearrangement of the Tyr128 side chain due to steric hindrance of the larger hydrophobic side chain of Val127. The extended α1 helix gets closer to the α2 and α3 helices. NMR dynamics analysis revealed that Tyr128, Gly131 and Tyr163 underwent significant conformational exchanges. Molecular dynamics simulations suggest that HuPrP(G127V) prevents the formation of stable ß-sheets and dimers. Unique structural and dynamic features potentially inhibit the conformational conversion of the G127V mutant. This work is beneficial for understanding the molecular mechanisms underlying the complete resistance of the G127V mutant to prion disease and for developing new therapeutics for prion disease.

Spindlin-1 recognizes methylations of K20 and R23 of histone H4 tail.

Using methods combining cross-linking, pull-down assays, and stable isotope labeling by amino acids in cell culture with mass spectrometry, we identified that the Tudor domain-containing protein Spindlin-1 recognizes trimethylation of histone H4 lysine 20 (H4K20me3). The binding affinity of Spindlin-1 to H4K20me3 is weaker than that to H3K4me3, indicating H4K20me3 as a secondary substrate for Spindlin-1. Structural studies of Spindlin-1 in complex with the H4K20me3 peptide indicate that Spindlin-1 attains a distinct binding mode for H4K20me3 recognition. Further biochemical analysis identified that Spindlin-1 also binds methylated R23 of H4, providing new clues for the function of Spindlin-1.

Structural plasticity of the TDRD3 Tudor domain probed by a fragment screening hit.

As a reader of di-methylated arginine on various proteins, such as histone, RNA polymerase II, PIWI and Fragile X mental retardation protein, the Tudor domain of Tudor domain-containing protein 3 (TDRD3) mediates transcriptional activation in nucleus and formation of stress granules in the cytoplasm. Despite the TDRD3 implication in cancer cell proliferation and invasion, warheads to block the di-methylated arginine recognition pocket of the TDRD3 Tudor domain have not yet been uncovered. Here we identified 14 small molecule hits against the TDRD3 Tudor domain through NMR fragment-based screening. These hits were further cross-validated by using competitive fluorescence polarization and isothermal titration calorimetry experiments. The crystal structure of the TDRD3 Tudor domain in complex with hit 1 reveals a distinct binding mode from the nature substrate. Hit 1 protrudes into the aromatic cage of the TDRD3 Tudor domain, where the aromatic residues are tilted to accommodate a sandwich-like π-π interaction. The side chain of the conserved residue N596 swings away 3.1 Å to form a direct hydrogen bond with hit 1. Moreover, this compound shows a decreased affinity against the single Tudor domain of survival motor neuron protein, but no detectable binding to neither the tandem Tudor domain of TP53-binding protein 1 nor the extended Tudor domain of staphylococcal nuclease domain-containing protein 1. Our work depicts the structural plasticity of the TDRD3 Tudor domain and paves the way for the subsequent structure-guided discovery of selective inhibitors targeting Tudor domains. DATABASE: Structural data are available in the PDB under the accession number 5YJ8.

Ligand Proton Pseudocontact Shifts Determined from Paramagnetic Relaxation Dispersion in the Limit of NMR Intermediate Exchange.

Delineation of protein-ligand interaction modes is key for rational drug discovery. The availability of complex crystal structures is often limited by the aqueous solubility of the compounds, while lead-like compounds with micromolar affinities normally fall into the NMR intermediate exchange regime, in which severe line broadening to beyond the detection of interfacial resonances limits NMR applications. Here, we developed a new method to retrieve low-populated bound-state 1H pseudocontact shifts (PCSs) using paramagnetic relaxation dispersion (RD). We evaluated using a 1H PCS-RD approach in a BRM bromodomain lead-like inhibitor to filter molecular docking poses using multiple intermolecular structural restraints. Considering the universal presence of proton atoms in druglike compounds, our work will have wide application in structure-guided drug discovery even under an extreme condition of NMR intermediate exchange and low aqueous solubility of ligands.

Dynamic Nature of CTCF Tandem 11 Zinc Fingers in Multivalent Recognition of DNA As Revealed by NMR Spectroscopy.

The 11 zinc fingers (ZFs) of the transcription factor CTCF play a versatile role in the regulation of gene expression. CTCF binds to numerous genomic sites to form chromatin loops and topologically associated domains and thus mediates the 3D architecture of chromatin. Although CTCF inter-ZF plasticity is essential for the recognition of multiple genomic sites, the dynamic nature of its 11 ZFs remains unknown. We assigned the chemical shifts of the CTCF ZFs 1-11 and solved the solution structures of each ZF. NMR backbone dynamics, residual dipolar couplings, and small-angle X-ray scattering experiments suggest a high inter-ZF plasticity of the free-form ZFs 1-11. As exemplified by two different protocadherin DNA sequences, the titration of DNAs to 15N-labeled CTCF ZFs 1-11 enabled systematic mapping of binding of CTCF ZFs to various chromatin sites. Our work paves the way for illustrating the molecular basis of the versatile DNA recognized by CTCF and has interesting implications for its conformational transition during DNA binding.

The polar warhead of a TRIM24 bromodomain inhibitor rearranges a water-mediated interaction network.

Tripartite motif-containing protein 24 (TRIM24) is closely correlated with multiple cancers, and a recent study demonstrated that the bromodomain of TRIM24 is essential for the proliferation of lethal castration-resistant prostate cancer. Here, we identify three new inhibitors of the TRIM24 bromodomain using NMR fragment-based screening. The crystal structures of two new inhibitors in complex with the TRIM24 bromodomain reveal that the water-bridged interaction network is conserved in the same fashion as those for known benzoimidazolone inhibitors. Interestingly, the polar substitution on the warhead of one new inhibitor pulls the whole ligand approximately 2 Å into the inner side pocket of the TRIM24 bromodomain, and thus exhibits a binding mode significantly different from other known bromodomain ligands. This mode provides a useful handle for further hit-to-lead evolution toward novel inhibitors of the TRIM24 bromodomain. DATABASE: Structural data are available in the PDB under the accession numbers 5H1T, 5H1U, and 5H1V.

Automated NMR fragment based screening identified a novel interface blocker to the LARG/RhoA complex.

The small GTPase cycles between the inactive GDP form and the activated GTP form, catalyzed by the upstream guanine exchange factors. The modulation of such process by small molecules has been proven to be a fruitful route for therapeutic intervention to prevent the over-activation of the small GTPase. The fragment based approach emerging in the past decade has demonstrated its paramount potential in the discovery of inhibitors targeting such novel and challenging protein-protein interactions. The details regarding the procedure of NMR fragment screening from scratch have been rarely disclosed comprehensively, thus restricts its wider applications. To achieve a consistent screening applicable to a number of targets, we developed a highly automated protocol to cover every aspect of NMR fragment screening as possible, including the construction of small but diverse libray, determination of the aqueous solubility by NMR, grouping compounds with mutual dispersity to a cocktail, and the automated processing and visualization of the ligand based screening spectra. We exemplified our streamlined screening in RhoA alone and the complex of the small GTPase RhoA and its upstream guanine exchange factor LARG. Two hits were confirmed from the primary screening in cocktail and secondary screening over individual hits for LARG/RhoA complex, while one of them was also identified from the screening for RhoA alone. HSQC titration of the two hits over RhoA and LARG alone, respectively, identified one compound binding to RhoA.GDP at a 0.11 mM affinity, and perturbed the residues at the switch II region of RhoA. This hit blocked the formation of the LARG/RhoA complex, validated by the native gel electrophoresis, and the titration of RhoA to ¹5N labeled LARG in the absence and presence the compound, respectively. It therefore provides us a starting point toward a more potent inhibitor to RhoA activation catalyzed by LARG.