Determination of the secondary structure of the DNA binding protein Ner from phage Mu using 1H homonuclear and 15N-1H heteronuclear NMR spectroscopy

The sequential resonance assignment of the 1H and 15N NMR spectra of the DNA binding protein Ner from phage Mu is presented. This is carried out by using a combination of 1H-1H and 1H-15N two-dimensional experiments. The availability of completely labeled 15N protein enabled us to record a variety of relayed heteronuclear multiple quantum coherence experiments, thereby enabling the correlation of proton-proton through-space and through-bond connectivities with the chemical shift of the directly bonded 15N atom. These heteronuclear experiments were crucial for the sequential assignment as the proton chemical shift dispersion of the Ner protein is limited and substantial overlap precluded unambiguous assignment of the homonuclear spectra in several cases. From a qualitative interpretation of the NOE data involving the NH, C alpha H, and C beta H protons, it is shown that Ner is composed of five helices extending from residues 11 to 22, 27 to 34, 38 to 45, 50 to 60, and 63 to 73.     

15N NMR assignments and chemical shift analysis of uniformly labeled 15N calbindin D9k in the apo, (Cd2+)1 and (Ca2+)2 states.

15N has been uniformly incorporated into the EF-hand Ca(2+)-binding protein calbindin D9k so that heteronuclear experiments can be used to further characterize the structure and dynamics of the apo, (Cd2+)1 and (Ca2+)2 states of the protein. The 15N NMR resonances were assigned by 2D 15N-resolved 1H experiments, which also allowed the identification of a number of sequential and medium-range 1H-1H contacts that are obscured by chemical shift degeneracy in homonuclear experiments. The 15N chemical shifts are analyzed with respect to correlations with protein secondary structure. In addition, the changes in 15N chemical shift found for the apo----(Cd2+)1----(Ca2+)2 binding sequence confirm that the effects on the protein are mainly associated with chelation of the first ion.     

NMR: prediction of protein flexibility

We present a protocol for predicting protein flexibility from NMR chemical shifts. The protocol consists of (i) ensuring that the chemical shift assignments are correctly referenced or, if not, performing a reference correction using information derived from the chemical shift index, (ii) calculating the random coil index (RCI), and (iii) predicting the expected root mean square fluctuations (RMSFs) and order parameters (S2) of the protein from the RCI. The key advantages of this protocol over existing methods for studying protein dynamics are that (i) it does not require prior knowledge of a protein's tertiary structure, (ii) it is not sensitive to the protein's overall tumbling and (iii) it does not require additional NMR measurements beyond the standard experiments for backbone assignments. When chemical shift assignments are available, protein flexibility parameters, such as S2 and RMSF, can be calculated within 1-2 h using a spreadsheet program.     

Demonstration of protein-protein interaction specificity by NMR chemical shift mapping.

Chemical shift mapping is becoming a popular method for studying protein-protein interactions in solution. The technique is used to identify putative sites of interaction on a protein surface by detecting chemical shift perturbations in simple (1H, 15N)-HSQC NMR spectra of a uniformly labeled protein as a function of added (unlabeled) target protein. The high concentrations required for these experiments raise questions concerning the possibility for non-specific interactions being detected, thereby compromising the information obtained. We demonstrate here that the simple chemical shift mapping approach faithfully reproduces the known functional specificities among pairs of closely related proteins from the phosphoenolpyruvate:sugar phosphotransferase systems of Escherichia coli and Bacillus subtilis.     

Prediction of peak overlap in NMR spectra

Peak overlap is one of the major factors complicating the analysis of biomolecular NMR spectra. We present a general method for predicting the extent of peak overlap in multidimensional NMR spectra and its validation using both, experimental data sets and Monte Carlo simulation. The method is based on knowledge of the magnetization transfer pathways of the NMR experiments and chemical shift statistics from the Biological Magnetic Resonance Data Bank. Assuming a normal distribution with characteristic mean value and standard deviation for the chemical shift of each observable atom, an analytic expression was derived for the expected overlap probability of the cross peaks. The analytical approach was verified to agree with the average peak overlap in a large number of individual peak lists simulated using the same chemical shift statistics. The method was applied to eight proteins, including an intrinsically disordered one, for which the prediction results could be compared with the actual overlap based on the experimentally measured chemical shifts. The extent of overlap predicted using only statistical chemical shift information was in good agreement with the overlap that was observed when the measured shifts were used in the virtual spectrum, except for the intrinsically disordered protein. Since the spectral complexity of a protein NMR spectrum is a crucial factor for protein structure determination, analytical overlap prediction can be used to identify potentially difficult proteins before conducting NMR experiments. Overlap predictions can be tailored to particular classes of proteins by preparing statistics from corresponding protein databases. The method is also suitable for optimizing recording parameters and labeling schemes for NMR experiments and improving the reliability of automated spectra analysis and protein structure determination.     

SimShift: identifying structural similarities from NMR chemical shifts

MOTIVATION: An important quantity that arises in NMR spectroscopy experiments is the chemical shift. The interpretation of these data is mostly done by human experts; to our knowledge there are no algorithms that predict protein structure from chemical shift sequences alone. One approach to facilitate this process could be to compare two such sequences, where the structure of one protein has already been resolved. Our claim is that similarity of chemical shifts thereby found implies structural similarity of the respective proteins. RESULTS: We present an algorithm to identify structural similarities of proteins by aligning their associated chemical shift sequences. To evaluate the correctness of our predictions, we propose a benchmark set of protein pairs that have high structural similarity, but low sequence similarity (because with high sequence similarity the structural similarities could easily be detected by a sequence alignment algorithm). We compare our results with those of HHsearch and SSEA and show that our method outperforms both in >50% of all cases.     

Resolution enhancement by homonuclear J-decoupling: application to three-dimensional solid-state magic angle spinning NMR spectroscopy

We describe a simple protocol to achieve homonuclear J-decoupling in the indirect dimensions of multidimensional experiments, and to enhance spectral resolution of the backbone Calpha carbons in the 3D NCACX experiment. In the proposed protocol, the refocusing of the Calpha-CO homonuclear J-couplings is achieved by applying an off-resonance selective pi pulse to the CO spectral region in the middle of Calpha chemical shift evolution. As is commonly used in solution NMR, a compensatory echo period is used to refocus the unwanted chemical shift evolution of Calpha spins, which takes place during the off-resonance selective pulse. The experiments were carried out on the beta1 immunoglobulin binding domain of protein G (GB1). In GB1, such implementation results in significantly reduced line widths, and leads to an overall sensitivity enhancement.     

Auto-FACE: An NMR Based Binding Site Mapping Program for Fast Chemical Exchange Protein-Ligand Systems

Background

Nuclear Magnetic Resonance (NMR) spectroscopy offers a variety of experiments to study protein-ligand interactions at atomic resolution. Among these experiments, N Heteronuclear Single Quantum Correlation (HSQC) experiment is simple, less time consuming and highly informative in mapping the binding site of the ligand. The interpretation of N HSQC becomes ambiguous when the chemical shift perturbations are caused by non-specific interactions like allosteric changes and local structural rearrangement. Under such cases, detailed chemical exchange analysis based on chemical shift perturbation will assist in locating the binding site accurately.

Methodology/Principal Findings

We have automated the mapping of binding sites for fast chemical exchange systems using information obtained from N HSQC spectra of protein serially titrated with ligand of increasing concentrations. The automated program Auto-FACE (Auto-FAst Chemical Exchange analyzer) determines the parameters, e.g. rate of change of perturbation, binding equilibrium constant and magnitude of chemical shift perturbation to map the binding site residues. Interestingly, the rate of change of perturbation at lower ligand concentration is highly sensitive in differentiating the binding site residues from the non-binding site residues. To validate this program, the interaction between the protein and the ligand BH3I-1 was studied. Residues in the hydrophobic BH3 binding groove of were easily identified to be crucial for interaction with BH3I-1 from other residues that also exhibited perturbation. The geometrically averaged equilibrium constant () calculated for the residues present at the identified binding site is consistent with the values obtained by other techniques like isothermal calorimetry and fluorescence polarization assays (). Adjacent to the primary site, an additional binding site was identified which had an affinity of 3.8 times weaker than the former one. Further NMR based model fitting for individual residues suggest single site model for residues present at these binding sites and two site model for residues present between these sites. This implies that chemical shift perturbation can represent the local binding event much more accurately than the global binding event.

Conclusion/Significance

Detail NMR chemical shift perturbation analysis enabled binding site residues to be distinguished from non-binding site residues for accurate mapping of interaction site in complex fast exchange system between small molecule and protein. The methodology is automated and implemented in a program called “Auto-FACE”, which also allowed quantitative information of each interaction site and elucidation of binding mechanism.     

Four-dimensional heteronuclear correlation experiments for chemical shift assignment of solid proteins

Chemical shift assignment is the first step in all established protocols for structure determination of uniformly labeled proteins by NMR. The explosive growth in recent years of magic-angle spinning (MAS) solid-state NMR (SSNMR) applications is largely attributable to improved methods for backbone and side-chain chemical shift correlation spectroscopy. However, the techniques developed so far have been applied primarily to proteins in the size range of 5–10 kDa, despite the fact that SSNMR has no inherent molecular weight limits. Rather, the degeneracy inherent to many 2D and 3D SSNMR spectra of larger proteins has prevented complete unambiguous chemical shift assignment. Here we demonstrate the implementation of 4D backbone chemical shift correlation experiments for assignment of solid proteins. The experiments greatly reduce spectral degeneracy at a modest cost in sensitivity, which is accurately described by theory. We consider several possible implementations and investigate the CANCOCX pulse sequence in detail. This experiment involves three cross polarization steps, from H to CA[i], CA[i] to N[i], and N[i] to C′[i−1], followed by a final homonuclear mixing period. With short homonuclear mixing times (<20 ms), backbone correlations are observed with high sensitivity; with longer mixing times (>200 ms), long-range correlations are revealed. For example, a single 4D experiment with 225 ms homonuclear mixing time reveals ∼200 uniquely resolved medium and long-range correlations in the 56-residue protein GB1. In addition to experimental demonstrations in the 56-residue protein GB1, we present a theoretical analysis of anticipated improvements in resolution for much larger proteins and compare these results in detail with the experiments, finding good agreement between experiment and theory under conditions of stable instrumental performance.     

The interaction between calcium- and integrin-binding protein 1 and the alphaIIb integrin cytoplasmic domain involves a novel C-terminal displacement mechanism

Calcium- and integrin-binding protein 1 (CIB1) regulates platelet aggregation in hemostasis through a specific interaction with the alphaIIb cytoplasmic domain of platelet integrin alphaIIbbeta3. In this work we report the structural characteristics of CIB1 in solution and the mechanistic details of its interaction with a synthetic peptide derived from the alphaIIb cytoplasmic domain. NMR spectroscopy experiments using perdeuterated CIB1 together with heteronuclear nuclear Overhauser effect experiments have revealed a well folded alpha-helical structure for both the ligand-free and alphaIIb-bound forms of the protein. Residual dipolar coupling experiments have shown that the N and C domains of CIB1 are positioned side by side, and chemical shift perturbation mapping has identified the alphaIIb-binding site as a hydrophobic channel spanning the entire C domain and part of the N domain. Data obtained with a truncated version of CIB1 suggest that the extreme C-terminal end of the protein weakly interacts with this channel in the absence of a biological target, but it is displaced by the alphaIIb cytoplasmic domain, suggesting a novel mechanism to increase binding specificity.     

3d Trosy-HncaCodedcb and Trosy-HncaCodedco Experiments: Triple Resonance nmr Experiments With two Sequential Connectivity Pathways and High Sensitivity

The concept of chemical shift-coding monitors chemical shifts in multi-dimensional NMR experiments without additional polarization transfer elements and evolution periods. The chemical shifts are coded in the line-shape of the cross-peak through an apparent scalar coupling dependent upon the chemical shift. This concept is applied to the three-dimensional triple-resonance experiment HNCA adding the information of (13)C(beta) or (13)C' chemical shifts. On average, the proposed TROSY-HNCA(coded)CB experiment is a factor of 2 less sensitive than the HNCA experiment. However, it contains correlations via the chemical shifts of both (13)C(alpha) and (13)C(beta), and provides up to three times higher resolution along the (13)C(alpha) chemical shift axis. Thus, it dramatically reduces ambiguities in linking the spin systems of adjacent residues in the protein sequence during the sequential assignment. The TROSY-HNCA(coded)CO experiment which is conceptually similar contains correlations via the chemical shifts of (13)C(alpha) and (13)C' without major signal losses. The proposed triple resonance experiments are applied to a approximately 70% (2)H, approximately 85% (13)C,(15)N labeled protein with a molecular weight of 25 kDa.     

The J-UNIO protocol for automated protein structure determination by NMR in solution

The J-UNIO (JCSG protocol using the software UNIO) procedure for automated protein structure determination by NMR in solution is introduced. In the present implementation, J-UNIO makes use of APSY-NMR spectroscopy, 3D heteronuclear-resolved [(1)H,(1)H]-NOESY experiments, and the software UNIO. Applications with proteins from the JCSG target list with sizes up to 150 residues showed that the procedure is highly robust and efficient. In all instances the correct polypeptide fold was obtained in the first round of automated data analysis and structure calculation. After interactive validation of the data obtained from the automated routine, the quality of the final structures was comparable to results from interactive structure determination. Special advantages are that the NMR data have been recorded with 6-10 days of instrument time per protein, that there is only a single step of chemical shift adjustments to relate the backbone signals in the APSY-NMR spectra with the corresponding backbone signals in the NOESY spectra, and that the NOE-based amino acid side chain chemical shift assignments are automatically focused on those residues that are heavily weighted in the structure calculation. The individual working steps of J-UNIO are illustrated with the structure determination of the protein YP_926445.1 from Shewanella amazonensis, and the results obtained with 17 JCSG targets are critically evaluated.     

Relationship between nuclear magnetic resonance chemical shift and protein secondary structure

An analysis of the 1H nuclear magnetic resonance chemical shift assignments and secondary structure designations for over 70 proteins has revealed some very strong and unexpected relationships. Similar studies, performed on smaller databases, for 13C and 15N chemical shifts reveal equally strong correlations to protein secondary structure. Among the more interesting results to emerge from this work is the finding that all 20 naturally occurring amino acids experience a mean alpha-1H upfield shift of 0.39 parts per million (from the random coil value) when placed in a helical configuration. In a like manner, the alpha-1H chemical shift is found to move downfield by an average of 0.37 parts per million when the residue is placed in a beta-strand or extended configuration. Similar changes are also found for amide 1H, carbonyl 13C, alpha-13C and amide 15N chemical shifts. Other relationships between chemical shift and protein conformation are also uncovered; in particular, a correlation between helix dipole effects and amide proton chemical shifts as well as a relationship between alpha-proton chemical shifts and main-chain flexibility. Additionally, useful relationships between alpha-proton chemical shifts and backbone dihedral angles as well as correlations between amide proton chemical shifts and hydrogen bond effects are demonstrated.     

NMR assignments of a hypothetical pseudo-knotted protein HP0242 from Helicobacter pylori

Many knotted proteins have been discovered recently, but the folding process of which remains elusive. HP0242 is a hypothetical protein from Helicobacter pylori, which is a model system for studying the folding pathway of a knotted protein. In this study, we report the ¹H, ¹³C, and ¹⁵N chemical shift assignments of HP0242. The results will enable us to further investigate HP0242 by NMR experiments.     

How do chemical denaturants affect the mechanical folding and unfolding of proteins?

We present the first single-molecule atomic force microscopy study on the effect of chemical denaturants on the mechanical folding/unfolding kinetics of a small protein GB1 (the B1 immunoglobulin-binding domain of protein G from Streptococcus). Upon increasing the concentration of the chemical denaturant guanidinium chloride (GdmCl), we observed a systematic decrease in the mechanical stability of GB1, indicating the softening effect of the chemical denaturant on the mechanical stability of proteins. This mechanical softening effect originates from the reduced free-energy barrier between the folded state and the unfolding transition state, which decreases linearly as a function of the denaturant concentration. Chemical denaturants, however, do not alter the mechanical unfolding pathway or shift the position of the transition state for mechanical unfolding. We also found that the folding rate constant of GB1 is slowed down by GdmCl in mechanical folding experiments. By combining the mechanical folding/unfolding kinetics of GB1 in GdmCl solution, we developed the “mechanical chevron plot” as a general tool to understand how chemical denaturants influence the mechanical folding/unfolding kinetics and free-energy diagram in a quantitative fashion. This study demonstrates great potential in combining chemical denaturation with single-molecule atomic force microscopy techniques to reveal invaluable information on the energy landscape underlying protein folding/unfolding reactions.     

Sugar-to-base correlation in nucleic acids with a 5D APSY-HCNCH or two 3D APSY-HCN experiments

A five-dimensional (5D) APSY (automated projection spectroscopy) HCNCH experiment is presented, which allows unambiguous correlation of sugar to base nuclei in nucleic acids. The pulse sequence uses multiple quantum (MQ) evolution which enables long constant-time evolution periods in all dimensions, an improvement that can also benefit non-APSY applications. Applied with an RNA with 23 nucleotides the 5D APSY-HCNCH experiment produced a complete and highly precise 5D chemical shift list within 1.5 h. Alternatively, and for molecules where the out-and-stay 5D experiment sensitivity is not sufficient, a set of out-and-back 3D APSY-HCN experiments is proposed: an intra-base (3D APSY-b-HCN) experiment in an MQ or in a TROSY version, and an MQ sugar-to-base (3D APSY-s-HCN) experiment. The two 3D peak lists require subsequent matching via the N1/9 chemical shift values to one 5D peak list. Optimization of the 3D APSY experiments for maximal precision in the N1/9 dimension allowed matching of all ¹⁵N chemical shift values contained in both 3D peak lists. The precise 5D chemical shift correlation lists resulting from the 5D experiment or a pair of 3D experiments also provide a valuable basis for subsequent connection to chemical shifts derived with other experiments.     

<Superscript>1</Superscript>H, <Superscript>15</Superscript>N, <Superscript>13</Superscript>C resonance assignment of plant dehydrin early response to dehydration 10 (ERD10)

Early response to dehydration 10 protein (ERD10) is an intrinsically disordered protein from Arabidopsis thaliana. The protein is upregulated during stress however its mechanism of action at atomic level is not well understood. In the present work multidimensional NMR methodologies are used in order to facilitate the process of chemical shift assignment. The information provided here supports further NMR spectroscopy experiments aimed at elucidation of ERD10 behaviour during molecular recognition events with other proteins.     

A mammalian type I fatty acid synthase acyl carrier protein domain does not sequester acyl chains

The synthases that produce fatty acids in mammals (FASs) are arranged as large multidomain polypeptides. The growing fatty acid chain is bound covalently during chain elongation and reduction to the acyl carrier protein (ACP) domain that is then able to access each catalytic site. In this work we report the high-resolution nuclear magnetic resonance (NMR) solution structure of the isolated rat fatty acid synthase apoACP domain. The final ensemble of NMR structures and backbone (15)N relaxation studies show that apoACP adopts a single, well defined fold. On conversion to the holo form, several small chemical shift changes are observed on the ACP for residues surrounding the phosphopantetheine attachment site (as monitored by backbone (1)H-(15)N correlation experiments). However, there are negligible chemical shift changes when the holo form is modified to either the hexanoyl or palmitoyl forms. For further NMR analysis, a (13)C,(15)N-labeled hexanoyl-ACP sample was prepared and full chemical shift assignments completed. Analysis of two-dimensional F(2)-filtered and three-dimensional (13)C-edited nuclear Overhauser effect spectroscopy experiments revealed no detectable NOEs to the acyl chain. These experiments demonstrate that unlike other FAS ACPs studied, this Type I ACP does not sequester a covalently linked acyl moiety, although transient interactions cannot be ruled out. This is an important mechanistic difference between the ACPs from Type I and Type II FASs and may be significant for the modulation and regulation of these important mega-synthases.     

Automatic <Superscript>13</Superscript>C chemical shift reference correction for unassigned protein NMR spectra

Poor chemical shift referencing, especially for ¹³C in protein Nuclear Magnetic Resonance (NMR) experiments, fundamentally limits and even prevents effective study of biomacromolecules via NMR, including protein structure determination and analysis of protein dynamics. To solve this problem, we constructed a Bayesian probabilistic framework that circumvents the limitations of previous reference correction methods that required protein resonance assignment and/or three-dimensional protein structure. Our algorithm named Bayesian Model Optimized Reference Correction (BaMORC) can detect and correct ¹³C chemical shift referencing errors before the protein resonance assignment step of analysis and without three-dimensional structure. By combining the BaMORC methodology with a new intra-peaklist grouping algorithm, we created a combined method called Unassigned BaMORC that utilizes only unassigned experimental peak lists and the amino acid sequence. Unassigned BaMORC kept all experimental three-dimensional HN(CO)CACB-type peak lists tested within ±?0.4 ppm of the correct ¹³C reference value. On a much larger unassigned chemical shift test set, the base method kept ¹³C chemical shift referencing errors to within ±?0.45 ppm at a 90% confidence interval. With chemical shift assignments, Assigned BaMORC can detect and correct ¹³C chemical shift referencing errors to within ±?0.22 at a 90% confidence interval. Therefore, Unassigned BaMORC can correct ¹³C chemical shift referencing errors when it will have the most impact, right before protein resonance assignment and other downstream analyses are started. After assignment, chemical shift reference correction can be further refined with Assigned BaMORC. These new methods will allow non-NMR experts to detect and correct ¹³C referencing error at critical early data analysis steps, lowering the bar of NMR expertise required for effective protein NMR analysis.     

Letter to the Editor: Backbone 1H, 15N and 13C Assignments for the Subunit a of the E. Coli ATP Synthase

The structure of the 30 KDa subunit a of the membrane component (F(0)) of E. coli ATP synthase is investigated in a mixture of chloroform, methanol and water, a solvent previously used for solving the structure of another integral membrane protein, subunit c. Near complete backbone chemical shift assignments were made from a set of TROSY experiments including HNCO, HNCA, HN(CA)CB, HN(CO)CACB and 4D HNCOCA and HNCACO. Secondary structure of subunit a was predicted from the backbone chemical shifts using TALOS program. The protein was found to consist of multiple elongated alpha-helical segments. This finding is generally consistent with previous predictions of multiple transmembrane alpha-helices in this polytopic protein.