Filtering and selection of structural models: combining docking and NMR

It is generally accepted that protein structures are more conserved than protein sequences, and 3D structure determination by computer simulations have become an important necessity in the postgenomic area. Despite major successes no robust, fast, and automated ab initio prediction algorithms for deriving accurate folds of single polypeptide chains or structures of intermolecular complexes exist at present. Here we present a methodology that uses selection and filtering of structural models generated by docking of known substructures such as individual proteins or domains through easily obtainable experimental NMR constraints. In particular, residual dipolar couplings and chemical shift mapping are used. Heuristic inclusion of chemical or biochemical knowledge about point-to-point interactions is combined in our selection strategy with the NMR data and commonly used contact potentials. We demonstrate the approach for the determination of protein-protein complexes using the EIN/HPr complex as an example and for establishing the domain-domain orientation in a chimeric protein, the recently determined hybrid human-Escherichia. coli thioredoxin.     

NMR approaches for monitoring domain orientations in calcium-binding proteins in solution using partial replacement of Ca by Tb

This work shows that the partial replacement of diamagnetic Ca²⁺ by paramagnetic Tb³⁺ in Ca²⁺/calmodulin systems in solution allows the measurement of interdomain NMR pseudocontact shifts and leads to magnetic alignment of the molecule such that significant residual dipolar couplings can be measured. Both these parameters can be used to provide structural information. Species in which Tb³⁺ ions are bound to only one domain of calmodulin (the N-domain) and Ca²⁺ ions to the other (the C-domain) provide convenient systems for measuring these parameters. The nuclei in the C-domain experience the local magnetic field induced by the paramagnetic Tb³⁺ ions bound to the other domain at distances of over 40 Å from the Tb³⁺ ion, shifting the resonances for these nuclei. In addition, the Tb³⁺ ions bound to the N-domain of calmodulin greatly enhance the magnetic susceptibility anisotropy of the molecule so that a certain degree of alignment is produced due to interaction with the external magnetic field. In this way, dipolar couplings between nuclear spins are not averaged to zero due to solution molecular tumbling and yield dipolar coupling contributions to, for example, the one-bond ¹⁵N-¹H splittings of up to 17 Hz in magnitude. The degree of alignment of the C-domain will also depend on the degree of orientational freedom of this domain with respect to the N-domain containing the Tb³⁺ ions. Pseudocontact shifts for NH groups and ¹H-¹⁵N residual dipolar couplings for the directly bonded atoms have been measured for calmodulin itself, where the domains have orientational freedom, and for the complex of calmodulin with a target peptide from skeletal muscle myosin light chain kinase, where the domains have fixed orientations with respect to each other. The simultaneous measurements of these parameters for systems with domains in fixed orientations show great potential for the determination of the relative orientation of the domains.     

A target function for quaternary structural refinement from small angle scattering and NMR orientational restraints

We present a novel target function based on atomic coordinates that permits quaternary structural refinement of multi-domain protein–protein or protein–RNA complexes. It requires that the high-resolution structures of the individual domains are known and that small angle scattering (SAS) data as well as NMR orientational restraints from residual dipolar couplings (RDCs) of the complex are available. We show that, when used in combination, the translational and rotational restraints contained in SAS intensities and RDCs, respectively, define a target potential function that permits to determine the overall topology of complexes made up of domains with low internal symmetry. We apply the target function on a modestly anisotropic model system, the Barnase/Barstar complex, and discuss factors that influence the structural refinement such as data errors and the geometrical properties of the individual domains.     

Perspectives of biomolecular NMR in drug discovery: the blessing and curse of versatility

The versatility of NMR and its broad applicability to several stages in the drug discovery process is well known and generally considered one of the major strengths of NMR (Pellecchia et al., Nature Rev Drug Discov 1:211–219, 2002; Stockman and Dalvit, Prog Nucl Magn Reson Spectrosc 41:187–231, 2002; Lepre et al., Comb Chem High throughput screen 5:583–590, 2002; Wyss et al., Curr Opin Drug Discov Devel 5:630–647, 2002; Jahnke and Widmer, Cell Mol Life Sci 61:580–599, 2004; Huth et al., Methods Enzymol 394:549–571, 2005b; Klages et al., Mol Biosyst 2:318–332, 2006; Takeuchi and Wagner, Curr Opin Struct Biol 16:109–117, 2006; Zartler and Shapiro, Curr Pharm Des 12:3963–3972, 2006). Indeed, NMR is the only biophysical technique which can detect and quantify molecular interactions, and at the same time provide detailed structural information with atomic level resolution. NMR should therefore be ideally suited and widely requested as a tool for drug discovery research, and numerous examples of drug discovery projects which have substantially benefited from NMR contributions or were even driven by NMR have been described in the literature. However, not all pharmaceutical companies have rigorously implemented NMR as integral tool of their research processes. Some companies invest with limited resources, and others do not use biomolecular NMR at all. This discrepancy in assessing the value of a technology is striking, and calls for clarification—under which circumstances can NMR provide added value to the drug discovery process? What kind of contributions can NMR make, and how is it implemented and integrated for maximum impact? This perspectives article suggests key areas of impact for NMR, and a model of integrating NMR with other technologies to realize synergies and maximize their value for drug discovery.     

Developing cell-free biology for industrial applications

Although cell-free protein synthesis has been practiced for decades as a research tool, only recently have advances suggested its feasibility for commercial protein production. This focused review, based on the 2005 Amgen Award lecture, summarizes the relevant progress from the Swartz laboratory. When our program began, projected costs were much too high, proteins with disulfide bonds could not be folded effectively, and no economical scale-up technologies were available. By focusing on basic biochemical reactions and by controlling cell-free metabolism, these limitations have been methodically addressed. Amino acid supply has been stabilized and central metabolism activated to dramatically reduce substrate costs. Control of the sulfhydral redox potential has been gained and a robust disulfide isomerase added to facilitate oxidative protein folding. Finally, simple scale-up technologies have been developed. These advances not only suggest production feasibility for pharmaceutical proteins, they also provide enabling technology for producing patient-specific vaccines, for evolving new enzymes to enable biological hydrogen production from sunlight, and for developing new and highly effective water filters. Although many challenges remain, this newly expanded ability to activate and control protein production holds much promise for both research and commercial applications.     

Quantitative NMR studies of high molecular weight proteins: application to domain orientation and ligand binding in the 723 residue enzyme malate synthase G

A high-resolution multidimensional NMR study of ligand-binding to Escherichia coli malate synthase G (MSG), a 723-residue monomeric enzyme (81.4 kDa), is presented. MSG catalyzes the condensation of glyoxylate with an acetyl group of acetyl-CoA, producing malate, an intermediate in the citric-acid cycle. We show that despite the size of the protein, important structural and dynamic information about the molecule can be obtained on a per-residue basis. 15N-1HN residual dipolar couplings and carbonyl chemical shift changes upon alignment in Pf1 phage establish that there are no significant domain reorientations in the molecule upon ligand binding, in contrast to what was anticipated on the basis of both the X-ray structure of the glyoxylate-bound form of the enzyme and structural studies of a related set of proteins. The chemical shift changes of 1HN, 15N and 13CO nuclei upon binding of pyruvate, a glyoxylate-mimicking inhibitor, and acetyl-CoA have been mapped onto the three-dimensional structure of the molecule. Binding constants of pyruvate, glyoxylate, and acetyl-CoA (in the presence of pyruvate) have been measured, along with the kinetic parameters for glyoxylate and pyruvate binding. The on-rates of pyruvate and glyoxalate binding, approximately 1.2 x 10(6)M(-1)s(-1) and approximately 2.7 x 10(6)M(-1)s(-1), respectively, are significantly lower than what is anticipated from a simple diffusion-controlled process. Some structural implications of the chemical shift perturbations upon binding and the estimated ligand on-rates are discussed.     

NMR methods for studying quadruplex nucleic acids

Solution NMR spectroscopy has traditionally played a central role in examining quadruplex structure, dynamics, and interactions. Here, an overview is given of the methods currently applied to structural, dynamics, thermodynamics, and kinetics studies of nucleic acid quadruplexes and associated cations.     

A large data set comparison of protein structures determined by crystallography and NMR: statistical test for structural differences and the effect of crystal packing

The existence of a large number of proteins for which both nuclear magnetic resonance (NMR) and X-ray crystallographic coordinates have been deposited into the Protein Data Bank (PDB) makes the statistical comparison of the corresponding crystal and NMR structural models over a large data set possible, and facilitates the study of the effect of the crystal environment and other factors on structure. We present an approach for detecting statistically significant structural differences between crystal and NMR structural models which is based on structural superposition and the analysis of the distributions of atomic positions relative to a mean structure. We apply this to a set of 148 protein structure pairs (crystal vs NMR), and analyze the results in terms of methodological and physical sources of structural difference. For every one of the 148 structure pairs, the backbone root-mean-square distance (RMSD) over core atoms of the crystal structure to the mean NMR structure is larger than the average RMSD of the members of the NMR ensemble to the mean, with 76% of the structure pairs having an RMSD of the crystal structure to the mean more than a factor of two larger than the average RMSD of the NMR ensemble. On average, the backbone RMSD over core atoms of crystal structure to the mean NMR is approximately 1 A. If non-core atoms are included, this increases to 1.4 A due to the presence of variability in loops and similar regions of the protein. The observed structural differences are only weakly correlated with the age and quality of the structural model and differences in conditions under which the models were determined. We examine steric clashes when a putative crystalline lattice is constructed using a representative NMR structure, and find that repulsive crystal packing plays a minor role in the observed differences between crystal and NMR structures. The observed structural differences likely have a combination of physical and methodological causes. Stabilizing attractive interactions arising from intermolecular crystal contacts which shift the equilibrium of the crystal structure relative to the NMR structure is a likely physical source which can account for some of the observed differences. Methodological sources of apparent structural difference include insufficient sampling or other issues which could give rise to errors in the estimates of the precision and/or accuracy.     

Unique polyamines produced by an extreme thermophile, <Emphasis Type="Italic">Thermus thermophilus</Emphasis>

Summary. Recent research progress on polyamines in extreme thermophiles is reviewed. Extreme thermophiles produce two types of unique polyamines; one is longer polyamines such as caldopentamine and caldohexamine, and the other is branched polyamines such as tetrakis(3-aminopropyl)ammonium. The protein synthesis catalyzed by a cell-free extract of Thermus thermophilus, an extreme thermophile, required the presence of a polyamine and the highest activity was found in the presence of tetrakis(3-aminopropyl)ammonium. In vitro experiments, longer polyamines efficiently stabilized double stranded nucleic acids and a branched polyamine, tetrakis(3-aminropyl)ammonium, stabilized stem-and-loop structures. In T. thermophilus, polyamines are synthesized from arginine by a new metabolic pathway; arginine is converted to agmatine and then agmatine is aminopropylated to N¹-aminopropylagmatine which is converted to spermidine by an enzyme coded by a gene homologous to speB (a gene for agmatinase). In this new pathway spermidine is not synthesized from putrescine. Reverse genetic studies indicated that the unique polyamines are synthesized from spermidine.     

Measurement of carbonyl chemical shifts of excited protein states by relaxation dispersion NMR spectroscopy: comparison between uniformly and selectively 13C labeled samples

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion nuclear magnetic resonance (NMR) spectroscopy has emerged as a powerful method for quantifying chemical shifts of excited protein states. For many applications of the technique that involve the measurement of relaxation rates of carbon magnetization it is necessary to prepare samples with isolated (13)C spins so that experiments do not suffer from magnetization transfer between coupled carbon spins that would otherwise occur during the CPMG pulse train. In the case of (13)CO experiments however the large separation between (13)CO and (13)C(alpha) chemical shifts offers hope that robust (13)CO dispersion profiles can be recorded on uniformly (13)C labeled samples, leading to the extraction of accurate (13)CO chemical shifts of the invisible, excited state. Here we compare such chemical shifts recorded on samples that are selectively labeled, prepared using [1-(13)C]-pyruvate and NaH(13)CO(3,) or uniformly labeled, generated from (13)C-glucose. Very similar (13)CO chemical shifts are obtained from analysis of CPMG experiments recorded on both samples, and comparison with chemical shifts measured using a second approach establishes that the shifts measured from relaxation dispersion are very accurate.     

NMR experiments for the measurement of proton-proton and carbon-carbon residual dipolar couplings in uniformly labelled oligosaccharides

A 2D-HSQC-carbon selective/proton selective-constant time COSY, 2D-HSQC-(sel C, sel H)-CT COSY experiment, which is applicable to uniformly ¹³C isotopically enriched samples (U-¹³C) of oligosaccharides or oligonucleotides is proposed for the measurement of proton–proton RDC in crowded regions of 2D-spectra. In addition, a heteronuclear constant time-COSY experiment, ¹³C-¹³C CT-COSY, is proposed for the measurement of one bond carbon–carbon RDC in these molecules. These two methods provide an extension, to U-¹³C molecules, of the original homonuclear constant time-COSY experiment proposed by Tian et al. (1999) for saccharides. The combination of a number of these RDC with NOE data may provide the method of choice to study oligosaccharide conformation in the free and receptor-bound state.     

The structural analysis of protein–protein interactions by NMR spectroscopy

A comprehensive understanding of protein–protein interactions is an important next step in our quest to understand how the information contained in a genome is put into action. Although a number of experimental techniques can report on the existence of a protein– protein interaction, very few can provide detailed structural information. NMR spectroscopy is one of these, and in recent years several complementary NMR approaches, including residual dipolar couplings and the use of paramagnetic effects, have been developed that can provide insight into the structure of protein–protein complexes. In this article, we review these approaches and comment on their strengths and weaknesses.     

Heme methyl 1H chemical shifts as structural parameters in some low-spin ferriheme proteins

 The different paramagnetic shifts of the four methyl groups in ferriheme proteins have been described as being due to the effect of the axial ligand nodal plane orientation. An equation, heuristically found and theoretically explained, describing the relation between contact and pseudocontact shifts and the position of the axial ligand(s) has been derived for bis-histidine ferriheme proteins and for cyanide-histidine ferriheme proteins. The values of the heuristic parameters contained in the equations were found by fitting the shifts of bovine cytochrome b ₅ and several bis-histidine cytochromes c ₃ and histidine-cyanide systems. The agreement between the observed and the calculated shifts was found to be good. Therefore, by taking advantage of this study, information on the position of the axial ligands, that can be used as a constraint for structure determination, can be obtained from the shifts of the methyl protons. Received: 13 April 1999 / Accepted: 4 June 1999     

Helix orientations in membrane-associated Bcl-X<Subscript>L</Subscript> determined by <Superscript>15</Superscript>N-solid-state NMR spectroscopy

Controlled cell death is fundamental to tissue hemostasis and apoptosis malfunctions can lead to a wide range of diseases. Bcl-x_L is an anti-apoptotic protein the function of which is linked to its reversible interaction with mitochondrial outer membranes. Its interfacial and intermittent bilayer association makes prediction of its bound structure difficult without using methods able to extract data from dynamic systems. Here we investigate Bcl-x_L associated with oriented lipid bilayers at physiological pH using solid-state NMR spectroscopy. The data are consistent with a C-terminal transmembrane anchoring sequence and an average alignment of the remaining helices, i.e. including helices 5 and 6, approximately parallel to the membrane surface. Data from several biophysical approaches confirm that after removal of the C-terminus from Bcl-x_L its membrane interactions are weak. In the presence of membranes Bcl-x_L can still interact with a Bak BH3 domain peptide suggesting a model where the hydrophobic C-terminus of the protein unfolds and inserts into the membrane. During this conformational change the Bcl-x_L hydrophobic binding pocket becomes accessible for protein–protein interactions whilst the structure of the N-terminal region remains intact. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Theoretical framework for NMR residual dipolar couplings in unfolded proteins

A theoretical framework for the prediction of nuclear magnetic resonance (NMR) residual dipolar couplings (RDCs) in unfolded proteins under weakly aligning conditions is presented. The unfolded polypeptide chain is modeled as a random flight chain while the alignment medium is represented by a set of regularly arranged obstacles. For the case of bicelles oriented perpendicular to the magnetic field, a closed-form analytical result is derived. With the obtained analytical expression the RDCs are readily accessible for any locus along the chain, for chains of differing length, and for varying bicelle concentrations. The two general features predicted by the model are (i) RDCs in the center segments of a polypeptide chain are larger than RDCs in the end segments, resulting in a bell-shaped sequential distribution of RDCs, and (ii) couplings are larger for shorter chains than for longer chains at a given bicelle concentration. Experimental data available from the literature confirm the first prediction of the model, providing a tool for recognizing fully unfolded polypeptide chains. With less certainty experimental data appear to support the second prediction as well. However, more systematic experimental studies are needed in order to validate or disprove the predictions of the model. The presented framework is an important step towards a solid theoretical foundation for the analysis of experimentally measured RDCs in unfolded proteins in the case of alignment media such as polyacrylamide gels and neutral bicelle systems which align biomacromolecules by a steric mechanism. Various improvements and generalizations are possible within the suggested approach.     

The dynamic duo: Combining NMR and small angle scattering in structural biology

Structural biology provides essential information for elucidating molecular mechanisms that underlie biological function. Advances in hardware, sample preparation, experimental methods, and computational approaches now enable structural analysis of protein complexes with increasing complexity that more closely represent biologically entities in the cellular environment. Integrated multidisciplinary approaches are required to overcome limitations of individual methods and take advantage of complementary aspects provided by different structural biology techniques. Although X‐ray crystallography remains the method of choice for structural analysis of large complexes, crystallization of flexible systems is often difficult and does typically not provide insights into conformational dynamics present in solution. Nuclear magnetic resonance spectroscopy (NMR) is well‐suited to study dynamics at picosecond to second time scales, and to map binding interfaces even of large systems at residue resolution but suffers from poor sensitivity with increasing molecular weight. Small angle scattering (SAS) methods provide low resolution information in solution and can characterize dynamics and conformational equilibria complementary to crystallography and NMR. The combination of NMR, crystallography, and SAS is, thus, very useful for analysis of the structure and conformational dynamics of (large) protein complexes in solution. In high molecular weight systems, where NMR data are often sparse, SAS provides additional structural information and can differentiate between NMR‐derived models. Scattering data can also validate the solution conformation of a crystal structure and indicate the presence of conformational equilibria. Here, we review current state‐of‐the‐art approaches for combining NMR, crystallography, and SAS data to characterize protein complexes in solution.     

Fast approximations of the rotational diffusion tensor and their application to structural assembly of molecular complexes

We present and evaluate a rigid-body, deterministic, molecular docking method, called ELMDOCK, that relies solely on the three-dimensional structure of the individual components and the overall rotational diffusion tensor of the complex, obtained from nuclear spin-relaxation measurements. We also introduce a docking method, called ELMPATIDOCK, derived from ELMDOCK and based on the new concept of combining the shape-related restraints from rotational diffusion with those from residual dipolar couplings, along with ambiguous contact/interface-related restraints obtained from chemical shift perturbations. ELMDOCK and ELMPATIDOCK use two novel approximations of the molecular rotational diffusion tensor that allow computationally efficient docking. We show that these approximations are accurate enough to properly dock the two components of a complex without the need to recompute the diffusion tensor at each iteration step. We analyze the accuracy, robustness, and efficiency of these methods using synthetic relaxation data for a large variety of protein-protein complexes. We also test our method on three protein systems for which the structure of the complex and experimental relaxation data are available, and analyze the effect of flexible unstructured tails on the outcome of docking. Additionally, we describe a method for integrating the new approximation methods into the existing docking approaches that use the rotational diffusion tensor as a restraint. The results show that the proposed docking method is robust against experimental errors in the relaxation data or structural rearrangements upon complex formation and is computationally more efficient than current methods. The developed approximations are accurate enough to be used in structure refinement protocols.     

Principal component method for assessing structural heterogeneity across multiple alignment media

The recent availability of residual dipolar coupling measurements in a variety of different alignment media raises the question to what extent biomolecular structure and dynamics are differentially affected by their presence. A computational method is presented that allows the sensitive assessment of such changes using dipolar couplings measured in six or more alignment media. The method is based on a principal component analysis of the covariance matrix of the dipolar couplings. It does not require a priori structural or dynamic information nor knowledge of the alignment tensors and their orientations. In the absence of experimental errors, the covariance matrix has at most five nonzero eigenvalues if the structure and dynamics of the biomolecule is the same in all media. In contrast, differential structural and dynamic changes lead to additional nonzero eigenvalues. Characteristic features of the eigenvalue distribution in the absence and presence of noise are discussed using dipolar coupling data calculated from conformational ensembles taken from a molecular dynamics trajectory of native ubiquitin.