Mapping the ligand binding site at protein side-chains in protein-ligand complexes through NOE difference spectroscopy

This report describes a novel NMR approach for mapping the interaction surface between an unlabeled ligand and a ¹³C,¹⁵N-labeled protein. The method relies on the spin inversion properties of the dipolar relaxation pathways and records the differential relaxation of two spin modes, where ligand and protein ¹H magnetizations are aligned either in a parallel or anti-parallel manner. Selective inversion of protein protons is achieved in a straightforward manner by exploiting the one-bond heteronuclear scalar couplings (¹J_CH, ¹J_NH). Suppression of indirect relaxation pathways mediated by bulk water or rapidly exchanging protons is achieved by selective inversion of the water signal in the middle of the NOESY mixing period. The method does not require deuteration of the protein or well separated spectral regions for the protein and the ligand, respectively. Additionally, in contrast to previous methods, the new experiment identifies side-chain enzyme ligand interactions along the intermolecular binding interface. The method is demonstrated with an application to the B₁₂-binding subunit of glutamate mutase from Clostridium tetanomorphum for which NMR chemical shift changes upon B₁₂-nucleotide loop binding and a high-resolution solution structure are available.     

NMR resonance assignment of selectively labeled proteins by the use of paramagnetic ligands

Selective isotopic labeling of larger proteins greatly simplifies protein NMR spectra and reduces signal overlap, but selectively labeled proteins cannot be easily assigned since the sequential assignment method is not applicable. Here we describe a strategy for resonance assignment in selectively labeled proteins. Our approach involves a spin-labeled analog of a ligand of which the three-dimensional structure in complex with the target protein is known. Other methods for introduction of the spin label are possible. The paramagnetic center causes faster relaxation of all neighboring nuclei in a distance-dependent manner. Measurement of this effect allows to deduce distances between isotopically labeled residues and the paramagnetic center which can be used for resonance assignment. The method is demonstrated for the catalytic domain of Abl kinase in complex with the inhibitor, STI571.     

An all atom energy based computational protocol for predicting binding affinities of protein-ligand complexes

We report here a computationally fast protocol for predicting binding affinities of non-metallo protein-ligand complexes. The protocol builds in an all atom energy based empirical scoring function comprising electrostatics, van der Waals, hydrophobicity and loss of conformational entropy of protein side chains upon ligand binding. The method is designed to ensure transferability across diverse systems and has been validated on a heterogenous dataset of 161 complexes consisting of 55 unique protein targets. The scoring function trained on a dataset of 61 complexes yielded a correlation of r=0.92 for the predicted binding free energies against the experimental binding affinities. Model validation and parameter analysis studies ensure the predictive ability of the scoring function. When tested on the remaining 100 protein-ligand complexes a correlation of r=0.92 was recovered. The high correlation obtained underscores the potential applicability of the methodology in drug design endeavors. The scoring function has been web enabled at as binding affinity prediction of protein-ligand (BAPPL) server.     

Analysing six types of protein-protein interfaces

Non-covalent residue side-chain interactions occur in many different types of proteins and facilitate many biological functions. Are these differences manifested in the sequence compositions and/or the residue-residue contact preferences of the interfaces? Previous studies analysed small data sets and gave contradictory answers. Here, we introduced a new data-mining method that yielded the largest high-resolution data set of interactions analysed. We introduced an information theory-based analysis method. On the basis of sequence features, we were able to differentiate six types of protein interfaces, each corresponding to a different functional or structural association between residues. Particularly, we found significant differences in amino acid composition and residue-residue preferences between interactions of residues within the same structural domain and between different domains, between permanent and transient interfaces, and between interactions associating homo-oligomers and hetero-oligomers. The differences between the six types were so substantial that, using amino acid composition alone, we could predict statistically to which of the six types of interfaces a pool of 1000 residues belongs at 63-100% accuracy. All interfaces differed significantly from the background of all residues in SWISS-PROT, from the group of surface residues, and from internal residues that were not involved in non-trivial interactions. Overall, our results suggest that the interface type could be predicted from sequence and that interface-type specific mean-field potentials may be adequate for certain applications.     

Development, validation, and application of adapted PEOE charges to estimate pKa values of functional groups in protein-ligand complexes

For routine pK(a) calculations of protein-ligand complexes in drug design, the PEOE method to compute partial charges was modified. The new method is applicable to a large scope of proteins and ligands. The adapted charges were parameterized using experimental free energies of solvation of amino acids and small organic ligands. For a data set of 80 small organic molecules, a correlation coefficient of r(2) = 0.78 between calculated and experimental solvation free energies was obtained. Continuum electrostatics pK(a) calculations based on the Poisson-Boltzmann equation were carried out on a validation set of nine proteins for which 132 experimental pK(a) values are known. In total, an overall RMSD of 0.88 log units between calculated and experimentally determined data is achieved. In particular, the predictions of significantly shifted pK(a) values are satisfactory, and reasonable estimates of protonation states in the active sites of lysozyme and xylanase could be obtained. Application of the charge-assignment and pK(a)-calculation procedure to protein-ligand complexes provides clear structural interpretations of experimentally observed changes of protonation states of functional groups upon complex formation. This information is essential for the interpretation of thermodynamic data of protein-ligand complex formation and provides the basis for the reliable factorization of the free energy of binding in enthalpic and entropic contributions. The modified charge-assignment procedure forms the basis for future automated pK(a) calculations of protein-ligand complexes.     

A sensitive and robust method for obtaining intermolecular NOEs between side chains in large protein complexes

A method for measuring intermolecular NOEs in protein complexes based on asymmetric sample deuteration is described. ¹³C/¹H-I,L,V-methyl, U-²H labeled protein is produced using the biosynthetic precursors [-¹³C]--ketobutyrate and [,-¹³C₂]--ketoisovalerate. The labeled protein is mixed with its unlabeled binding partner and a 3D ¹³C-HMQC-NOESY is recorded, yielding unambiguous intermolecular aromatic/methyl NOEs. A simple synthesis of the biosynthetic precursors via reaction of diethyl oxalate with alkyl Grignard compounds is reported. The method is demonstrated for a 35 kDa heterodimeric protein complex dissolved in a CHAPS micelle. This approach will facilitate the solution structure determination of protein/protein, protein/ligand or protein/nucleic acid complexes.These authors contributed equallyThese authors contributed equally     

Characterizing hydration sites in protein-ligand complexes towards the design of novel ligands

Water is an essential part of protein binding sites and mediates interactions to ligands. Its displacement by ligand parts affects the free binding energy of resulting protein-ligand complexes. Therefore the characterization of solvation properties is important for design. Of particular interest is the propensity of localized water to be favorably displaced by a ligand. This review discusses two popular computational approaches addressing these questions, namely WaterMap based on statistical mechanics analysis of MD simulations and 3D RISM based on integral equation theory of liquids. The theoretical background and recent applications in structure-based design will be presented.     

Similarities between intra- and intermolecular hydrogen bonds in RNA kissing complexes found by means of cross-correlated relaxation

The bond lengths and dynamics of intra- and intermolecular hydrogen bonds in an RNA kissing complex have been characterized by determining the NMR relaxation rates of various double- and triple-quantum coherences that involve an imino proton and two neighboring nitrogen-15 nuclei belonging to opposite bases. New experiments allow one to determine the chemical shift anisotropy of the imino protons. The bond lengths derived from dipolar relaxation and the lack of modulations of the nitrogen chemical shifts indicate that the intermolecular hydrogen bonds which hold the kissing complex together are very similar to the intramolecular hydrogen bonds in the double-stranded stem of the RNA.     

Syntheses, characterisation, infrared and Mo NMR spectroscopy of some coordinated oxo—molybdenum(VI) complexes

Adducts with MoO₄²⁻ tetrahedra coordinated to Cr(III) or Co(III) complexes have been synthesized and studied by IR and high resolution ⁹⁵Mo NMR spectroscopy. The ⁹⁵Mo chemical shifts of the adducts with cobalt(III) lie in the range −33.2 to + 49.4 ppm. This may be compared with an overall known chemical shift range in excess of 7000 ppm and implies a similarity in the molybdenum environment in all cases. For adducts with chelated cobalt(III) complexes several rather broad ⁹⁵Mo singnals are obtained with linewidths up to 260 Hz.     

Solvation study of the non-specific lipid transfer protein from wheat by intermolecular NOEs with water and small organic molecules

Intermolecular nuclear Overhauser effects (NOEs) were measured between the protons of various small solvent or gas molecules and the non-specific lipid transfer protein (ns-LTP) from wheat. Intermolecular NOEs were observed with the hydrophobic pocket in the interior of wheat ns-LTP, which grew in intensity in the order cyclopropane (saturated solution) < methane (140 bar) < ethane (40 bar) < acetonitrile (5% in water) < cyclohexane (saturated solution) < benzene (saturated solution). No intermolecular NOEs were observed with dioxane (5% in water). The intermolecular NOEs were negative for all of the organic molecules tested. Intermolecular NOEs between wheat ns-LTP and water were weak or could not be distinguished from exchange-relayed NOEs. As illustrated by the NOEs with cyclohexane versus dioxane, the hydrophobic pocket in wheat ns-LTP preferably binds non-polar molecules. Yet, polar molecules like acetonitrile can also be accommodated. The pressure dependence of the NOEs between methane and wheat ns-LTP indicated incomplete occupancy, even at 190 bar methane pressure. In general, NOE intensities increased with the size of the ligand molecule and its vapor pressure. NMR of the vapor phase showed excellent resolution between the signals from the gas phase and those from the liquid phase. The vapor concentration of cyclohexane was fivefold higher than that of the dioxane solution, supporting the binding of cyclohexane versus uptake of dioxane.     

NMR structure of a complex between MDM2 and a small molecule inhibitor

MDM2 is a regulator of cell growth processes that acts by binding to the tumor suppressor protein p53 and ultimately restraining its activity. While inactivation of p53 by mutation is commonly observed in human cancers, a substantial percentage of tumors express wild type p53. In many of these cases, MDM2 is overexpressed, and it is believed that suppression of MDM2 activity could yield therapeutic benefits. Therefore, we have been focusing on the p53-MDM2 interaction as the basis of a drug discovery program and have been able to develop a series of small molecule inhibitors. We herein report a high resolution NMR structure of a complex between the p53-binding domain of MDM2 and one of these inhibitors. The form of MDM2 utilized was an engineered hybrid between the human and Xenopus sequences, which provided a favorable combination of relevancy and stability. The inhibitor is found to bind in the same site as does a highly potent peptide fragment of p53. The inhibitor is able to successfully mimic the peptide by duplicating interactions in three subpockets normally made by amino acid sidechains, and by utilizing a scaffold that presents substituents with rigidity and spatial orientation comparable to that provided by the alpha helical backbone of the peptide. The structure also suggests opportunities for modifying the inhibitor to increase its potency.     

Formation of cyclo- and polystannanes by dehydrogenative stannane coupling catalyzed by platinum(II) complexes

Reaction of (PhMe₂P)₂PtMe₂ or [(κ²-P,N)-Ph₂PC₂H₄NMe₂]PtMe₂ with an excess of H₂SnBu₂ or H₂SnPh₂ resulted in the catalytic formation of cyclo-, oligo- and/or polystannanes. In the reaction of (PhMe₂P)₂PtMe₂ with H₂SnBu₂, linear oligomeric species H(SnBu₂)_nH were observed in the initial stage of the reaction, which eventually converted into cyclostannanes. Only polystannanes were observed in the reaction of [(κ²-P,N)-Ph₂PC₂H₄NMe₂]PtMe₂ with H₂SnBu₂. The reactions of H₂SnPh₂ were similar, but more difficult to analyze due to redistribution reactions and the formation of insoluble products. The mechanism of the reactions is clearly different to that previously observed for HSnR₃ because metal complexes indicative of oxidative addition/reductive elimination reactions were only observed as minor products.     

Syntheses and structural characterization of mononuclear Rh-Cp* and Ir-Cp* complexes with η-phenanthrene, η-pyrene and η-triphenylene

Four novel mononuclear Rh-Cp* and Ir-Cp* complexes with polycyclic aromatic hydrocarbons (PAHs), [M(Cp*)(η⁶-PAHs)](BF₄)₂ (M = Rh and Ir; Cp* = η⁵-C₅Me₅; PAHs = phenanthrene (phn), pyrene (pyr) and triphenylene (triph)), were prepared by the reactions of the intermediate [M(Cp*)(Me₂CO)₃]²⁺ with appreciable PAHs. Their structures were characterized by a single crystal X-ray analysis, ¹H, ¹³C {¹H} NMR and 2D NMR techniques. The X-ray crystallographic studies showed that the [M(Cp*)]²⁺ fragment is η⁶-coordinated to one terminal benzene ring in each PAH. In particular, it is interesting to note that the partial π/π/π/π interaction was formed in the Ir-pyr complex [Ir(Cp*)(η⁶-pyr)](BF₄)₂. The 1D and 2D NMR studies described that the Rh-Cp* and Ir-Cp* complexes with PAHs gave unique ¹H and ¹³C {¹H} NMR spectra with positive coordination shifts (Δδ(¹H, ¹³C)) in (CD₃)₂CO at 23 °C, which are likely induced by the local effect and the non-local effect on the coordination of the [M(Cp*)]²⁺ fragment to PAHs. The decreasing of the coupling constants (³J_H-H) in the η⁶-coordinated benzene ring is also induced, with no changes in the uncoordinated benzene rings. The time-course of ¹H NMR spectra showed that Rh-Cp* and Ir-Cp* complexes with PAHs are partially dissociated to [M(Cp*)(Me₂CO)₃]²⁺ and metal-free PAHs in (CD₃)₂CO at 23 °C. It was demonstrated that their stabilities are in the order of Ir-triph, Ir-phn, Ir-pyr and Rh-triph complexes in (CD₃)₂CO.     

Measurement of amide proton exchange rates and NOEs with water in 13C/15N-enriched calcineurin B

Summary A rapid and sensitive 2D approach is presented for measuring amide proton exchange rates and the NOE interaction between amide protons and water. The approach is applicable to uniformly ¹³C/¹⁵N-enriched proteins and can measure magnetization exchange rates in the 0.02 to >20s^–1 range. The experiments rely on selective excitation of the water resonance, coupled with purging of underlying H resonances, followed by NOESY-or ROESY-type transfer to amide protons, which are dispersed by the amide ¹⁵N frequencies in an HSQC-type experiment. Two separate but interleaved experiments, with and without selective inversion of the H₂O resonance, yield quantitative results. The method is demonstrated for a sample of the calcium-binding protein calcineurin B. Results indicate rapid amide exchange for the five calcineurin B residues that are analogous to the five rapidly exchanging residues in the central helix of the homologous protein calmodulin.     

Hydrogen-deuterium exchange strategy for delineation of contact sites in protein complexes

We use NMR spectra to determine protein-protein contact sites by observing differences in amide proton hydrogen-deuterium exchange in the complex compared to the free protein in solution. Aprotic organic solvents are used to preserve H/D labeling patterns that would be scrambled in water solutions. The binding site between the mammalian co-chaperone Aha1 with the middle domain of the chaperone Hsp90 obtained by our H/D exchange method corresponds well with that in the X-ray crystal structure of the homologous complex from yeast, even to the observation of a secondary binding site. This method can potentially provide data for complexes with unknown structure and for large or dynamic complexes inaccessible via NMR and X-ray methods.     

Anatomy of protein pockets and cavities: measurement of binding site geometry and implications for ligand design.

Identification and size characterization of surface pockets and occluded cavities are initial steps in protein structure-based ligand design. A new program, CAST, for automatically locating and measuring protein pockets and cavities, is based on precise computational geometry methods, including alpha shape and discrete flow theory. CAST identifies and measures pockets and pocket mouth openings, as well as cavities. The program specifies the atoms lining pockets, pocket openings, and buried cavities; the volume and area of pockets and cavities; and the area and circumference of mouth openings. CAST analysis of over 100 proteins has been carried out; proteins examined include a set of 51 monomeric enzyme-ligand structures, several elastase-inhibitor complexes, the FK506 binding protein, 30 HIV-1 protease-inhibitor complexes, and a number of small and large protein inhibitors. Medium-sized globular proteins typically have 10-20 pockets/cavities. Most often, binding sites are pockets with 1-2 mouth openings; much less frequently they are cavities. Ligand binding pockets vary widely in size, most within the range 10(2)-10(3)A3. Statistical analysis reveals that the number of pockets and cavities is correlated with protein size, but there is no correlation between the size of the protein and the size of binding sites. Most frequently, the largest pocket/cavity is the active site, but there are a number of instructive exceptions. Ligand volume and binding site volume are somewhat correlated when binding site volume is < or =700 A3, but the ligand seldom occupies the entire site. Auxiliary pockets near the active site have been suggested as additional binding surface for designed ligands (Mattos C et al., 1994, Nat Struct Biol 1:55-58). Analysis of elastase-inhibitor complexes suggests that CAST can identify ancillary pockets suitable for recruitment in ligand design strategies. Analysis of the FK506 binding protein, and of compounds developed in SAR by NMR (Shuker SB et al., 1996, Science 274:1531-1534), indicates that CAST pocket computation may provide a priori identification of target proteins for linked-fragment design. CAST analysis of 30 HIV-1 protease-inhibitor complexes shows that the flexible active site pocket can vary over a range of 853-1,566 A3, and that there are two pockets near or adjoining the active site that may be recruited for ligand design.     

Stereochemical structure determination of p-cymene Ru(II) complexes containing the PPh2Py ligand with 2-D NOESY and HMQC NMR experiments

Complexes [Ru(p-cymene)Cl₂(PPh₂Py)] (1) and [Ru(p-cymene)Cl(PPh₂Py)]BF₄ (2) were studied by means of ¹H, ¹³C{¹H} 2-D NOESY and HMQC NMR spectral methods. NMR data agree with C₁ and C_s symmetries for complexes 1 and 2, respectively. NOESY cross-peaks allowed the assignment of signals to CH arene protons and in the case of complex 2 the determination of the molecular stereochemistry. These results are in agreement with the X-ray molecular structures of both complexes.     

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.