Binding site identification and structure determination of protein-ligand complexes by NMR a semiautomated approach

Equivariant Flexible Modeling of the Protein-Ligand Binding Pose with Geometric Deep Learning

Flexible modeling of the protein-ligand complex structure is a fundamental challenge for in silico drug development. Recent studies have improved commonly used docking tools by incorporating extra-deep learning-based steps. However, such strategies limit their accuracy and efficiency because they retain massive sampling pressure and lack consideration for flexible biomolecular changes. In this study, we propose FlexPose, a geometric graph network capable of direct flexible modeling of complex structures in Euclidean space without the following conventional sampling and scoring strategies. Our model adopts two key designs: scalar-vector dual feature representation and SE(3)-equivariant network, to manage dynamic structural changes, as well as two strategies: conformation-aware pretraining and weakly supervised learning, to boost model generalizability in unseen chemical space. Benefiting from these paradigms, our model dramatically outperforms all tested popular docking tools and recently advanced deep learning methods, especially in tasks involving protein conformation changes. We further investigate the impact of protein and ligand similarity on the model performance with two conformation-aware strategies. Moreover, FlexPose provides an affinity estimation and model confidence for postanalysis.

Fragment-based screening by protein-detected NMR spectroscopy

Fragment-based drug discovery (FBDD) identifies low molecular weight compounds that can be developed into ligands with high affinity and selectivity for therapeutic targets. Screening fragment libraries (<10,000 molecules) with biophysical techniques against macromolecules provides information about novel chemical spaces that bind the macromolecule and scaffolds that can be modified to increase potency. A fragment-screening pipeline requires a standardized protocol for target selection, library assembly and maintenance, library screening, and hit validation to ensure hit integrity. Herein, the fundamental aspects of a fragment screening pipeline-focusing on protein-detected NMR data collection and analysis-are discussed in detail for researchers to use as a resource in their FBDD projects. Selected screening targets must undergo rigorous stability and buffer testing by NMR spectroscopy to ensure the protein structure is stable for the entire screen. Biophysical instrumentation that rapidly measures protein thermostability is helpful in buffer screening. Molecules in fragment libraries are analyzed computationally and physically, stored at appropriate temperatures, and multiplexed in well plates for library conservation. The screening protocol is streamlined using liquid handling robotics for sample preparation and customized Python scripts for protein-detected NMR data analysis. Molecules identified from the screen are titrated to determine their binding site(s) and Kd values and confirmed with an orthogonal biophysical assay. This detailed FBDD screening pipeline developed by the Program in Chemical Biology at the Medical College of Wisconsin has successfully screened many unrelated target proteins to identified novel molecules that selectively bind to these target proteins.

Fragment-based drug discovery of small molecule ligands for the human chemokine CCL28

The mucosal chemokine CCL28 is a promising target for immunotherapy drug development due to its elevated expression level in epithelial cells and critical role in creating and maintaining an immunosuppressive tumor microenvironment. Using sulfotyrosine as a probe, NMR chemical shift mapping identified a potential receptor-binding hotspot on the human CCL28 surface. CCL28 was screened against 2,678 commercially available chemical fragments by 2D NMR, yielding thirteen verified hits. Computational docking predicted that two fragments could occupy adjoining subsites within the sulfotyrosine recognition cleft. Dual NMR titrations confirmed their ability to bind CCL28 simultaneously, thereby validating an initial fragment pair for linking and merging strategies to design high-potency CCL28 inhibitors.

Protocol to identify drug-binding sites in proteins using solution NMR spectroscopy

Dusquetide is a next-generation IDR (innate defense regulator) targeting the major autophagy receptor protein SQSTM1/p62 and modulating the innate immune response. Here, we describe a protocol for determining dusquetide-binding sites of p62 by solution NMR spectroscopy. Step-by-step technique details were provided, including sample preparation, NMR experiment setup, data processing, and binding site analysis. This protocol could be applied to characterize other small molecules targeting the ZZ domain of p62 (9 kDa) or other proteins containing ZZ domains. For complete details on the use and execution of this protocol, please refer to Zhang et al. (2022).

Binding Adaptation of GS-441524 Diversifies Macro Domains and Downregulates SARS-CoV-2 de-MARylation Capacity

Viral infection in cells triggers a cascade of molecular defense mechanisms to maintain host-cell homoeostasis. One of these mechanisms is ADP-ribosylation, a fundamental post-translational modification (PTM) characterized by the addition of ADP-ribose (ADPr) on substrates. Poly(ADP-ribose) polymerases (PARPs) are implicated in this process and they perform ADP-ribosylation on host and pathogen proteins. Some viral families contain structural motifs that can reverse this PTM. These motifs known as macro domains (MDs) are evolutionarily conserved protein domains found in all kingdoms of life. They are divided in different classes with the viral belonging to Macro-D-type class because of their properties to recognize and revert the ADP-ribosylation. Viral MDs are potential pharmaceutical targets, capable to counteract host immune response. Sequence and structural homology between viral and human MDs are an impediment for the development of new active compounds against their function. Remdesivir, is a drug administrated in viral infections inhibiting viral replication through RNA-dependent RNA polymerase (RdRp). Herein, GS-441524, the active metabolite of the remdesivir, is tested as a hydrolase inhibitor for several viral MDs and for its binding to human homologs found in PARPs. This study presents biochemical and biophysical studies, which indicate that GS-441524 selectively modifies SARS-CoV-2 MD de-MARylation activity, while it does not interact with hPARP14 MD2 and hPARP15 MD2. The structural investigation of MD•GS-441524 complexes, using solution NMR and X-ray crystallography, discloses the impact of certain amino acids in ADPr binding cavity suggesting that F360 and its adjacent residues tune the selective binding of the inhibitor to SARS-CoV-2 MD.

Implications of critical node-dependent unidirectional cross-talk of Plasmodium SUMO pathway proteins

The endoparasitic pathogen, Plasmodium falciparum (Pf), modulates protein-protein interactions to employ post-translational modifications like SUMOylation to establish successful infections. The interaction between E1 and E2 (Ubc9) enzymes governs species specificity in the Plasmodium SUMOylation pathway. Here, we demonstrate that a unidirectional cross-species interaction exists between Pf-SUMO and human E2, whereas Hs-SUMO1 failed to interact with Pf-E2. Biochemical and biophysical analyses revealed that surface-accessible aspartates of Pf-SUMO determine the efficacy and specificity of SUMO-Ubc9 interactions. Furthermore, we demonstrate that critical residues of the Pf-Ubc9 N terminus are responsible for diminished Hs-SUMO1 and Pf-Ubc9 interaction. Mutating these residues to corresponding Hs-Ubc9 residues restores electrostatic, π-π, and hydrophobic interactions and allows efficient cross-species interactions. We suggest that, in comparison with human counterparts, Plasmodium SUMO and Ubc9 proteins have acquired critical changes on their surfaces as nodes, which Plasmodium can use to exploit the host SUMOylation machinery.

Parallel detection and spatial mapping of large nuclear spin clusters

Nuclear magnetic resonance imaging (MRI) at the atomic scale offers exciting prospects for determining the structure and function of individual molecules and proteins. Quantum defects in diamond have recently emerged as a promising platform towards reaching this goal, and allowed for the detection and localization of single nuclear spins under ambient conditions. Here, we present an efficient strategy for extending imaging to large nuclear spin clusters, fulfilling an important requirement towards a single-molecule MRI technique. Our method combines the concepts of weak quantum measurements, phase encoding and simulated annealing to detect three-dimensional positions from many nuclei in parallel. Detection is spatially selective, allowing us to probe nuclei at a chosen target radius while avoiding interference from strongly-coupled proximal nuclei. We demonstrate our strategy by imaging clusters containing more than 20 carbon-13 nuclear spins within a radius of 2.4 nm from single, near-surface nitrogen-vacancy centers at room temperature. The radius extrapolates to 5-6 nm for ¹H. Beside taking an important step in nanoscale MRI, our experiment also provides an efficient tool for the characterization of large nuclear spin registers in the context of quantum simulators and quantum network nodes.

Structure and Gating Behavior of the Human Integral Membrane Protein VDAC1 in a Lipid Bilayer

The voltage-dependent anion channel (VDAC), the most abundant protein in the outer mitochondrial membrane, is responsible for the transport of all ions and metabolites into and out of mitochondria. Larger than any of the β-barrel structures determined to date by magic-angle spinning (MAS) NMR, but smaller than the size limit of cryo-electron microscopy (cryo-EM), VDAC1’s 31 kDa size has long been a bottleneck in determining its structure in a near-native lipid bilayer environment. Using a single two-dimensional (2D) crystalline sample of human VDAC1 in lipids, we applied proton-detected fast magic-angle spinning NMR spectroscopy to determine the arrangement of β strands. Combining these data with long-range restraints from a spin-labeled sample, chemical shift-based secondary structure prediction, and previous MAS NMR and atomic force microscopy (AFM) data, we determined the channel’s structure at a 2.2 Å root-mean-square deviation (RMSD). The structure, a 19-stranded β-barrel, with an N-terminal α-helix in the pore is in agreement with previous data in detergent, which was questioned due to the potential for the detergent to perturb the protein’s functional structure. Using a quintuple mutant implementing the channel’s closed state, we found that dynamics are a key element in the protein’s gating behavior, as channel closure leads to the destabilization of not only the C-terminal barrel residues but also the α2 helix. We showed that cholesterol, previously shown to reduce the frequency of channel closure, stabilizes the barrel relative to the N-terminal helix. Furthermore, we observed channel closure through steric blockage by a drug shown to selectively bind to the channel, the Bcl2-antisense oligonucleotide G3139.

Selective binding of small molecules to Vibrio cholerae DsbA offers a starting point for the design of novel antibacterials

DsbA enzymes catalyze oxidative folding of proteins that are secreted into the periplasm of Gram‐negative bacteria, and they are indispensable for the virulence of human pathogens such as Vibrio cholerae and Escherichia coli. Therefore, targeting DsbA represents an attractive approach to control bacterial virulence. X‐ray crystal structures reveal that DsbA enzymes share a similar fold, however, the hydrophobic groove adjacent to the active site, which is implicated in substrate binding, is shorter and flatter in the structure of V. cholerae DsbA (VcDsbA) compared to E. coli DsbA (EcDsbA). The flat and largely featureless nature of this hydrophobic groove is challenging for the development of small molecule inhibitors. Using fragment‐based screening approaches, we have identified a novel small molecule, based on the benzimidazole scaffold, that binds to the hydrophobic groove of oxidized VcDsbA with a K_D of 446±10 μM. The same benzimidazole compound has ∼8‐fold selectivity for VcDsbA over EcDsbA and binds to oxidized EcDsbA, with K_D>3.5 mM. We generated a model of the benzimidazole complex with VcDsbA using NMR data but were unable to determine the structure of the benzimidazole bound EcDsbA using either NMR or X‐ray crystallography. Therefore, a structural basis for the observed selectivity is unclear. To better understand ligand binding to these two enzymes we crystallized each of them in complex with a known ligand, the bile salt sodium taurocholate. The crystal structures show that taurocholate adopts different binding poses in complex with VcDsbA and EcDsbA, and reveal the protein‐ligand interactions that stabilize the different modes of binding. This work highlights the capacity of fragment‐based drug discovery to identify inhibitors of challenging protein targets. In addition, it provides a starting point for development of more potent and specific VcDsbA inhibitors that act through a novel anti‐virulence mechanism.

Identification and characterization of two drug-like fragments that bind to the same cryptic binding pocket of Burkholderia pseudomallei DsbA

Disulfide-bond-forming proteins (Dsbs) play a crucial role in the pathogenicity of many Gram-negative bacteria. Disulfide-bond-forming protein A (DsbA) catalyzes the formation of the disulfide bonds necessary for the activity and stability of multiple substrate proteins, including many virulence factors. Hence, DsbA is an attractive target for the development of new drugs to combat bacterial infections. Here, two fragments, bromophenoxy propanamide (1) and 4-methoxy-N-phenylbenzenesulfonamide (2), were identified that bind to DsbA from the pathogenic bacterium Burkholderia pseudomallei, the causative agent of melioidosis. The crystal structures of oxidized B. pseudomallei DsbA (termed BpsDsbA) co-crystallized with 1 or 2 show that both fragments bind to a hydrophobic pocket that is formed by a change in the side-chain orientation of Tyr110. This conformational change opens a `cryptic' pocket that is not evident in the apoprotein structure. This binding location was supported by 2D-NMR studies, which identified a chemical shift perturbation of the Tyr110 backbone amide resonance of more than 0.05 p.p.m. upon the addition of 2 mM fragment 1 and of more than 0.04 p.p.m. upon the addition of 1 mM fragment 2. Although binding was detected by both X-ray crystallography and NMR, the binding affinity (Kd) for both fragments was low (above 2 mM), suggesting weak interactions with BpsDsbA. This conclusion is also supported by the crystal structure models, which ascribe partial occupancy to the ligands in the cryptic binding pocket. Small fragments such as 1 and 2 are not expected to have a high energetic binding affinity due to their relatively small surface area and the few functional groups that are available for intermolecular interactions. However, their simplicity makes them ideal for functionalization and optimization. The identification of the binding sites of 1 and 2 to BpsDsbA could provide a starting point for the development of more potent novel antimicrobial compounds that target DsbA and bacterial virulence.

Studying the Structures of Relaxed and Fuzzy Interactions: The Diverse World of S100 Complexes

S100 proteins are small, dimeric, Ca²⁺-binding proteins of considerable interest due to their associations with cancer and rheumatic and neurodegenerative diseases. They control the functions of numerous proteins by forming protein–protein complexes with them. Several of these complexes were found to display “fuzzy” properties. Examining these highly flexible interactions, however, is a difficult task, especially from a structural biology point of view. Here, we summarize the available in vitro techniques that can be deployed to obtain structural information about these dynamic complexes. We also review the current state of knowledge about the structures of S100 complexes, focusing on their often-asymmetric nature.

Elaboration of a benzofuran scaffold and evaluation of binding affinity and inhibition of Escherichia coli DsbA: A fragment-based drug design approach to novel antivirulence compounds

Bacterial thiol-disulfide oxidoreductase DsbA is essential for bacterial virulence factor assembly and has been identified as a viable antivirulence target. Herein, we report a structure-based elaboration of a benzofuran hit that bound to the active site groove of Escherichia coli DsbA. Substituted phenyl groups were installed at the 5- and 6-position of the benzofuran using Suzuki-Miyaura coupling. HSQC NMR titration experiments showed dissociation constants of this series in the high µM to low mM range and X-ray crystallography produced three co-structures, showing binding in the hydrophobic groove, comparable with that of the previously reported benzofurans. The 6-(m-methoxy)phenyl analogue (2b), which showed a promising binding pose, was chosen for elaboration from the C-2 position. The 2,6-disubstituted analogues bound to the hydrophobic region of the binding groove and the C-2 groups extended into the more polar, previously un-probed, region of the binding groove. Biochemical analysis of the 2,6-disubsituted analogues showed they inhibited DsbA oxidation activity in vitro. The results indicate the potential to develop the elaborated benzofuran series into a novel class of antivirulence compounds.

Overcoming challenges in developing small molecule inhibitors for GPVI and CLEC-2

GPVI and CLEC-2 have emerged as promising targets for long-term prevention of both arterial thrombosis and thrombo-inflammation with a decreased bleeding risk relative to current drugs. However, while there are potent blocking antibodies of both receptors, their protein nature comes with decreased bioavailability, making formulation for oral medication challenging. Small molecules are able to overcome these limitations, but there are many challenges in developing antagonists of nanomolar potency, which is necessary when considering the structural features that underlie the interaction of CLEC-2 and GPVI with their protein ligands. In this review, we describe current small-molecule inhibitors for both receptors and strategies to overcome such limitations, including considerations when it comes to in silico drug design and the importance of complex compound library selection.

Chemo-informatics guided study of natural inhibitors targeting rho GTPase: a lead for treatment of glaucoma

Glaucoma, the most perilous disease leading to blindness is a result of optical neuropathy. Accumulation of aqueous humor in the posterior chamber due to a large difference in the rate of formation and its drainage in the anterior chamber causes an increase in intraocular pressure (IOP) leading to damage of nerve cells. A literature survey has revealed that inhibition of the Rho guanosine triphosphatases (rho GTPase) pathway by specific inhibitors leads to the relaxation of contractile cells involved in the aqueous outflow pathway. Relaxation of the strained contractile cells results in increased outflow thereby releasing IOP. In the present study molecular docking has been used to screen twenty seven bioactive (17 natural compounds and 10 conventional drugs) compounds that may play a significant role in relaxing contractile cells by inhibiting rho-GTPase protein. Docking results showed that among all-natural bioactive compounds Cyanidin and Delphinidine have a good binding affinity (- 8.4 kcal/mol) than the top screened conventional drug molecule Mitomycin, (- 6.3 kcal/mol) when docked with rho-GTPase protein. Cyanidin and Delphinidin belong to anthocyanidin, a glycoside form of anthocyanins from Vaccinium myrtillus L. and Punica granatum. The resembling potential of Cyanidin and Delphinidin concerning the drug Mitomycin was confirmed through simulation analysis. Molecular dynamics study (MDS) for 100 ns, showed that the rho GTPase-Delphinidine complex structure was energetically more stable than rho GTPase-Cyaniding complex in comparison to rho GTPase-Mitomycin complex. The comparative study of both the selected hits (Cyanidin and Delphinidin) was assessed by RMSD, RMSF, Rg, SASA, H-bond, PCA MM/PBSA analysis. The analysis revealed that Delphinidine is more potent to inhibit the rho GTPase as compare to Cyaniding and available conventional drugs in terms of stability and binding free energy. Based on the results, these molecules have good pharmacokinetic and pharmacodynamics properties and will prove to be a promising lead compound as a future drug for Glaucoma.

Functional and Structural Aspects of La Protein Overexpression in Lung Cancer

La is an abundant phosphoprotein that protects polymerase III transcripts from 3'-5' exonucleolytic degradation and facilitates their folding. Consisting of the evolutionary conserved La motif (LAM) and two consecutive RNA Recognition Motifs (RRMs), La was also found to bind additional RNA transcripts or RNA domains like internal ribosome entry site (IRES), through sequence-independent binding modes which are poorly understood. Although it has been reported overexpressed in certain cancer types and depletion of its expression sensitizes cancer cells to certain chemotherapeutic agents, its role in cancer remains essentially uncharacterized. Herein, we study the effects of La overexpression in A549 lung adenocarcinoma cells, which leads to increased cell proliferation and motility. Expression profiling of several transcription and translation factors indicated that La overexpression leads to downregulation of global translation through hypophosphorylation of 4E-BPs and upregulation of IRES-mediated translation. Moreover, analysis of La localization after nutrition deprivation of the transfected cells showed a normal distribution in the nucleus and nucleoli. Although the RNA binding capacity of La has been primarily linked to the synergy between the conserved LAM and RRM1 domains which act as a module, we show that recombinant stand-alone LAM can specifically bind a pre-tRNA ligand, based on binding experiments combined with NMR analysis. We propose that LAM RNA binding properties could support the expanding and diverse RNA ligand repertoire of La, thus promoting its modulatory role, both under normal and pathogenic conditions like cancer.

NMR fragment screening reveals a novel small molecule binding site near the catalytic surface of the disulfide-dithiol oxidoreductase enzyme DsbA from Burkholderia pseudomallei

The presence of suitable cavities or pockets on protein structures is a general criterion for a therapeutic target protein to be classified as 'druggable'. Many disease-related proteins that function solely through protein-protein interactions lack such pockets, making development of inhibitors by traditional small-molecule structure-based design methods much more challenging. The 22 kDa bacterial thiol oxidoreductase enzyme, DsbA, from the gram-negative bacterium Burkholderia pseudomallei (BpsDsbA) is an example of one such target. The crystal structure of oxidized BpsDsbA lacks well-defined surface pockets. BpsDsbA is required for the correct folding of numerous virulence factors in B. pseudomallei, and genetic deletion of dsbA significantly attenuates B. pseudomallei virulence in murine infection models. Therefore, BpsDsbA is potentially an attractive drug target. Herein we report the identification of a small molecule binding site adjacent to the catalytic site of oxidized BpsDsbA. ¹H^N CPMG relaxation dispersion NMR measurements suggest that the binding site is formed transiently through protein dynamics. Using fragment-based screening, we identified a small molecule that binds at this site with an estimated affinity of K_D ~ 500 µM. This fragment inhibits BpsDsbA enzymatic activity in vitro. The binding mode of this molecule has been characterized by NMR data-driven docking using HADDOCK. These data provide a starting point towards the design of more potent small molecule inhibitors of BpsDsbA.

NMR as a "Gold Standard" Method in Drug Design and Discovery

Studying disease models at the molecular level is vital for drug development in order to improve treatment and prevent a wide range of human pathologies. Microbial infections are still a major challenge because pathogens rapidly and continually evolve developing drug resistance. Cancer cells also change genetically, and current therapeutic techniques may be (or may become) ineffective in many cases. The pathology of many neurological diseases remains an enigma, and the exact etiology and underlying mechanisms are still largely unknown. Viral infections spread and develop much more quickly than does the corresponding research needed to prevent and combat these infections; the present and most relevant outbreak of SARS-CoV-2, which originated in Wuhan, China, illustrates the critical and immediate need to improve drug design and development techniques. Modern day drug discovery is a time-consuming, expensive process. Each new drug takes in excess of 10 years to develop and costs on average more than a billion US dollars. This demonstrates the need of a complete redesign or novel strategies. Nuclear Magnetic Resonance (NMR) has played a critical role in drug discovery ever since its introduction several decades ago. In just three decades, NMR has become a "gold standard" platform technology in medical and pharmacology studies. In this review, we present the major applications of NMR spectroscopy in medical drug discovery and development. The basic concepts, theories, and applications of the most commonly used NMR techniques are presented. We also summarize the advantages and limitations of the primary NMR methods in drug development.

Solution structure of the Hop TPR2A domain and investigation of target druggability by NMR, biochemical and in silico approaches

Heat shock protein 90 (Hsp90) is a molecular chaperone that plays an important role in tumour biology by promoting the stabilisation and activity of oncogenic 'client' proteins. Inhibition of Hsp90 by small-molecule drugs, acting via its ATP hydrolysis site, has shown promise as a molecularly targeted cancer therapy. Owing to the importance of Hop and other tetratricopeptide repeat (TPR)-containing cochaperones in regulating Hsp90 activity, the Hsp90-TPR domain interface is an alternative site for inhibitors, which could result in effects distinct from ATP site binders. The TPR binding site of Hsp90 cochaperones includes a shallow, positively charged groove that poses a significant challenge for druggability. Herein, we report the apo, solution-state structure of Hop TPR2A which enables this target for NMR-based screening approaches. We have designed prototype TPR ligands that mimic key native 'carboxylate clamp' interactions between Hsp90 and its TPR cochaperones and show that they block binding between Hop TPR2A and the Hsp90 C-terminal MEEVD peptide. We confirm direct TPR-binding of these ligands by mapping 1H-15N HSQC chemical shift perturbations to our new NMR structure. Our work provides a novel structure, a thorough assessment of druggability and robust screening approaches that may offer a potential route, albeit difficult, to address the chemically challenging nature of the Hop TPR2A target, with relevance to other TPR domain interactors.

Identifying Ortholog Selective Fragment Molecules for Bacterial Glutaredoxins by NMR and Affinity Enhancement by Modification with an Acrylamide Warhead

Illustrated here is the development of a new class of antibiotic lead molecules targeted at Pseudomonas aeruginosa glutaredoxin (PaGRX). This lead was produced to (a) circumvent efflux-mediated resistance mechanisms via covalent inhibition while (b) taking advantage of species selectivity to target a fundamental metabolic pathway. This work involved four components: a novel workflow for generating protein specific fragment hits via independent nuclear magnetic resonance (NMR) measurements, NMR-based modeling of the target protein structure, NMR guided docking of hits, and synthetic modification of the fragment hit with a vinyl cysteine trap moiety, i.e., acrylamide warhead, to generate the chimeric lead. Reactivity of the top warhead-fragment lead suggests that the ortholog selectivity observed for a fragment hit can translate into a substantial kinetic advantage in the mature warhead lead, which bodes well for future work to identify potent, species specific drug molecules targeted against proteins heretofore deemed undruggable.

The Fragment-Based Development of a Benzofuran Hit as a New Class of Escherichia coli DsbA Inhibitors

A fragment-based drug discovery approach was taken to target the thiol-disulfide oxidoreductase enzyme DsbA from Escherichia coli (EcDsbA). This enzyme is critical for the correct folding of virulence factors in many pathogenic Gram-negative bacteria, and small molecule inhibitors can potentially be developed as anti-virulence compounds. Biophysical screening of a library of fragments identified several classes of fragments with affinity to EcDsbA. One hit with high mM affinity, 2-(6-bromobenzofuran-3-yl)acetic acid (6), was chemically elaborated at several positions around the scaffold. X-ray crystal structures of the elaborated analogues showed binding in the hydrophobic binding groove adjacent to the catalytic disulfide bond of EcDsbA. Binding affinity was calculated based on NMR studies and compounds 25 and 28 were identified as the highest affinity binders with dissociation constants (K_D) of 326 ± 25 and 341 ± 57 µM respectively. This work suggests the potential to develop benzofuran fragments into a novel class of EcDsbA inhibitors.

Deciphering the Nucleotide and RNA Binding Selectivity of the Mayaro Virus Macro Domain

Mayaro virus (MAYV) is a member of Togaviridae family, which also includes Chikungunya virus as a notorious member. MAYV recently emerged in urban areas of the Americas, and this emergence emphasized the current paucity of knowledge about its replication cycle. The macro domain (MD) of MAYV belongs to the N-terminal region of its non-structural protein 3, part of the replication complex. Here, we report the first structural and dynamical characterization of a previously unexplored Alphavirus MD investigated through high-resolution NMR spectroscopy, along with data on its ligand selectivity and binding properties. The structural analysis of MAYV MD reveals a typical "macro" (ββαββαβαβα) fold for this polypeptide, while NMR-driven interaction studies provide in-depth insights into MAYV MD-ligand adducts. NMR data in concert with thermodynamics and biochemical studies provide convincing experimental evidence for preferential binding of adenosine diphosphate ribose (ADP-r) and adenine-rich RNAs to MAYV MD, thus shedding light on the structure-function relationship of a previously unexplored viral MD. The emerging differences with any other related MD are expected to enlighten distinct functions.

A Cysteine-Reactive Small Photo-Crosslinker Possessing Caged-Fluorescence Properties: Binding-Site Determination of a Combinatorially-Selected Peptide by Fluorescence Imaging/Tandem Mass Spectrometry

To determine the binding-site of a combinatorially-selected peptide possessing a fluoroprobe, a novel cysteine reactive small photo-crosslinker that can be excited by a conventional long-wavelength ultraviolet handlamp (365 nm) was synthesized via Suzuki coupling with three steps. The crosslinker is rationally designed, not only as a bioisostere of the fluoroprobe, but as a caged-fluorophore, and the photo-crosslinked target protein became fluorescent with a large Stokes-shift. By introducing the crosslinker to a designated sulfhydryl (SH) group of a combinatorially-selected peptide, the protein-binding site of the targeted peptide was deduced by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE)/fluorescence imaging followed by matrix-assisted laser desorption ionization-time of flight tandem mass spectrometry (MALDI-TOF-MS/MS) analysis.

Exploring a new ligand binding site of G protein-coupled receptors

Identifying a target ligand binding site is an important step for structure-based rational drug design as shown here for G protein-coupled receptors (GPCRs), which are among the most popular drug targets. We applied long-time scale molecular dynamics simulations, coupled with mutagenesis studies, to two prototypical GPCRs, the M3 and M4 muscarinic acetylcholine receptors. Our results indicate that unlike synthetic antagonists, which bind to the classic orthosteric site, the endogenous agonist acetylcholine is able to diffuse into a much deeper binding pocket. We also discovered that the most recently resolved crystal structure of the LTB4 receptor comprised a bound inverse agonist, which extended its benzamidine moiety to the same binding pocket discovered in this work. Analysis on all resolved GPCR crystal structures indicated that this new pocket could exist in most receptors. Our findings provide new opportunities for GPCR drug discovery.

Fragment-Based Discovery of Inhibitors of the Bacterial DnaG-SSB Interaction

In bacteria, the DnaG primase is responsible for synthesis of short RNA primers used to initiate chain extension by replicative DNA polymerase(s) during chromosomal replication. Among the proteins with which Escherichia coli DnaG interacts is the single-stranded DNA-binding protein, SSB. The C-terminal hexapeptide motif of SSB (DDDIPF; SSB-Ct) is highly conserved and is known to engage in essential interactions with many proteins in nucleic acid metabolism, including primase. Here, fragment-based screening by saturation-transfer difference nuclear magnetic resonance (STD-NMR) and surface plasmon resonance assays identified inhibitors of the primase/SSB-Ct interaction. Hits were shown to bind to the SSB-Ct-binding site using 15N-¹H HSQC spectra. STD-NMR was used to demonstrate binding of one hit to other SSB-Ct binding partners, confirming the possibility of simultaneous inhibition of multiple protein/SSB interactions. The fragment molecules represent promising scaffolds on which to build to discover new antibacterial compounds.

The NMR2 Method to Determine Rapidly the Structure of the Binding Pocket of a Protein–Ligand Complex with High Accuracy

Structural characterization of complexes is crucial for a better understanding of biological processes and structure-based drug design. However, many protein–ligand structures are not solvable by X-ray crystallography, for example those with low affinity binders or dynamic binding sites. Such complexes are usually targeted by solution-state NMR spectroscopy. Unfortunately, structure calculation by NMR is very time consuming since all atoms in the complex need to be assigned to their respective chemical shifts. To circumvent this problem, we recently developed the Nuclear Magnetic Resonance Molecular Replacement (NMR2) method. NMR2 very quickly provides the complex structure of a binding pocket as measured by solution-state NMR. NMR2 circumvents the assignment of the protein by using previously determined structures and therefore speeds up the whole process from a couple of months to a couple of days. Here, we recall the main aspects of the method, show how to apply it, discuss its advantages over other methods and outline its limitations and future directions.

Evidence for Novel Action at the Cell-Binding Site of Human Angiogenin Revealed by Heteronuclear NMR Spectroscopy, in silico and in vivo Studies

A member of the ribonuclease A superfamily, human angiogenin (hAng) is a potent angiogenic factor. Heteronuclear NMR spectroscopy combined with induced-fit docking revealed a dual binding mode for the most antiangiogenic compound of a series of ribofuranosyl pyrimidine nucleosides that strongly inhibit hAng's angiogenic activity in vivo. While modeling suggests the potential for simultaneous binding of the inhibitors at the active and cell-binding sites, NMR studies indicate greater affinity for the cell-binding site than for the active site. Additionally, molecular dynamics simulations at 100 ns confirmed the stability of binding at the cell-binding site with the predicted protein-ligand interactions, in excellent agreement with the NMR data. This is the first time that a nucleoside inhibitor is reported to completely inhibit the angiogenic activity of hAng in vivo by exerting dual inhibitory activity on hAng, blocking both the entrance of hAng into the cell and its ribonucleolytic activity.

An efficient combination of BEST and NUS methods in multidimensional NMR spectroscopy for high throughput analysis of proteins

Application of Non Uniform Sampling (NUS) along with Band-selective Excitation Short-Transient (BEST) NMR experiments has been demonstrated for obtaining the important residue-specific atomic level backbone chemical shift values in short durations of time. This application has been demonstrated with both well-folded (ubiquitin) and unfolded (α-synuclein) proteins alike. With this strategy, the experiments required for determining backbone chemical shifts can be performed very rapidly, i.e., in ∼2 hours of spectrometer time, and this data can be used to calculate the backbone folds of proteins using well established algorithms. This will be of great value for structural proteomic investigations on one hand, where the speed of structure determination is a limiting factor and for application in the study of slow kinetic processes involving proteins, such as fibrillization, on the other hand.

Application of NUS along with BEST NMR experiments has been demonstrated for obtaining the important residue-specific atomic level backbone chemical shift values in short durations of time.

Discovery of stimulator binding to a conserved pocket in the heme domain of soluble guanylyl cyclase

Soluble guanylyl cyclase (sGC) is the receptor for nitric oxide and a highly sought-after therapeutic target for the management of cardiovascular diseases. New compounds that stimulate sGC show clinical promise, but where these stimulator compounds bind and how they function remains unknown. Here, using a photolyzable diazirine derivative of a novel stimulator compound, IWP-051, and MS analysis, we localized drug binding to the β1 heme domain of sGC proteins from the hawkmoth Manduca sexta and from human. Covalent attachments to the stimulator were also identified in bacterial homologs of the sGC heme domain, referred to as H-NOX domains, including those from Nostoc sp. PCC 7120, Shewanella oneidensis, Shewanella woodyi, and Clostridium botulinum, indicating that the binding site is highly conserved. The identification of photoaffinity-labeled peptides was aided by a signature MS fragmentation pattern of general applicability for unequivocal identification of covalently attached compounds. Using NMR, we also examined stimulator binding to sGC from M. sexta and bacterial H-NOX homologs. These data indicated that stimulators bind to a conserved cleft between two subdomains in the sGC heme domain. L12W/T48W substitutions within the binding pocket resulted in a 9-fold decrease in drug response, suggesting that the bulkier tryptophan residues directly block stimulator binding. The localization of stimulator binding to the sGC heme domain reported here resolves the longstanding question of where stimulators bind and provides a path forward for drug discovery.

Emerging Diagnostic and Therapeutic Potentials of Human Hair Proteomics

The use of noninvasive human substrates to interrogate pathophysiological conditions has become essential in the post- Human Genome Project era. Due to its high turnover rate, and its long term capability to incorporate exogenous and endogenous substances from the circulation, hair testing is emerging as a key player in monitoring long term drug compliance, chronic alcohol abuse, forensic toxicology, and biomarker discovery, among other things. Novel high-throughput 'omics based approaches like proteomics have been underutilized globally in comprehending human hair morphology and its evolving use as a diagnostic testing substrate in the era of precision medicine. There is paucity of scientific evidence that evaluates the difference in drug incorporation into hair based on lipid content, and very few studies have addressed hair growth rates, hair forms, and the biological consequences of hair grooming or bleaching. It is apparent that protein-based identification using the human hair proteome would play a major role in understanding these parameters akin to DNA single nucleotide polymorphism profiling, up to single amino acid polymorphism resolution. Hence, this work seeks to identify and discuss the progress made thus far in the field of molecular hair testing using proteomic approaches, and identify ways in which proteomics would improve the field of hair research, considering that the human hair is mostly composed of proteins. Gaps in hair proteomics research are identified and the potential of hair proteomics in establishing a historic medical repository of normal and disease-specific proteome is also discussed.

HN, N, Cα and Cβ assignments of the two periplasmic domains of Neisseria meningitidis DsbD

DsbD is a disulfide bond reductase present in the inner membrane of many Gamma-Proteobacteria. In the human pathogen Neisseria meningitidis, DsbD is required for viability and represents a potential target for the development of antibiotics. Here we report the chemical shift assignments (H^N, N, C^α and C^β) for the reduced and oxidized forms of the two periplasmic domains of N. meningitidis DsbD, n-NmDsbD and c-NmDsbD. The backbone amide resonances in all four forms were completely assigned, and the secondary structures for the core regions of the proteins were calculated using ¹³C^αβ shifts. The reduced and oxidized forms of each domain have similar secondary shifts suggesting they retain the same fold. We anticipate that these data will provide an important basis for studying the interaction between n-NmDsbD and c-NmDsbD, which is required for electron transfer across the bacterial cytoplasmic membrane.

Comparative functional analysis of ribonuclease 1 homologs: molecular insights into evolving vertebrate physiology

Pancreatic-type ribonucleases (ptRNases) comprise a class of highly conserved secretory endoribonucleases in vertebrates. The prototype of this enzyme family is ribonuclease 1 (RNase 1). Understanding the physiological roles of RNase 1 is becoming increasingly important, as engineered forms of the enzyme progress through clinical trials as chemotherapeutic agents for cancer. Here, we present an in-depth biochemical characterization of RNase 1 homologs from a broad range of mammals (human, bat, squirrel, horse, cat, mouse, and cow) and nonmammalian species (chicken, lizard, and frog). We discover that the human homolog of RNase 1 has a pH optimum for catalysis, ability to degrade double-stranded RNA, and affinity for cell-surface glycans that are distinctly higher than those of its homologs. These attributes have relevance for human health. Moreover, the functional diversification of the 10 RNase 1 homologs illuminates the regulation of extracellular RNA and other aspects of vertebrate evolution.

Fragment library screening identifies hits that bind to the non-catalytic surface of Pseudomonas aeruginosa DsbA1

At a time when the antibiotic drug discovery pipeline has stalled, antibiotic resistance is accelerating with catastrophic implications for our ability to treat bacterial infections. Globally we face the prospect of a future when common infections can once again kill. Anti-virulence approaches that target the capacity of the bacterium to cause disease rather than the growth or survival of the bacterium itself offer a tantalizing prospect of novel antimicrobials. They may also reduce the propensity to induce resistance by removing the strong selection pressure imparted by bactericidal or bacteriostatic agents. In the human pathogen Pseudomonas aeruginosa, disulfide bond protein A (PaDsbA1) plays a central role in the oxidative folding of virulence factors and is therefore an attractive target for the development of new anti-virulence antimicrobials. Using a fragment-based approach we have identified small molecules that bind to PaDsbA1. The fragment hits show selective binding to PaDsbA1 over the DsbA protein from Escherichia coli, suggesting that developing species-specific narrow-spectrum inhibitors of DsbA enzymes may be feasible. Structures of a co-complex of PaDsbA1 with the highest affinity fragment identified in the screen reveal that the fragment binds on the non-catalytic surface of the protein at a domain interface. This biophysical and structural data represent a starting point in the development of higher affinity compounds, which will be assessed for their potential as selective PaDsbA1 inhibitors.

Polyoxometalates as sialidase mimics: selective and non-destructive removal of sialic acid from a glycoprotein promoted by phosphotungstic acid

The selective hydrolysis of the glycosidic bond between the terminal sialic acid and the penultimate sugar has been achieved in the alpha-2-HS-glycoprotein (Fetuin-A) in the presence of H₃PW₁₂O₄₀, a Keggin type polyoxometalate.

Determination of ligand binding modes in weak protein-ligand complexes using sparse NMR data

We describe a general approach to determine the binding pose of small molecules in weakly bound protein-ligand complexes by deriving distance constraints between the ligand and methyl groups from all methyl-containing residues of the protein. We demonstrate that using a single sample, which can be prepared without the use of expensive precursors, it is possible to generate high-resolution data rapidly and obtain the resonance assignments of Ile, Leu, Val, Ala and Thr methyl groups using triple resonance scalar correlation data. The same sample may be used to obtain Met ^εCH₃ assignments using NOESY-based methods, although the superior sensitivity of NOESY using [U-¹³C,¹⁵N]-labeled protein makes the use of this second sample more efficient. We describe a structural model for a weakly binding ligand bound to its target protein, DsbA, derived from intermolecular methyl-to-ligand nuclear Overhauser enhancements, and demonstrate that the ability to assign all methyl resonances in the spectrum is essential to derive an accurate model of the structure. Once the methyl assignments have been obtained, this approach provides a rapid means to generate structural models for weakly bound protein-ligand complexes. Such weak complexes are often found at the beginning of programs of fragment based drug design and can be challenging to characterize using X-ray crystallography.

A Sequence-Based Dynamic Ensemble Learning System for Protein Ligand-Binding Site Prediction

BACKGROUND

Proteins have the fundamental ability to selectively bind to other molecules and perform specific functions through such interactions, such as protein-ligand binding. Accurate prediction of protein residues that physically bind to ligands is important for drug design and protein docking studies. Most of the successful protein-ligand binding predictions were based on known structures. However, structural information is not largely available in practice due to the huge gap between the number of known protein sequences and that of experimentally solved structures.

RESULTS

This paper proposes a dynamic ensemble approach to identify protein-ligand binding residues by using sequence information only. To avoid problems resulting from highly imbalanced samples between the ligand-binding sites and non ligand-binding sites, we constructed several balanced data sets and we trained a random forest classifier for each of them. We dynamically selected a subset of classifiers according to the similarity between the target protein and the proteins in the training data set. The combination of the predictions of the classifier subset to each query protein target yielded the final predictions. The ensemble of these classifiers formed a sequence-based predictor to identify protein-ligand binding sites.

CONCLUSIONS

Experimental results on two Critical Assessment of protein Structure Prediction datasets and the ccPDB dataset demonstrated that of our proposed method compared favorably with the state-of-the-art.

AVAILABILITY

http://www2.ahu.edu.cn/pchen/web/LigandDSES.htm.

Secondary Structure Analysis of a Functional Construct of Caveolin-1 Reveals a Long C-Terminal Helix

Caveolin-1 is an integral membrane protein that is the primary component of cell membrane invaginations called caveolae. While caveolin-1 is known to participate in a myriad of vital cellular processes, structural data on caveolin-1 of any kind is severely limited. In order to rectify this dearth, secondary structure analysis of a functional construct of caveolin-1, containing the intact C-terminal domain, was performed using NMR spectroscopy in lyso-myristoylphosphatidylglycerol micelles. Complete backbone assignments of caveolin-1 (residues 62-178) were made, and it was determined that residues 62-79 were dynamic; residues 89-107, 111-128, and 132-175 were helical; and residues 80-88, 108-110, and 129-131 represent unstructured breaks between the helices.

The Involvement of His50 during Protein Disulfide Isomerase Binding Is Essential for Inhibiting α-Syn Fibril Formation

An increased level of protein disulfide isomerase (PDI) is a protective response to various neurodegenerative disorders, including Parkinson's disease. Interaction of PDI with α-synuclein (α-Syn) has been shown to inhibit its aggregation. Here, we report the residue-specific mapping of binding of PDI to α-Syn. We demonstrate that α-Syn N-terminal residues V3-S9 and L38-V40 bind more strongly to PDI than residues V49-V52 do, as do C-terminal residues E123-M127 and D135-E137. In addition, we show that residue H50 is key in preventing aggregation. These findings improve our understanding of PDI-protected aggregation of wild-type α-Syn and its H50Q familial mutant.

Orphan Nuclear Receptor NR4A1 Binds a Novel Protein Interaction Site on Anti-apoptotic B Cell Lymphoma Gene 2 Family Proteins

B cell lymphoma gene 2 (Bcl-2) family proteins are key regulators of programmed cell death and important targets for drug discovery. Pro-apoptotic and anti-apoptotic Bcl-2 family proteins reciprocally modulate their activities in large part through protein interactions involving a motif known as BH3 (Bcl-2 homology 3). Nur77 is an orphan member of the nuclear receptor family that lacks a BH3 domain but nevertheless binds certain anti-apoptotic Bcl-2 family proteins (Bcl-2, Bfl-1, and Bcl-B), modulating their effects on apoptosis and autophagy. We used a combination of NMR spectroscopy-based methods, mutagenesis, and functional studies to define the interaction site of a Nur77 peptide on anti-apoptotic Bcl-2 family proteins and reveal a novel interaction surface. Nur77 binds adjacent to the BH3 peptide-binding crevice, suggesting the possibility of cross-talk between these discrete binding sites. Mutagenesis of residues lining the identified interaction site on Bcl-B negated the interaction with Nur77 protein in cells and prevented Nur77-mediated modulation of apoptosis and autophagy. The findings establish a new protein interaction site with the potential to modulate the apoptosis and autophagy mechanisms governed by Bcl-2 family proteins.

Structure-Based Identification of Novel Ligands Targeting Multiple Sites within a Chemokine-G-Protein-Coupled-Receptor Interface

CXCL12 is a human chemokine that recognizes the CXCR4 receptor and is involved in immune responses and metastatic cancer. Interactions between CXCL12 and CXCR4 are an important drug target but, like other elongated protein-protein interfaces, present challenges for small molecule ligand discovery due to the relatively shallow and featureless binding surfaces. Calculations using an NMR complex structure revealed a binding hot spot on CXCL12 that normally interacts with the I4/I6 residues from CXCR4. Virtual screening was performed against the NMR model, and subsequent testing has verified the specific binding of multiple docking hits to this site. Together with our previous results targeting two other binding pockets that recognize sulfotyrosine residues (sY12 and sY21) of CXCR4, including a new analog against the sY12 binding site reported herein, we demonstrate that protein-protein interfaces can often possess multiple sites for engineering specific small molecule ligands that provide lead compounds for subsequent optimization by fragment based approaches.

Molecular Basis for the Interaction Between AP4 β4 and its Accessory Protein, Tepsin

The adaptor protein 4 (AP4) complex (ϵ/β4/μ4/σ4 subunits) forms a non-clathrin coat on vesicles departing the trans-Golgi network. AP4 biology remains poorly understood, in stark contrast to the wealth of molecular data available for the related clathrin adaptors AP1 and AP2. AP4 is important for human health because mutations in any AP4 subunit cause severe neurological problems, including intellectual disability and progressive spastic para- or tetraplegias. We have used a range of structural, biochemical and biophysical approaches to determine the molecular basis for how the AP4 β4 C-terminal appendage domain interacts with tepsin, the only known AP4 accessory protein. We show that tepsin harbors a hydrophobic sequence, LFxG[M/L]x[L/V], in its unstructured C-terminus, which binds directly and specifically to the C-terminal β4 appendage domain. Using nuclear magnetic resonance chemical shift mapping, we define the binding site on the β4 appendage by identifying residues on the surface whose signals are perturbed upon titration with tepsin. Point mutations in either the tepsin LFxG[M/L]x[L/V] sequence or in its cognate binding site on β4 abolish in vitro binding. In cells, the same point mutations greatly reduce the amount of tepsin that interacts with AP4. However, they do not abolish the binding between tepsin and AP4 completely, suggesting the existence of additional interaction sites between AP4 and tepsin. These data provide one of the first detailed mechanistic glimpses at AP4 coat assembly and should provide an entry point for probing the role of AP4-coated vesicles in cell biology, and especially in neuronal function.

NMR reveals structural rearrangements associated to substrate insertion in nucleotide-adding enzymes

The protein NP_344798.1 from Streptococcus pneumoniae TIGR4 exhibits a head and base-interacting neck domain architecture, as observed in class II nucleotide-adding enzymes. Although it has less than 20% overall sequence identity with any member of this enzyme family, the residues involved in substrate-recognition and catalysis are highly conserved in NP_344798.1. NMR studies showed binding affinity of NP_344798.1 for nucleotides and revealed μs to ms time scale rate processes involving residues constituting the active site. The results thus obtained indicate that large-amplitude rearrangements of regular secondary structures facilitate the penetration of the substrate into the occluded nucleotide-binding site of NP_344798.1 and, by inference based on sequence and structural homology, probably a wide range of other nucleotide-adding enzymes.

Mechanism of N-Acylthiourea-mediated activation of human histone deacetylase 8 (HDAC8) at molecular and cellular levels

We reported previously that an N-acylthiourea derivative (TM-2-51) serves as a potent and isozyme-selective activator for human histone deacetylase 8 (HDAC8). To probe the molecular mechanism of the enzyme activation, we performed a detailed account of the steady-state kinetics, thermodynamics, molecular modeling, and cell biology studies. The steady-state kinetic data revealed that TM-2-51 binds to HDAC8 at two sites in a positive cooperative manner. Isothermal titration calorimetric and molecular modeling data conformed to the two-site binding model of the enzyme-activator complex. We evaluated the efficacy of TM-2-51 on SH-SY5Y and BE(2)-C neuroblastoma cells, wherein the HDAC8 expression has been correlated with cellular malignancy. Whereas TM-2-51 selectively induced cell growth inhibition and apoptosis in SH-SY5Y cells, it showed no such effects in BE(2)-C cells, and this discriminatory feature appears to be encoded in the p53 genotype of the above cells. Our mechanistic and cellular studies on HDAC8 activation have the potential to provide insight into the development of novel anticancer drugs.

Protein-Inhibitor Interaction Studies Using NMR

Solution-state NMR has been widely applied to determine the three-dimensional structure, dynamics, and molecular interactions of proteins. The designs of experiments used in protein NMR differ from those used for small-molecule NMR, primarily because the information available prior to an experiment, such as molecular mass and knowledge of the primary structure, is unique for proteins compared to small molecules. In this review article, protein NMR for structural biology is introduced with comparisons to small-molecule NMR, such as descriptions of labeling strategies and the effects of molecular dynamics on relaxation. Next, applications for protein NMR are reviewed, especially practical aspects for protein-observed ligand-protein interaction studies. Overall, the following topics are described: (1) characteristics of protein NMR, (2) methods to detect protein-ligand interactions by NMR, and (3) practical aspects of carrying out protein-observed inhibitor-protein interaction studies.

LigandRFs: random forest ensemble to identify ligand-binding residues from sequence information alone

Background

Protein-ligand binding is important for some proteins to perform their functions. Protein-ligand binding sites are the residues of proteins that physically bind to ligands. Despite of the recent advances in computational prediction for protein-ligand binding sites, the state-of-the-art methods search for similar, known structures of the query and predict the binding sites based on the solved structures. However, such structural information is not commonly available.

Results

In this paper, we propose a sequence-based approach to identify protein-ligand binding residues. We propose a combination technique to reduce the effects of different sliding residue windows in the process of encoding input feature vectors. Moreover, due to the highly imbalanced samples between the ligand-binding sites and non ligand-binding sites, we construct several balanced data sets, for each of which a random forest (RF)-based classifier is trained. The ensemble of these RF classifiers forms a sequence-based protein-ligand binding site predictor.

Conclusions

Experimental results on CASP9 and CASP8 data sets demonstrate that our method compares favorably with the state-of-the-art protein-ligand binding site prediction methods.

Interaction studies of the human and Arabidopsis thaliana Med25-ACID proteins with the herpes simplex virus VP16- and plant-specific Dreb2a transcription factors

Mediator is an evolutionary conserved multi-protein complex present in all eukaryotes. It functions as a transcriptional co-regulator by conveying signals from activators and repressors to the RNA polymerase II transcription machinery. The Arabidopsis thaliana Med25 (aMed25) ACtivation Interaction Domain (ACID) interacts with the Dreb2a activator which is involved in plant stress response pathways, while Human Med25-ACID (hMed25) interacts with the herpes simplex virus VP16 activator. Despite low sequence similarity, hMed25-ACID also interacts with the plant-specific Dreb2a transcriptional activator protein. We have used GST pull-down-, surface plasmon resonance-, isothermal titration calorimetry and NMR chemical shift experiments to characterize interactions between Dreb2a and VP16, with the hMed25 and aMed25-ACIDs. We found that VP16 interacts with aMed25-ACID with similar affinity as with hMed25-ACID and that the binding surface on aMed25-ACID overlaps with the binding site for Dreb2a. We also show that the Dreb2a interaction region in hMed25-ACID overlaps with the earlier reported VP16 binding site. In addition, we show that hMed25-ACID/Dreb2a and aMed25-ACID/Dreb2a display similar binding affinities but different binding energetics. Our results therefore indicate that interaction between transcriptional regulators and their target proteins in Mediator are less dependent on the primary sequences in the interaction domains but that these domains fold into similar structures upon interaction.

Sulfopeptide probes of the CXCR4/CXCL12 interface reveal oligomer-specific contacts and chemokine allostery

Tyrosine sulfation is a post-translational modification that enhances protein-protein interactions and may identify druggable sites in the extracellular space. The G protein-coupled receptor CXCR4 is a prototypical example with three potential sulfation sites at positions 7, 12, and 21. Each receptor sulfotyrosine participates in specific contacts with its chemokine ligand in the structure of a soluble, dimeric CXCL12:CXCR4(1-38) complex, but their relative importance for CXCR4 binding and activation by the monomeric chemokine remains undefined. NMR titrations with short sulfopeptides showed that the tyrosine motifs of CXCR4 varied widely in their contributions to CXCL12 binding affinity and site specificity. Whereas the Tyr21 sulfopeptide bound the same site as in previously solved structures, the Tyr7 and Tyr12 sulfopeptides interacted nonspecifically. Surprisingly, the unsulfated Tyr7 peptide occupied a hydrophobic site on the CXCL12 monomer that is inaccessible in the CXCL12 dimer. Functional analysis of CXCR4 mutants validated the relative importance of individual CXCR4 sulfotyrosine modifications (Tyr21 > Tyr12 > Tyr7) for CXCL12 binding and receptor activation. Biophysical measurements also revealed a cooperative relationship between sulfopeptide binding at the Tyr21 site and CXCL12 dimerization, the first example of allosteric behavior in a chemokine. Future ligands that occupy the sTyr21 recognition site may act as both competitive inhibitors of receptor binding and allosteric modulators of chemokine function. Together, our data suggests that sulfation does not ubiquitously enhance complex affinity and that distinct patterns of tyrosine sulfation could encode oligomer selectivity, implying another layer of regulation for chemokine signaling.