Communication breakdown: protein dynamics and drug design

Environmentally Ultrasensitive Fluorine Probe to Resolve Protein Conformational Ensembles by 19F NMR and Cryo-EM

Limited chemical shift dispersion represents a significant barrier to studying multistate equilibria of large membrane proteins by 19F NMR. We describe a novel monofluoroethyl 19F probe that dramatically increases the chemical shift dispersion. The improved conformational sensitivity and line shape enable the detection of previously unresolved states in one-dimensional (1D) 19F NMR spectra of a 134 kDa membrane transporter. Changes in the populations of these states in response to ligand binding, mutations, and temperature correlate with population changes of distinct conformations in structural ensembles determined by single-particle cryo-electron microscopy (cryo-EM). Thus, 19F NMR can guide sample preparation to discover and visualize novel conformational states and facilitate image analysis and three-dimensional (3D) classification.

Exploring the Potential Mechanisms of Melilotus officinalis (L.) Pall. in Chronic Muscle Repair Patterns Using Single Cell Receptor-Ligand Marker Analysis and Molecular Dynamics Simulations

Information regarding the function of Melilotus officinalis (L.) Pall. in skeletal muscles is still unknown. In this study, we explored the possible regulatory targets of M. (L.) Pall. that affects the repair patterns in chronic muscle injury. We analyzed the potential target genes and chemical composition of M. (L.) Pall. and constructed a "drug-component-disease target genes" network analysis. Five active ingredients and 87 corresponding targets were obtained. Muscle-tendon junction (MTJ) cells were used to perform receptor-ligand marker analysis using the CellphoneDB algorithm. Targets of M. (L.) Pall. were screened further for the cellular ligand-receptor protein action on MTJs. Enrichment analysis suggests that those protein-associated ligand receptors may be associated with a range of intercellular signaling pathways. Molecular docking validation was then performed. Five proteins (CCL2, VEGFA, MMP2, MET, and EGFR) may be regulated by the active ingredient luteolin and scoparone. Finally, molecular dynamics simulations revealed that luteolin can stably target binding to MMP2. M. (L.) Pall. influences skeletal muscle repair patterns by affecting the fibroblast interactions in the muscle-tendon junctions through the active ingredients luteolin and scoparone.

Distal Regions Regulate Dihydrofolate Reductase-Ligand Interactions

Protein motions play a fundamental role in enzyme catalysis and ligand binding. The relationship between protein motion and function has been extensively investigated in the model enzyme dihydrofolate reductase (DHFR). DHFR is an essential enzyme that catalyzes the reduction of dihydrofolate to tetrahydrofolate. Numerous experimental and computational methods have been used to probe the motions of DHFR through the catalytic cycle and to investigate the effect of distal mutations on DHFR motions and ligand binding. These experimental investigations have pushed forward the study of protein motions and their role in protein-ligand interactions. The introduction of mutations distal to the active site has been shown to have profound effects on ligand binding, hydride transfer rates and catalytic efficacy and these changes are captured by enzyme kinetics measurements. Distal mutations have been shown to exert their effects through a network of correlated amino acids and these effects have been investigated by NMR, protein dynamics, and analysis of coupled amino acids. The experimental methods and the findings that are reviewed here have broad implications for our understanding of enzyme mechanisms, ligand binding and for the future design and discovery of enzyme inhibitors.

What is the structural chemistry of the living organism at its temperature and pressure?

The three probes of the structure of matter (X-rays, neutrons and electrons) in biology have complementary properties and strengths. The balance between these three probes within their strengths and weaknesses is perceived to change, even dramatically so at times. For the study of combined states of order and disorder, NMR crystallography is also applicable. Of course, to understand biological systems the required perspectives are surely physiologically relevant temperatures and relevant chemical conditions, as well as a minimal perturbation owing to the needs of the probe itself. These remain very tough challenges because, for example, cryoEM by its very nature will never be performed at room temperature, crystallization often requires nonphysiological chemical conditions, and X-rays and electrons cause beam damage. However, integrated structural biology techniques and functional assays provide a package towards physiological relevance of any given study. Reporting of protein crystal structures, and their associated database entries, could usefully indicate how close to the biological situation they are, as discussed in detail in this feature article.

Perspectives on electrostatics and conformational motions in enzyme catalysis

Enzymes are essential for all living organisms, and their effectiveness as chemical catalysts has driven more than a half century of research seeking to understand the enormous rate enhancements they provide. Nevertheless, a complete understanding of the factors that govern the rate enhancements and selectivities of enzymes remains elusive, due to the extraordinary complexity and cooperativity that are the hallmarks of these biomolecules. We have used a combination of site-directed mutagenesis, pre-steady-state kinetics, X-ray crystallography, nuclear magnetic resonance (NMR), vibrational and fluorescence spectroscopies, resonance energy transfer, and computer simulations to study the implications of conformational motions and electrostatic interactions on enzyme catalysis in the enzyme dihydrofolate reductase (DHFR).

We have demonstrated that modest equilibrium conformational changes are functionally related to the hydride transfer reaction. Results obtained for mutant DHFRs illustrated that reductions in hydride transfer rates are correlated with altered conformational motions, and analysis of the evolutionary history of DHFR indicated that mutations appear to have occurred to preserve both the hydride transfer rate and the associated conformational changes. More recent results suggested that differences in local electrostatic environments contribute to finely tuning the substrate pK_a in the initial protonation step. Using a combination of primary and solvent kinetic isotope effects, we demonstrated that the reaction mechanism is consistent across a broad pH range, and computer simulations suggested that deprotonation of the active site Tyr100 may play a crucial role in substrate protonation at high pH.

Site-specific incorporation of vibrational thiocyanate probes into the ecDHFR active site provided an experimental tool for interrogating these microenvironments and for investigating changes in electrostatics along the DHFR catalytic cycle. Complementary molecular dynamics simulations in conjunction with mixed quantum mechanical/molecular mechanical calculations accurately reproduced the vibrational frequency shifts in these probes and provided atomic-level insight into the residues influencing these changes. Our findings indicate that conformational and electrostatic changes are intimately related and functionally essential. This approach can be readily extended to the study of other enzyme systems to identify more general trends in the relationship between conformational fluctuations and electrostatic interactions. These results are relevant to researchers seeking to design novel enzymes as well as those seeking to develop therapeutic agents that function as enzyme inhibitors.

Conformation-based restrictions and scaffold replacements in the design of hepatitis C virus polymerase inhibitors: discovery of deleobuvir (BI 207127)

Conformational restrictions of flexible torsion angles were used to guide the identification of new chemotypes of HCV NS5B inhibitors. Sites for rigidification were based on an acquired conformational understanding of compound binding requirements and the roles of substituents in the free and bound states. Chemical bioisosteres of amide bonds were explored to improve cell-based potency. Examples are shown, including the design concept that led to the discovery of the phase III clinical candidate deleobuvir (BI 207127). The structure-based strategies employed have general utility in drug design.

Ligand bioactive conformation plays a critical role in the design of drugs that target the hepatitis C virus NS3 protease

A ligand-focused strategy employed NMR, X-ray, modeling, and medicinal chemistry to expose the critical role that bioactive conformation played in the design of a variety of drugs that target the HCV protease. The bioactive conformation (bound states) were determined for key inhibitors identified along our drug discovery pathway from the hit to clinical compounds. All adopt similar bioactive conformations for the common core derived from the hit peptide DDIVPC. A carefully designed SAR analysis, based on the advanced inhibitor 1 in which the P1 to P3 side chains and the N-terminal Boc were sequentially truncated, revealed a correlation between affinity and the relative predominance of the bioactive conformation in the free state. Interestingly, synergistic conformation effects on potency were also noted. Comparisons with clinical and recently marketed drugs from the pharmaceutical industry showed that all have the same core and similar bioactive conformations. This suggested that the variety of appendages discovered for these compounds also properly satisfy the bioactive conformation requirements and allowed for a large variety of HCV protease drug candidates to be designed.

Ligand and receptor dynamics contribute to the mechanism of graded PPARγ agonism

Ligand binding to proteins is not a static process, but rather involves a number of complex dynamic transitions. A flexible ligand can change conformation upon binding its target. The conformation and dynamics of a protein can change to facilitate ligand binding. The conformation of the ligand, however, is generally presumed to have one primary binding mode, shifting the protein conformational ensemble from one state to another. We report solution nuclear magnetic resonance (NMR) studies that reveal peroxisome proliferator-activated receptor γ (PPARγ) modulators can sample multiple binding modes manifesting in multiple receptor conformations in slow conformational exchange. Our NMR, hydrogen/deuterium exchange and docking studies reveal that ligand-induced receptor stabilization and binding mode occupancy correlate with the graded agonist response of the ligand. Our results suggest that ligand and receptor dynamics affect the graded transcriptional output of PPARγ modulators.

Direct and allosteric inhibition of the FGF2/HSPGs/FGFR1 ternary complex formation by an antiangiogenic, thrombospondin-1-mimic small molecule

Fibroblast growth factors (FGFs) are recognized targets for the development of therapies against angiogenesis-driven diseases, including cancer. The formation of a ternary complex with the transmembrane tyrosine kinase receptors (FGFRs), and heparan sulphate proteoglycans (HSPGs) is required for FGF2 pro-angiogenic activity. Here by using a combination of techniques including Nuclear Magnetic Resonance, Molecular Dynamics, Surface Plasmon Resonance and cell-based binding assays we clarify the molecular mechanism of inhibition of an angiostatic small molecule, sm27, mimicking the endogenous inhibitor of angiogenesis, thrombospondin-1. NMR and MD data demonstrate that sm27 engages the heparin-binding site of FGF2 and induces long-range dynamics perturbations along FGF2/FGFR1 interface regions. The functional consequence of the inhibitor binding is an impaired FGF2 interaction with both its receptors, as demonstrated by SPR and cell-based binding assays. We propose that sm27 antiangiogenic activity is based on a twofold-direct and allosteric-mechanism, inhibiting FGF2 binding to both its receptors.

Rational prediction with molecular dynamics for hit identification

Although the motions of proteins are fundamental for their function, for pragmatic reasons, the consideration of protein elasticity has traditionally been neglected in drug discovery and design. This review details protein motion, its relevance to biomolecular interactions and how it can be sampled using molecular dynamics simulations. Within this context, two major areas of research in structure-based prediction that can benefit from considering protein flexibility, binding site detection and molecular docking, are discussed. Basic classification metrics and statistical analysis techniques, which can facilitate performance analysis, are also reviewed. With hardware and software advances, molecular dynamics in combination with traditional structure-based prediction methods can potentially reduce the time and costs involved in the hit identification pipeline.

Automated NMR relaxation dispersion data analysis using NESSY

Background

Proteins are dynamic molecules with motions ranging from picoseconds to longer than seconds. Many protein functions, however, appear to occur on the micro to millisecond timescale and therefore there has been intense research of the importance of these motions in catalysis and molecular interactions. Nuclear Magnetic Resonance (NMR) relaxation dispersion experiments are used to measure motion of discrete nuclei within the micro to millisecond timescale. Information about conformational/chemical exchange, populations of exchanging states and chemical shift differences are extracted from these experiments. To ensure these parameters are correctly extracted, accurate and careful analysis of these experiments is necessary.

Results

The software introduced in this article is designed for the automatic analysis of relaxation dispersion data and the extraction of the parameters mentioned above. It is written in Python for multi platform use and highest performance. Experimental data can be fitted to different models using the Levenberg-Marquardt minimization algorithm and different statistical tests can be used to select the best model. To demonstrate the functionality of this program, synthetic data as well as NMR data were analyzed. Analysis of these data including the generation of plots and color coded structures can be performed with minimal user intervention and using standard procedures that are included in the program.

Conclusions

NESSY is easy to use open source software to analyze NMR relaxation data. The robustness and standard procedures are demonstrated in this article.

Direct detection of structurally resolved dynamics in a multiconformation receptor-ligand complex

Structure-based drug design relies on static protein structures despite significant evidence for the need to include protein dynamics as a serious consideration. In practice, dynamic motions are neglected because they are not understood well enough to model, a situation resulting from a lack of explicit experimental examples of dynamic receptor-ligand complexes. Here, we report high-resolution details of pronounced ~1 ms time scale motions of a receptor-small molecule complex using a combination of NMR and X-ray crystallography. Large conformational dynamics in Escherichia coli dihydrofolate reductase are driven by internal switching motions of the drug-like, nanomolar-affinity inhibitor. Carr-Purcell-Meiboom-Gill relaxation dispersion experiments and NOEs revealed the crystal structure to contain critical elements of the high energy protein-ligand conformation. The availability of accurate, structurally resolved dynamics in a protein-ligand complex should serve as a valuable benchmark for modeling dynamics in other receptor-ligand complexes and prediction of binding affinities.

Importance of ligand bioactive conformation in the discovery of potent indole-diamide inhibitors of the hepatitis C virus NS5B

Significant advances have led to receptor induced-fit and conformational selection models for describing bimolecular recognition, but a more comprehensive view must evolve to also include ligand shape and conformational changes. Here, we describe an example where a ligand's "structural hinge" influences potency by inducing an "L-shape" bioactive conformation, and due to its solvent exposure in the complex, reasonable conformation-activity-relationships can be qualitatively attributed. From a ligand design perspective, this feature was exploited by successful linker hopping to an alternate "structural hinge" that led to a new and promising chemical series which matched the ligand bioactive conformation and the pocket bioactive space. Using a combination of X-ray crystallography, NMR and modeling with support from binding-site resistance mutant studies and photoaffinity labeling experiments, we were able to derive inhibitor-polymerase complexes for various chemical series.

Application of NMR and molecular docking in structure-based drug discovery

Drug discovery is a complex and costly endeavor, where few drugs that reach the clinical testing phase make it to market. High-throughput screening (HTS) is the primary method used by the pharmaceutical industry to identify initial lead compounds. Unfortunately, HTS has a high failure rate and is not particularly efficient at identifying viable drug leads. These shortcomings have encouraged the development of alternative methods to drive the drug discovery process. Specifically, nuclear magnetic resonance (NMR) spectroscopy and molecular docking are routinely being employed as important components of drug discovery research. Molecular docking provides an extremely rapid way to evaluate likely binders from a large chemical library with minimal cost. NMR ligand-affinity screens can directly detect a protein-ligand interaction, can measure a corresponding dissociation constant, and can reliably identify the ligand binding site and generate a co-structure. Furthermore, NMR ligand affinity screens and molecular docking are perfectly complementary techniques, where the combination of the two has the potential to improve the efficiency and success rate of drug discovery. This review will highlight the use of NMR ligand affinity screens and molecular docking in drug discovery and describe recent examples where the two techniques were combined to identify new and effective therapeutic drugs.

An introduction to NMR-based approaches for measuring protein dynamics

Proteins are inherently flexible at ambient temperature. At equilibrium, they are characterized by a set of conformations that undergo continuous exchange within a hierarchy of spatial and temporal scales ranging from nanometers to micrometers and femtoseconds to hours. Dynamic properties of proteins are essential for describing the structural bases of their biological functions including catalysis, binding, regulation and cellular structure. Nuclear magnetic resonance (NMR) spectroscopy represents a powerful technique for measuring these essential features of proteins. Here we provide an introduction to NMR-based approaches for studying protein dynamics, highlighting eight distinct methods with recent examples, contextualized within a common experimental and analytical framework. The selected methods are (1) Real-time NMR, (2) Exchange spectroscopy, (3) Lineshape analysis, (4) CPMG relaxation dispersion, (5) Rotating frame relaxation dispersion, (6) Nuclear spin relaxation, (7) Residual dipolar coupling, (8) Paramagnetic relaxation enhancement. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Dynameomics: a comprehensive database of protein dynamics

The dynamic behavior of proteins is important for an understanding of their function and folding. We have performed molecular dynamics simulations of the native state and unfolding pathways of over 2000 protein/peptide systems (approximately 11,000 independent simulations) representing the majority of folds in globular proteins. These data are stored and organized using an innovative database approach, which can be mined to obtain both general and specific information about the dynamics and folding/unfolding of proteins, relevant subsets thereof, and individual proteins. Here we describe the project in general terms and the type of information contained in the database. Then we provide examples of mining the database for information relevant to protein folding, structure building, the effect of single-nucleotide polymorphisms, and drug design. The native state simulation data and corresponding analyses for the 100 most populated metafolds, together with related resources, are publicly accessible through http://www.dynameomics.org.

Advances in Nuclear Magnetic Resonance for Drug Discovery

BACKGROUND: Drug discovery is a complex and unpredictable endeavor with a high failure rate. Current trends in the pharmaceutical industry have exasperated these challenges and are contributing to the dramatic decline in productivity observed over the last decade. The industrialization of science by forcing the drug discovery process to adhere to assembly-line protocols is imposing unnecessary restrictions, such as short project time-lines. Recent advances in nuclear magnetic resonance are responding to these self-imposed limitations and are providing opportunities to increase the success rate of drug discovery. OBJECTIVE/METHOD: A review of recent advancements in NMR technology that have the potential of significantly impacting and benefiting the drug discovery process will be presented. These include fast NMR data collection protocols and high-throughput protein structure determination, rapid protein-ligand co-structure determination, lead discovery using fragment-based NMR affinity screens, NMR metabolomics to monitor in vivo efficacy and toxicity for lead compounds, and the identification of new therapeutic targets through the functional annotation of proteins by FAST-NMR. CONCLUSION: NMR is a critical component of the drug discovery process, where the versatility of the technique enables it to continually expand and evolve its role. NMR is expected to maintain this growth over the next decade with advancements in automation, speed of structure calculation, in-cell imaging techniques, and the expansion of NMR amenable targets.