Structures of protein complexes by multidimensional heteronuclear magnetic resonance spectroscopy

Protein stability governed by its structural plasticity is inferred by physicochemical factors and salt bridges

Several organisms, specifically microorganisms survive in a wide range of harsh environments including extreme temperature, pH, and salt concentration. We analyzed systematically a large number of protein sequences with their structures to understand their stability and to discriminate extremophilic proteins from their non-extremophilic orthologs. Our results highlighted that the strategy for the packing of the protein core was influenced by the environmental stresses through substitutive structural events through better ionic interaction. Statistical analysis showed that a significant difference in number and composition of amino acid exist among them. The negative correlation of pairwise sequence alignments and structural alignments indicated that most of the extremophile and non-extremophile proteins didn’t contain any association for maintaining their functional stability. A significant numbers of salt bridges were noticed on the surface of the extremostable proteins. The Ramachandran plot data represented more occurrences of amino acids being present in helix and sheet regions of extremostable proteins. We also found that a significant number of small nonpolar amino acids and moderate number of charged amino acids like Arginine and Aspartic acid represented more nonplanar Omega angles in their peptide bond. Thus, extreme conditions may predispose amino acid composition including geometric variability for molecular adaptation of extremostable proteins against atmospheric variations and associated changes under natural selection pressure. The variation of amino acid composition and structural diversifications in proteins play a major role in evolutionary adaptation to mitigate climate change.

Characterization of Dynamic IDP Complexes by NMR Spectroscopy

NMR spectroscopy has proven to be a key method for studying intrinsically disordered proteins (IDPs). Nonetheless, traditional NMR methods developed for solving structures of ordered protein complexes are insufficient for the full characterization of dynamic IDP complexes, where the energy landscape is broader and more rugged. Furthermore, due to their high sensitivity to environmental changes, NMR studies of IDP complexes must be conducted with extra care and the observed NMR parameters thoroughly evaluated to enable disentanglement of binding events from ensemble distribution changes. In this chapter, written for the non-NMR expert, we start out by outlining sample preparation for IDP complexes, guide through the recording and evaluation of diagnostic ¹H,¹⁵N-HSQC spectra, and delineate more sophisticated NMR strategies to follow for the particular type of complex. The most relevant experiments are then described in terms of aims, needs, pitfalls, analysis, and expected outcomes, with references to recent examples.

Chemical shift-based methods in NMR structure determination

Chemical shifts are highly sensitive probes harnessed by NMR spectroscopists and structural biologists as conformational parameters to characterize a range of biological molecules. Traditionally, assignment of chemical shifts has been a labor-intensive process requiring numerous samples and a suite of multidimensional experiments. Over the past two decades, the development of complementary computational approaches has bolstered the analysis, interpretation and utilization of chemical shifts for elucidation of high resolution protein and nucleic acid structures. Here, we review the development and application of chemical shift-based methods for structure determination with a focus on ab initio fragment assembly, comparative modeling, oligomeric systems, and automated assignment methods. Throughout our discussion, we point out practical uses, as well as advantages and caveats, of using chemical shifts in structure modeling. We additionally highlight (i) hybrid methods that employ chemical shifts with other types of NMR restraints (residual dipolar couplings, paramagnetic relaxation enhancements and pseudocontact shifts) that allow for improved accuracy and resolution of generated 3D structures, (ii) the utilization of chemical shifts to model the structures of sparsely populated excited states, and (iii) modeling of sidechain conformations. Finally, we briefly discuss the advantages of contemporary methods that employ sparse NMR data recorded using site-specific isotope labeling schemes for chemical shift-driven structure determination of larger molecules. With this review, we aim to emphasize the accessibility and versatility of chemical shifts for structure determination of challenging biological systems, and to point out emerging areas of development that lead us towards the next generation of tools.

High quality NMR structures: a new force field with implicit water and membrane solvation for Xplor-NIH

Structure determination of proteins by NMR is unique in its ability to measure restraints, very accurately, in environments and under conditions that closely mimic those encountered in vivo. For example, advances in solid-state NMR methods enable structure determination of membrane proteins in detergent-free lipid bilayers, and of large soluble proteins prepared by sedimentation, while parallel advances in solution NMR methods and optimization of detergent-free lipid nanodiscs are rapidly pushing the envelope of the size limit for both soluble and membrane proteins. These experimental advantages, however, are partially squandered during structure calculation, because the commonly used force fields are purely repulsive and neglect solvation, Van der Waals forces and electrostatic energy. Here we describe a new force field, and updated energy functions, for protein structure calculations with EEFx implicit solvation, electrostatics, and Van der Waals Lennard-Jones forces, in the widely used program Xplor-NIH. The new force field is based primarily on CHARMM22, facilitating calculations with a wider range of biomolecules. The new EEFx energy function has been rewritten to enable OpenMP parallelism, and optimized to enhance computation efficiency. It implements solvation, electrostatics, and Van der Waals energy terms together, thus ensuring more consistent and efficient computation of the complete nonbonded energy lists. Updates in the related python module allow detailed analysis of the interaction energies and associated parameters. The new force field and energy function work with both soluble proteins and membrane proteins, including those with cofactors or engineered tags, and are very effective in situations where there are sparse experimental restraints. Results obtained for NMR-restrained calculations with a set of five soluble proteins and five membrane proteins show that structures calculated with EEFx have significant improvements in accuracy, precision, and conformation, and that structure refinement can be obtained by short relaxation with EEFx to obtain improvements in these key metrics. These developments broaden the range of biomolecular structures that can be calculated with high fidelity from NMR restraints.

A practical implicit solvent potential for NMR structure calculation

The benefits of protein structure refinement in water are well documented. However, performing structure refinement with explicit atomic representation of the solvent molecules is computationally expensive and impractical for NMR-restrained structure calculations that start from completely extended polypeptide templates. Here we describe a new implicit solvation potential, EEFx (Effective Energy Function for XPLOR-NIH), for NMR-restrained structure calculations of proteins in XPLOR-NIH. The key components of EEFx are an energy term for solvation energy that works together with other nonbonded energy functions, and a dedicated force field for conformational and nonbonded protein interaction parameters. The initial results obtained with EEFx show that significant improvements in structural quality can be obtained. EEFx is computationally efficient and can be used both to fold and refine structures. Overall, EEFx improves the quality of protein conformation and nonbonded atomic interactions. Moreover, such benefits are accompanied by enhanced structural precision and enhanced structural accuracy, reflected in improved agreement with the cross-validated dipolar coupling data. Finally, implementation of EEFx calculations is straightforward and computationally efficient. Overall, EEFx provides a useful method for the practical calculation of experimental protein structures in a physically realistic environment.

Chirally deuterated benzyl chlorides from benzyl alcohols via hexachloroacetone/polymer-supported triphenylphosphine: synthesis of protected (2S, 3S)-[3-(2)H, (15)N]-tyrosine

Chirally deuterated benzyl chlorides were prepared using novel, general hexachloroacetone/polymer-supported triphenylphosphine treatment of chirally deuterated benzyl alcohols. Doubly labeled protected tyrosine was obtained in 62% yield with 86% de at the α-carbon and 82% de at the β-carbon. Key in the synthesis was the alkylation of (15)N-labeled (-)-8-phenylmenthylhippurate with R-(-)-4-triisopropylsilyloxybenzyl-α-d chloride.

The structure of a truncated phosphoribosylanthranilate isomerase suggests a unified model for evolution of the (βα)8 barrel fold

The (βα)(8) barrel is one of the most common protein folds, and enzymes with this architecture display a remarkable range of catalytic activities. Many of these functions are associated with ancient metabolic pathways, and phylogenetic reconstructions suggest that the (βα)(8) barrel was one of the very first protein folds to emerge. Consequently, there is considerable interest in understanding the evolutionary processes that gave rise to this fold. In particular, much attention has been focused on the plausibility of (βα)(8) barrel evolution from homodimers of half barrels. However, we previously isolated a three-quarter-barrel-sized fragment of a (βα)(8) barrel, termed truncated phosphoribosylanthranilate isomerase (trPRAI), that is soluble and almost as thermostable as full-length N-(5'-phosphoribosyl)anthranilate isomerase (PRAI). Here, we report the NMR-derived structure of trPRAI. The subdomain is monomeric, is well ordered and adopts a native-like structure in solution. Side chains from strands β(1) (Glu3 and Lys5), β(2) (Tyr25) and β(6) (Lys122) of trPRAI repack to shield the hydrophobic core from the solvent. This result demonstrates that three-quarter barrels were viable intermediates in the evolution of the (βα)(8) barrel fold. We propose a unified model for (βα)(8) barrel evolution that combines our data, previously published work and plausible scenarios for the emergence of (initially error-prone) genetic systems. In this model, the earliest proto-cells contained diverse pools of part-barrel subdomains. Combinatorial assembly of these subdomains gave rise to many distinct lineages of (βα)(8) barrel proteins, that is, our model excludes the possibility that there was a single (βα)(8) barrel from which all present examples are descended.

Structural modification of acyl carrier protein by butyryl group

Fatty acid synthesis in bacteria is catalyzed by a set of individual enzymes known as the type II fatty acid synthase. Acyl carrier protein (ACP) shuttles the acyl intermediates between individual pathway enzymes. In this study, we determined the solution structures of three different forms of ACP, apo-ACP, ACP, and butyryl-ACP under identical experimental conditions. The structural studies revealed that attachment of butyryl acyl intermediate to ACP alters the conformation of ACP. This finding supports the more general notion that the attachment of different acyl intermediates alters the ACP structure to facilitate their recognition and turnover by the appropriate target enzymes.

The spliceosomal phosphopeptide P140 controls the lupus disease by interacting with the HSC70 protein and via a mechanism mediated by gammadelta T cells

The phosphopeptide P140 issued from the spliceosomal U1-70K snRNP protein is recognized by lupus CD4(+) T cells, transiently abolishes T cell reactivity to other spliceosomal peptides in P140-treated MRL/lpr mice, and ameliorates their clinical features. P140 modulates lupus patients' T cell response ex vivo and is currently included in phase IIb clinical trials. Its underlying mechanism of action remains elusive. Here we show that P140 peptide binds a unique cell-surface receptor, the constitutively-expressed chaperone HSC70 protein, known as a presenting-protein. P140 induces apoptosis of activated MRL/lpr CD4(+) T cells. In P140-treated mice, it increases peripheral blood lymphocyte apoptosis and decreases B cell, activated T cell, and CD4(-)CD8(-)B220(+) T cell counts via a specific mechanism strictly depending on gammadelta T cells. Expression of inflammation-linked genes is rapidly regulated in CD4(+) T cells. This work led us to identify a powerful pathway taken by a newly-designed therapeutic peptide to immunomodulate lupus autoimmunity.

Discrimination between biological interfaces and crystal-packing contacts

A discrimination method between biologically relevant interfaces and artificial crystal-packing contacts in crystal structures was constructed. The method evaluates protein-protein interfaces in terms of complementarities for hydrophobicity, electrostatic potential and shape on the protein surfaces, and chooses the most probable biological interfaces among all possible contacts in the crystal. The method uses a discriminator named as "COMP", which is a linear combination of the complementarities for the above three surface features and does not correlate with the contact area. The discrimination of homo-dimer interfaces from symmetry-related crystal-packing contacts based on the COMP value achieved the modest success rate. Subsequent detailed review of the discrimination results raised the success rate to about 88.8%. In addition, our discrimination method yielded some clues for understanding the interaction patterns in several examples in the PDB. Thus, the COMP discriminator can also be used as an indicator of the "biological-ness" of protein-protein interfaces.

The point mutation A34F causes dimerization of GB1

The immunoglobulin-binding domain B1 of streptococcal protein G (GB1), a very stable, small, single-domain protein, is one of the most extensively used models in the area of protein folding and design. Variants derived from a library of randomized hydrophobic core residues previously revealed alternative folds, namely a completely intertwined tetramer (Frank et al., Nat Struct Biol 2002;9:877-885) and a domain-swapped dimer (Byeon et al., J Mol Biol 2003;333:141-152). Here, we report the NMR structure of the single amino acid mutant Ala-34-Phe which exists as side-by-side dimer. The dimer dissociation constant is 27 +/- 4 microM. The dimer interface comprises two structural elements: First, the beta-sheets of the two monomers pair in an antiparallel arrangement, thereby forming an eight-stranded beta-sheet. Second, the alpha-helix is shortened, ending in a loop that engages in intermolecular contacts. The largest difference between the monomer unit in the A34F dimer and the monomeric wild-type GB1 is the dissolution of the C-terminal half of the alpha-helix associated with a pronounced slow conformational motion of the interface loop. This involves a large movement of the Tyr-33 side chain that swings out from the monomer to engage in dimer contacts.

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.

NMR and X-ray analysis of structural additivity in metal binding site-swapped hybrids of rubredoxin

Background

Chimeric hybrids derived from the rubredoxins of Pyrococcus furiosus (Pf) and Clostridium pasteurianum (Cp) provide a robust system for the characterization of protein conformational stability and dynamics in a differential mode. Interchange of the seven nonconserved residues of the metal binding site between the Pf and Cp rubredoxins yields a complementary pair of hybrids, for which the sum of the thermodynamic stabilities is equal to the sum for the parental proteins. Furthermore, the increase in amide hydrogen exchange rates for the hyperthermophile-derived metal binding site hybrid is faithfully mirrored by a corresponding decrease for the complementary hybrid that is derived from the less thermostable rubredoxin, indicating a degree of additivity in the conformational fluctuations that underlie these exchange reactions.

Results

Initial NMR studies indicated that the structures of the two complementary hybrids closely resemble "cut-and-paste" models derived from the parental Pf and Cp rubredoxins. This protein system offers a robust opportunity to characterize differences in solution structure, permitting the quantitative NMR chemical shift and NOE peak intensity data to be analyzed without recourse to the conventional conversion of experimental NOE peak intensities into distance restraints. The intensities for 1573 of the 1652 well-resolved NOE crosspeaks from the hybrid rubredoxins were statistically indistinguishable from the intensities of the corresponding parental crosspeaks, to within the baseplane noise level of these high sensitivity data sets. The differences in intensity for the remaining 79 NOE crosspeaks were directly ascribable to localized dynamical processes. Subsequent X-ray analysis of the metal binding site-swapped hybrids, to resolution limits of 0.79 Å and 1.04 Å, demonstrated that the backbone and sidechain heavy atoms in the NMR-derived structures lie within the range of structural variability exhibited among the individual molecules in the crystallographic asymmetric unit (~0.3 Å), indicating consistency with the "cut-and-paste" structuring of the hybrid rubredoxins in both crystal and solution.

Conclusion

Each of the significant energetic interactions in the metal binding site-swapped hybrids appears to exhibit a 1-to-1 correspondence with the interactions present in the corresponding parental rubredoxin structure, thus providing a structural basis for the observed additivity in conformational stability and dynamics. The congruence of these X-ray and NMR experimental data offers additional support for the interpretation that the conventional treatment of NOE distance restraints contributes substantially to the systematic differences that are commonly reported between NMR- and X-ray-derived protein structures.

Hydrophobic potential of mean force as a solvation function for protein structure prediction

We have developed a solvation function that combines a Generalized Born model for polarization of protein charge by the high dielectric solvent, with a hydrophobic potential of mean force (HPMF) as a model for hydrophobic interaction, to aid in the discrimination of native structures from other misfolded states in protein structure prediction. We find that our energy function outperforms other reported scoring functions in terms of correct native ranking for 91% of proteins and low Z scores for a variety of decoy sets, including the challenging Rosetta decoys. This work shows that the stabilizing effect of hydrophobic exposure to aqueous solvent that defines the HPMF hydration physics is an apparent improvement over solvent-accessible surface area models that penalize hydrophobic exposure. Decoys generated by thermal sampling around the native-state basin reveal a potentially important role for side-chain entropy in the future development of even more accurate free energy surfaces.

On side-chain conformational entropy of proteins

The role of side-chain entropy (SCE) in protein folding has long been speculated about but is still not fully understood. Utilizing a newly developed Monte Carlo method, we conducted a systematic investigation of how the SCE relates to the size of the protein and how it differs among a protein's X-ray, NMR, and decoy structures. We estimated the SCE for a set of 675 nonhomologous proteins, and observed that there is a significant SCE for both exposed and buried residues for all these proteins-the contribution of buried residues approaches approximately 40% of the overall SCE. Furthermore, the SCE can be quite different for structures with similar compactness or even similar conformations. As a striking example, we found that proteins' X-ray structures appear to pack more "cleverly" than their NMR or decoy counterparts in the sense of retaining higher SCE while achieving comparable compactness, which suggests that the SCE plays an important role in favouring native protein structures. By including a SCE term in a simple free energy function, we can significantly improve the discrimination of native protein structures from decoys.

A simple and reliable approach to docking protein-protein complexes from very sparse NOE-derived intermolecular distance restraints

A simple and reliable approach for docking protein-protein complexes from very sparse NOE-derived intermolecular distance restraints (as few as three from a single point) in combination with a novel representation for an attractive potential between mapped interaction surfaces is described. Unambiguous assignments of very sparse intermolecular NOEs are obtained using a reverse labeling strategy in which one the components is fully deuterated with the exception of selective protonation of the delta-methyl groups of isoleucine, while the other component is uniformly (13)C-labeled. This labeling strategy can be readily extended to selective protonation of Ala, Leu, Val or Met. The attractive potential is described by a 'reduced' radius of gyration potential applied specifically to a subset of interfacial residues (those with an accessible surface area > or = 50% in the free proteins) that have been delineated by chemical shift perturbation. Docking is achieved by rigid body minimization on the basis of a target function comprising the sparse NOE distance restraints, a van der Waals repulsion potential and the 'reduced' radius of gyration potential. The method is demonstrated for two protein-protein complexes (EIN-HPr and IIA(Glc)-HPr) from the bacterial phosphotransferase system. In both cases, starting from 100 different random orientations of the X-ray structures of the free proteins, 100% convergence is achieved to a single cluster (with near identical atomic positions) with an overall backbone accuracy of approximately 2 A. The approach described is not limited to NMR, since interfaces can also be mapped by alanine scanning mutagenesis, and sparse intermolecular distance restraints can be derived from double cycle mutagenesis, cross-linking combined with mass spectrometry, or fluorescence energy transfer.

Error distribution derived NOE distance restraints

Errors and imprecisions in distance restraints derived from NOESY peak volumes are usually accounted for by generous lower and upper bounds on the distances. In this paper, we propose a new form of distance restraints, replacing the subjective bounds by a potential function obtained from the error distribution of the distances. We derived the shape of the potential from molecular dynamics calculations and by comparison of NMR data with X-ray crystal structures. We used complete cross-validation to derive the optimal weight for the data in the calculation. In a model system with synthetic restraints, the accuracy of the structures improved significantly compared to calculations with the usual form of restraints. For experimental data sets, the structures systematically approach the X-ray crystal structures of the same protein. Also standard quality indicators improve compared to standard calculations. The results did not depend critically on the exact shape of the potential. The new approach is less subjective and uses fewer assumptions in the interpretation of NOESY peak volumes as distance restraints than the usual approach. Figures of merit for the structures, such as the RMS difference from the average structure or the RMS difference from the data, are therefore less biased and more meaningful measures of structure quality than with the usual form of restraints.

Phosphorylated and nonphosphorylated epitopes of the La/SSB autoantigen: comparison of their antigenic and conformational characteristics

La/SSB phosphoprotein is the target antigen of autoantibodies in sera of patients with Sjögren's syndrome (SS) and systemic lupus erythematosus (SLE). Among other structural and function motifs, four phosphorylation sites are encompassed in the primary sequence of La/SSB. Two of them (Thr-362 and Ser-366) are located within GSGKGKVQFQGKKTKFASDD (346-368) and one (Thr-302) within VTWEVLEGEVEKEALKKI (301-318), which are main B-cell epitopes of La/SSB. With the aim to investigate how phosphorylation, one of the most common posttranslational protein modifications, affects the antigenic and conformational characteristics of the La/SSB epitopes, we synthesized and studied the phosphorylated epitopes La/SSB(346-368)-P, La/SSB(359-368)-P, and La/SSB(301-318)-P with respect to their nonphosphorylated counterparts. Anti-La/SSB positive sera from SS and SLE patients are better recognized by the phosphorylated epitopes compared to their nonphosphorylated counterparts. Conformational analysis by (1)H nuclear magnetic resonance spectroscopy and molecular dynamics showed that the phosphorylated epitopes adopt different structural characteristics from those of the corresponding nonphosphorylated epitopes. It is concluded that phosphorylation can create neoepitopes with altered functions, compared to the nonphosphorylated epitopes, which might be seen from the immune system as "foreign."

Structure of bacteriocin AS-48: from soluble state to membrane bound state

The bacteriocin AS-48 is a membrane-interacting peptide, which displays a broad anti-microbial spectrum against Gram-positive and Gram-negative bacteria. The NMR structure of AS-48 at pH 3 has been solved. The analysis of this structure suggests that the mechanism of AS-48 anti-bacterial activity involves the accumulation of positively charged molecules at the membrane surface leading to a disruption of the membrane potential. Here, we report the high-resolution crystal structure of AS-48 and sedimentation equilibrium experiments showing that this bacteriocin is able to adopt different oligomeric structures according to the physicochemical environment. The analysis of these structures suggests a mechanism for molecular function of AS-48 involving a transition from a water-soluble form to a membrane-bound state upon membrane binding.

Novel p-arylthio cinnamides as antagonists of leukocyte function-associated antigen-1/intracellular adhesion molecule-1 interaction. 2. Mechanism of inhibition and structure-based improvement of pharmaceutical properties

The interaction between leukocyte function-associated antigen-1 (LFA-1) and intracellular adhesion molecule-1 (ICAM-1) has been implicated in inflammatory and immune diseases. Recently, a novel series of p-arylthio cinnamides has been described as potent antagonists of the LFA-1/ICAM-1 interaction. These compounds were found to bind to the I domain of LFA-1 using two-dimensional NMR spectroscopy of 15N-labeled LFA-1 I domain. On the basis of NOE studies between compound 1 and the I domain of LFA-1, a model of the complex was constructed. This model revealed that compound 1 does not directly inhibit ICAM-1 binding by interacting with the metal ion dependent adhesion site (MIDAS). Instead, it binds to the previously proposed I domain allosteric site (IDAS) of LFA-1 and likely modulates the activation of LFA-1 through its interaction with this regulatory site. A fragment-based NMR screening strategy was applied to identify small, more water-soluble ligands that bind to a specific region of the IDAS. When incorporated into the parent cinnamide template, the resulting analogues exhibited increased aqueous solubility and improved pharmacokinetic profiles in rats, demonstrating the power of this NMR-based screening approach for rapidly modifying high-affinity ligands.

NMR-restrained docking of a peptidic inhibitor to the N-terminal domain of the phosphoenolpyruvate:sugar phosphotransferase enzyme I

Starting from the NMR structure of the binary complex between the N-terminal domain of the unphosphorylated enzyme I (EIN) of the phosphoenolpyruvate:sugar phosphotransferase (PTS) and the histidine-containing phosphocarrier protein (HPr), a molecular model of the phosphorylated transition state of the related complex was established using constrained simulated annealing. The coordinates of the phosphorylated EIN enzyme were then used in a second step for flexible docking of a decapeptide inhibitor of EIN whose enzyme-bound conformation itself was determined by NMR using transferred nuclear Overhauser effects. Two phosphorylation models of the peptide inhibitor were investigated and shown to be both functional. Interestingly, one model is very similar to that of the complex between EIN and its natural substrate HPr. The present study demonstrates that NMR-guided flexible docking constitutes an interesting tool for docking highly flexible peptide ligands and facilitates the upcoming protein-based design of nonpeptide EIN inhibitors for discovering new antibiotics.

Structure determination and binding kinetics of a DNA aptamer-argininamide complex

The structure of a DNA aptamer, which was selected for specific binding to arginine, was determined using NMR spectroscopy. The sequence forms a hairpin loop, with residues important for binding occurring in the loop region. Binding of argininamide induces formation of one Watson-Crick and two non-Watson-Crick base pairs, which facilitate generation of a binding pocket. The specificity for arginine seems to arise from contacts between the guanidino end of the arginine and phosphates, with atoms positioned by the shape of the pocket. Complex binding kinetics are observed suggesting that there is a slow interconversion of two forms of the DNA, which have different binding affinities. These data provide information on the process of adaptive recognition of a ligand by an aptamer.

Exploring the dynamic information content of a protein NMR structure: comparison of a molecular dynamics simulation with the NMR and X-ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7-ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high-resolution (1.48 A) crystal structure than with either the lower-resolution (2.05 A) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side-chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn-pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean-square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non-rotameric/syn-pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X-ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed.

NMR solution structure of the catalytic fragment of human fibroblast collagenase complexed with a sulfonamide derivative of a hydroxamic acid compound

The solution structure of the catalytic fragment of human fibroblast collagenase (MMP-1) complexed with a sulfonamide derivative of a hydroxamic acid compound (CGS-27023A) has been determined using two-dimensional and three-dimensional heteronuclear NMR spectroscopy. The solution structure of the complex was calculated by means of hybrid distance geometry-simulated annealing using a combination of experimental NMR restraints obtained from the previous refinement of the inhibitor-free MMP-1 (1) and recent restraints for the MMP-1:CGS-27023A complex. The hydroxamic acid moiety of CGS-27023A was found to chelate to the "right" of the catalytic zinc where the p-methoxyphenyl sits in the S1' active-site pocket, the isopropyl group is in contact with H83 and N80, and the pyridine ring is solvent exposed. The sulfonyl oxygens are in hydrogen-bonding distance to the backbone NHs of L81 and A82. This is similar to the conformation determined by NMR of the inhibitor bound to stromelysin (2, 3). A total of 48 distance restraints were observed between MMP-1 and CGS-27023A from 3D 13C-edited/12C-filtered NOESY and 3D 15N-edited NOESY experiments. An additional 18 intramolecular restraints were observed for CGS-27023A from a 2D 12C-filtered NOESY experiment. A minimal set of NMR experiments in combination with the free MMP-1 assignments were used to assign the MMP-1 (1)H, 13C, and 15N resonances in the MMP-1:CGS-27023A complex. The assignments of CGS-27023A in the complex were obtained from 2D 12C-filtered NOESY and 2D 12C-filtered TOCSY experiments.

New general approach for determining the solution structure of a ligand bound weakly to a receptor: structure of a fibrinogen Aalpha-like peptide bound to thrombin (S195A) obtained using NOE distance constraints and an ECEPP/3 flexible docking program

A new approach incorporating flexible docking simulations and NMR data is presented for calculating the bound conformation of a ligand that interacts weakly with an enzyme. This approach consists of sampling directly the conformation of a flexible ligand inside a receptor active site containing surrounding flexible loops. To make this sampling efficient, a ligand-growing procedure has been adopted. Optimization of the ECEPP/3-plus-NOE constraint function is carried out by using a collective variable Monte Carlo minimization technique. Numerous energy minimizations are made possible for such a large system by using a Bezier splines energy grid technique. This new flexible docking approach was applied to determine the structure of a fibrinogen Aalpha-like peptide (7DFLAEGGGVRGPRV20) bound to an active site mutant of thrombin [thrombin(S195A)]. Structure calculations of the bound ligand, using 2D-transferred NOESY distance constraints in the DIANA program, showed that the N-terminal portion of the peptide (D7-R16) involves a chain reversal, whereas the C-terminal portion (G17-V20) adopts a fold that exists in several different orientations. In addition, the ECEPP/3 flexible docking package was used to assess the conformational variability of the ligand and surrounding 60D-insertion loop of thrombin. Amino acid residues (17-20) of the peptide interact with a region of the enzyme that exhibits broad specificity, with a preferred direction between the 60D-insertion loop and Pro37 of thrombin.

NMR study on the binding of d(GGAAATTTCC)2 with a positively charged pentacosapeptide

To obtain a better understanding of the electrostatic nature of protein-nucleic acid interactions, we have investigated the interaction of a double-stranded decamer d(GGAAATTTCC)2 with a synthetic arginine and lysine-rich pentacosapeptide (Pep25), using NMR and optical spectroscopy. The chemical shift data of the decamer under various experimental conditions show that the binding of Pep25 changes the conformation of the decamer in a different way, as compared to the conformational changes induced by a variation in temperature or ionic strength. The chemical shift results are interpreted in terms of ring current effects that emerge into a model for the conformational change, in which the double-stranded helix of the decamer undergoes a decrease of twist and rise to accommodate Pep25. The binding results indicate that the positively charged arginine and lysine side chains of Pep25 not only have a stabilising electrostatic interaction with the negatively charged backbone phosphates of d(GGAAATTTCC)2, but also that a stabilisation of the base pairs of d(GGAAATTTCC)2 by Pep25 takes place.

Binding of peptides in solution by the Escherichia coli chaperone PapD as revealed using an inhibition ELISA and NMR spectroscopy

PapD is the prototype member of a family of periplasmic chaperones which are required for assembly of virulence associated pili in pathogenic, gram-negative bacteria. In the present investigation, an ELISA has been developed for evaluation of compounds as inhibitors of PapD. Synthetic peptides, including an octamer, derived from the C-terminus of the pilus adhesin PapG were able to inhibit PapD in the ELISA. Evaluation of a panel of octapeptides in the ELISA, in combination with NMR studies, showed that the peptides were bound as extended beta-strands by PapD in aqueous solution. The PapD-peptide complex was stabilized by backbone to backbone hydrogen bonds and interactions involving three hydrophobic peptide side chains. This structural information, together with previous crystal structure data, provides a starting point in efforts to design and synthesize compounds which bind to chaperones and interfere with pilus assembly in pathogenic bacteria.

New methods of structure refinement for macromolecular structure determination by NMR

Recent advances in multidimensional NMR methodology have permitted solution structures of proteins in excess of 250 residues to be solved. In this paper, we discuss several methods of structure refinement that promise to increase the accuracy of macromolecular structures determined by NMR. These methods include the use of a conformational database potential and direct refinement against three-bond coupling constants, secondary 13C shifts, 1H shifts, T1/T2 ratios, and residual dipolar couplings. The latter two measurements provide long range restraints that are not accessible by other solution NMR parameters.

The solution structure of a fungal AREA protein-DNA complex: an alternative binding mode for the basic carboxyl tail of GATA factors

The solution structure of a complex between the DNA binding domain of a fungal GATA factor and a 13 base-pair oligonucleotide containing its physiologically relevant CGATAG target sequence has been determined by multidimensional nuclear magnetic resonance spectroscopy. The AREA DNA binding domain, from Aspergillus nidulans, possesses a single Cys2-Cys2 zinc finger module and a basic C-terminal tail, which recognize the CGATAG element via an extensive network of hydrophobic interactions with the bases in the major groove and numerous non-specific contacts along the sugar-phosphate backbone. The zinc finger core of the AREA DNA binding domain has the same global fold as that of the C-terminal DNA binding domain of chicken GATA-1. In contrast to the complex with the DNA binding domain of GATA-1 in which the basic C-terminal tail wraps around the DNA and lies in the minor groove, the structure of complex with the AREA DNA binding domain reveals that the C-terminal tail of the fungal domain runs parallel with the sugar phosphate backbone along the edge of the minor groove. This difference is principally attributed to amino acid substitutions at two positions of the AREA DNA binding domain (Val55, Asn62) relative to that of GATA-1 (Gly55, Lys62). The impact of the different C-terminal tail binding modes on the affinity and specificity of GATA factors is discussed.

The solution structure of the Leu22-->Val mutant AREA DNA binding domain complexed with a TGATAG core element defines a role for hydrophobic packing in the determination of specificity

The seemingly innocuous leucine-to-valine mutation at position 22 of the AREA DNA binding domain results in dramatic changes in the in vivo expression profile of genes controlled by this GATA transcription factor. This is associated with a preference of the Leu22-->Val mutant for TGATAG sites over (A/C)GATAG sites. Quantitative gel retardation assays confirm this observation and show that the Leu22-->Val mutant AREA DNA binding domain has a approximately 30-fold lower affinity than the wild-type domain for a 13 base-pair oligonucleotide containing the wild-type CGATAG target. To gain insight into the measured affinity data and further explore sequence specificity of the AREA protein, the solution structure of a complex between the Leu22-->Val mutant AREA DNA binding domain and a 13 base-pair oligonucleotide containing its physiologically relevant TGATAG target sequence has been determined by multidimensional nuclear magnetic resonance spectroscopy. Comparison of this structure with that of the wild-type AREA DNA binding domain complexed to its cognate CGATAG target site shows how subtle changes in amino acid side-chain length and hydrophobic packing can affect affinity and specificity for GATA-containing sequences, and how changes in DNA sequence can be compensated for by changes in protein sequence.

The helix-hinge-helix structural motif in human apolipoprotein A-I determined by NMR spectroscopy

The conformation of a synthetic peptide of 46 residues from apoA-I was investigated by fluorescence, CD, and 2D NMR spectroscopies in lipid-mimetic environments. ApoA-I(142-187) is mainly unstructured in water but helical in SDS or dodecylphosphocholine (DPC), although the peptide only associates with DPC at approximately the critical micellar concentration. Solution structures of apoA-I(142-187) were determined by distance geometry calculations based on 450 (in DPC-d38) or 397 (in SDS-d25) NOE-derived distance restraints, respectively. Backbone RMSDs for superimposing the two helical regions 146-162 and 168-182 are 0.98 +/- 0.22 (2.38 +/- 0.20) and 1.99 +/- 0.42 (2.02 +/- 0.21) A in DPC (SDS), respectively. No interhelical NOE was found, suggesting that helix-helix interactions between the two helical domains in apoA-I(142-187) are unlikely. Similar average, curved helix-hinge-helix structures were found in both SDS and DPC micelles with the hydrophobic residues occupying the concave face, indicating that hydrophobic interactions dominate. Intermolecular NOESY experiments, performed in the presence of 50% protonated SDS, confirm that the two amphipathic helices and Y166 in the hinge all interact with the micelle. The involvement of Y166 in lipid binding is supported by fluorescence spectroscopy as well. On the basis of all the data above, we propose a model for the peptide-lipid complexes wherein the curved amphipathic helix-hinge-helix structural motif straddles the micelle. The peptide-aided signal assignment achieved for apoA-I(122-187) (66mer) and apoA-I suggests that such a structural motif is retained in the longer peptide and most likely in the intact protein.

Protein NMR extends into new fields of structural biology

Advances in protein NMR have opened new doors for the understanding of macromolecular structure and interactions. Isotope-labelling approaches have extended the size limit for structure determination, the mapping of protein-ligand interactions is now widely used and has led to a new drug discovery approach, pathways of protein folding processes can be followed, and the structures of molecules from membranes can be determined.

Determination of the secondary structure and global topology of the 44 kDa ectodomain of gp41 of the simian immunodeficiency virus by multidimensional nuclear magnetic resonance spectroscopy

The gp41 protein of the human (HIV) and simian (SIV) immunodeficiency viruses is part of the envelope glycoprotein complex gp41/gp120 which plays an essential role in viral infection. We present a multidimensional NMR study on the trimeric 44 kDa soluble ectodomain of SIV gp41 (e-gp41), comprising residues 27 to 149. Despite the large molecular weight and very limited spectral dispersion, complete backbone 1H, 13C, 13CO and 15N assignments have been made using a combination of triple resonance experiments on uniformly 13C/15N and 2H/13C/15N-labeled samples. The secondary structure of SIV e-gp41, derived on the basis of 13C chemical shifts, NH exchange rates, medium range nuclear Overhauser enhancements (NOE), and 3JHNalpha coupling constants, consists of a 49 residue helix at the N terminus (residues 29 to 77) and a 40 residue helix at the C terminus (residues 108 to 147), connected by a 30 residue loop which does not display any of the characteristics of regular secondary structure. The cross-peak intensities of the loop region in scalar correlation experiments suggests that it is more mobile than the core helical regions. The presence, however, of numerous long range NOEs, both intra and inter-subunit, within the loop indicates that it adopts a well-defined structure in which the loops from the three subunits interact with each other. Based on a number of long range intra and inter-subunit NOEs, a topological model is presented for the symmetric SIV e-gp41 trimer in which the N-terminal helices are packed within the protein interior in a parallel trimeric coiled-coil arrangement, while the C-terminal helices are located on the protein exterior, oriented antiparallel to the N-terminal helices.

The solution structure of an HMG-I(Y)-DNA complex defines a new architectural minor groove binding motif

The solution structure of a complex between a truncated form of HMG-I(Y), consisting of the second and third DNA binding domains (residues 51-90), and a DNA dodecamer containing the PRDII site of the interferon-beta promoter has been solved by multidimensional nuclear magnetic resonance spectroscopy. The stoichiometry of the complex is one molecule of HMG-I(Y) to two molecules of DNA. The structure reveals a new architectural minor groove binding motif which stabilizes B-DNA, thereby facilitating the binding of other transcription factors in the opposing major groove. The interactions involve a central Arg-Gly-Arg motif together with two other modules that participate in extensive hydrophobic and polar contracts. The absence of one of these modules in the third DNA binding domain accounts for its-100 fold reduced affinity relative to the second one.

Solution structure of the N-terminal zinc binding domain of HIV-1 integrase

The solution structure of the N-terminal zinc binding domain (residues 1-55; IN1-55) of HIV-1 integrase has been solved by NMR spectroscopy. IN1-55 is dimeric, and each monomer comprises four helices with the zinc tetrahedrally coordinated to His 12, His 16, Cys 40 and Cys 43. IN1-55 exists in two interconverting conformational states that differ with regard to the coordination of the two histidine side chains to zinc. The different histidine arrangements are associated with large conformational differences in the polypeptide backbone (residues 9-18) around the coordinating histidines. The dimer interface is predominantly hydrophobic and is formed by the packing of the N-terminal end of helix 1, and helices 3 and 4. The monomer fold is remarkably similar to that of a number of helical DNA binding proteins containing a helix-turn-helix (HTH) motif with helices 2 and 3 of IN1-55 corresponding to the HTH motif. In contrast to the DNA binding proteins where the second helix of the HTH motif is employed for DNA recognition, IN1-55 uses this helix for dimerization.

NMR spectroscopy of peptides and proteins. Practical considerations

High resolution nuclear magnetic resonance (NMR) spectroscopy is the only method available for determining the three-dimensional structures of peptides and proteins in solution at atomic resolution. This article deals with a range of practical considerations associated with such studies, including sample preparation, instrumental setup, one- and two-dimensional NMR methods, interpretation of spectral data, and structure calculations.

Leptin is a four-helix bundle: secondary structure by NMR

Leptin is a signaling protein that in its mutant forms has been associated with obesity and Type II diabetes. The lack of sequence similarity has precluded analogies based on structural resemblance to known systems. Backbone NMR signals for mouse leptin (13C/15N -labeled) have been assigned and its secondary structure reveals it to be a four-helix bundle cytokine. Helix lengths and disulfide pattern are in agreement with leptin as a member of the short-helix cytokine family. A three-dimensional model was built verifying the mechanical consistency of the identified elements with a short-helix cytokine core.

The solution structure of a specific GAGA factor-DNA complex reveals a modular binding mode

The structure of a complex between the DNA binding domain of the GAGA factor (GAGA-DBD) and an oligonucleotide containing its GAGAG consensus binding site has been determined by nuclear magnetic resonance spectroscopy. The GAGA-DBD comprises a single classical Cys2-His2 zinc finger core, and an N-terminal extension containing two highly basic regions, BR1 and BR2. The zinc finger core binds in the major groove and recognizes the first three GAG bases of the consensus in a manner similar to that seen in other classical zinc finger-DNA complexes. Unlike the latter, which require tandem zinc finger repeats with a minimum of two units for high affinity binding, the GAGA-DBD makes use of only a single finger complemented by BR1 and BR2. BR2 forms a helix that interacts in the major groove recognizing the last G of the consensus, while BR1 wraps around the DNA in the minor groove and recognizes the A in the fourth position of the consensus. The implications of the structure of the GAGA-DBD-DNA complex for chromatin remodelling are discussed.

Conformations of human apolipoprotein E(263-286) and E(267-289) in aqueous solutions of sodium dodecyl sulfate by CD and 1H NMR

Structures of apoE(263-286) and apoE(267-289) have been determined in aqueous solution containing 90-fold molar excess of perdeuterated sodium dodecyl sulfate by CD and 1H NMR. Conformations were calculated by distance geometry based on 370 and 276 NOE distance restraints, respectively. RMSD for superimposing the region 265-284 from an ensemble of 41 structures for apoE(263-286) is 0.64 +/- 0.17 A for backbone atoms (N, C alpha, C = O) and 1.51 +/- 0.13 A for all atoms. The backbone RMSD for an ensemble of 37 structures for apoE(267-289) is 0.74 +/- 0.21 A for the region 268-275 and 0.34 +/- 0.10 A for the region 276-286. A two-domain structure was found for apoE(267-289) with the C-terminal half adopting a very well defined helix and the N-terminal segment 268-275 a less well defined helix, suggesting that the N-terminus may weakly bind to SDS. For apoE(263-286), an amphipathic helix-bend-helix structural motif was found with all hydrophobic side chains on the concave face. The existence of a bend around residues Q273 to G278 is consistent with their temperature coefficients of amide protons as well as secondary shifts of alpha-protons. Comparison of the structures of the two peptides revealed that the enhanced binding of apoE(263-286) to lipid could be attributed to the formation of a hydrophobic cluster consisting of residues W264, F265, L268, and V269. Aromatic side chains are proposed to be especially important in anchoring apolipoprotein fragments to micelles.

Improving the quality of NMR and crystallographic protein structures by means of a conformational database potential derived from structure databases

A new conformational database potential involving dihedral angle relationships in databases of high-resolution highly refined protein crystal structures is presented as a method for improving the quality of structures generated from NMR data. The rationale for this procedure is based on the observation that uncertainties in the description of the nonbonded contacts present a key limiting factor in the attainable accuracy of protein NMR structures and that the nonbonded interaction terms presently used have poor discriminatory power between high- and low-probability local conformations. The idea behind the conformational database potential is to restrict sampling during simulated annealing refinement to conformations that are likely to be energetically possible by effectively limiting the choices of dihedral angles to those that are known to be physically realizable. In this manner, the variability in the structures produced by this method is primarily a function of the experimental restraints, rather than an artifact of a poor nonbonded interaction model. We tested this approach with the experimental NMR data (comprising an average of about 30 restraints per residue and consisting of interproton distances, torsion angles, 3JHN alpha coupling constants, and 13C chemical shifts) used previously to calculate the solution structure of reduced human thioredoxin (Qin J, Clore GM, Gronenborn AM, 1994, Structure 2:503-522). Incorporation of the conformational database potential into the target function used for refinement (which also includes terms for the experimental restraints, covalent geometry, and nonbonded interactions in the form of either a repulsive, repulsive-attractive, or 6-12 Lennard-Jones potential) results in a significant improvement in various quantitative measures of quality (Ramachandran plot, side-chain torsion angles, overall packing). This is achieved without compromising the agreement with the experimental restraints and the deviations from idealized covalent geometry that remain within experimental error, and the agreement between calculated and observed 1H chemical shifts that provides an independent NMR parameter of accuracy. The method is equally applicable to crystallographic refinement, and should be particular useful during the early stages of either an NMR or crystallographic structure determination and in cases where relatively few experimental restraints can be derived from the measured data (due, for example, to broad lines in the NMR spectra or to poorly diffracting crystals).