Chemical Shifts Of Backbone Atoms Research Articles

Negative thermal expansion of a disordered native protein

Despite immense interest in thermal expansion of solid materials, negative thermal expansion in particular, expansivity of proteins near ambient temperature is minimally reported. One reason for the paucity of protein expansion studies could be the susceptibility of proteins to conformational changes often within subdenaturing temperatures, given that the expansivity phenomenon should include only vibrational fluctuations in atom positions in the absence of a conformational change. This study encounters a disordered plant protein called At PP16-1 ( Arabidopsis thaliana phloem protein type 16-1) which shows negative thermal expansion in the 292–309(±1) K range of temperature at no expense of conformational transition. Temperature dependence of NMR chemical shifts of backbone atoms, global tumbling time derived from 15 N relaxations, and volume compressibility measurements yield a negative linear expansion coefficient 〈α〉 ∼ −5.3 × 10 −3 K −1 averaged over a sizable set of hydrogen bonds, and a decrease in the root mean square volume fluctuation by ∼36%. A frequency spectrum of normal modes in which one or a fewer low frequency vibrational modes of large fluctuation amplitudes are excited at the lower end, and many high-frequency vibrational modes of lower fluctuation amplitudes are excited at the higher end of the temperature range can explain negative thermal expansion.

Statistical Model To Decipher Protein Folding/Unfolding at a Local Scale.

Protein folding/unfolding can be analyzed experimentally at a local scale by monitoring the physical properties of local probes as a function of the temperature, for example, the distance between fluorophores or the values of chemical shifts of backbone atoms. Here, the analytical Lifson-Roig model for the helix-coil transition is modified to analyze local thermal unfolding of the fast-folder W protein of bacteriophage lambda (gpW) simulated by all-atom molecular dynamics (MD) simulations in explicit solvent at 15 different temperatures. The protein structure is described by the coarse-grained dihedral angles (γ) and bond angles (θ) built between successive Cα-Cα virtual bonds. Each (γ,θ) pair serves as a local probe of protein unfolding. Local native/non-native states are defined for each pair of (γ,θ) angles by analyzing the free-energy landscapes Δ G(γ,θ) computed from MD trajectories. The three local elementary equilibrium constants of the model are extracted for each (γ,θ) pair along the sequence from MD simulations, and the model predictions are compared to the MD data. Using only the local equilibrium constants as an input, we show that the local denaturation curves computed from the model partition function fit their MD simulated counterparts in a satisfying manner without any adjustment. In the model and MD simulations, gpW unfolds gradually between 320 and 340 K, with an average native percentage decreasing from 0.8 (320 K) to 0.2 (340 K). In the prism of the model, there is no stable structure at the local scale in this 20 K unfolding temperature range. The enthalpy change upon local unfolding computed from the model and from MD trajectories suggests that the unfolded state between 320 and 340 K corresponds to a dynamical equilibrium between a large ensemble of constantly changing structures. The present results confirm the downhill unfolding of gpW, which does not obey a two-state global folding/unfolding model, and shed light on the interpretation of local denaturation curves.

Refolded recombinant Siglec5 for NMR investigation of complex carbohydrate binding

Sialic-acid-binding immunoglobulin-like lectin (Siglec5) is a carbohydrate-binding surface receptor expressed on neutrophils, monocytes and B cells in human lymphoid and myeloid cell lineages. Existing structural and functional data fail to define the clear ligand specificity of Siglec5, though like other Siglec family members, it binds a variety of complex carbohydrates containing a sialic acid at the non-reducing terminus. Prokaryotic expression of this protein has proven challenging due to disulfide bonds and Asn-linked glycosylation. We developed an expression and purification protocol that uses an on-column strategy to refold Escherichia coli expressed protein that produced a high yield (2mg/L) of the single N-terminal Siglec5 carbohydrate recognition domain (CRD). A 2D heteronuclear single-quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectrum showed this material was folded, and a secondary structure prediction based on the assigned chemical shifts of backbone atoms was consistent with a previously determined X-ray model. NMR chemical shift mapping of Siglec5 binding to three carbohydrate ligands revealed similarities in binding interfaces and affinities. In addition, the role of alternate protein conformations identified by NMR in ligand binding is discussed. These studies demonstrate the Siglec5 CRD alone is sufficient for binding sialylated carbohydrates and provide a foundation for further investigation of Siglec5 structure and function.

NMR studies of the C-terminus of alpha4 reveal possible mechanism of its interaction with MID1 and protein phosphatase 2A.

Alpha4 is a regulatory subunit of the protein phosphatase family of enzymes and plays an essential role in regulating the catalytic subunit of PP2A (PP2Ac) within the rapamycin-sensitive signaling pathway. Alpha4 also interacts with MID1, a microtubule-associated ubiquitin E3 ligase that appears to regulate the function of PP2A. The C-terminal region of alpha4 plays a key role in the binding interaction of PP2Ac and MID1. Here we report on the solution structure of a 45-amino acid region derived from the C-terminus of alpha4 (alpha45) that binds tightly to MID1. In aqueous solution, alpha45 has properties of an intrinsically unstructured peptide although chemical shift index and dihedral angle estimation based on chemical shifts of backbone atoms indicate the presence of a transient α-helix. Alpha45 adopts a helix-turn-helix HEAT-like structure in 1% SDS micelles, which may mimic a negatively charged surface for which alpha45 could bind. Alpha45 binds tightly to the Bbox1 domain of MID1 in aqueous solution and adopts a structure consistent with the helix-turn-helix structure observed in 1% SDS. The structure of alpha45 reveals two distinct surfaces, one that can interact with a negatively charged surface, which is present on PP2A, and one that interacts with the Bbox1 domain of MID1.

Recognition and Specificity Determinants of the Human Cbx Chromodomains

The eight mammalian Cbx proteins are chromodomain-containing proteins involved in regulation of heterochromatin, gene expression, and developmental programs. They are evolutionarily related to the Drosophila HP1 (dHP1) and Pc (dPc) proteins that are key components of chromatin-associated complexes capable of recognizing repressive marks such as trimethylated Lys-9 and Lys-27, respectively, on histone H3. However, the binding specificity and function of the human homologs, Cbx1–8, remain unclear. To this end we employed structural, biophysical, and mutagenic approaches to characterize the molecular determinants of sequence contextual methyllysine binding to human Cbx1–8 proteins. Although all three human HP1 homologs (Cbx1, -3, -5) replicate the structural and binding features of their dHP counterparts, the five Pc homologs (Cbx2, -4, -6, -7, -8) bind with lower affinity to H3K9me3 or H3K27me3 peptides and are unable to distinguish between these two marks. Additionally, peptide permutation arrays revealed a greater sequence tolerance within the Pc family and suggest alternative nonhistone sequences as potential binding targets for this class of chromodomains. Our structures explain the divergence of peptide binding selectivity in the Pc subfamily and highlight previously unrecognized features of the chromodomain that influence binding and specificity.

A Bayesian-probability-based method for assigning protein backbone dihedral angles based on chemical shifts and local sequences

Chemical shifts contain substantial information about protein local conformations. We present a method to assign individual protein backbone dihedral angles into specific regions on the Ramachandran map based on the amino acid sequences and the chemical shifts of backbone atoms of tripeptide segments. The method uses a scoring function derived from the Bayesian probability for the central residue of a query tripeptide segment to have a particular conformation. The Ramachandran map is partitioned into representative regions at two levels of resolution. The lower resolution partitioning is equivalent to the conventional definitions of different secondary structure regions on the map. At the higher resolution level, the alpha and beta regions are further divided into subregions. Predictions are attempted at both levels of resolution. We compared our method with TALOS using the original TALOS database, and obtained comparable results. Although TALOS may produce the best results with currently available databases which are much enlarged, the Bayesian-probability-based approach can provide a quantitative measure for the reliability of predictions.

Protein backbone angle restraints from searching a database for chemical shift and sequence homology.

Chemical shifts of backbone atoms in proteins are exquisitely sensitive to local conformation, and homologous proteins show quite similar patterns of secondary chemical shifts. The inverse of this relation is used to search a database for triplets of adjacent residues with secondary chemical shifts and sequence similarity which provide the best match to the query triplet of interest. The database contains 13C alpha, 13C beta, 13C', 1H alpha and 15N chemical shifts for 20 proteins for which a high resolution X-ray structure is available. The computer program TALOS was developed to search this database for strings of residues with chemical shift and residue type homology. The relative importance of the weighting factors attached to the secondary chemical shifts of the five types of resonances relative to that of sequence similarity was optimized empirically. TALOS yields the 10 triplets which have the closest similarity in secondary chemical shift and amino acid sequence to those of the query sequence. If the central residues in these 10 triplets exhibit similar phi and psi backbone angles, their averages can reliably be used as angular restraints for the protein whose structure is being studied. Tests carried out for proteins of known structure indicate that the root-mean-square difference (rmsd) between the output of TALOS and the X-ray derived backbone angles is about 15 degrees. Approximately 3% of the predictions made by TALOS are found to be in error.

Lipid Binding Ridge on Loops 2 and 3 of the C2A Domain of Synaptotagmin I as Revealed by NMR Spectroscopy

The C2A domain of synaptotagmin I, which binds Ca2+ and anionic phospholipids, serves as a Ca2+ sensor during excitation-secretion coupling. We have used multidimensional NMR to locate the region of C2A from rat synaptotagmin I that interacts, in the presence of Ca2+, with phosphatidylserine. Untagged, recombinant C2A was double-labeled with 13C and 15N, and triple-resonance NMR data were collected from C2A samples containing either Ca2+ alone or Ca2+ plus 6:0 phosphatidylserine. Phospholipid binding led to changes in chemical shifts of backbone atoms in residues Arg233 and Phe234 of loop 3 (a loop that also binds Ca2+) and His198, Val205, and Phe206 of loop 2. These residues lie along a straight line on a surface ridge of the C2A domain. The only other residue that exhibited appreciable chemical shift changes upon adding lipid was His254; however, because His254 is located on the other side of the molecule from the phospholipid docking site defined by the other residues, its shifts may result from nonspecific interactions. The results show that the "docking ridge" responsible for Ca2+-dependent membrane association is localized on the opposite side of the C2A domain from the transmembrane and C2B domains of synaptotagmin.

NMR-based structural studies of the pNR-2/pS2 single domain trefoil peptide. Similarities to porcine spasmolytic peptide and evidence for a monomeric structure.

NMR spectroscopy measurements have been used to obtain structural information about the pNR-2/pS2 single-domain trefoil peptide. NMR data from 2D (two dimensional) double-quantum-filtered correlation spectroscopy (DQF-COSY), total correlation spectroscopy (TOCSY), NOE spectroscopy (NOESY), rotating frame NOE spectroscopy (ROESY) and 2D 13C-1H heteronuclear single-quantum coherence (HSQC) and 13C-1H HSQC-TOCSY spectra have been analysed to provide essentially complete 1H and 13C sequence-specific assignments for the pNR-2/pS2 protein. From a consideration of the NOE intensities, 3J(NH-alpha CH) coupling constants, 1H and 13C chemical shifts of backbone atoms and amide-proton exchange rates, the pNR-2/pS2 was found to contain two short antiparallel beta-strands (32-35 and 43-46), a short helix (25-30) and a type I beta-turn (11-15). These elements of secondary structure are very similar to those found in the two trefoil domains of pSP for which detailed structural information is already available. Similar 1H chemical shifts were noted for several conserved residues in pNR-2/pS2 and pSP and a characteristic Phe residue with a slowly flipping ring was found in the pNR-2/pS2 variant and in both domains of pSP. The tertiary structures of the domains therefore appear to be very similar in the two proteins and it is likely that the pNR-2/pS2 has the same pattern of disulphide bonds (1-5, 2-4, 3-6) as pSP. Correlation time measurements derived from 1H-1H NOE measurements indicate that the Cys58-->Ser form of the pNR-2/pS2 protein used in this study is monomeric in solution at approximately 2 mM.