Mapping Membrane Protein Backbone Dynamics: A Comparison of Site-Directed Spin Labeling with NMR N-Relaxation Measurements

The ability to detect nanosecond backbone dynamics with site-directed spin labeling (SDSL) in soluble proteins has been well established. However, for membrane proteins, the nitroxide appears to have more interactions with the protein surface, potentially hindering the sensitivity to backbone motions. To determine whether membrane protein backbone dynamics could be mapped with SDSL, a nitroxide was introduced at 55 independent sites in a model polytopic membrane protein, TM0026. Electron paramagnetic resonance spectral parameters were compared with NMR ¹⁵N-relaxation data. Sequential scans revealed backbone dynamics with the same trends observed for the R₁ relaxation rate, suggesting that nitroxide dynamics remain coupled to the backbone on membrane proteins.     

Mapping Membrane Protein Backbone Dynamics: A Comparison of Site-Directed Spin Labeling with NMR 15N-Relaxation Measurements

The ability to detect nanosecond backbone dynamics with site-directed spin labeling (SDSL) in soluble proteins has been well established. However, for membrane proteins, the nitroxide appears to have more interactions with the protein surface, potentially hindering the sensitivity to backbone motions. To determine whether membrane protein backbone dynamics could be mapped with SDSL, a nitroxide was introduced at 55 independent sites in a model polytopic membrane protein, TM0026. Electron paramagnetic resonance spectral parameters were compared with NMR ¹⁵N-relaxation data. Sequential scans revealed backbone dynamics with the same trends observed for the R₁ relaxation rate, suggesting that nitroxide dynamics remain coupled to the backbone on membrane proteins.     

Retrieving backbone string neighbors provides insights into structural modeling of membrane proteins

Identification of protein structural neighbors to a query is fundamental in structure and function prediction. Here we present BS-align, a systematic method to retrieve backbone string neighbors from primary sequences as templates for protein modeling. The backbone conformation of a protein is represented by the backbone string, as defined in Ramachandran space. The backbone string of a query can be accurately predicted by two innovative technologies: a knowledge-driven sequence alignment and encoding of a backbone string element profile. Then, the predicted backbone string is employed to align against a backbone string database and retrieve a set of backbone string neighbors. The backbone string neighbors were shown to be close to native structures of query proteins. BS-align was successfully employed to predict models of 10 membrane proteins with lengths ranging between 229 and 595 residues, and whose high-resolution structural determinations were difficult to elucidate both by experiment and prediction. The obtained TM-scores and root mean square deviations of the models confirmed that the models based on the backbone string neighbors retrieved by the BS-align were very close to the native membrane structures although the query and the neighbor shared a very low sequence identity. The backbone string system represents a new road for the prediction of protein structure from sequence, and suggests that the similarity of the backbone string would be more informative than describing a protein as belonging to a fold.     

Computational design of high-affinity epitope scaffolds by backbone grafting of a linear epitope

Computational grafting of functional motifs onto scaffold proteins is a promising way to engineer novel proteins with pre-specified functionalities. Typically, protein grafting involves the transplantation of protein side chains from a functional motif onto structurally homologous regions of scaffold proteins. Using this approach, we previously transplanted the human immunodeficiency virus 2F5 and 4E10 epitopes onto heterologous proteins to design novel "epitope-scaffold" antigens. However, side-chain grafting is limited by the availability of scaffolds with compatible backbone for a given epitope structure and offers no route to modify backbone structure to improve mimicry or binding affinity. To address this, we report here a new and more aggressive computational method-backbone grafting of linear motifs-that transplants the backbone and side chains of linear functional motifs onto scaffold proteins. To test this method, we first used side-chain grafting to design new 2F5 epitope scaffolds with improved biophysical characteristics. We then independently transplanted the 2F5 epitope onto three of the same parent scaffolds using the newly developed backbone grafting procedure. Crystal structures of side-chain and backbone grafting designs showed close agreement with both the computational models and the desired epitope structure. In two cases, backbone grafting scaffolds bound antibody 2F5 with 30- and 9-fold higher affinity than corresponding side-chain grafting designs. These results demonstrate that flexible backbone methods for epitope grafting can significantly improve binding affinities over those achieved by fixed backbone methods alone. Backbone grafting of linear motifs is a general method to transplant functional motifs when backbone remodeling of the target scaffold is necessary.     

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.     

Recapitulation of protein family divergence using flexible backbone protein design

We use flexible backbone protein design to explore the sequence and structure neighborhoods of naturally occurring proteins. The method samples sequence and structure space in the vicinity of a known sequence and structure by alternately optimizing the sequence for a fixed protein backbone using rotamer based sequence search, and optimizing the backbone for a fixed amino acid sequence using atomic-resolution structure prediction. We find that such a flexible backbone design method better recapitulates protein family sequence variation than sequence optimization on fixed backbones or randomly perturbed backbone ensembles for ten diverse protein structures. For the SH3 domain, the backbone structure variation in the family is also better recapitulated than in randomly perturbed backbones. The potential application of this method as a model of protein family evolution is highlighted by a concerted transition to the amino acid sequence in the structural core of one SH3 domain starting from the backbone coordinates of an homologous structure.     

Fast Assignment of <Superscript>15</Superscript>N-HSQC Peaks using High-Resolution 3D HNcocaNH Experiments with Non-Uniform Sampling

We describe an efficient NMR triple resonance approach for fast assignment of backbone amide resonance peaks in the ¹⁵N-HSQC spectrum. The exceptionally high resolutions achieved in the 3D HncocaNH and hNcocaNH experiments together with non-uniform sampling facilitate error-free sequential connection of backbone amides. Data required for the complete backbone amide assignment of the 56-residue protein GB1 domain were obtained in 14 h. Data analysis was vastly streamlined using a ‘backbone NH walk’ method to determine sequential connectivities without the need for ¹³C chemical shifts comparison. Amino acid residues in the sequentially connected NH chains are classified into two groups by a simple variation of the NMR pulse sequence, and the resulting ‘ZeBra’ stripe patterns are useful for mapping these chains to the protein sequence. In addition to resolving ambiguous assignments derived from conventional backbone experiments, this approach can be employed to rapidly assign small proteins or flexible regions in larger proteins, and to transfer assignments to mutant proteins or proteins in different ligand-binding states.     

Backbone-driven collapse in unfolded protein chains

Collapse of unfolded protein chains is an early event in folding. It affects structural properties of intrinsically disordered proteins, which take a considerable fraction of the human proteome. Collapse is generally believed to be driven by hydrophobic forces imposed by the presence of nonpolar amino acid side chains. Contributions from backbone hydrogen bonds to protein folding and stability, however, are controversial. To date, the experimental dissection of side-chain and backbone contributions has not yet been achieved because both types of interactions are integral parts of protein structure. Here, we realized this goal by applying mutagenesis and chemical modification on a set of disordered peptides and proteins. We measured the protein dimensions and kinetics of intra-chain diffusion of modified polypeptides at the level of individual molecules using fluorescence correlation spectroscopy, thereby avoiding artifacts commonly caused by aggregation of unfolded protein material in bulk. We found no contributions from side chains to collapse but, instead, identified backbone interactions as a source sufficient to form globules of native-like dimensions. The presence of backbone hydrogen bonds decreased polypeptide water solubility dramatically and accelerated the nanosecond kinetics of loop closure, in agreement with recent predictions from computer simulation. The presence of side chains, instead, slowed loop closure and modulated the dimensions of intrinsically disordered domains. It appeared that the transient formation of backbone interactions facilitates the diffusive search for productive conformations at the early stage of folding and within intrinsically disordered proteins.     

Analysis of protein folding and function using backbone modified proteins

With the recent development of chemical and biological methods to introduce backbone modifications into the polypeptide chains of proteins, there have been a growing number of site-directed mutagenesis experiments focused on understanding the role of the polypeptide backbone in protein folding and function. The substitution of a main chain amide bond with an ester bond is now a popular mutation to investigate the role of the polypeptide backbone in ligand, binding, enzyme catalysis, and protein folding. Here we review the results of studies on some 25 ester-bond containing analogues from nine different protein systems. The structural, thermodynamic, and functional consequences of introducing backbone amide- to ester-bond mutations into these protein systems are discussed.     

Protein folding and wring resonances

The polypeptide chain of a protein is shown to obey topological constraints which enable long range excitations in the form of wring modes of the protein backbone. Wring modes of proteins of specific lengths can therefore resonate with molecular modes present in the cell. It is suggested that protein folding takes place when the amplitude of a wring excitation becomes so large that it is energetically favorable to bend the protein backbone. The condition under which such structural transformations can occur is found, and it is shown that both cold and hot denaturation (the unfolding of proteins) are natural consequences of the suggested wring mode model. Native (folded) proteins are found to possess an intrinsic standing wring mode.     

Comparison of Rosetta flexible-backbone computational protein design methods on binding interactions

Computational design of binding sites in proteins remains difficult, in part due to limitations in our current ability to sample backbone conformations that enable precise and accurate geometric positioning of side chains during sequence design. Here we present a benchmark framework for comparison between flexible-backbone design methods applied to binding interactions. We quantify the ability of different flexible backbone design methods in the widely used protein design software Rosetta to recapitulate observed protein sequence profiles assumed to represent functional protein/protein and protein/small molecule binding interactions. The CoupledMoves method, which combines backbone flexibility and sequence exploration into a single acceptance step during the sampling trajectory, better recapitulates observed sequence profiles than the BackrubEnsemble and FastDesign methods, which separate backbone flexibility and sequence design into separate acceptance steps during the sampling trajectory. Flexible-backbone design with the CoupledMoves method is a powerful strategy for reducing sequence space to generate targeted libraries for experimental screening and selection.     

A comprehensive library of blocked dipeptides reveals intrinsic backbone conformational propensities of unfolded proteins

Despite prolonged scientific efforts to elucidate the intrinsic peptide backbone preferences of amino-acids based on understanding of intermolecular forces, many open questions remain, particularly concerning neighboring peptide interaction effects on the backbone conformational distribution of short peptides and unfolded proteins. Here, we show that spectroscopic studies of a complete library of 400 dipeptides reveal that, irrespective of side-chain properties, the backbone conformation distribution is narrow and they adopt polyproline II and β-strand, indicating the importance of backbone peptide solvation and electronic effects. By directly comparing the dipeptide circular dichroism and NMR results with those of unfolded proteins, the comprehensive dipeptides form a complete set of structural motifs of unfolded proteins. We thus anticipate that the present dipeptide library with spectroscopic data can serve as a useful database for understanding the nature of unfolded protein structures and for further refinements of molecular mechanical parameters.     

Deep-UV resonance Raman analysis of the Rhodobacter capsulatus cytochrome bc₁complex reveals a potential marker for the transmembrane peptide backbone

Classical strategies for structure analysis of proteins interacting with a lipid phase typically correlate ensemble secondary structure content measurements with changes in the spectroscopic responses of localized aromatic residues or reporter molecules to map regional solvent environments. Deep-UV resonance Raman (DUVRR) spectroscopy probes the vibrational modes of the peptide backbone itself, is very sensitive to the ensemble secondary structures of a protein, and has been shown to be sensitive to the extent of solvent interaction with the peptide backbone [ Wang , Y. , Purrello , R. , Georgiou , S. , and Spiro , T. G. ( 1991 ) J. Am. Chem. Soc. 113 , 6368 - 6377 ]. Here we show that a large detergent solubilized membrane protein, the Rhodobacter capsulatus cytochrome bc(1) complex, has a distinct DUVRR spectrum versus that of an aqueous soluble protein with similar overall secondary structure content. Cross-section calculations of the amide vibrational modes indicate that the peptide backbone carbonyl stretching modes differ dramatically between these two proteins. Deuterium exchange experiments probing solvent accessibility confirm that the contribution of the backbone vibrational mode differences are derived from the lipid solubilized or transmembrane α-helical portion of the protein complex. These findings indicate that DUVRR is sensitive to both the hydration status of a protein's peptide backbone, regardless of primary sequence, and its secondary structure content. Therefore, DUVRR may be capable of simultaneously measuring protein dynamics and relative water/lipid solvation of the protein.     

Disorder predictors also predict backbone dynamics for a family of disordered proteins

Several algorithms have been developed that use amino acid sequences to predict whether or not a protein or a region of a protein is disordered. These algorithms make accurate predictions for disordered regions that are 30 amino acids or longer, but it is unclear whether the predictions can be directly related to the backbone dynamics of individual amino acid residues. The nuclear Overhauser effect between the amide nitrogen and hydrogen (NHNOE) provides an unambiguous measure of backbone dynamics at single residue resolution and is an excellent tool for characterizing the dynamic behavior of disordered proteins. In this report, we show that the NHNOE values for several members of a family of disordered proteins are highly correlated with the output from three popular algorithms used to predict disordered regions from amino acid sequence. This is the first test between an experimental measure of residue specific backbone dynamics and disorder predictions. The results suggest that some disorder predictors can accurately estimate the backbone dynamics of individual amino acids in a long disordered region.     

Covariation of backbone motion throughout a small protein domain

The synchronization (correlation) of conformational fluctuations in folded proteins may influence the rates of enzyme catalysis and ligand binding as well as the stabilities of native proteins and their complexes. However, experimental detection of correlated motions remains difficult. Herein, we present an analysis of the covariation of NMR-derived backbone dynamical parameters among a family of ten mutants of a small protein. Both the spatial restriction and the time scales of backbone motions exhibit a higher degree of covariation than would be expected if the internal motions of each group were independent, providing experimental support for correlated dynamics. Application of this approach to other proteins may reveal dynamical correlations that influence catalysis, ligand-binding and/or protein stability.     

Backrub-like backbone simulation recapitulates natural protein conformational variability and improves mutant side-chain prediction

Incorporation of effective backbone sampling into protein simulation and design is an important step in increasing the accuracy of computational protein modeling. Recent analysis of high-resolution crystal structures has suggested a new model, termed backrub, to describe localized, hinge-like alternative backbone and side-chain conformations observed in the crystal lattice. The model involves internal backbone rotations about axes between C-alpha atoms. Based on this observation, we have implemented a backrub-inspired sampling method in the Rosetta structure prediction and design program. We evaluate this model of backbone flexibility using three different tests. First, we show that Rosetta backrub simulations recapitulate the correlation between backbone and side-chain conformations in the high-resolution crystal structures upon which the model was based. As a second test of backrub sampling, we show that backbone flexibility improves the accuracy of predicting point-mutant side-chain conformations over fixed backbone rotameric sampling alone. Finally, we show that backrub sampling of triosephosphate isomerase loop 6 can capture the millisecond/microsecond oscillation between the open and closed states observed in solution. Our results suggest that backrub sampling captures a sizable fraction of localized conformational changes that occur in natural proteins. Application of this simple model of backbone motions may significantly improve both protein design and atomistic simulations of localized protein flexibility.     

Backbone assignment of proteins with known structure using residual dipolar couplings

A prerequisite for NMR studies of protein-ligand interactions or protein dynamics is the assignment of backbone resonances. Here we demonstrate that protein assignment can significantly be enhanced when experimental dipolar couplings (RDCs) are matched to values back-calculated from a known three-dimensional structure. In case of small proteins, the program MARS allows assignment of more than 90% of backbone resonances without the need for sequential connectivity information. For bigger proteins, we show that the combination of sequential connectivity information with RDC-matching enables more residues to be assigned reliably and backbone assignment to be more robust against missing data. Structural or dynamic deviations from the employed 3D coordinates do not lead to an increased error rate in RDC-supported assignment. RDC-enhanced assignment is particularly useful when chemical shifts and sequential connectivity only provide a few reliable assignments.