Targeting the Main Protease (Mpro, nsp5) by Growth of Fragment Scaffolds Exploiting Structure-Based Methodologies.

The main protease Mpro, nsp5, of SARS-CoV-2 (SCoV2) is one of its most attractive drug targets. Here, we report primary screening data using nuclear magnetic resonance spectroscopy (NMR) of four different libraries and detailed follow-up synthesis on the promising uracil-containing fragment Z604 derived from these libraries. Z604 shows time-dependent binding. Its inhibitory effect is sensitive to reducing conditions. Starting with Z604, we synthesized and characterized 13 compounds designed by fragment growth strategies. Each compound was characterized by NMR and/or activity assays to investigate their interaction with Mpro. These investigations resulted in the four-armed compound 35b that binds directly to Mpro. 35b could be cocrystallized with Mpro revealing its noncovalent binding mode, which fills all four active site subpockets. Herein, we describe the NMR-derived fragment-to-hit pipeline and its application for the development of promising starting points for inhibitors of the main protease of SCoV2.

Survivin prevents the polycomb repressor complex 2 from methylating histone 3 lysine 27.

This study investigates the role of survivin in epigenetic control of gene transcription through interaction with the polycomb repressive complex 2 (PRC2). PRC2 is responsible for silencing gene expression by trimethylating lysine 27 on histone 3. We observed differential expression of PRC2 subunits in CD4+ T cells with varying levels of survivin expression, and ChIP-seq results indicated that survivin colocalizes with PRC2 along DNA. Inhibition of survivin resulted in a significant increase in H3K27 trimethylation, implying that survivin prevents PRC2 from functioning. Peptide microarray showed that survivin interacts with peptides from PRC2 subunits, and machine learning revealed that amino acid composition contains relevant information for predicting survivin interaction. NMR and BLI experiments supported the interaction of survivin with PRC2 subunit EZH2. Finally, protein-protein docking revealed that the survivin-EZH2 interaction interface overlaps with catalytic residues of EZH2, potentially inhibiting its H3K27 methylation activity. These findings suggest that survivin inhibits PRC2 function.

Entropy-Entropy Compensation between the Protein, Ligand, and Solvent Degrees of Freedom Fine-Tunes Affinity in Ligand Binding to Galectin-3C.

Molecular recognition is fundamental to biological signaling. A central question is how individual interactions between molecular moieties affect the thermodynamics of ligand binding to proteins and how these effects might propagate beyond the immediate neighborhood of the binding site. Here, we investigate this question by introducing minor changes in ligand structure and characterizing the effects of these on ligand affinity to the carbohydrate recognition domain of galectin-3, using a combination of isothermal titration calorimetry, X-ray crystallography, NMR relaxation, and computational approaches including molecular dynamics (MD) simulations and grid inhomogeneous solvation theory (GIST). We studied a congeneric series of ligands with a fluorophenyl-triazole moiety, where the fluorine substituent varies between the ortho, meta, and para positions (denoted O, M, and P). The M and P ligands have similar affinities, whereas the O ligand has 3-fold lower affinity, reflecting differences in binding enthalpy and entropy. The results reveal surprising differences in conformational and solvation entropy among the three complexes. NMR backbone order parameters show that the O-bound protein has reduced conformational entropy compared to the M and P complexes. By contrast, the bound ligand is more flexible in the O complex, as determined by 19F NMR relaxation, ensemble-refined X-ray diffraction data, and MD simulations. Furthermore, GIST calculations indicate that the O-bound complex has less unfavorable solvation entropy compared to the other two complexes. Thus, the results indicate compensatory effects from ligand conformational entropy and water entropy, on the one hand, and protein conformational entropy, on the other hand. Taken together, these different contributions amount to entropy-entropy compensation among the system components involved in ligand binding to a target protein.

Signaling Mechanism of Phytochromes in Solution.

Phytochrome proteins guide the red/far-red photoresponse of plants, fungi, and bacteria. Crystal structures suggest that the mechanism of signal transduction from the chromophore to the output domains involves refolding of the so-called PHY tongue. It is currently not clear how the two other notable structural features of the phytochrome superfamily, the so-called helical spine and a knot in the peptide chain, are involved in photoconversion. Here, we present solution NMR data of the complete photosensory core module from Deinococcus radiodurans. Photoswitching between the resting and the active states induces changes in amide chemical shifts, residual dipolar couplings, and relaxation dynamics. All observables indicate a photoinduced structural change in the knot region and lower part of the helical spine. This implies that a conformational signal is transduced from the chromophore to the helical spine through the PAS and GAF domains. The discovered pathway underpins functional studies of plant phytochromes and may explain photosensing by phytochromes under biological conditions.

Modulation of Structural Heterogeneity Controls Phytochrome Photoswitching.

Phytochromes sense red/far-red light and control many biological processes in plants, fungi, and bacteria. Although the crystal structures of dark- and light-adapted states have been determined, the molecular mechanisms underlying photoactivation remain elusive. Here, we demonstrate that the conserved tongue region of the PHY domain of a 57-kDa photosensory module of Deinococcus radiodurans phytochrome changes from a structurally heterogeneous dark state to an ordered, light-activated state. The results were obtained in solution by utilizing a laser-triggered activation approach detected on the atomic level with high-resolution protein NMR spectroscopy. The data suggest that photosignaling of phytochromes relies on careful modulation of structural heterogeneity of the PHY tongue.

Assessing the Ability of Spectroscopic Methods to Determine the Difference in the Folding Propensities of Highly Similar ß-Hairpins.

We have evaluated the ability of nuclear magnetic resonance (NMR) and circular dichroism (CD) spectroscopies to describe the difference in the folding propensities of two structurally highly similar cyclic ß-hairpins, comparing the outcome to that of molecular dynamics simulations. NAMFIS-type NMR ensemble analysis and CD spectroscopy were observed to accurately describe the consequence of altering a single interaction site, whereas a single-site 13C NMR chemical shift melting curve-based technique was not.

Paramagnetic Ligand Tagging To Identify Protein Binding Sites.

Transient biomolecular interactions are the cornerstones of the cellular machinery. The identification of the binding sites for low affinity molecular encounters is essential for the development of high affinity pharmaceuticals from weakly binding leads but is hindered by the lack of robust methodologies for characterization of weakly binding complexes. We introduce a paramagnetic ligand tagging approach that enables localization of low affinity protein-ligand binding clefts by detection and analysis of intermolecular protein NMR pseudocontact shifts, which are invoked by the covalent attachment of a paramagnetic lanthanoid chelating tag to the ligand of interest. The methodology is corroborated by identification of the low millimolar volatile anesthetic interaction site of the calcium sensor protein calmodulin. It presents an efficient route to binding site localization for low affinity complexes and is applicable to rapid screening of protein-ligand systems with varying binding affinity.

The impact of interchain hydrogen bonding on ß-hairpin stability is readily predicted by molecular dynamics simulation.

Peptides are frequently used model systems for protein folding. They are also gaining increased importance as therapeutics. Here, the ability of molecular dynamics (MD) simulation for describing the structure and dynamics of ß-hairpin peptides was investigated, with special attention given to the impact of a single interstrand sidechain to sidechain interaction. The MD trajectories were compared to structural information gained from solution NMR. By assigning frames from restraint-free MD simulations to an intuitive hydrogen bond on/off pattern, folding ratios and folding pathways were predicted. The computed molecular model successfully reproduces the folding ratios determined by NMR, indicating that MD simulation may be straightforwardly used as a screening tool in ß-hairpin design.

Mapping the sevoflurane-binding sites of calmodulin.

General anesthetics, with sevoflurane (SF) being the first choice inhalational anesthetic agent, provide reversible, broad depressor effects on the nervous system yet have a narrow margin of safety. As characterization of low-affinity binding interactions of volatile substances is exceptionally challenging with the existing methods, none of the numerous cellular targets proposed as chief protagonists in anesthesia could yet be confirmed. The recognition that most critical functions modulated by volatile anesthetics are under the control of intracellular Ca(2+) concentration, which in turn is primarily regulated by calmodulin (CaM), motivated us for characterization of the SF-CaM interaction. Solution NMR (Nuclear Magnetic Resonance) spectroscopy was used to identify SF-binding sites using chemical shift displacement, NOESY and heteronuclear Overhauser enhancement spectroscopy (HOESY) experiments. Binding affinities were measured using ITC (isothermal titration calorimetry). SF binds to both lobes of (Ca(2+))4-CaM with low mmol/L affinity whereas no interaction was observed in the absence of Ca(2+). SF does not affect the calcium binding of CaM. The structurally closely related SF and isoflurane are shown to bind to the same clefts. The SF-binding clefts overlap with the binding sites of physiologically relevant ion channels and bioactive small molecules, but the binding affinity suggests it could only interfere with very weak CaM targets.

Off-resonance rotating-frame relaxation dispersion experiment for 13C in aromatic side chains using L-optimized TROSY-selection.

Protein dynamics on the microsecond-millisecond time scales often play a critical role in biological function. NMR relaxation dispersion experiments are powerful approaches for investigating biologically relevant dynamics with site-specific resolution, as shown by a growing number of publications on enzyme catalysis, protein folding, ligand binding, and allostery. To date, the majority of studies has probed the backbone amides or side-chain methyl groups, while experiments targeting other sites have been used more sparingly. Aromatic side chains are useful probes of protein dynamics, because they are over-represented in protein binding interfaces, have important catalytic roles in enzymes, and form a sizable part of the protein interior. Here we present an off-resonance R 1ρ experiment for measuring microsecond to millisecond conformational exchange of aromatic side chains in selectively (13)C labeled proteins by means of longitudinal- and transverse-relaxation optimization. Using selective excitation and inversion of the narrow component of the (13)C doublet, the experiment achieves significant sensitivity enhancement in terms of both signal intensity and the fractional contribution from exchange to transverse relaxation; additional signal enhancement is achieved by optimizing the longitudinal relaxation recovery of the covalently attached (1)H spins. We validated the L-TROSY-selected R 1ρ experiment by measuring exchange parameters for Y23 in bovine pancreatic trypsin inhibitor at a temperature of 328 K, where the ring flip is in the fast exchange regime with a mean waiting time between flips of 320 µs. The determined chemical shift difference matches perfectly with that measured from the NMR spectrum at lower temperatures, where separate peaks are observed for the two sites. We further show that potentially complicating effects of strong scalar coupling between protons (Weininger et al. in J Phys Chem B 117: 9241-9247, 2013b) can be accounted for using a simple expression, and provide recommendations for data acquisition when the studied system exhibits this behavior. The present method extends the repertoire of relaxation methods tailored for aromatic side chains by enabling studies of faster processes and improved control over artifacts due to strong coupling.

Symmetric halogen bonding is preferred in solution.

Halogen bonding is a recently rediscovered secondary interaction that shows potential to become a complementary molecular tool to hydrogen bonding in rational drug design and in material sciences. Whereas hydrogen bond symmetry has been the subject of systematic studies for decades, the understanding of the analogous three-center halogen bonds is yet in its infancy. The isotopic perturbation of equilibrium (IPE) technique with (13)C NMR detection was applied to regioselectively deuterated pyridine complexes to investigate the symmetry of [N-I-N](+) and [N-Br-N](+) halogen bonding in solution. Preference for a symmetric arrangement was observed for both a freely adjustable and for a conformationally restricted [N-X-N](+) model system, as also confirmed by computation on the DFT level. A closely attached counterion is shown to be compatible with the preferred symmetric arrangement. The experimental observations and computational predictions reveal a high energetic gain upon formation of symmetric, three-center four-electron halogen bonding. Whereas hydrogen bonds are generally asymmetric in solution and symmetric in the crystalline state, the analogous bromine and iodine centered halogen bonds prefer symmetric arrangement in solution.

Protein loop compaction and the origin of the effect of arginine and glutamic acid mixtures on solubility, stability and transient oligomerization of proteins.

Addition of a 50 mM mixture of L: -arginine and L: -glutamic acid (RE) is extensively used to improve protein solubility and stability, although the origin of the effect is not well understood. We present Small Angle X-ray Scattering (SAXS) and Nuclear Magnetic Resonance (NMR) results showing that RE induces protein compaction by collapsing flexible loops on the protein core. This is suggested to be a general mechanism preventing aggregation and improving resistance to proteases and to originate from the polyelectrolyte nature of RE. Molecular polyelectrolyte mixtures are expected to display long range correlation effects according to dressed interaction site theory. We hypothesize that perturbation of the RE solution by dissolved proteins is proportional to the volume occupied by the protein. As a consequence, loop collapse, minimizing the effective protein volume, is favored in the presence of RE.

Differential responses of the backbone and side-chain conformational dynamics in FKBP12 upon binding the transition-state analog FK506: implications for transition-state stabilization and target protein recognition.

FKBP12 serves a dual role as a peptidyl-prolyl cis-trans isomerase and as a modulator of several cell signaling pathways. The macrolide FK506 is a transition-state analog of the catalyzed reaction and displaces FKBP12 from its natural target proteins. We compared the conformational exchange dynamics of the backbone and methyl-bearing side chains of FKBP12 in the free and FK506-bound states using NMR relaxation-dispersion experiments. Our results show that the free enzyme exchanges between the ground state and an excited state that resembles the ligand-bound state or Michaelis complex. In FK506-bound FKBP12, the backbone is confined to a single conformation, while conformational exchange prevails for many methyl groups. The residual side-chain dynamics in the transition-state analog-bound state suggests that the transition-state ensemble involves multiple conformations, a finding that challenges the long-standing concept of conformational restriction in the transition-state complex. Furthermore, exchange between alternative conformations is observed in the bound state for an extended network of methyl groups that includes locations remote from the active site. Several of these locations are known to be important for interactions with cellular target proteins, including calcineurin and the ryanodine receptor, suggesting that the conformational heterogeneity might play a role in the promiscuous binding of FKBP12 to different targets.

Functional dynamics of human FKBP12 revealed by methyl 13C rotating frame relaxation dispersion NMR spectroscopy.

Transverse relaxation dispersion NMR spectroscopy can provide atom-specific information about time scales, populations, and the extent of structural reorganization in proteins under equilibrium conditions. A method is described that uses side-chain methyl groups as local reporters for conformational transitions taking place in the microsecond regime. The experiment measures carbon nuclear spin relaxation rates in the presence of continuous wave off-resonance irradiation, in proteins uniformly enriched with 13C, and partially randomly labeled with 2H. The method was applied to human FK-506 binding protein (FKBP12), which uses a common surface for binding substrates in its dual role as both an immunophilin and folding assistant. Conformational dynamics on a time scale of approximately 130 micros were detected for methyl groups located in the substrate binding pocket, demonstrating its plasticity in the absence of substrate. The spatial arrangement of affected side-chain atoms suggests that substrate recognition involves the rapid relative movement of the subdomain comprising residues Ala81-Thr96 and that the observed dynamics play an important role in facilitating the interaction of this protein with its many partners, including calcineurin.

Biosynthetic 13C labeling of aromatic side chains in proteins for NMR relaxation measurements.

Site-specific 13C labeling offers a desirable means of eliminating unwanted relaxation pathways and coherent magnetization transfer in NMR relaxation experiments. Here we use [1-13C]-glucose as the sole carbon source in the growth media for protein overexpression in Escherichia coli. The approach results in specific incorporation of 13C at isolated positions in the side chains of aromatic amino acids, which greatly simplifies the measurements and interpretation of 13C relaxation rates in these spin systems. The method is well suited for characterization of chemical exchange by CPMG or spin-lock relaxation methods. We validated the method by acquiring 13C rotating-frame relaxation dispersion data on the E140Q mutant of the C-terminal domain of calmodulin, which reveal conformational exchange dynamics with a time constant of 71 mus for Y138.