HADDOCK: a protein-protein docking approach based on biochemical or biophysical information

The structure determination of protein-protein complexes is a rather tedious and lengthy process, by both NMR and X-ray crystallography. Several methods based on docking to study protein complexes have also been well developed over the past few years. Most of these approaches are not driven by experimental data but are based on a combination of energetics and shape complementarity. Here, we present an approach called HADDOCK (High Ambiguity Driven protein-protein Docking) that makes use of biochemical and/or biophysical interaction data such as chemical shift perturbation data resulting from NMR titration experiments or mutagenesis data. This information is introduced as Ambiguous Interaction Restraints (AIRs) to drive the docking process. An AIR is defined as an ambiguous distance between all residues shown to be involved in the interaction. The accuracy of our approach is demonstrated with three molecular complexes. For two of these complexes, for which both the complex and the free protein structures have been solved, NMR titration data were available. Mutagenesis data were used in the last example. In all cases, the best structures generated by HADDOCK, that is, the structures with the lowest intermolecular energies, were the closest to the published structure of the respective complexes (within 2.0 A backbone RMSD).     

Characterization of the conformational equilibrium between the two major substates of RNase A using NMR chemical shifts

Following the recognition that NMR chemical shifts can be used for protein structure determination, rapid advances have recently been made in methods for extending this strategy for proteins and protein complexes of increasing size and complexity. A remaining major challenge is to develop approaches to exploit the information contained in the chemical shifts about conformational fluctuations in native states of proteins. In this work we show that it is possible to determine an ensemble of conformations representing the free energy surface of RNase A using chemical shifts as replica-averaged restraints in molecular dynamics simulations. Analysis of this surface indicates that chemical shifts can be used to characterize the conformational equilibrium between the two major substates of this protein.     

CABS‐NMR—De novo tool for rapid global fold determination from chemical shifts,residual dipolar couplings and sparse methyl‐methyl noes

Recent development of nuclear magnetic resonance (NMR) techniques provided new types of structural restraints that can be successfully used in fast and low‐cost global protein fold determination. Here, we present CABS‐NMR, an efficient protein modeling tool, which takes advantage of such structural restraints. The restraints are converted from original NMR data to fit the coarse grained protein representation of the C‐Alpha‐Beta‐Side‐group (CABS) algorithm. CABS is a Monte Carlo search algorithm that uses a knowledge‐based force field. Its versatile structure enables a variety of protein‐modeling protocols, including purely de novo folding, folding guided by restraints derived from template structures or, structure assembly based on experimental data. In particular, CABS‐NMR uses the distance and angular restraints set derived from various NMR experiments. This new modeling technique was successfully tested in structure determination of 10 globular proteins of size up to 216 residues, for which sparse NMR data were available. Additional detailed analysis was performed for a S100A1 protein. Namely, we successfully predicted Nuclear Overhauser Effect signals on the basis of low‐energy structures obtained from chemical shifts by CABS‐NMR. It has been observed that utility of chemical shifts and other types of experimental data (i.e. residual dipolar couplings and methyl‐methyl Nuclear Overhauser Effect signals) in the presented modeling pipeline depends mainly on size of a protein and complexity of its topology. In this work, we have provided tools for either post‐experiment processing of various kinds of NMR data or fast and low‐cost structural analysis in the still challenging field of new fold predictions. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Platinum-modified adenines: unprecedented protonation behavior revealed by NMR spectroscopy and relativistic density-functional theory calculations

Two novel Pt(IV) complexes of aromatic cytokinins with possible antitumor properties were prepared by reaction of selected aminopurines with K(2)PtCl(6). The structures of both complexes, 9-[6-(benzylamino)purine] pentachloroplatinate (IV) and 9-[6-(furfurylamino)purine] pentachloroplatinate (IV), were characterized in detail by using two-dimensional NMR spectroscopy ((1)H, (13)C, (15)N, and (195)Pt) in solution and CP/MAS NMR techniques in the solid state. We report for the first time the X-ray structure of a nucleobase adenine derivative coordinated to Pt(IV) via the N9 atom. The protonation equilibria for the complexes in solution were characterized by using NMR spectroscopy (isotropic chemical shifts and indirect nuclear spin-spin coupling constants) and the structural conclusions drawn from the NMR analysis are supported by relativistic density-functional theory (DFT) calculations. Because of the presence of the Pt atom, hybrid GGA functionals and scalar-relativistic and spin-orbit corrections were employed for both the DFT calculations of the molecular structure and particularly for the NMR chemical shifts. In particular, the populations of the N7-protonated and neutral forms of the complexes in solution were characterized by correlating the experimental and the DFT-calculated NMR chemical shifts. In contrast to the chemical exchange process involving the N7-H group, the hydrogen atom at N3 was determined to be unexpectedly rigid, probably because of the presence of the stabilizing intramolecular interaction N3-H···Cl. The described methodology combining the NMR spectroscopy and relativistic DFT calculations can be employed for characterizing the tautomeric and protonation equilibria in a large family of transition-metal-modified purine bases.     

3D structure determination of the Crh protein from highly ambiguous solid-state NMR restraints

In a wide variety of proteins, insolubility presents a challenge to structural biology, as X-ray crystallography and liquid-state NMR are unsuitable. Indeed, no general approach is available as of today for studying the three-dimensional structures of membrane proteins and protein fibrils. We here demonstrate, at the example of the microcrystalline model protein Crh, how high-resolution 3D structures can be derived from magic-angle spinning solid-state NMR distance restraints for fully labeled protein samples. First, we show that proton-mediated rare-spin correlation spectra, as well as carbon-13 spin diffusion experiments, provide enough short, medium, and long-range structural restraints to obtain high-resolution structures of this 2 x 10.4 kDa dimeric protein. Nevertheless, the large number of 13C/15N spins present in this protein, combined with solid-state NMR line widths of about 0.5-1 ppm, induces substantial ambiguities in resonance assignments, preventing 3D structure determination by using distance restraints uniquely assigned on the basis of their chemical shifts. In the second part, we thus demonstrate that an automated iterative assignment algorithm implemented in a dedicated solid-state NMR version of the program ARIA permits to resolve the majority of ambiguities and to calculate a de novo 3D structure from highly ambiguous solid-state NMR data, using a unique fully labeled protein sample. We present, using distance restraints obtained through the iterative assignment process, as well as dihedral angle restraints predicted from chemical shifts, the 3D structure of the fully labeled Crh dimer refined at a root-mean-square deviation of 1.33 A.     

Chemical shifts in nucleic acids studied by density functional theory calculations and comparison with experiment

NMR chemical shifts are highly sensitive probes of local molecular conformation and environment and form an important source of structural information. In this study, the relationship between the NMR chemical shifts of nucleic acids and the glycosidic torsion angle, χ, has been investigated for the two commonly occurring sugar conformations. We have calculated by means of DFT the chemical shifts of all atoms in the eight DNA and RNA mono-nucleosides as a function of these two variables. From the DFT calculations, structures and potential energy surfaces were determined by using constrained geometry optimizations at the BP86/TZ2P level of theory. The NMR parameters were subsequently calculated by single-point calculations at the SAOP/TZ2P level of theory. Comparison of the (1) H and (13) C?NMR shifts calculated for the mono-nucleosides with the shifts determined by NMR spectroscopy for nucleic acids demonstrates that the theoretical shifts are valuable for the characterization of nucleic acid conformation. For example, a clear distinction can be made between χ angles in the anti and syn domains. Furthermore, a quantitative determination of the χ angle in the syn domain is possible, in particular when (13) C and (1) H chemical shift data are combined. The approximate linear dependence of the C1' shift on the χ angle in the anti domain provides a good estimate of the angle in this region. It is also possible to derive the sugar conformation from the chemical shift information. The DFT calculations reported herein were performed on mono-nucleosides, but examples are also provided to estimate intramolecularly induced shifts as a result of hydrogen bonding, polarization effects, or ring-current effects.     

A novel approach for characterizing protein ligand complexes: molecular basis for specificity of small-molecule Bcl-2 inhibitors

The increasing diversity of small molecule libraries has been an important source for the development of new drugs and, more recently, for unraveling the mechanisms of cellular events-a process termed chemical genetics.(1) Unfortunately, the majority of currently available compounds are mechanism-based enzyme inhibitors, whereas most of cellular activity regulation proceeds on the level of protein-protein interactions. Hence, the development of small molecule inhibitors of protein-protein interactions is important. When screening compound libraries, low-micromolar inhibitors of protein interactions can be routinely found. The enhancement of affinities and rationalization of the binding mechanism require structural information about the protein-ligand complexes. Crystallization of low-affinity complexes is difficult, and their NMR analysis suffers from exchange broadening, which limits the number of obtainable intermolecular constraints. Here we present a novel method of ligand validation and optimization, which is based on the combination of structural and computational approaches. We successfully used this method to analyze the basis for structure-activity relationships of previously selected (2) small molecule inhibitors of the antiapoptotic protein Bcl-xL and identified new members of this inhibitor family.     

Design, synthesis, and evaluation of a lanthanide chelating protein probe: CLaNP-5 yields predictable paramagnetic effects independent of environment

Immobilized lanthanide ions offer the opportunity to refine structures of proteins and the complexes they form by using restraints obtained from paramagnetic NMR experiments. We report the design, synthesis, and spectroscopic evaluation of the lanthanide chelator, Caged Lanthanide NMR Probe 5 (CLaNP-5) readily attachable to a protein surface via two cysteine residues. The probe causes tunable pseudocontact shifts, alignment, paramagnetic relaxation enhancement, and luminescence, by chelating it to the appropriate lanthanide ion. The observation of single shifts and the finding that the magnetic susceptibility tensors obtained from shifts and alignment analyses are highly similar strongly indicate that the probe is rigid with respect to the protein backbone. By placing the probe at various positions on a model protein it is demonstrated that the size and orientation of the magnetic susceptibility tensor of the probe are independent of the local protein environment. Consequently, the effects of the probe are readily predictable using a protein structure only. These findings designate CLaNP-5 as a protein probe to deliver unambiguous high quality structural restraints in studies on protein-protein and protein-ligand interactions.     

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.     

Mass spectrometry-based structural proteomics

Structural proteomics is the application of protein chemistry and modern mass spectrometric techniques to problems such as the characterization of protein structures and assemblies and the detailed determination of protein-protein interactions. The techniques used in structural proteomics include crosslinking, photoaffinity labeling, limited proteolysis, chemical protein modification and hydrogen/deuterium exchange, all followed by mass spectrometric analysis. None of these methods alone can provide complete structural information, but a "combination" of these complementary approaches can be used to provide enough information for answering important biological questions. Structural proteomics can help to determine, for example, the detailed structure of the interfaces between proteins that may be important drug targets and the interactions between proteins and ligands. In this review, we have tried to provide a brief overview of structural proteomics methodologies, illustrated with examples from our laboratory and from the literature.     

Effects of amino acid phi,psi propensities and secondary structure interactions in modulating H alpha chemical shifts in peptide and protein beta-sheet.

H alpha chemical shifts are often used as indicators of secondary structure formation in protein structural analysis and peptide folding studies. On the basis of NMR analysis of model beta-sheet and alpha-helical peptides, together with a statistical analysis of protein structures for which NMR data are available, we show that although the gross pattern of H alpha chemical shifts reflects backbone torsion angles, longer range effects from distant amino acids are the dominant factor determining experimental chemical shifts in beta-sheets of peptides and proteins. These show context-dependent variations that aid structural assignment and highlight anomalous shifts that may be of structural significance and provide insights into beta-sheet stability.     

PI by NMR: Probing CH–π Interactions in Protein–Ligand Complexes by NMR Spectroscopy

While CH–π interactions with target proteins are crucial determinants for the affinity of arguably every drug molecule, no method exists to directly measure the strength of individual CH–π interactions in drug–protein complexes. Herein, we present a fast and reliable methodology called PI (π interactions) by NMR, which can differentiate the strength of protein–ligand CH–π interactions in solution. By combining selective amino‐acid side‐chain labeling with ¹H‐¹³C NMR, we are able to identify specific protein protons of side‐chains engaged in CH–π interactions with aromatic ring systems of a ligand, based solely on ¹H chemical‐shift values of the interacting protein aromatic ring protons. The information encoded in the chemical shifts induced by such interactions serves as a proxy for the strength of each individual CH–π interaction. PI by NMR changes the paradigm by which chemists can optimize the potency of drug candidates: direct determination of individual π interactions rather than averaged measures of all interactions.     

Characterization of metal centers in bioinorganic complexes using ab initio calculations of 113Cd chemical shifts

Determination of (113)Cd chemical shift is of significant interest in NMR characterization of metal porphyrins, metal-histidine interactions, and other metal-ligand interactions in many bioinorganic complexes and metalloproteins. In this study, we present a detailed account of a number of quantum chemical investigations aimed at relating isotropic and anisotropic (113)Cd chemical shifts to the structure of several biologically relevant complexes with discrete and polymeric structures. Calculated and experimentally determined chemical shift values are compared to correlate the variation of the chemical shift values with the structural changes around the metal center. Our results infer that the density functional theory using the Sadlej basis set on the cadmium atom is a suitable method for estimating cadmium shielding values to a reasonable accuracy.     

SOS-NMR: a saturation transfer NMR-based method for determining the structures of protein-ligand complexes

An NMR-based alternative to traditional X-ray crystallography and NMR methods for structure-based drug design is described that enables the structure determination of ligands complexed to virtually any biomolecular target regardless of size, composition, or oligomeric state. The method utilizes saturation transfer difference (STD) NMR spectroscopy performed on a ligand complexed to a series of target samples that have been deuterated everywhere except for specific amino acid types. In this way, the amino acid composition of the ligand-binding site can be defined, and, given the three-dimensional structure of the protein target, the three-dimensional structure of the protein-ligand complex can be determined. Unlike earlier NMR methods for solving the structures of protein-ligand complexes, no protein resonance assignments are necessary. Thus, the approach has broad potential applications--especially in cases where X-ray crystallography and traditional NMR methods have failed to produce structural data. The method is called SOS-NMR for structural information using Overhauser effects and selective labeling and is validated on two protein-ligand complexes: FKBP complexed to 2-(3'-pyridyl)-benzimidazole and MurA complexed to uridine diphosphate N-acetylglucosamine.     

Structural studies of a protein using the assigned back-bone carbonyl carbon-13 NMR resonances

Several new methods to observe and assign the back-bone peptide carbonyl carbon NMR resonances of a protein have been developed and applied for Streptomyces subtilisin inhibitor. With the assigned carbonyl resonances various site-specific structural information can be obtained for virtually any site of the peptide chain. By putting those site-specific NMR parameters, including chemical shifts, rotational structures around peptide bond, and solvent accessibilities, together a new algorithm to predict the protein structures in solution has emerged.     

Characterization of noninnocent metal complexes using solid-state NMR spectroscopy: o-dioxolene vanadium complexes

(51)V solid-state NMR (SSNMR) studies of a series of noninnocent vanadium(V) catechol complexes have been conducted to evaluate the possibility that (51)V NMR observables, quadrupolar and chemical shift anisotropies, and electronic structures of such compounds can be used to characterize these compounds. The vanadium(V) catechol complexes described in these studies have relatively small quadrupolar coupling constants, which cover a surprisingly small range from 3.4 to 4.2 MHz. On the other hand, isotropic (51)V NMR chemical shifts cover a wide range from -200 to 400 ppm in solution and from -219 to 530 ppm in the solid state. A linear correlation of (51)V NMR isotropic solution and solid-state chemical shifts of complexes containing noninnocent ligands is observed. These experimental results provide the information needed for the application of (51)V SSNMR spectroscopy in characterizing the electronic properties of a wide variety of vanadium-containing systems and, in particular, those containing noninnocent ligands and that have chemical shifts outside the populated range of -300 to -700 ppm. The studies presented in this report demonstrate that the small quadrupolar couplings covering a narrow range of values reflect the symmetric electronic charge distribution, which is also similar across these complexes. These quadrupolar interaction parameters alone are not sufficient to capture the rich electronic structure of these complexes. In contrast, the chemical shift anisotropy tensor elements accessible from (51)V SSNMR experiments are a highly sensitive probe of subtle differences in electronic distribution and orbital occupancy in these compounds. Quantum chemical (density functional theory) calculations of NMR parameters for [VO(hshed)(Cat)] yield a (51)V chemical shift anisotropy tensor in reasonable agreement with the experimental results, but surprisingly the calculated quadrupolar coupling constant is significantly greater than the experimental value. The studies demonstrate that substitution of the catechol ligand with electron-donating groups results in an increase in the HOMO-LUMO gap and can be directly followed by an upfield shift for the vanadium catechol complex. In contrast, substitution of the catechol ligand with electron-withdrawing groups results in a decrease in the HOMO-LUMO gap and can directly be followed by a downfield shift for the complex. The vanadium catechol complexes were used in this work because (51)V is a half-integer quadrupolar nucleus whose NMR observables are highly sensitive to the local environment. However, the results are general and could be extended to other redox-active complexes that exhibit coordination chemistry similar to that of the vanadium catechol complexes.     

Predicting 9Be nuclear magnetic resonance chemical shielding tensors utilizing density functional theory

The structures of a series of beryllium containing complexes have been optimized at the B3LYP/6-31G(d) level and their (9)Be magnetic shielding values have been determined using B3LYP/6-311G+g(2d,p) and the gauge-including atomic orbital (GIAO) method. The calculated chemical shifts are in excellent agreement with experimental values. The performance of a variety of NMR methods (SGO, IGAIM, CSGT) were also examined but were found to be inferior to the GIAO method at the chosen level of theory employed. The theoretical method has been utilized to predict the beryllium chemical shifts of structurally characterized complexes for which no measured (9)Be NMR spectrum exists, and to investigate a literature complex with an unusual (9)Be NMR chemical shift. A new standard for beryllium NMR in nonaqueous solvents has been suggested.     

Structure of a Protein–RNA Complex by Solid‐State NMR Spectroscopy

Solid‐state NMR (ssNMR) is applicable to high molecular‐weight (MW) protein assemblies in a non‐amorphous precipitate. The technique yields atomic resolution structural information on both soluble and insoluble particles without limitations of MW or requirement of crystals. Herein, we propose and demonstrate an approach that yields the structure of protein–RNA complexes (RNP) solely from ssNMR data. Instead of using low‐sensitivity magnetization transfer steps between heteronuclei of the protein and the RNA, we measure paramagnetic relaxation enhancement effects elicited on the RNA by a paramagnetic tag coupled to the protein. We demonstrate that this data, together with chemical‐shift‐perturbation data, yields an accurate structure of an RNP complex, starting from the bound structures of its components. The possibility of characterizing protein–RNA interactions by ssNMR may enable applications to large RNP complexes, whose structures are not accessible by other methods.     

PI by NMR: Probing CH–π Interactions in Protein–Ligand Complexes by NMR Spectroscopy

While CH–π interactions with target proteins are crucial determinants for the affinity of arguably every drug molecule, no method exists to directly measure the strength of individual CH–π interactions in drug–protein complexes. Herein, we present a fast and reliable methodology called PI (π interactions) by NMR, which can differentiate the strength of protein–ligand CH–π interactions in solution. By combining selective amino-acid side-chain labeling with ¹H-¹³C NMR, we are able to identify specific protein protons of side-chains engaged in CH–π interactions with aromatic ring systems of a ligand, based solely on ¹H chemical-shift values of the interacting protein aromatic ring protons. The information encoded in the chemical shifts induced by such interactions serves as a proxy for the strength of each individual CH–π interaction. PI by NMR changes the paradigm by which chemists can optimize the potency of drug candidates: direct determination of individual π interactions rather than averaged measures of all interactions.     

Heteronuclear NMR Studies of Cobalt Corrinoids. 18. Correlation of Structure and Magnetic Resonance Parameters in Base-On Cobalamins(1)

Recent X-ray crystal structure determinations (including a new X-ray determination of the structure of cyano-13-epicobalamin reported herein) create a series of seven base-on cobalamins structurally characterized by modern crystallographic techniques in which the intramolecular equilibrium constant for coordination of the axial benzimidazole ligand (Bzm) varies from 76.6 to 4.90 x 10(7). For the five normal, unepimerized cobalamins, the free energy change for this equilibrium correlates linearly with the axial Co-N bond length (r(2) = 0.99). Absolute assignment of the (1)H and (13)C NMR spectra of two of these structurally characterized cobalamins (CH(3)Cbl and CN-13-epiCbl) together with literature assignments for the other complexes now provides reliable (13)C NMR assignments and chemical shifts for all seven complexes. The magnetic anisotropies of the central cobalt atom of all seven complexes, estimated by a method described earlier, are well correlated with the axial Co-N bond distance (r(2) = 0.97) and the free energy of coordination of the Bzm ligand (r(2) = 0.95). The (31)P NMR chemical shift of the phosphodiester moiety of the nucleotide loop is excellently correlated to the axial Co-N bond length (r(2) = 0.996) of the unepimerized cobalamins and provides a reliable method of estimating this bond length. The (15)N chemical shifts of the axially coordinated Bzm nitrogen vary strongly with the axial Co-N bond distance and correlate linearly with this structural parameter (r(2) = 0.991) except for the case of H(2)OCbl(+), which deviates substantially. However, there is a good linear correlation (r(2) = 0.98) of this (15)N chemical shift with the free energy of Bzm coordination for the five unepimerized cobalamins. Attempts to correlate (13)C NMR chemical shifts with structural, thermodynamic, and corrin ring conformational parameters are discussed.