Efficiently driving protein-based fragment screening and lead discovery using two-dimensional NMR.

Fragment-based drug discovery (FBDD) and validation of small molecule binders using NMR spectroscopy is an established and widely used method in the early stages of drug discovery. Starting from a library of small compounds, ligand- or protein-observed NMR methods are employed to detect binders, typically weak, that become the starting points for structure-activity relationships (SAR) by NMR. Unlike the more frequently used ligand-observed 1D NMR techniques, protein-observed 2D 1H-15N or 1H-13C heteronuclear correlation (HSQC or HMQC) methods offer insights that include the mechanism of ligand engagement on the target and direct binding affinity measurements in addition to routine screening. We hereby present the development of a set of software tools within the MestReNova (Mnova) package for analyzing 2D NMR for FBDD and hit validation purposes. The package covers three main tasks: (1) unsupervised profiling of raw data to identify outlier data points to exclude in subsequent analyses; (2) batch processing of single-point spectra to identify and rank binders based on chemical shift perturbations or spectral peak intensity changes; and (3) batch processing of multiple titration series to derive binding affinities (KD) by tracing the changes in peak locations or measuring global spectral changes. Toward this end, we implemented and evaluated a set of algorithms for automated peak tracing, spectral binning, and variance analysis by PCA, and a new tool for spectral data intensity comparison using ECHOS. The accuracy and speed of the tools are demonstrated on 2D NMR binding data collected on ligands used in the development of potential inhibitors of the anti-apoptotic MCL-1 protein.

A cryptic pocket in Ebola VP35 allosterically controls RNA binding.

Protein-protein and protein-nucleic acid interactions are often considered difficult drug targets because the surfaces involved lack obvious druggable pockets. Cryptic pockets could present opportunities for targeting these interactions, but identifying and exploiting these pockets remains challenging. Here, we apply a general pipeline for identifying cryptic pockets to the interferon inhibitory domain (IID) of Ebola virus viral protein 35 (VP35). VP35 plays multiple essential roles in Ebola's replication cycle but lacks pockets that present obvious utility for drug design. Using adaptive sampling simulations and machine learning algorithms, we predict VP35 harbors a cryptic pocket that is allosterically coupled to a key dsRNA-binding interface. Thiol labeling experiments corroborate the predicted pocket and mutating the predicted allosteric network supports our model of allostery. Finally, covalent modifications that mimic drug binding allosterically disrupt dsRNA binding that is essential for immune evasion. Based on these results, we expect this pipeline will be applicable to other proteins.

Opening of a cryptic pocket in ß-lactamase increases penicillinase activity.

Understanding the functional role of protein-excited states has important implications in protein design and drug discovery. However, because these states are difficult to find and study, it is still unclear if excited states simply result from thermal fluctuations and generally detract from function or if these states can actually enhance protein function. To investigate this question, we consider excited states in ß-lactamases and particularly a subset of states containing a cryptic pocket which forms under the Ω-loop. Given the known importance of the Ω-loop and the presence of this pocket in at least two homologs, we hypothesized that these excited states enhance enzyme activity. Using thiol-labeling assays to probe Ω-loop pocket dynamics and kinetic assays to probe activity, we find that while this pocket is not completely conserved across ß-lactamase homologs, those with the Ω-loop pocket have a higher activity against the substrate benzylpenicillin. We also find that this is true for TEM ß-lactamase variants with greater open Ω-loop pocket populations. We further investigate the open population using a combination of NMR chemical exchange saturation transfer experiments and molecular dynamics simulations. To test our understanding of the Ω-loop pocket's functional role, we designed mutations to enhance/suppress pocket opening and observed that benzylpenicillin activity is proportional to the probability of pocket opening in our designed variants. The work described here suggests that excited states containing cryptic pockets can be advantageous for function and may be favored by natural selection, increasing the potential utility of such cryptic pockets as drug targets.

Coupled intra- and interdomain dynamics support domain cross-talk in Pin1.

The functional mechanisms of multidomain proteins often exploit interdomain interactions, or "cross-talk." An example is human Pin1, an essential mitotic regulator consisting of a Trp-Trp (WW) domain flexibly tethered to a peptidyl-prolyl isomerase (PPIase) domain, resulting in interdomain interactions important for Pin1 function. Substrate binding to the WW domain alters its transient contacts with the PPIase domain via means that are only partially understood. Accordingly, we have investigated Pin1 interdomain interactions using NMR paramagnetic relaxation enhancement (PRE) and molecular dynamics (MD) simulations. The PREs show that apo-Pin1 samples interdomain contacts beyond the range suggested by previous structural studies. They further show that substrate binding to the WW domain simultaneously alters interdomain separation and the internal conformation of the WW domain. A 4.5-µs all-atom MD simulation of apo-Pin1 suggests that the fluctuations of interdomain distances are correlated with fluctuations of WW domain interresidue contacts involved in substrate binding. Thus, the interdomain/WW domain conformations sampled by apo-Pin1 may already include a range of conformations appropriate for binding Pin1's numerous substrates. The proposed coupling between intra-/interdomain conformational fluctuations is a consequence of the dynamic modular architecture of Pin1. Such modular architecture is common among cell-cycle proteins; thus, the WW-PPIase domain cross-talk mechanisms of Pin1 may be relevant for their mechanisms as well.

The Cap-Snatching SFTSV Endonuclease Domain Is an Antiviral Target.

Severe fever with thrombocytopenia syndrome virus (SFTSV) is a tick-borne virus with 12%-30% case mortality rates and is related to the Heartland virus (HRTV) identified in the United States. Together, SFTSV and HRTV are emerging segmented, negative-sense RNA viral (sNSV) pathogens with potential global health impact. Here, we characterize the amino-terminal cap-snatching endonuclease domain of SFTSV polymerase (L) and solve a 2.4-Å X-ray crystal structure. While the overall structure is similar to those of other cap-snatching sNSV endonucleases, differences near the C terminus of the SFTSV endonuclease suggest divergence in regulation. Influenza virus endonuclease inhibitors, including the US Food and Drug Administration (FDA) approved Baloxavir (BXA), inhibit the endonuclease activity in in vitro enzymatic assays and in cell-based studies. BXA displays potent activity with a half maximal inhibitory concentration (IC50) of ∼100 nM in enzyme inhibition and an EC50 value of ∼250 nM against SFTSV and HRTV in plaque assays. Together, our data support sNSV endonucleases as an antiviral target.

A gratuitous ß-Lactamase inducer uncovers hidden active site dynamics of the Staphylococcus aureus BlaR1 sensor domain.

Increasing evidence shows that active sites of proteins have non-trivial conformational dynamics. These dynamics include active site residues sampling different local conformations that allow for multiple, and possibly novel, inhibitor binding poses. Yet, active site dynamics garner only marginal attention in most inhibitor design efforts and exert little influence on synthesis strategies. This is partly because synthesis requires a level of atomic structural detail that is frequently missing in current characterizations of conformational dynamics. In particular, while the identity of the mobile protein residues may be clear, the specific conformations they sample remain obscure. Here, we show how an appropriate choice of ligand can significantly sharpen our abilities to describe the interconverting binding poses (conformations) of protein active sites. Specifically, we show how 2-(2'-carboxyphenyl)-benzoyl-6-aminopenicillanic acid (CBAP) exposes otherwise hidden dynamics of a protein active site that binds ß-lactam antibiotics. When CBAP acylates (binds) the active site serine of the ß-lactam sensor domain of BlaR1 (BlaRS), it shifts the time scale of the active site dynamics to the slow exchange regime. Slow exchange enables direct characterization of inter-converting protein and bound ligand conformations using NMR methods. These methods include chemical shift analysis, 2-d exchange spectroscopy, off-resonance ROESY of the bound ligand, and reduced spectral density mapping. The active site architecture of BlaRS is shared by many ß-lactamases of therapeutic interest, suggesting CBAP could expose functional motions in other ß-lactam binding proteins. More broadly, CBAP highlights the utility of identifying chemical probes common to structurally homologous proteins to better expose functional motions of active sites.

Electron Cryo-microscopy Structure of Ebola Virus Nucleoprotein Reveals a Mechanism for Nucleocapsid-like Assembly.

Ebola virus nucleoprotein (eNP) assembles into higher-ordered structures that form the viral nucleocapsid (NC) and serve as the scaffold for viral RNA synthesis. However, molecular insights into the NC assembly process are lacking. Using a hybrid approach, we characterized the NC-like assembly of eNP, identified novel regulatory elements, and described how these elements impact function. We generated a three-dimensional structure of the eNP NC-like assembly at 5.8 Å using electron cryo-microscopy and identified a new regulatory role for eNP helices α22-α23. Biochemical, biophysical, and mutational analyses revealed that inter-eNP contacts within α22-α23 are critical for viral NC assembly and regulate viral RNA synthesis. These observations suggest that the N terminus and α22-α23 of eNP function as context-dependent regulatory modules (CDRMs). Our current study provides a framework for a structural mechanism for NC-like assembly and a new therapeutic target.

Endogenous retinoid X receptor ligands in mouse hematopoietic cells.

The retinoid X receptor α (RXRA) has been implicated in diverse hematological processes. To identify natural ligands of RXRA that are present in hematopoietic cells, we adapted an upstream activation sequence-green fluorescent protein (UAS-GFP) reporter mouse to detect natural RXRA ligands in vivo. We observed reporter activity in diverse types of hematopoietic cells in vivo. Reporter activity increased during granulocyte colony-stimulating factor (G-CSF)-induced granulopoiesis and after phenylhydrazine (PHZ)-induced anemia, suggesting the presence of dynamically regulated natural RXRA ligands in hematopoietic cells. Mouse plasma activated Gal4-UAS reporter cells in vitro, and plasma from mice treated with G-CSF or PHZ recapitulated the patterns of reporter activation that we observed in vivo. Plasma from mice with dietary vitamin A deficiency only mildly reduced RXRA reporter activity, whereas plasma from mice on a fatty acid restriction diet reduced reporter activity, implicating fatty acids as plasma RXRA ligands. Through differential extraction coupled with mass spectrometry, we identified the long-chain fatty acid C24:5 as a natural RXRA ligand that was greatly increased in abundance in response to hematopoietic stress. Together, these data suggest that natural RXRA ligands are present and dynamically increased in abundance in mouse hematopoietic cells in vivo.

Repurposing Hsp90 inhibitors as antibiotics targeting histidine kinases.

To address the growing need for new antimicrobial agents, we explored whether inhibition of bacterial signaling machinery could inhibit bacterial growth. Because bacteria rely on two-component signaling systems to respond to environmental changes, and because these systems are both highly conserved and mediated by histidine kinases, inhibiting histidine kinases may provide broad spectrum antimicrobial activity. The histidine kinase ATP binding domain is conserved with the ATPase domain of eukaryotic Hsp90 molecular chaperones. To find a chemical scaffold for compounds that target histidine kinases, we leveraged this conservation. We screened ATP competitive Hsp90 inhibitors against CckA, an essential histidine kinase in Caulobacter crescentus that controls cell growth, and showed that the diaryl pyrazole is a promising scaffold for histidine kinase inhibition. We synthesized a panel of derivatives and found that they inhibit the histidine kinases C. crescentus CckA and Salmonella PhoQ but not C. crescentus DivJ; and they inhibit bacterial growth in both Gram-negative and Gram-positive bacterial strains.

Designing small molecules to target cryptic pockets yields both positive and negative allosteric modulators.

Allosteric drugs, which bind to proteins in regions other than their main ligand-binding or active sites, make it possible to target proteins considered "undruggable" and to develop new therapies that circumvent existing resistance. Despite growing interest in allosteric drug discovery, rational design is limited by a lack of sufficient structural information about alternative binding sites in proteins. Previously, we used Markov State Models (MSMs) to identify such "cryptic pockets," and here we describe a method for identifying compounds that bind in these cryptic pockets and modulate enzyme activity. Experimental tests validate our approach by revealing both an inhibitor and two activators of TEM ß-lactamase (TEM). To identify hits, a library of compounds is first virtually screened against either the crystal structure of a known cryptic pocket or an ensemble of structures containing the same cryptic pocket that is extracted from an MSM. Hit compounds are then screened experimentally and characterized kinetically in individual assays. We identify three hits, one inhibitor and two activators, demonstrating that screening for binding to allosteric sites can result in both positive and negative modulation. The hit compounds have modest effects on TEM activity, but all have higher affinities than previously identified inhibitors, which bind the same cryptic pocket but were found, by chance, via a computational screen targeting the active site. Site-directed mutagenesis of key contact residues predicted by the docking models is used to confirm that the compounds bind in the cryptic pocket as intended. Because hit compounds are identified from docking against both the crystal structure and structures from the MSM, this platform should prove suitable for many proteins, particularly targets whose crystal structures lack obvious druggable pockets, and for identifying both inhibitory and activating small-molecule modulators.

Prediction of New Stabilizing Mutations Based on Mechanistic Insights from Markov State Models.

Protein stabilization is fundamental to enzyme function and evolution, yet understanding the determinants of a protein's stability remains a challenge. This is largely due to a shortage of atomically detailed models for the ensemble of relevant protein conformations and their relative populations. For example, the M182T substitution in TEM ß-lactamase, an enzyme that confers antibiotic resistance to bacteria, is stabilizing but the precise mechanism remains unclear. Here, we employ Markov state models (MSMs) to uncover how M182T shifts the distribution of different structures that TEM adopts. We find that M182T stabilizes a helix that is a key component of a domain interface. We then predict the effects of other mutations, including a novel stabilizing mutation, and experimentally test our predictions using a combination of stability measurements, crystallography, NMR, and in vivo measurements of bacterial fitness. We expect our insights and methodology to provide a valuable foundation for protein design.

Investigation of signal transduction routes within the sensor/transducer protein BlaR1 of Staphylococcus aureus.

The transmembrane antibiotic sensor/signal transducer protein BlaR1 is part of a cohort of proteins that confer ß-lactam antibiotic resistance in methicillin-resistant Staphylococcus aureus (MRSA) [Fisher, J. F., Meroueh, S. O., and Mobashery, S. (2005) Chem. Rev. 105, 395-424; Llarrull, L. I., Fisher, J. F., and Mobashery, S. (2009) Antimicrob. Agents Chemother. 53, 4051-4063; Llarrull, L. I., Toth, M., Champion, M. M., and Mobashery, S. (2011) J. Biol. Chem. 286, 38148-38158]. Specifically, BlaR1 regulates the inducible expression of ß-lactamases that hydrolytically destroy ß-lactam antibiotics. The resistance phenotype starts with ß-lactam antibiotic acylation of the BlaR1 extracellular domain (BlaRS). The acylation activates the cytoplasmic protease domain through an obscure signal transduction mechanism. Here, we compare protein dynamics of apo versus antibiotic-acylated BlaRS using nuclear magnetic resonance. Our analyses reveal inter-residue interactions that relay acylation-induced perturbations within the antibiotic-binding site to the transmembrane helix regions near the membrane surface. These are the first insights into the process of signal transduction by BlaR1.

Revealing cell-surface intramolecular interactions in the BlaR1 protein of methicillin-resistant Staphylococcus aureus by NMR spectroscopy.

In methicillin-resistant Staphylococcus aureus, ß-lactam antibiotic resistance is mediated by the transmembrane protein BlaR1. The antibiotic sensor domain BlaR(S) and the L2 loop of BlaR1 are on the membrane surface. We used NMR to investigate interactions between BlaR(S) and a water-soluble peptide from L2. This peptide binds BlaR(S) proximal to the antibiotic acylation site as an amphipathic helix. Acylation of BlaR(S) by penicillin G does not disrupt binding. These results suggest a signal transduction mechanism whereby the L2 helix, partially embedded in the membrane, propagates conformational changes caused by BlaR(S) acylation through the membrane via transmembrane segments, leading to antibiotic resistance.

Effects of the endosomal lipid bis(monoacylglycero)phosphate on the thermotropic properties of DPPC: A 2H NMR and spin label EPR study.

Bis(monoacylglycero)phosphate (BMP) is an endosomal lipid with a unique structure that is implicated in the formation of intraendosomal vesicular bodies. Here we have characterized the effects of dioleoyl-BMP (BMP(18:1)) at concentrations of 5, 10, 15 and 20mol% on the thermotropic behavior of dipalmitoyl phosphatidylcholine (DPPC) vesicles, and compared them to those of equimolar concentrations of dioleoyl phosphatidylglycerol (DOPG), a structural isoform of BMP(18:1). Because BMP is found in the acidic environments of the late endosome and intralysosomal vesicles, samples were prepared at pH 4.2 to mimic the pH of the lysosome. Both (2)H NMR of perdeuterated DPPC and spin-labeled EPR with 16-doxyl phosphatidylcholine were utilized in these investigations. NMR and EPR results show that BMP(18:1) induces a lowering in the main phase transition temperature of DPPC similar to that of DOPG. The EPR studies reveal that BMP(18:1) induced more disorder in the L(beta) phase when compared to equimolar concentrations of DOPG. Analysis from dePaked (2)H NMR spectra in the L(alpha) phase reveals that BMP(18:1) induces less disorder than equal concentrations of DOPG. Additionally, the results demonstrate that BMP mixes with other phospholipids as a phospholipid and not as a detergent molecule as once speculated.

Bis(monoacylglycero)phosphate and ganglioside GM1 spontaneously form small homogeneous vesicles at specific concentrations.

The morphology and size of hydrated lipid dispersions of bis(monoacylglycero)phosphate (BMP) mixed with varying mole percentages of the ganglioside GM1 were investigated by dynamic light scattering (DLS) and transmission electron microscopy (TEM). Electron paramagnetic resonance (EPR) spectroscopy of these same mixtures, doped at 0.5 mol% with doxyl labeled lipids, was used to investigate acyl-chain packing. Results show that for 20-30% GM1, hydrated BMP:GM1 mixtures spontaneously form small spherical vesicles with diameters approximately 100 nm and a narrow size distribution profile. For other concentrations of GM1, hydrated dispersions with BMP have non-spherical shapes and heterogeneous size profiles, with average vesicle diameters>400 nm. All samples were prepared at pH 5.5 to mimic the lumen acidity of the late endosome where BMP is an essential component of intraendosomal vesicle budding, lipid sorting and trafficking. These findings indicate that GM1 and BMP under a limited concentration range spontaneously form small vesicles of homogeneous size in an energy independent manner without the need of protein templating. Because BMP is essential for intraendosomal vesicle formation, these results imply that lipid-lipid interactions may play a critical role in the endosomal process of lipid sorting and trafficking.

Bis(monoacylglycero)phosphate forms stable small lamellar vesicle structures: insights into vesicular body formation in endosomes.

Bis(monoacylglycero)phosphate (BMP) is an unusually shaped lipid found in relatively high percentage in the late endosome. Here, we report the characterization of the morphology and molecular organization of dioleoyl-BMP (DOBMP) with dynamic light scattering, transmission electron microscopy, nuclear magnetic resonance (NMR) spectroscopy, and electron paramagnetic resonance spectroscopy. The morphology of hydrated DOBMP dispersions varies with pH and ionic strength, and DOBMP vesicles are significantly smaller in diameter than phosphatidylcholine dispersions. At neutral pH, DOBMP forms highly structured, clustered dispersions 500 nm in size. On the other hand, at acidic pH, spherically shaped vesicles are formed. NMR and spin-labeled electron paramagnetic resonance demonstrate that DOBMP forms a lamellar mesophase with acyl-chain packing similar to that of other unsaturated phospholipids. (31)P NMR reveals an orientation of the phosphate group in DOBMP that differs significantly from that of other phospholipids. These macroscopic and microscopic structural characterizations suggest that the biosynthesis of BMP on the inner luminal membrane of maturing endosomes may possibly produce budded vesicles high in BMP content, which form small vesicular structures stabilized by the physical properties of the BMP lipid.