Alignment of weakly interacting molecules to protein surfaces using simulations of chemical shift perturbations

A Step toward NRF2-DNA Interaction Inhibitors by Fragment-Based NMR Methods

The NRF2 transcription factor is a key regulator in cellular oxidative stress response, and acts as a tumor suppressor. Aberrant activation of NRF2 has been implicated in promoting chemo-resistance, tumor growth, and metastasis by activating its downstream target genes. Hence, inhibition of NRF2 promises to be an attractive therapeutic strategy to suppress cell proliferation and enhance cell apoptosis in cancer. Direct targeting of NRF2 with small-molecules to discover protein-DNA interaction inhibitors is challenging as it is a largely intrinsically disordered protein. To discover molecules that bind to NRF2 at the DNA binding interface, we performed an NMR-based fragment screen against its DNA-binding domain. We discovered several weakly binding fragment hits that bind to a region overlapping with the DNA binding site. Using SAR by catalogue we developed an initial structure-activity relationship for the most interesting initial hit series. By combining NMR chemical shift perturbations and data-driven docking, binding poses which agreed with NMR information and the observed SAR were elucidated. The herein discovered NRF2 hits and proposed binding modes form the basis for future structure-based optimization campaigns on this important but to date 'undrugged' cancer driver.

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.

Automation of peak-tracking analysis of stepwise perturbed NMR spectra

We describe a new algorithmic approach able to automatically pick and track the NMR resonances of a large number of 2D NMR spectra acquired during a stepwise variation of a physical parameter. The method has been named Trace in Track (TINT), referring to the idea that a gaussian decomposition traces peaks within the tracks recognised through 3D mathematical morphology. It is capable of determining the evolution of the chemical shifts, intensity and linewidths of each tracked peak.The performances obtained in term of track reconstruction and correct assignment on realistic synthetic spectra were high above 90% when a noise level similar to that of experimental data were considered. TINT was applied successfully to several protein systems during a temperature ramp in isotope exchange experiments. A comparison with a state-of-the-art algorithm showed promising results for great numbers of spectra and low signal to noise ratios, when the graduality of the perturbation is appropriate. TINT can be applied to different kinds of high throughput chemical shift mapping experiments, with quasi-continuous variations, in which a quantitative automated recognition is crucial.

An NMR-Guided Screening Method for Selective Fragment Docking and Synthesis of a Warhead Inhibitor

Selective hits for the glutaredoxin ortholog of Brucella melitensis are determined using STD NMR and verified by trNOE and (15)N-HSQC titration. The most promising hit, RK207, was docked into the target molecule using a scoring function to compare simulated poses to experimental data. After elucidating possible poses, the hit was further optimized into the lead compound by extension with an electrophilic acrylamide warhead. We believe that focusing on selectivity in this early stage of drug discovery will limit cross-reactivity that might occur with the human ortholog as the lead compound is optimized. Kinetics studies revealed that lead compound 5 modified with an ester group results in higher reactivity than an acrylamide control; however, after modification this compound shows little selectivity for bacterial protein versus the human ortholog. In contrast, hydrolysis of compound 5 to the acid form results in a decrease in the activity of the compound. Together these results suggest that more optimization is warranted for this simple chemical scaffold, and opens the door for discovery of drugs targeted against glutaredoxin proteins-a heretofore untapped reservoir for antibiotic agents.

Protein-ligand structure guided by backbone and side-chain proton chemical shift perturbations

The fragment-based drug design approach consists of screening libraries of fragment-like ligands, to identify hits that typically bind the protein target with weak affinity (100 μM-5 mM). The determination of the protein-fragment complex 3D structure constitutes a crucial step for uncovering the key interactions responsible for the protein-ligand recognition, and for growing the initial fragment into potent active compounds. The vast majority of fragments are aromatic compounds that induce chemical shift perturbations (CSP) on protein NMR spectra. These experimental CSPs can be quantitatively used to guide the ligand docking, through the comparison between experimental CSPs and CSP back-calculation based on the ring current effect. Here we implemented the CSP back-calculation into the scoring function of the program PLANTS. We compare the results obtained with CSPs measured either on amide or aliphatic protons of the human peroxiredoxin 5. We show that the different kinds of protons lead to different results for resolving the 3D structures of protein-fragment complexes, with the best results obtained with the Hα protons.

Performance of protein-ligand docking with simulated chemical shift perturbations

Protein chemical shift perturbations (CSPs) that result from the binding of a ligand to the protein contain structural information about the complex. Therefore, the CSP data, typically obtained during library screening from two-dimensional (2D) nuclear magnetic resonance (NMR) spectra, are often available before attempts to solve the experimental structure of the complex are started, and can be used to solve the complex structure with CSP-based docking. Here, we compare the performance of the post-docking filter and the guided-docking approaches using either amide or α-proton CSPs with 10 protein-ligand complexes. We show that the comparison of experimental CSPs with CSPs simulated for virtual ligand positions can be used to evidence protein conformational change upon binding and possibly improve the CSP-based docking.

Comparing binding modes of analogous fragments using NMR in fragment-based drug design: application to PRDX5

Fragment-based drug design is one of the most promising approaches for discovering novel and potent inhibitors against therapeutic targets. The first step of the process consists of identifying fragments that bind the protein target. The determination of the fragment binding mode plays a major role in the selection of the fragment hits that will be processed into drug-like compounds. Comparing the binding modes of analogous fragments is a critical task, not only to identify specific interactions between the protein target and the fragment, but also to verify whether the binding mode is conserved or differs according to the fragment modification. While X-ray crystallography is the technique of choice, NMR methods are helpful when this fails. We show here how the ligand-observed saturation transfer difference (STD) experiment and the protein-observed 15N-HSQC experiment, two popular NMR screening experiments, can be used to compare the binding modes of analogous fragments. We discuss the application and limitations of these approaches based on STD-epitope mapping, chemical shift perturbation (CSP) calculation and comparative CSP sign analysis, using the human peroxiredoxin 5 as a protein model.

Using chemical shift perturbation to characterise ligand binding

Chemical shift perturbation (CSP, chemical shift mapping or complexation-induced changes in chemical shift, CIS) follows changes in the chemical shifts of a protein when a ligand is added, and uses these to determine the location of the binding site, the affinity of the ligand, and/or possibly the structure of the complex. A key factor in determining the appearance of spectra during a titration is the exchange rate between free and bound, or more specifically the off-rate koff. When koff is greater than the chemical shift difference between free and bound, which typically equates to an affinity Kd weaker than about 3μM, then exchange is fast on the chemical shift timescale. Under these circumstances, the observed shift is the population-weighted average of free and bound, which allows Kd to be determined from measurement of peak positions, provided the measurements are made appropriately. (1)H shifts are influenced to a large extent by through-space interactions, whereas (13)Cα and (13)Cβ shifts are influenced more by through-bond effects. (15)N and (13)C' shifts are influenced both by through-bond and by through-space (hydrogen bonding) interactions. For determining the location of a bound ligand on the basis of shift change, the most appropriate method is therefore usually to measure (15)N HSQC spectra, calculate the geometrical distance moved by the peak, weighting (15)N shifts by a factor of about 0.14 compared to (1)H shifts, and select those residues for which the weighted shift change is larger than the standard deviation of the shift for all residues. Other methods are discussed, in particular the measurement of (13)CH3 signals. Slow to intermediate exchange rates lead to line broadening, and make Kd values very difficult to obtain. There is no good way to distinguish changes in chemical shift due to direct binding of the ligand from changes in chemical shift due to allosteric change. Ligand binding at multiple sites can often be characterised, by simultaneous fitting of many measured shift changes, or more simply by adding substoichiometric amounts of ligand. The chemical shift changes can be used as restraints for docking ligand onto protein. By use of quantitative calculations of ligand-induced chemical shift changes, it is becoming possible to determine not just the position but also the orientation of ligands.

BcL-xL conformational changes upon fragment binding revealed by NMR

Protein-protein interactions represent difficult but increasingly important targets for the design of therapeutic compounds able to interfere with biological processes. Recently, fragment-based strategies have been proposed as attractive approaches for the elaboration of protein-protein surface inhibitors from fragment-like molecules. One major challenge in targeting protein-protein interactions is related to the structural adaptation of the protein surface upon molecular recognition. Methods capable of identifying subtle conformational changes of proteins upon fragment binding are therefore required at the early steps of the drug design process. In this report we present a fast NMR method able to probe subtle conformational changes upon fragment binding. The approach relies on the comparison of experimental fragment-induced Chemical Shift Perturbation (CSP) of amine protons to CSP simulated for a set of docked fragment poses, considering the ring-current effect from fragment binding. We illustrate the method by the retrospective analysis of the complex between the anti-apoptotic Bcl-xL protein and the fragment 4′-fluoro-[1,1′-biphenyl]-4-carboxylic acid that was previously shown to bind one of the Bcl-xL hot spots. The CSP-based approach shows that the protein undergoes a subtle conformational rearrangement upon interaction, for residues located in helices 2, 3 and the very beginning of 5. Our observations are corroborated by residual dipolar coupling measurements performed on the free and fragment-bound forms of the Bcl-xL protein. These NMR-based results are in total agreement with previous molecular dynamic calculations that evidenced a high flexibility of Bcl-xL around the binding site. Here we show that CSP of protein amine protons are useful and reliable structural probes. Therefore, we propose to use CSP simulation to assess protein conformational changes upon ligand binding in the fragment-based drug design approach.

Quantitative use of chemical shifts for the modeling of protein complexes

Despite recent advances in the modeling of protein-protein complexes by docking, additional information is often required to identify the best solutions. For this purpose, NMR data deliver valuable restraints that can be used in the sampling and/or the scoring stage, like in the data-driven docking approach HADDOCK that can make use of NMR chemical shift perturbation (CSP) data to define the binding site on each protein and drive the docking. We show here that a quantitative use of chemical shifts (CS) in the scoring stage can help to resolve ambiguities. A quantitative CS-RMSD score based on (1) H(α) ,(13) C(α) , and (15) N chemical shifts ranks the best solutions always at the top, as demonstrated on a small benchmark of four complexes. It is implemented in a new docking protocol, CS-HADDOCK, which combines CSP data as ambiguous interaction restraints in the sampling stage with the CS-RMSD score in the final scoring stage. This combination of qualitative and quantitative use of chemical shifts increases the reliability of data-driven docking for the structure determination of complexes from limited NMR data.

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.

Discovery of fragment molecules that bind the human peroxiredoxin 5 active site

The search for protein ligands is a crucial step in the inhibitor design process. Fragment screening represents an interesting method to rapidly find lead molecules, as it enables the exploration of a larger portion of the chemical space with a smaller number of compounds as compared to screening based on drug-sized molecules. Moreover, fragment screening usually leads to hit molecules that form few but optimal interactions with the target, thus displaying high ligand efficiencies. Here we report the screening of a homemade library composed of 200 highly diverse fragments against the human Peroxiredoxin 5 protein. Peroxiredoxins compose a family of peroxidases that share the ability to reduce peroxides through a conserved cysteine. The three-dimensional structures of these enzymes ubiquitously found throughout evolution have been extensively studied, however, their biological functions are still not well understood and to date few inhibitors have been discovered against these enzymes. Six fragments from the library were shown to bind to the Peroxiredoxin 5 active site and ligand-induced chemical shift changes were used to drive the docking of these small molecules into the protein structure. The orientation of the fragments in the binding pocket was confirmed by the study of fragment homologues, highlighting the role of hydroxyl functions that hang the ligands to the Peroxiredoxin 5 protein. Among the hit fragments, the small catechol molecule was shown to significantly inhibit Peroxiredoxin 5 activity in a thioredoxin peroxidase assay. This study reports novel data about the ligand-Peroxiredoxin interactions that will help considerably the development of potential Peroxiredoxin inhibitors.

Influence of conformational flexibility on complexation-induced changes in chemical shift in a neocarzinostatin protein-ligand complex

In this paper is described an analysis of the effects of protein flexibility on the observed CIS values and the impact on the accuracy of 3D structures determined using a (1)H NMR CIS approach. The effects of protein conformational mobility have been investigated by using a set of different protein structures as starting points for the calculation: the unbound X-ray crystal structure, the unbound NMR solution structure, and the bound NMR solution structure of the protein. The results indicated that loop movement does have a significant impact on the quality of the structure generated by the CIS structure determination methodology. The implementation of methods to treat loop flexibility within our protocol, however, did not improve the results for calculations based on the unbound protein frame.

Determination of protein-ligand binding modes using complexation-induced changes in (1)h NMR chemical shift

A new method for determining three-dimensional solution structures of protein-ligand complexes using experimentally determined complexation-induced changes in (1)H NMR chemical shift (CIS) is introduced. The method has been validated using the complex formed between the protein antitumor antibiotic neocarzinostatin (NCS) and a synthetic chromophore analogue. The X-ray crystal structure of the unbound protein and the backbone amide proton CIS were the input data used in the determination of the three-dimensional structure of the complex. The experimental CIS values were used in a continuous direct structure refinement process based on genetic algorithms to sample conformational space. The calculated structure of the complex agrees well with the NMR solution structure, indicating the potential of this approach for structure determination.

NMR-Based Strategies to Elucidate Bioactive Conformations of Weakly Binding Ligands

Key processes in molecular biology are regulated by interactions between biomolecules. Protein-proteinand protein-ligand interactions, e.g., in signal transduction pathways, rely on the subtle interactionsbetween atoms at the binding interface of the involved molecules. Because biomolecules often havemany interacting partners, these interactions are not necessarily strong. The study of molecularrecognition gives insight into the complex network of signaling in life and is the basis of structure-baseddrug design.In the situation where the interaction is weak, one of the traditional methods that can be appliedto obtain structural information (internuclear distances) of the bound ligand is the so-called transferredNOE (trNOE) method. Recently, it became possible to use transferred cross-correlated relaxation (trCCR)to directly measure dihedral angles. The combined use of these two techniques significantly improvesthe precision of the structure determination of ligands weakly bound to macromolecules.The application of these techniques will be discussed in detail for a peptide derived fromIKKβ bound to the protein NEMO that plays an important rolein the NFκB pathway.

Protein-ligand NOE matching: a high-throughput method for binding pose evaluation that does not require protein NMR resonance assignments

Given the three-dimensional (3D) structure of a protein, the binding pose of a ligand can be determined using distance restraints derived from assigned intra-ligand and protein-ligand nuclear Overhauser effects (NOEs). A primary limitation of this approach is the need for resonance assignments of the ligand-bound protein. We have developed an approach that utilizes data from 3D 13C-edited, 13C/15N-filtered HSQC-NOESY spectra for evaluating ligand binding poses without requiring protein NMR resonance assignments. Only the 1H NMR assignments of the bound ligand are essential. Trial ligand binding poses are generated by any suitable method (e.g., computational docking). For each trial binding pose, the 3D 13C-edited, 13C/15N-filtered HSQC-NOESY spectrum is predicted, and the predicted and observed patterns of protein-ligand NOEs are matched and scored using a fast, deterministic bipartite graph matching algorithm. The best scoring (lowest "cost") poses are identified. Our method can incorporate any explicit restraints or protein assignment data that are available, and many extensions of the basic procedure are feasible. Only a single sample is required, and the method can be applied to both slowly and rapidly exchanging ligands. The method was applied to three test cases: one complex involving muscle fatty acid-binding protein (mFABP) and two complexes involving the leukocyte function-associated antigen 1 (LFA-1) I-domain. Without using experimental protein NMR assignments, the method identified the known binding poses with good accuracy. The addition of experimental protein NMR assignments improves the results. Our "NOE matching" approach is expected to be widely applicable; i.e., it does not appear to depend on a fortuitous distribution of binding pocket residues.

Discovery of ligands for Nurr1 by combined use of NMR screening with different isotopic and spin-labeling strategies

A comprehensive approach to target screening, hit validation, and binding site determination by nuclear magnetic resonance (NMR) spectroscopy is presented. NMR (19)F signal perturbation was used to screen a small compound library and identify candidate ligands to the target of interest. Ligand dissociation constants were measured using a pegylated form of the protein, which resulted in a 2-fold increase in the strength of the saturation transfer difference signal. The initial small-molecule hits were further optimized by combining a residue-specific labeling strategy, to identify the specific sites of interaction with the protein, with a second site screening approach based on relaxation enhancement using a paramagnetic probe. The advantages of this combination strategy in the identification and optimization of weak binding chemical entities early in a program are illustrated with the discovery of a low micromolar ligand (K(d) = 20 microM) for Nurr1 and identification of the binding site location through residue-specific (15)N isotope labeling and derivatization of Cys residues with 2-mercaptoethanol-1-(13)C.

A Profile of the Residues in the First Intracellular Loop Critical for Gs-Mediated Signaling of Human Prostacyclin Receptor Characterized by an Integrative Approach of NMR-Experiment and Mutagenesis

The first intracellular loop (iLP1, residues 39-51) of human prostacyclin receptor (IP) was proposed to be involved in signaling via its interaction with the Galphas protein. First, evidence of the IP iLP1 interaction with the C-terminus of the Galphas protein was observed by the fluorescence and NMR spectroscopy using the synthetic peptide (Galphas-Ct) mimicking the C-terminal 11 residues of the Galphas protein in the presence of a constrained synthetic peptide mimicking the IP iLP1. Then, the residues (Arg42, Ala44, and Arg45) in the IP iLP1 peptide possibly involved in contacting the Galphas-Ct peptide were initially assigned by observation of the significant proton resonance shifts of the side chains of the constrained IP iLP1 peptide using 2D (1)H NMR spectroscopy. The results of the NMR studies were used as a guide for further identification of the residues in the IP important to the receptor signaling using a recombinant protein approach. A profile of the residues in the IP iLP1, including the residues observed from the NMR studies involved in the Galphas mediated signaling, was mapped out by mutagenesis. According to our results, it can be predicted that the seven residues (Arg42-Ala48) with the conserved Arg45 at the center will form an epitope with a specific conformation involved in the Galphas mediated signaling. The conservation of the basic residues (Arg45 in the IP) in all of the prostanoid receptors suggests that the iLP1 regions of the other prostanoid receptors may also contain the epitopes important to their signaling.

Utilization of NMR-derived fragment leads in drug design

The advent of large-scale NMR-based screening has enabled new strategies for the design of novel, potent inhibitors of therapeutic targets. In particular, fragment-based strategies, in which molecular portions of the final high-affinity ligand are experimentally identified prior to chemical synthesis, have found widespread utility. This chapter will discuss some of the practical considerations for identifying and utilizing these fragment leads in drug design, with special emphasis on some of the lessons learned from more than a decade of industry experience.

Validation of the Binding Site Structure of the Cellular Retinol-Binding Protein (CRBP) by Ligand NMR Chemical Shift Perturbations

We have calculated proton chemical shift perturbations (CSPs) of retinol in the cellular retinol-binding protein (CRBP) through the use of a recently developed computational approach (Wang et al. J. Chem. Phys. 2004, 120, 11392-11400). Excellent agreement with experimental values was obtained for the X-ray structure, whereas the lack of a key hydrogen bond and the distorted isoprene tail of retinol for some NMR models lead to large CSP RMSDs. Therefore, a comparison of computed CSPs of retinol with experiment offers a convenient way to validate the structure of retinol and its orientation in the binding site for the NMR structures.

Pose scoring by NMR

Recently, we have developed a fast approach to calculate NMR chemical shifts using the divide and conquer method at the semiempirical level. To demonstrate the utility of this approach for characterizing protein-ligand interactions, we used the deviation of calculated chemical shift perturbations from experiment to determine the orientation of a ligand (GPI-1046) in the binding pocket of the FK506 binding protein (FKBP12). Moreover, we were able to select the native state of the ligand from a collection of decoy poses. A key hydrogen bond between O1 and HN in Ile56 was also identified. Our results suggest that ligand-induced chemical shift perturbations can be used to refine protein/ligand structures.

Non-Peptidic Small-Molecule Inhibitors of the Single-Chain Hepatitis C Virus NS3 Protease/NS4A Cofactor Complex Discovered by Structure-Based NMR Screening

NMR-based screening of a customized fragment library identified 16 small-molecule hits that bind weakly (K(D) approximately 100 microM to 10 mM) to substrate binding sites of the NS4A-bound NS3 protease of the hepatitis C virus (HCV). Analogues for five classes of NMR hits were evaluated by a combination of NMR and biochemical data yielding SAR and, in most cases, optimized hits with improved potencies (K(D) approximately K(I) approximately 40 microM to 1 mM). NMR chemical shift perturbation data were used to establish the binding location and orientation of the active site directed scaffolds in these five analogue series. Two of these scaffolds, which bind the enzyme at the proximal S1-S3 and S2' substrate binding sites, were linked together producing competitive inhibitors of the HCV NS3 protease with potencies in the micromolar range. This example illustrates that the low molecular weight scaffolds discovered from structure-based NMR screening can be optimized with focused structure-guided chemistry to produce potent nonpeptidic small-molecule inhibitors of the HCV NS3 protease.

NMR Studies of the Interaction of Tryparedoxin with Redox-Inactive Substrate Homologues

Tryparedoxins (TXNs) are trypanothione-dependent peroxiredoxin oxidoreductases involved in hydroperoxide detoxification that have been shown to determine virulence in trypanosomatids. The structure of (15)N,(13)C-doubly-labeled, C-terminally-His-tagged tryparedoxin 1 from Crithidia fasciculata (Cf TXN1) was elucidated by three-dimensional NMR spectroscopy. Global folding was found to be similar to the crystal structure, but regions near the active site, especially the onset of helix alpha1 with the redox-active Cys 43 and helix alpha2 relevant to substrate binding, were less well defined in solution. The redox-inactive inhibitory substrate analogue N(1),N(8)-bis(ophthalmyl)spermidine was used to study the substrate/TXN interaction by two-dimensional (1)H,(15)N NMR spectroscopy. The NMR data complemented by molecular modeling revealed several alternative modes of ligand binding. The results confirm and extend the concept of TXN action and specificity derived from X-ray analysis and site-directed mutagenesis and thus improve the rational basis for inhibitor design.

Consequences of Single-Site Mutations in the Intestinal Fatty Acid Binding Protein

The intestinal fatty acid binding protein (IFABP) is a small (15 kDa) protein consisting mostly of 10 antiparallel beta-strands (A-J) and a small helical region that serves as a portal for the ligand. Two beta-sheet structures (strands A-E and F-J) surround a cavity into which the ligand binds. In this work, we investigated how changes in the side chains of specific residues are propagated through the structure. To determine what these changes were and how they relate to changes in stability, (15)N chemical shift perturbations were measured and compared to those of the wild-type protein. Seven mutations, five of which change either valine or leucine to glycine, have been examined. All these mutants were less stable than wild-type IFABP, suggesting some structural changes. For five of the mutants, the data suggest that destabilization of a small region of the protein propagates throughout the structure, resulting in an overall decrease in stability. In two (Leu38Gly and Leu89Gly), the loss of cooperativity in the equilibrium denaturation curves suggests that the destabilization of one region may not be transmitted to other regions in a cooperative manner. It is shown that the effect of mutating hydrophobic residues is much greater than that observed upon mutation of a solvent-exposed polar residue.

Spatial localization of ligand binding sites from electron current density surfaces calculated from NMR chemical shift perturbations

Rapid, accurate structure determination of protein-ligand complexes is an essential component in structure-based drug design. We have developed a method that uses NMR protein chemical shift perturbations to spatially localize a ligand when it is complexed with a protein. Chemical shift perturbations on the protein arise primarily from the close proximity of electron current density from the ligand. In our approach the location of the center of the electron current density for a ligand aromatic ring was approximated by a point-dipole, and dot densities were used to represent ligand positions that are allowed by the experimental data. The dot density is increased in the region of space that is consistent for the most data. A surface can be formed in regions of the highest dot density that correlates to the center of the ligand aromatic ring. These surfaces allow for the rapid evaluation of ligand binding, which is demonstrated on a model system and on real data from HCV NS3 protease and HCV NS3 helicase, where the location of ligand binding can be compared to that obtained from difference electron density from X-ray crystallography.

Structures of protein-protein complexes are docked using only NMR restraints from residual dipolar coupling and chemical shift perturbations

NMR structures of protein-protein and protein-ligand complexes rely heavily on intermolecular NOEs. Recent work has shown that if no significant conformational changes occur upon complex formation residual dipolar coupling can replace most of the NOE restraints in protein-protein complexes, while restraints derived from chemical shift perturbations can largely replace intermolecular NOEs in protein-ligand structures. By combining restraints from chemical shift perturbations with orientation restraints derived from measurements of residual dipolar couplings, we show that the structure of the EIN-HPr complex can be calculated without NOE restraints. The final structure, built from the crystal structures of EIN and HPr in their uncomplexed form and docked only with NMR restraints, places HPr within 2.5 A of the position determined from the mean NMR structure of the complex.

Mapping protein-protein interactions in solution by NMR spectroscopy

NMR is very well suited to the study of especially weak protein-protein interactions, as no crystallization is required. The available NMR methods to this end are reviewed and illustrated with applications from the recent biochemical literature: intermolecular NOEs, cross-saturation, chemical shift perturbation, dynamics and exchange perturbation, paramagnetic methods, and dipolar orientation. Most of these methods are now routinely applied for complexes with total molecular mass of 60 kDa and can likely be applied to systems up to 1000 kDa. A substantial fraction of complexes studied show distinct effects of induced fit affecting structural and dynamical properties beyond the contact interface.