Three-dimensional structure of the weakly associated protein homodimer SeR13 using RDCs and paramagnetic surface mapping

The traditional NMR‐based method for determining oligomeric protein structure usually involves distinguishing and assigning intra‐ and intersubunit NOEs. This task becomes challenging when determining symmetric homo‐dimer structures because NOE cross‐peaks from a given pair of protons occur at the same position whether intra‐ or intersubunit in origin. While there are isotope‐filtering strategies for distinguishing intra from intermolecular NOE interactions in these cases, they are laborious and often prove ineffectual in cases of weak dimers, where observation of intermolecular NOEs is rare. Here, we present an efficient procedure for weak dimer structure determination based on residual dipolar couplings (RDCs), chemical shift changes upon dilution, and paramagnetic surface perturbations. This procedure is applied to the Northeast Structural Genomics Consortium protein target, SeR13, a negatively charged Staphylococcus epidermidis dimeric protein (K_d 3.4 ± 1.4 mM) composed of 86 amino acids. A structure determination for the monomeric form using traditional NMR methods is presented, followed by a dimer structure determination using docking under orientation constraints from RDCs data, and scoring under residue pair potentials and shape‐based predictions of RDCs. Validation using paramagnetic surface perturbation and chemical shift perturbation data acquired on sample dilution is also presented. The general utility of the dimer structure determination procedure and the possible relevance of SeR13 dimer formation are discussed.     

Refined solution structure of the 82-kDa enzyme malate synthase G from joint NMR and synchrotron SAXS restraints

Determination of the accurate three-dimensional structure of large proteins by NMR remains challenging due to a loss in the density of experimental restraints resulting from the often prerequisite perdeuteration. Solution small-angle scattering, which carries long-range translational information, presents an opportunity to enhance the structural accuracy of derived models when used in combination with global orientational NMR restraints such as residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs). We have quantified the improvements in accuracy that can be obtained using this strategy for the 82 kDa enzyme Malate Synthase G (MSG), currently the largest single chain protein solved by solution NMR. Joint refinement against NMR and scattering data leads to an improvement in structural accuracy as evidenced by a decrease from approximately 4.5 to approximately 3.3 A of the backbone rmsd between the derived model and the high-resolution X-ray structure, PDB code 1D8C. This improvement results primarily from medium-angle scattering data, which encode the overall molecular shape, rather than the lowest angle data that principally determine the radius of gyration and the maximum particle dimension. The effect of the higher angle data, which are dominated by internal density fluctuations, while beneficial, is also found to be relatively small. Our results demonstrate that joint NMR/SAXS refinement can yield significantly improved accuracy in solution structure determination and will be especially well suited for the study of systems with limited NMR restraints such as large proteins, oligonucleotides, or their complexes.     

Combining NMR and small angle X-ray and neutron scattering in the structural analysis of a ternary protein-RNA complex

Many processes in the regulation of gene expression and signaling involve the formation of protein complexes involving multi-domain proteins. Individual domains that mediate protein-protein and protein-nucleic acid interactions are typically connected by flexible linkers, which contribute to conformational dynamics and enable the formation of complexes with distinct binding partners. Solution techniques are therefore required for structural analysis and to characterize potential conformational dynamics. Nuclear magnetic resonance spectroscopy (NMR) provides such information but often only sparse data are obtained with increasing molecular weight of the complexes. It is therefore beneficial to combine NMR data with additional structural restraints from complementary solution techniques. Small angle X-ray/neutron scattering (SAXS/SANS) data can be efficiently combined with NMR-derived information, either for validation or by providing additional restraints for structural analysis. Here, we show that the combination of SAXS and SANS data can help to refine structural models obtained from data-driven docking using HADDOCK based on sparse NMR data. The approach is demonstrated with the ternary protein-protein-RNA complex involving two RNA recognition motif (RRM) domains of Sex-lethal, the N-terminal cold shock domain of Upstream-to-N-Ras, and msl-2 mRNA. Based on chemical shift perturbations we have mapped protein-protein and protein-RNA interfaces and complemented this NMR-derived information with SAXS data, as well as SANS measurements on subunit-selectively deuterated samples of the ternary complex. Our results show that, while the use of SAXS data is beneficial, the additional combination with contrast variation in SANS data resolves remaining ambiguities and improves the docking based on chemical shift perturbations of the ternary protein-RNA complex.     

Cross-linking and mass spectrometry methodologies to facilitate structural biology: finding a path through the maze

Multiprotein complexes, rather than individual proteins, make up a large part of the biological macromolecular machinery of a cell. Understanding the structure and organization of these complexes is critical to understanding cellular function. Chemical cross-linking coupled with mass spectrometry is emerging as a complementary technique to traditional structural biology methods and can provide low-resolution structural information for a multitude of purposes, such as distance constraints in computational modeling of protein complexes. In this review, we discuss the experimental considerations for successful application of chemical cross-linking-mass spectrometry in biological studies and highlight three examples of such studies from the recent literature. These examples (as well as many others) illustrate the utility of a chemical cross-linking-mass spectrometry approach in facilitating structural analysis of large and challenging complexes.     

Mapping protein-protein interactions in solution by NMR spectroscopy.

NMR is very well suited to the study of especially weak protein-protein interactions, as no crystallization is required. The available NMR methods to this end are reviewed and illustrated with applications from the recent biochemical literature: intermolecular NOEs, cross-saturation, chemical shift perturbation, dynamics and exchange perturbation, paramagnetic methods, and dipolar orientation. Most of these methods are now routinely applied for complexes with total molecular mass of 60 kDa and can likely be applied to systems up to 1000 kDa. A substantial fraction of complexes studied show distinct effects of induced fit affecting structural and dynamical properties beyond the contact interface.     

Expression, purification, and characterization of Thermotoga maritima membrane proteins for structure determination

Structural studies of integral membrane proteins typically rely upon detergent micelles as faithful mimics of the native lipid bilayer. Therefore, membrane protein structure determination would be greatly facilitated by biophysical techniques that are capable of evaluating and assessing the fold and oligomeric state of these proteins solubilized in detergent micelles. In this study, an approach to the characterization of detergent-solubilized integral membrane proteins is presented. Eight Thermotoga maritima membrane proteins were screened for solubility in 11 detergents, and the resulting soluble protein-detergent complexes were characterized with small angle X-ray scattering (SAXS), nuclear magnetic resonance (NMR) spectroscopy, circular dichroism (CD) spectroscopy, and chemical cross-linking to evaluate the homogeneity, oligomeric state, radius of gyration, and overall fold. A new application of SAXS is presented, which does not require density matching, and NMR methods, typically used to evaluate soluble proteins, are successfully applied to detergent-solubilized membrane proteins. Although detergents with longer alkyl chains solubilized the most proteins, further characterization indicates that some of these protein-detergent complexes are not well suited for NMR structure determination due to conformational exchange and protein oligomerization. These results emphasize the need to screen several different detergents and to characterize the protein-detergent complex in order to pursue structural studies. Finally, the physical characterization of the protein-detergent complexes indicates optimal solution conditions for further structural studies for three of the eight overexpressed membrane proteins.     

Determination of the interface of a large protein complex by transferred cross-saturation measurements

In an earlier paper, it was shown that the cross-saturation method enables us to identify the contact residues of large protein complexes in a more rigorous manner than is possible using chemical shift perturbation and hydrogen-deuterium exchange experiments. However, there are limitations within the determination of the contact residues by the cross-saturation method, in that the method is difficult to apply to protein complexes with a molecular mass over 150 kDa and/or with weak binding, since the resonances originating from the complexes should be observed directly in the method. In the present work, to overcome these limitations, we carried out the cross-saturation measurements under conditions of a fast exchange between free and bound states on the NMR time-scale, and determined the contact residues of the complex of the B domain of protein A and intact IgG, which has a molecular mass of 164 kDa and shows weak binding.     

Quantitative use of chemical shifts for the modeling of protein complexes

Despite recent advances in the modeling of protein-protein complexes by docking, additional information is often required to identify the best solutions. For this purpose, NMR data deliver valuable restraints that can be used in the sampling and/or the scoring stage, like in the data-driven docking approach HADDOCK that can make use of NMR chemical shift perturbation (CSP) data to define the binding site on each protein and drive the docking. We show here that a quantitative use of chemical shifts (CS) in the scoring stage can help to resolve ambiguities. A quantitative CS-RMSD score based on (1) H(α) ,(13) C(α) , and (15) N chemical shifts ranks the best solutions always at the top, as demonstrated on a small benchmark of four complexes. It is implemented in a new docking protocol, CS-HADDOCK, which combines CSP data as ambiguous interaction restraints in the sampling stage with the CS-RMSD score in the final scoring stage. This combination of qualitative and quantitative use of chemical shifts increases the reliability of data-driven docking for the structure determination of complexes from limited NMR data.     

A device for the measurement of residual chemical shift anisotropy and residual dipolar coupling in soluble and membrane-associated proteins

Residual dipolar coupling (RDC) and residual chemical shift anisotropy (RCSA) report on orientational properties of a dipolar bond vector and a chemical shift anisotropy principal axis system, respectively. They can be highly complementary in the analysis of backbone structure and dynamics in proteins as RCSAs generally include a report on vectors out of a peptide plane while RDCs usually report on in-plane vectors. Both RDC and RCSA average to zero in isotropic solutions and require partial orientation in a magnetic field to become observable. While the alignment and measurement of RDC has become routine, that of RCSA is less common. This is partly due to difficulties in providing a suitable isotopic reference spectrum for the measurement of the small chemical shift offsets coming from RCSA. Here we introduce a device (modified NMR tube) specifically designed for accurate measurement of reference and aligned spectra for RCSA measurements, but with a capacity for RDC measurements as well. Applications to both soluble and membrane anchored proteins are illustrated.     

Residual dipolar coupling derived orientational constraints on ligand geometry in a 53 kDa protein-ligand complex.

The geometric relationships between ligands and the functional groups that bind ligands in soluble ligand-protein complexes have traditionally been deduced from distance constraints between pairs of NMR active nuclei spanning the ligand-protein interface. Frequently, the steep inverse distance dependence of the nuclear Overhauser effect (NOE), from which the distance constraints are derived, makes identification of sufficient numbers of constraints difficult. In these cases the ability to supplement NOE-derived information with distance-independent angular information can be very important. Here, the observation of residual dipolar couplings from alpha-methyl mannose bound to mannose binding-protein in a dilute liquid crystalline medium has allowed the determination of a bound ligand's average orientation. The 3-fold rotational symmetry of mannose-binding protein defines its orientational tensor and obviates the need to determine experimentally the protein's average orientation. Through superimposition of ligand and protein orientational tensors we describe the binding geometry of alpha-methyl mannose bound to mannose-binding protein. This new method is of general applicability to the study of ligands bound to proteins, and it is of particular interest when neither X-ray crystallography nor NOE techniques can provide sufficient information to describe binding geometries.     

Rapid analysis of large protein-protein complexes using NMR-derived orientational constraints: the 95 kDa complex of LpxA with acyl carrier protein

Characterization of protein-protein interactions that are critical to the specific function of many biological systems has become a primary goal of structural biology research. Analysis of these interactions by structural techniques is, however, challenging due to inherent limitations of the techniques and because many of the interactions are transient, and suitable complexes are difficult to isolate. In particular, structural studies of large protein complexes by traditional solution NMR methods are difficult due to a priori requirement of extensive assignments and a large number of intermolecular restraints for the complex. An approach overcoming some of these challenges by utilizing orientational restraints from residual dipolar couplings collected on solution NMR samples is presented. The approach exploits existing structures of individual components, including the symmetry properties of some of these structures, to assemble rapidly models for relatively large protein-protein complexes. An application is illustrated with a 95 kDa homotrimeric complex of the acyltransferase protein, LpxA (UDP-N-acetylglucosamine acyltransferase), and acyl carrier protein. LpxA catalyzes the first step in the biosynthesis of the lipid A component of lipopolysaccharide in Gram-negative bacteria. The structural model generated for this complex can be useful in the design of new anti-bacterial agents that inhibit the biosynthesis of lipid A.     

Deciphering the molecular organization of GET pathway chaperones through native mass spectrometry

Get3/4/5 chaperone complex is responsible for targeting C-terminal tail-anchored membrane proteins to the endoplasmic reticulum. Despite the availability of several crystal structures of independent proteins and partial structures of subcomplexes, different models of oligomeric states and structural organization have been proposed for the protein complexes involved. Here, using native mass spectrometry (Native-MS), coupled with intact dissociation, we show that Get4/5 exclusively forms a tetramer using both Get5/5 and a novel Get4/4 dimerization interface. Addition of Get3 to this leads to a hexameric (Get3)₂-(Get4)₂-(Get5)₂ complex with closed-ring cyclic architecture. We further validate our claims through molecular modeling and mutational abrogation of the proposed interfaces. Native-MS has become a principal tool to determine the state of oligomeric organization of proteins. The work demonstrates that for multiprotein complexes, native-MS, coupled with molecular modeling and mutational perturbation, can provide an alternative route to render a detailed view of both the oligomeric states as well as the molecular interfaces involved. This is especially useful for large multiprotein complexes with large unstructured domains that make it recalcitrant to conventional structure determination approaches.     

Filtering and selection of structural models: combining docking and NMR

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

Gramicidin cation channel: an experimental determination of the right-handed helix sense and verification of beta-type hydrogen bonding

Due to the difficulty of obtaining protein/lipid cocrystals for diffraction studies, structural research on intrinsic membrane proteins and polypeptides has been largely restricted to indirect experimental techniques. Hence, many fundamental questions associated with peptide/lipid systems remain unanswered. In particular, the handedness of the gramicidin A transmembrane ion channel incorporated into lipid bilayers has been an open question for nearly two decades. In this study, solid-state 15N NMR spectroscopy is employed to probe directly the secondary structure of the polypeptide backbone. Recent determinations of the 15N chemical shift anisotropy tensor with respect to the molecular frame enable the quantitative evaluation of the 15N chemical shift resonances obtained from oriented dimyristoylphosphatidylcholine (DMPC) bilayer samples containing specific site 15N labeled gramicidin. This direct structural approach verifies the beta-sheet hydrogen-bonding pattern proposed by Urry [Urry, D. W. (1971) Proc. Natl. Acad. Sci. U.S.A. 68, 672-676] and determines that in our DMPC bilayer preparations the gramicidin channel is right-handed. Additional structural information is provided by the 15N chemical shift data in the form of orientational constraints on the C alpha-C alpha axis orientation of individual peptides relative to the helix axis. The significance of these solid-state NMR results lies in the direct determination of the helix sense and the verification of the beta-type hydrogen bonding, in the development of the solid-state NMR methods for obtaining such information, and in emphasizing the importance of having direct structural data at atomic resolution.     

NMR analysis of protein interactions

Recent technological advances in NMR spectroscopy have alleviated the size limitations for the determination of biomolecular structures in solution. At the same time, novel NMR parameters such as residual dipolar couplings are providing greater accuracy. As this review shows, the structures of protein-protein and protein-nucleic acid complexes up to 50 kDa can now be accurately determined. Although de novo structure determination still requires considerable effort, information on interaction surfaces from chemical shift perturbations is much easier to obtain. Advances in modelling and data-driven docking procedures allow this information to be used for determining approximate structures of biomolecular complexes. As a result, a wealth of information has become available on the way in which proteins interact with other biomolecules. Of particular interest is the fact that these NMR-based methods can be applied to weak and transient protein-protein complexes that are difficult to study by other structural methods.     

A graphical method for analyzing distance restraints using residual dipolar couplings for structure determination of symmetric protein homo-oligomers

High-resolution structure determination of homo-oligomeric protein complexes remains a daunting task for NMR spectroscopists. Although isotope-filtered experiments allow separation of intermolecular NOEs from intramolecular NOEs and determination of the structure of each subunit within the oligomeric state, degenerate chemical shifts of equivalent nuclei from different subunits make it difficult to assign intermolecular NOEs to nuclei from specific pairs of subunits with certainty, hindering structural analysis of the oligomeric state. Here, we introduce a graphical method, DISCO, for the analysis of intermolecular distance restraints and structure determination of symmetric homo-oligomers using residual dipolar couplings. Based on knowledge that the symmetry axis of an oligomeric complex must be parallel to an eigenvector of the alignment tensor of residual dipolar couplings, we can represent distance restraints as annuli in a plane encoding the parameters of the symmetry axis. Oligomeric protein structures with the best restraint satisfaction correspond to regions of this plane with the greatest number of overlapping annuli. This graphical analysis yields a technique to characterize the complete set of oligomeric structures satisfying the distance restraints and to quantitatively evaluate the contribution of each distance restraint. We demonstrate our method for the trimeric E. coli diacylglycerol kinase, addressing the challenges in obtaining subunit assignments for distance restraints. We also demonstrate our method on a dimeric mutant of the immunoglobulin-binding domain B1 of streptococcal protein G to show the resilience of our method to ambiguous atom assignments. In both studies, DISCO computed oligomer structures with high accuracy despite using ambiguously assigned distance restraints.     

High-resolution structure determination of the CylR2 homodimer using paramagnetic relaxation enhancement and structure-based prediction of molecular alignment

Structure determination of homooligomeric proteins by NMR spectroscopy is difficult due to the lack of chemical shift perturbation data, which is very effective in restricting the binding interface in heterooligomeric systems, and the difficulty of obtaining a sufficient number of intermonomer distance restraints. Here we solved the high-resolution solution structure of the 15.4 kDa homodimer CylR2, the regulator of cytolysin production from Enterococcus faecalis, which deviates by 1.1 angstroms from the previously determined X-ray structure. We studied the influence of different experimental information such as long-range distances derived from paramagnetic relaxation enhancement, residual dipolar couplings, symmetry restraints and intermonomer Nuclear Overhauser Effect restraints on the accuracy of the derived structure. In addition, we show that it is useful to combine experimental information with methods of ab initio docking when the available experimental data are not sufficient to obtain convergence to the correct homodimeric structure. In particular, intermonomer distances may not be required when residual dipolar couplings are compared to values predicted on the basis of the charge distribution and the shape of ab initio docking solutions.     

A geometric arrangement algorithm for structure determination of symmetric protein homo-oligomers from NOEs and RDCs

Nuclear magnetic resonance (NMR) spectroscopy is a primary tool to perform structural studies of proteins in physiologically-relevant solution conditions. Restraints on distances between pairs of nuclei in the protein, derived from the nuclear Overhauser effect (NOE), provide information about the structure of the protein in its folded state. NMR studies of symmetric protein homo-oligomers present a unique challenge. Using X-filtered NOESY experiments, it is possible to determine whether an NOE restrains a pair of protons across different subunits or within a single subunit, but current experimental techniques are unable to determine in which subunits the restrained protons lie. Consequently, it is difficult to assign NOEs to particular pairs of subunits with certainty, thus hindering the structural analysis of the oligomeric state. Computational approaches are needed to address this subunit ambiguity, but traditional solutions often rely on stochastic search coupled with simulated annealing and simulations of simplified molecular dynamics, which have many tunable parameters that must be chosen carefully and can also fail to report structures consistent with the experimental restraints. In addition, these traditional approaches rarely provide guarantees on running time or solution quality. We reduce the structure determination of homo-oligomers with cyclic symmetry to computing geometric arrangements of unions of annuli in a plane. Our algorithm, disco, runs in expected O(n2) time, where n is the number of distance restraints, potentially assigned ambiguously. disco is guaranteed to report the exact set of oligomer structures consistent with the distance restraints and also with orientational restraints from residual dipolar couplings (RDCs). We demonstrate our method using two symmetric protein complexes: the trimeric E. coli diacylglycerol kinase (DAGK) and a dimeric mutant of the immunoglobulin-binding domain B1 of streptococcal protein G (GB1). In both cases, disco computes oligomer structures with high precision and also finds distance restraints that are either mutually inconsistent or inconsistent with the RDCs. The entire protocol DISCO has been completely automated in a software package that is freely available and open-source at www.cs.duke.edu/donaldlab/software.php.     

A novel NMR method for determining the interfaces of large protein-protein complexes

Identification of the interfaces of large (Mr > 50,000) protein-protein complexes in solution by high resolution NMR has typically been achieved using experiments involving chemical shift perturbation and/or hydrogen-deuterium exchange of the main chain amide groups of the proteins. Interfaces identified using these techniques, however, are not always identical to those revealed using X-ray crystallography. In order to identify the contact residues in a large protein-protein complex more accurately, we developed a novel NMR method that uses cross-saturation phenomena in combination with TROSY detection in an optimally deuterium labeled system.