Measurement of one-bond 13Calpha-1Halpha residual dipolar coupling constants in proteins by selective manipulation of CalphaHalpha spins

Diversity in the structures and ligand-binding sites of nematode fatty acid and retinol-binding proteins revealed by Na-FAR-1 from Necator americanus

Fatty acid and retinol-binding proteins (FARs) comprise a family of unusual α-helix rich lipid-binding proteins found exclusively in nematodes. They are secreted into host tissues by parasites of plants, animals and humans. The structure of a FAR protein from the free-living nematode Caenorhabditis elegans is available, but this protein [C. elegans FAR-7 (Ce-FAR-7)] is from a subfamily of FARs that does not appear to be important at the host/parasite interface. We have therefore examined [Necator americanus FAR-1 (Na-FAR-1)] from the blood-feeding intestinal parasite of humans, N. americanus. The 3D structure of Na-FAR-1 in its ligand-free and ligand-bound forms, determined by NMR (nuclear magnetic resonance) spectroscopy and X-ray crystallography respectively, reveals an α-helical fold similar to Ce-FAR-7, but Na-FAR-1 possesses a larger and more complex internal ligand-binding cavity and an additional C-terminal α-helix. Titration of apo-Na-FAR-1 with oleic acid, analysed by NMR chemical shift perturbation, reveals that at least four distinct protein-ligand complexes can be formed. Na-FAR-1 and possibly other FARs may have a wider repertoire for hydrophobic ligand binding, as confirmed in the present study by our finding that a range of neutral and polar lipids co-purify with the bacterially expressed recombinant protein. Finally, we show by immunohistochemistry that Na-FAR-1 is present in adult worms with a tissue distribution indicative of possible roles in nutrient acquisition by the parasite and in reproduction in the male.

Protein loop closure using orientational restraints from NMR data

Protein loops often play important roles in biological functions. Modeling loops accurately is crucial to determining the functional specificity of a protein. Despite the recent progress in loop prediction approaches, which led to a number of algorithms over the past decade, few rigorous algorithmic approaches exist to model protein loops using global orientational restraints, such as those obtained from residual dipolar coupling (RDC) data in solution nuclear magnetic resonance (NMR) spectroscopy. In this article, we present a novel, sparse data, RDC-based algorithm, which exploits the mathematical interplay between RDC-derived sphero-conics and protein kinematics, and formulates the loop structure determination problem as a system of low-degree polynomial equations that can be solved exactly, in closed-form. The polynomial roots, which encode the candidate conformations, are searched systematically, using provable pruning strategies that triage the vast majority of conformations, to enumerate or prune all possible loop conformations consistent with the data; therefore, completeness is ensured. Results on experimental RDC datasets for four proteins, including human ubiquitin, FF2, DinI, and GB3, demonstrate that our algorithm can compute loops with higher accuracy, a three- to six-fold improvement in backbone RMSD, versus those obtained by traditional structure determination protocols on the same data. Excellent results were also obtained on synthetic RDC datasets for protein loops of length 4, 8, and 12 used in previous studies. These results suggest that our algorithm can be successfully applied to determine protein loop conformations, and hence, will be useful in high-resolution protein backbone structure determination, including loops, from sparse NMR data. Proteins 2012. © 2011 Wiley Periodicals, Inc.

A Bayesian approach for determining protein side-chain rotamer conformations using unassigned NOE data

A major bottleneck in protein structure determination via nuclear magnetic resonance (NMR) is the lengthy and laborious process of assigning resonances and nuclear Overhauser effect (NOE) cross peaks. Recent studies have shown that accurate backbone folds can be determined using sparse NMR data, such as residual dipolar couplings (RDCs) or backbone chemical shifts. This opens a question of whether we can also determine the accurate protein side-chain conformations using sparse or unassigned NMR data. We attack this question by using unassigned nuclear Overhauser effect spectroscopy (NOESY) data, which records the through-space dipolar interactions between protons nearby in three-dimensional (3D) space. We propose a Bayesian approach with a Markov random field (MRF) model to integrate the likelihood function derived from observed experimental data, with prior information (i.e., empirical molecular mechanics energies) about the protein structures. We unify the side-chain structure prediction problem with the side-chain structure determination problem using unassigned NMR data, and apply the deterministic dead-end elimination (DEE) and A* search algorithms to provably find the global optimum solution that maximizes the posterior probability. We employ a Hausdorff-based measure to derive the likelihood of a rotamer or a pairwise rotamer interaction from unassigned NOESY data. In addition, we apply a systematic and rigorous approach to estimate the experimental noise in NMR data, which also determines the weighting factor of the data term in the scoring function derived from the Bayesian framework. We tested our approach on real NMR data of three proteins: the FF Domain 2 of human transcription elongation factor CA150 (FF2), the B1 domain of Protein G (GB1), and human ubiquitin. The promising results indicate that our algorithm can be applied in high-resolution protein structure determination. Since our approach does not require any NOE assignment, it can accelerate the NMR structure determination process.

Protein side-chain resonance assignment and NOE assignment using RDC-defined backbones without TOCSY data

One bottleneck in NMR structure determination lies in the laborious and time-consuming process of side-chain resonance and NOE assignments. Compared to the well-studied backbone resonance assignment problem, automated side-chain resonance and NOE assignments are relatively less explored. Most NOE assignment algorithms require nearly complete side-chain resonance assignments from a series of through-bond experiments such as HCCH-TOCSY or HCCCONH. Unfortunately, these TOCSY experiments perform poorly on large proteins. To overcome this deficiency, we present a novel algorithm, called NASCA: (NOE Assignment and Side-Chain Assignment), to automate both side-chain resonance and NOE assignments and to perform high-resolution protein structure determination in the absence of any explicit through-bond experiment to facilitate side-chain resonance assignment, such as HCCH-TOCSY. After casting the assignment problem into a Markov Random Field (MRF), NASCA: extends and applies combinatorial protein design algorithms to compute optimal assignments that best interpret the NMR data. The MRF captures the contact map information of the protein derived from NOESY spectra, exploits the backbone structural information determined by RDCs, and considers all possible side-chain rotamers. The complexity of the combinatorial search is reduced by using a dead-end elimination (DEE) algorithm, which prunes side-chain resonance assignments that are provably not part of the optimal solution. Then an A* search algorithm is employed to find a set of optimal side-chain resonance assignments that best fit the NMR data. These side-chain resonance assignments are then used to resolve the NOE assignment ambiguity and compute high-resolution protein structures. Tests on five proteins show that NASCA: assigns resonances for more than 90% of side-chain protons, and achieves about 80% correct assignments. The final structures computed using the NOE distance restraints assigned by NASCA: have backbone RMSD 0.8-1.5 Å from the reference structures determined by traditional NMR approaches.

One-bond and two-bond J couplings help annotate protein secondary-structure motifs: J-coupling indexing applied to human endoplasmic reticulum protein ERp18

NMR coupling constants, both direct one-bond ((1)J) and geminal two-bond ((2)J), are employed to analyze the protein secondary structure of human oxidized ERp18. Coupling constants collected and evaluated for the 18 kDa protein comprise 1268 values of (1)J(CαHα), (1)J(CαCβ), (1)J(CαC'), (1)J(C'N'), (1)J(N'Cα), (1)J(N') (HN), (2)J(CαN'), (2)J(HNCα), (2)J(C'HN), and (2)J(HαC'). Comparison with (1)J and (2)J data from reference proteins and pattern analysis on a per-residue basis permitted main-chain ϕ,ψ torsion-angle combinations of many of the 149 amino-acid residues in ERp18 to be narrowed to particular secondary-structure motifs. J-coupling indexing is here being developed on statistical criteria and used to devise a ternary grid for interpreting patterns of relative values of J. To account for the influence of the varying substituent pattern in different amino-acid sidechains, a table of residue-type specific threshold values was compiled for discriminating small, medium, and large categories of J. For the 15-residue insertion that distinguishes the ERp18 fold from that of thioredoxin, the J-coupling data hint at a succession of five isolated Type-I β turns at progressively shorter sequence intervals, in agreement with the crystal structure.

Unconventional ubiquitin recognition by the ubiquitin-binding motif within the Y family DNA polymerases iota and Rev1

Translesion synthesis is an essential cell survival strategy to promote replication after DNA damage. The accumulation of Y family polymerases (pol) iota and Rev1 at the stalled replication machinery is mediated by the ubiquitin-binding motifs (UBMs) of the polymerases and enhanced by PCNA monoubiquitination. We report the solution structures of the C-terminal UBM of human pol iota and its complex with ubiquitin. Distinct from other ubiquitin-binding domains, the UBM binds to the hydrophobic surface of ubiquitin centered at L8. Accordingly, mutation of L8A, but not I44A, of ubiquitin abolishes UBM binding. Human pol iota contains two functional UBMs, both contributing to replication foci formation. In contrast, only the second UBM of Saccharomyces cerevisiae Rev1 binds to ubiquitin and is essential for Rev1-dependent cell survival and mutagenesis. Point mutations disrupting the UBM-ubiquitin interaction also impair the accumulation of pol iota in replication foci and Rev1-mediated DNA damage tolerance in vivo.

High-resolution protein structure determination starting with a global fold calculated from exact solutions to the RDC equations

We present a novel structure determination approach that exploits the global orientational restraints from RDCs to resolve ambiguous NOE assignments. Unlike traditional approaches that bootstrap the initial fold from ambiguous NOE assignments, we start by using RDCs to compute accurate secondary structure element (SSE) backbones at the beginning of structure calculation. Our structure determination package, called RDC-PANDA: (RDC-based SSE PAcking with NOEs for Structure Determination and NOE Assignment), consists of three modules: (1) RDC-EXACT: ; (2) PACKER: ; and (3) HANA: (HAusdorff-based NOE Assignment). RDC-EXACT: computes the global optimal solution of backbone dihedral angles for each secondary structure element by exactly solving a system of quartic RDC equations derived by Wang and Donald (Proceedings of the IEEE computational systems bioinformatics conference (CSB), Stanford, CA, 2004a; J Biomol NMR 29(3):223-242, 2004b), and systematically searching over the roots, each of which is a backbone dihedral varphi- or psi-angle consistent with the RDC data. Using a small number of unambiguous inter-SSE NOEs extracted using only chemical shift information, PACKER: performs a systematic search for the core structure, including all SSE backbone conformations. HANA: uses a Hausdorff-based scoring function to measure the similarity between the experimental spectra and the back-computed NOE pattern for each side-chain from a statistically-diverse rotamer library, and drives the selection of optimal position-specific rotamers for filtering ambiguous NOE assignments. Finally, a local minimization approach is used to compute the loops and refine side-chain conformations by fixing the core structure as a rigid body while allowing movement of loops and side-chains. RDC-PANDA: was applied to NMR data for the FF Domain 2 of human transcription elongation factor CA150 (RNA polymerase II C-terminal domain interacting protein), human ubiquitin, the ubiquitin-binding zinc finger domain of the human Y-family DNA polymerase Eta (pol eta UBZ), and the human Set2-Rpb1 interacting domain (hSRI). These results demonstrated the efficiency and accuracy of our algorithm, and show that RDC-PANDA: can be successfully applied for high-resolution protein structure determination using only a limited set of NMR data by first computing RDC-defined backbones.

Variation in protein C(alpha)-related one-bond J couplings

Four types of polypeptide (1)J(C alpha X) couplings are examined, involving the main-chain carbon C(alpha) and either of four possible substituents. A total 3105 values of (1)J(C alpha H alpha), (1)J(C alpha C beta), (1)J(C alpha C'), and (1)J(C alpha N') were collected from six proteins, averaging 143.4 +/- 3.3, 34.9 +/- 2.5, 52.6 +/- 0.9, and 10.7 +/- 1.2 Hz, respectively. Analysis of variances (ANOVA) reveals a variety of factors impacting on (1)J and ranks their relative statistical significance and importance to biomolecular NMR structure refinement. Accordingly, the spread in the (1)J values is attributed, in equal proportions, to amino-acid specific substituent patterns and to polypeptide-chain geometry, specifically torsions phi, psi, and chi(1) circumjacent to C(alpha). The (1)J coupling constants correlate with protein secondary structure. For alpha-helical phi, psi combinations, (1)J(C alpha H alpha) is elevated by more than one standard deviation (147.8 Hz), while both (1)J(C alpha N') and (1)J(C alpha C beta) fall short of their grand means (9.5 and 33.7 Hz). Rare positive phi torsion angles in proteins exhibit concomitant small (1)J(C alpha H alpha) and (1)J(C alpha N') (138.4 and 9.6 Hz) and large (1)J(C alpha C beta) (39.9 Hz) values. The (1)J(C alpha N') coupling varies monotonously over the phi torsion range typical of beta-sheet secondary structure and is largest (13.3 Hz) for phi around -160 degrees. All four coupling types depend on psi and thus help determine a torsion that is notoriously difficult to assess by traditional approaches using (3)J. Influences on (1)J stemming from protein secondary structure and other factors, such as amino-acid composition, are largely independent.

Structure of the N-terminal region of complement factor H and conformational implications of disease-linked sequence variations

Factor H is a regulatory glycoprotein of the complement system. We expressed the three N-terminal complement control protein modules of human factor H (FH1-3) and confirmed FH1-3 to be the minimal unit with cofactor activity for C3b proteolysis by factor I. We reconstructed FH1-3 from NMR-derived structures of FH1-2 and FH2-3 revealing an ∼105-Å-long rod-like arrangement of the modules. In structural comparisons with other C3b-engaging proteins, factor H module 3 most closely resembles factor B module 3, consistent with factor H competing with factor B for binding C3b. Factor H modules 1, 2, and 3 each has a similar backbone structure to first, second, and third modules, respectively, of functional sites in decay accelerating factor and complement receptor type 1; the equivalent intermodular tilt and twist angles are also broadly similar. Resemblance between molecular surfaces is closest for first modules but absent in the case of second modules. Substitution of buried Val-62 with Ile (a factor H single nucleotide polymorphism potentially protective for age-related macular degeneration and dense deposit disease) causes rearrangements within the module 1 core and increases thermal stability but does not disturb the interface with module 2. Replacement of partially exposed (in module 1) Arg-53 by His (an atypical hemolytic uremic syndrome-linked mutation) did not impair structural integrity at 37 °C, but this FH1-2 mutant was less stable at higher temperatures; furthermore, chemical shift differences indicated potential for small structural changes at the module 1-2 interface.