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

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

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

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

A novel strategy to determine protein structures using exclusively residual dipolar coupling

Recent advances in NMR techniques to measure anisotropic spin interactions such as residual dipolar coupling (RDC) have provided better insights into protein structure as well as dynamics. Exploitation of RDC, however, still remains challenging because its successful application requires a reasonable starting model. Using the singular value decomposition method, we have recently developed an RDC restraint potential to optimally extract orientational information from RDC without the prerequisite of any structural information. In the present study, its efficacy is further illustrated by folding a beta-hairpin and alpha-helix of protein G from extended conformations with RDC restraints alone by employing the replica exchange torsion angle molecular dynamics (REX-TAMD) technique. Subsequently, the entire structure of protein G has been determined accurately using the developed fragment superposition method (FRAGSUM). In FRAGSUM, each overlapping fragments (10 amino acids long) is first folded individually by REX-TAMD, and then the common amino acids are superimposed to determine the entire structure. Because FRAGSUM does not require any additional information besides RDC, it offers a new strategy for de novo structure determination using exclusively RDC.     

Refinement of NMR structures using implicit solvent and advanced sampling techniques

NMR biomolecular structure calculations exploit simulated annealing methods for conformational sampling and require a relatively high level of redundancy in the experimental restraints to determine quality three-dimensional structures. Recent advances in generalized Born (GB) implicit solvent models should make it possible to combine information from both experimental measurements and accurate empirical force fields to improve the quality of NMR-derived structures. In this paper, we study the influence of implicit solvent on the refinement of protein NMR structures and identify an optimal protocol of utilizing these improved force fields. To do so, we carry out structure refinement experiments for model proteins with published NMR structures using full NMR restraints and subsets of them. We also investigate the application of advanced sampling techniques to NMR structure refinement. Similar to the observations of Xia et al. (J.Biomol. NMR 2002, 22, 317-331), we find that the impact of implicit solvent is rather small when there is a sufficient number of experimental restraints (such as in the final stage of NMR structure determination), whether implicit solvent is used throughout the calculation or only in the final refinement step. The application of advanced sampling techniques also seems to have minimal impact in this case. However, when the experimental data are limited, we demonstrate that refinement with implicit solvent can substantially improve the quality of the structures. In particular, when combined with an advanced sampling technique, the replica exchange (REX) method, near-native structures can be rapidly moved toward the native basin. The REX method provides both enhanced sampling and automatic selection of the most native-like (lowest energy) structures. An optimal protocol based on our studies first generates an ensemble of initial structures that maximally satisfy the available experimental data with conventional NMR software using a simplified force field and then refines these structures with implicit solvent using the REX method. We systematically examine the reliability and efficacy of this protocol using four proteins of various sizes ranging from the 56-residue B1 domain of Streptococcal protein G to the 370-residue Maltose-binding protein. Significant improvement in the structures was observed in all cases when refinement was based on low-redundancy restraint data. The proposed protocol is anticipated to be particularly useful in early stages of NMR structure determination where a reliable estimate of the native fold from limited data can significantly expedite the overall process. This refinement procedure is also expected to be useful when redundant experimental data are not readily available, such as for large multidomain biomolecules and in solid-state NMR structure determination.     

Refinement of multidomain protein structures by combination of solution small-angle X-ray scattering and NMR data

Determination of the 3D structures of multidomain proteins by solution NMR methods presents a number of unique challenges related to their larger molecular size and the usual scarcity of constraints at the interdomain interface, often resulting in a decrease in structural accuracy. In this respect, experimental information from small-angle scattering of X-ray radiation in solution (SAXS) presents a suitable complement to the NMR data, as it provides an independent constraint on the overall molecular shape. A computational procedure is described that allows incorporation of such SAXS data into the mainstream high-resolution macromolecular structure refinement. The method is illustrated for a two-domain 177-amino-acid protein, gammaS crystallin, using an experimental SAXS data set fitted at resolutions from approximately 200 A to approximately 30 A. Inclusion of these data during structure refinement decreases the backbone coordinate root-mean-square difference between the derived model and the high-resolution crystal structure of a 54% homologous gammaB crystallin from 1.96 +/- 0.07 A to 1.31 +/- 0.04 A. Combining SAXS data with NMR restraints can be accomplished at a moderate computational expense and is expected to become useful for multidomain proteins, multimeric assemblies, and tight macromolecular complexes.     

Residual dipolar coupling measurements of transmembrane proteins using aligned low-q bicelles and high-resolution magic angle spinning NMR spectroscopy

Bicelles are a major medium form to produce weak alignment of soluble proteins for residual dipolar coupling (RDC) measurements. The obstacle to using the same type of bicelles for transmembrane proteins with solution-state NMR spectroscopy is the loss of signals due to the adhesion or penetration of the proteins into large bicelles, resulting in slow protein tumbling. In this study, weak alignment of the second and third transmembrane domains (TM23) of the human glycine receptor (GlyR) was achieved in low-q bicelles (q = DMPC/DHPC). Although protein-free bicelles with such low q would likely show isotropic properties, the insertion of TM23 induced weakly preferred orientations so that the RDC of the embedded protein can be measured. The extent of the alignment increased but the TM23 signal intensity decreased when q was varied from 0.19 to 0.60. A q of 0.50 was found to be an optimal compromise between alignment and the signal-to-noise ratio. In each pair of NMR experiments for RDC measurements, the same sample and pulse sequence were used, with one being performed at high-resolution magic-angle spinning to obtain pure J-couplings without RDC. A meaningful structure refinement in bicelles was possible by iteratively fitting the experimental RDCs to the back-calculated RDCs using the high-resolution NMR structure of GlyR TM23 in trifluoroethanol as the starting template. Combination of this method with the conventional high-resolution NMR in membrane mimicking mixtures of water and organic solvents offers an attractive way to derive structural information for membrane proteins in their native environment.     

Protein backbone structure determination using RDC: An inverse kinematics approach with fast and exact solutions

Residual dipolar couplings (RDC) of proteins dissolved in anisotropic media promise to speed up the determination of protein structures. We consider the backbone as a robotic mechanism and formulate inverse kinematics problems using RDC restraints from two media. The φ, ψ of each secondary structure element (SSE) are computed from oriented vectors in consecutive peptide planes. We search for the optimum conformation joining the solutions of two independent backbone halves. The matrix transforming the vector Z of a global frame from one SSE into the other determines their orientation. Three distance constraints between two oriented SSE determine their relative position by solving nine polynomial equations. The benefit of this method is that complete and accurate solutions are obtained overcoming the local minima problems of heuristic procedures. The algorithm is implemented on MAPLE using the least number of experimental data; the runtimes take an order of seconds on a common PC. © 2013 Wiley Periodicals, Inc.     

Structure calculation from unambiguous long-range amide and methyl 1H-1H distance restraints for a microcrystalline protein with MAS solid-state NMR spectroscopy

Magic-angle spinning (MAS) solid-state NMR becomes an increasingly important tool for the determination of structures of membrane proteins and amyloid fibrils. Extensive deuteration of the protein allows multidimensional experiments with exceptionally high sensitivity and resolution to be obtained. Here we present an experimental strategy to measure highly unambiguous spatial correlations for distances up to 13 ?. Two complementary three-dimensional experiments, or alternatively a four-dimensional experiment, yield highly unambiguous cross-peak assignments, which rely on four encoded chemical shift dimensions. Correlations to residual aliphatic protons are accessible via synchronous evolution of the (15)N and (13)C chemical shifts, which encode valuable amide-methyl distance restraints. On average, we obtain six restraints per residue. Importantly, 50% of all restraints correspond to long-range distances between residues i and j with |i - j| > 5, which are of particular importance in structure calculations. Using ARIA, we calculate a high-resolution structure for the microcrystalline 7.2 kDa α-spectrin SH3 domain with a backbone precision of ～1.1 ?.     

Structure determination of protein-protein complexes using NMR chemical shifts: case of an endonuclease colicin-immunity protein complex

Nuclear magnetic resonance (NMR) spectroscopy provides a range of powerful techniques for determining the structures and the dynamics of proteins. The high-resolution determination of the structures of protein-protein complexes, however, is still a challenging problem for this approach, since it can normally provide only a limited amount of structural information at protein-protein interfaces. We present here the determination using NMR chemical shifts of the structure (PDB code 2K5X) of the cytotoxic endonuclease domain from bacterial toxin colicin (E9) in complex with its cognate immunity protein (Im9). In order to achieve this result, we introduce the CamDock method, which combines a flexible docking procedure with a refinement that exploits the structural information provided by chemical shifts. The results that we report thus indicate that chemical shifts can be used as structural restraints for the determination of the conformations of protein complexes that are difficult to obtain by more standard NMR approaches.     

Determination of the structures of symmetric protein oligomers from NMR chemical shifts and residual dipolar couplings

Symmetric protein dimers, trimers, and higher-order cyclic oligomers play key roles in many biological processes. However, structural studies of oligomeric systems by solution NMR can be difficult due to slow tumbling of the system and the difficulty in identifying NOE interactions across protein interfaces. Here, we present an automated method (RosettaOligomers) for determining the solution structures of oligomeric systems using only chemical shifts, sparse NOEs, and domain orientation restraints from residual dipolar couplings (RDCs) without a need for a previously determined structure of the monomeric subunit. The method integrates previously developed Rosetta protocols for solving the structures of monomeric proteins using sparse NMR data and for predicting the structures of both nonintertwined and intertwined symmetric oligomers. We illustrated the performance of the method using a benchmark set of nine protein dimers, one trimer, and one tetramer with available experimental data and various interface topologies. The final converged structures are found to be in good agreement with both experimental data and previously published high-resolution structures. The new approach is more readily applicable to large oligomeric systems than conventional structure-determination protocols, which often require a large number of NOEs, and will likely become increasingly relevant as more high-molecular weight systems are studied by NMR.     

Paramagnetic-based NMR restraints lift residual dipolar coupling degeneracy in multidomain detergent-solubilized membrane proteins

Residual dipolar couplings (RDCs) are widely used as orientation-dependent NMR restraints to improve the resolution of the NMR conformational ensemble of biomacromolecules and define the relative orientation of multidomain proteins and protein complexes. However, the interpretation of RDCs is complicated by the intrinsic degeneracy of analytical solutions and protein dynamics that lead to ill-defined orientations of the structural domains (ghost orientations). Here, we illustrate how restraints from paramagnetic relaxation enhancement (PRE) experiments lift the orientational ambiguity of multidomain membrane proteins solubilized in detergent micelles. We tested this approach on monomeric phospholamban (PLN), a 52-residue membrane protein, which is composed of two helical domains connected by a flexible loop. We show that the combination of classical solution NMR restraints (NOEs and dihedral angles) with RDC and PRE constraints resolves topological ambiguities, improving the convergence of the PLN structural ensemble and giving the depth of insertion of the protein within the micelle. The combination of RDCs with PREs will be necessary for improving the accuracy and precision of membrane protein conformational ensembles, where three-dimensional structures are dictated by interactions with the membrane-mimicking environment rather than compact tertiary folds common in globular proteins.     

Direct structure determination using residual dipolar couplings: reaction-site conformation of methionine sulfoxide reductase in solution

Residual dipolar couplings (RDC) from partially aligned molecules provide long-range structural data and are thus particularly well adapted to rapid structure validation or protein fold recognition. Extensive measurements in two alignment media can also provide precise de novo structure from RDC alone. We have applied a novel combination of these approaches to the study of methionine sulfoxide reductase (MsrA) from Erwinia chrysanthemi, a 27 kDa enzyme essential for repairing oxidative stress damage. The tertiary fold was initially validated by comparing backbone RDC to expected values from the crystal structure of the homologous MsrA from Escherichia coli. Good agreement was found throughout the chain, verifying the overall topology of the molecule, with the exception of the catalytically important peptide P196-L202, where strong and systematic RDC violation was observed. No evidence for local differential mobility in this region was detected, implying that the structure of the strand differs in the two molecules. We have therefore applied the de novo approach meccano to determine the conformation of this peptide using only RDC. A single conformation is found that is in agreement with all measured data. The aligned peptide can be docked onto the expected covalence of the rest of the template molecule while respecting its strictly defined relative orientation. In contrast to the structure of MsrA from E. coli, the reactive side chain of Cys200 is oriented toward the interior of the molecule and therefore closer to the catalytic Cys53, obviating the need for previously proposed conformational reorganization prior to formation of this disulfide intermediate. This analysis requires only backbone assignment and uses unambiguously assigned and readily measurable structural data, thereby greatly economizing investigation time compared to established nuclear Overhauser effect- (nOe-) based structure calculation methods.     

Explicit treatment of force contribution from alignment tensor using overdetermined linear equations and its application in NMR structure determination

Residual dipolar coupling (RDC) provides valuable information about the orientation of each internuclear vector in a macromolecule with respect to the static magnetic field. However, structure determination utilizing RDC still remains challenging without additional restraints such as NOE. In this context, a novel approach has been developed to efficiently extract structural information from RDC by successive application of singular value decomposition (SVD) method in the course of NMR structure determination. Force contribution from the alignment tensor is rigorously formulated in the context of SVD, and assessments have been made to verify its numerical accuracy. The efficacy of this approach is illustrated by showing that RDC restraints alone can restore a distorted beta-hairpin to native-like structure using the replica-exchange molecular dynamics simulations.     

Completely automated, highly error-tolerant macromolecular structure determination from multidimensional nuclear overhauser enhancement spectra and chemical shift assignments

The major rate-limiting step in high-throughput NMR protein structure determination involves the calculation of a reliable initial fold, the elimination of incorrect nuclear Overhauser enhancement (NOE) assignments, and the resolution of NOE assignment ambiguities. We present a robust approach to automatically calculate structures with a backbone coordinate accuracy of 1.0-1.5 A from datasets in which as much as 80% of the long-range NOE information (i.e., between residues separated by more than five positions in the sequence) is incorrect. The current algorithm differs from previously published methods in that it has been expressly designed to ensure that the results from successive cycles are not biased by the global fold of structures generated in preceding cycles. Consequently, the method is highly error tolerant and is not easily funnelled down an incorrect path in either three-dimensional structure or NOE assignment space. The algorithm incorporates three main features: a linear energy function representation of the NOE restraints to allow maximization of the number of simultaneously satisfied restraints during the course of simulated annealing; a method for handling the presence of multiple possible assignments for each NOE cross-peak which avoids local minima by treating each possible assignment as if it were an independent restraint; and a probabilistic method to permit both inactivation and reactivation of all NOE restraints on the fly during the course of simulated annealing. NOE restraints are never removed permanently, thereby significantly reducing the likelihood of becoming trapped in a false minimum of NOE assignment space. The effectiveness of the algorithm is demonstrated using completely automatically peak-picked experimental NOE data from two proteins: interleukin-4 (136 residues) and cyanovirin-N (101 residues). The limits of the method are explored using simulated data on the 56-residue B1 domain of Streptococcal protein G.     

De novo determination of bond orientations and order parameters from residual dipolar couplings with high accuracy

We report the de novo determination of 15N-1H bond orientations and motional order parameters for the protein ubiquitin with high accuracy based solely on NMR residual dipolar coupling measurements made in six distinct alignment media. The resulting bond orientations are in agreement with RDC-refined orientations of either solid or solution state coordinates to within approximately 2 degrees , which is also the estimated precision of the resulting orientations. The squared generalized order parameters, which reflect amplitudes of motion spanning the picosecond to millisecond time scales, exhibit a correlation with picosecond time scale order parameters derived from conventional NMR 15N spin relaxation methods. Provided that RDC measurements can be obtained using many different alignment media, this approach (called direct interpretation of dipolar couplings) may significantly impact the attainable accuracy and the molecular weight range accessible to NMR structure determination in the solution state, as well as provide a route for the study of motions occurring on the nanosecond to microsecond time scales, which have been traditionally difficult to study at atomic resolution.     

Increasing the accuracy of solution NMR structures of membrane proteins by application of residual dipolar couplings. High-resolution structure of outer membrane protein A

The structure determination of membrane proteins is one of the most challenging applications of solution NMR spectroscopy. The paucity of distance information available from the highly deuterated proteins employed requires new approaches in structure determination. Here we demonstrate that significant improvement in the structure accuracy of the membrane protein OmpA can be achieved by refinement with residual dipolar couplings (RDCs). The application of charged polyacrylamide gels allowed us to obtain two alignments and accurately measure numerous heteronuclear dipolar couplings. Furthermore, we have demonstrated that using a large set of RDCs in the refinement can yield a structure with 1 A rms deviation to the backbone of the high-resolution crystal structure. Our simulations with various data sets indicate that dipolar couplings will be critical for obtaining accurate structures of membrane proteins.     

Protein Folding with a Reduced Model and Inaccurate Short‐Range Restraints

Summary: A reduced high‐coordination lattice protein model and the Replica Exchange Monte Carlo sampling were employed in de novo folding simulations of a set of representative small proteins. Three distinct situations were analyzed. In the first series of simulations, the folding was controlled purely by the generic force field of the model. In the second, a bias was introduced towards the theoretically predicted secondary structure. Finally, we superimposed soft restraints towards the native‐like local conformation of the backbone. The short‐range restraints used in these simulations are based on approximate values of ϕ and ψ dihedral angles, which may simulate restraints derived from inaccurate experimental measurements. Incorporating such data into the reduced model required developing a procedure, which transforms the ϕ and ψ coordinates into coordinates of the protein alpha carbon trace. It has been shown that such limited data are sufficient for de novo determination of three‐dimensional structures of small and topologically not too complex proteins.

Protein folding based on secondary structure prediction and simulated torsion angles data.