Protein Backbone Dynamics Research Articles

Effects of Phosphorylation on Protein Backbone Dynamics and Conformational Preferences.

Phosphorylations are the most common and extensively studied post-translational modification (PTM) of proteins in eukaryotes. They constitute a major regulatory mechanism, modulating protein function, protein-protein interactions, as well as subcellular localization. Phosphorylation sites are preferably located in intrinsically disordered regions and have been shown to trigger structural rearrangements and order-to-disorder transitions. They can therefore have a significant effect on protein backbone dynamics or conformation, but only sparse experimental data are available. To obtain a more general description of how and when phosphorylations have a significant effect on protein behavior, molecular dynamics (MD) currently provides the only suitable framework to study these effects at a large scale in atomistic detail. This study develops a systematic MD simulation framework to explore the influence of phosphorylations on the local backbone dynamics and conformational propensities of proteins. Through a series of glycine-backbone peptides, we studied the effects of amino acid residues including the three most common phosphorylations (Ser, Thr, and Tyr), on local backbone dynamics and conformational propensities. We further extended our study to investigate the interactions of all such residues between position i to positions i + 1, i + 2, i + 3, and i + 4 in such peptides. The final data set comprises structural ensembles for 3393 sequences with more than 1 μs of sampling for each ensemble. To validate the relevance of the results, the structural and conformational properties extracted from the MD simulations are compared to NMR data from the Biological Magnetic Resonance Data Bank. The systematic nature of this study enables the projection of the gained knowledge onto any phosphorylation site in the proteome and provides a general framework for the study of further PTMs. The full data set is publicly available, as a training and reference set.

ARCHE-NOAH: NMR supersequence with five different CEST experiments for studying protein conformational dynamics.

An NMR NOAH-supersequence is presented consisting of five CEST experiments for studying protein backbone and side-chain dynamics by 15N-CEST, carbonyl-13CO-CEST, aromatic-13Car-CEST, 13Cα-CEST, and methyl-13Cmet-CEST. The new sequence acquires the data for these experiments in a fraction of the time required for the individual experiments, saving over four days of NMR time per sample.

Genome sequencing and molecular characterisation of XDR Acinetobacter baumannii reveal complexities in resistance: Novel combination of sulbactam-durlobactam holds promise for therapeutic intervention.

Emerging nosocomial strains of Acinetobacter baumannii are of recent concern as they are expressing extensive drug resistance (XDR). Using whole-genome sequencing and molecular characterisation analysis, the current study reveals the presence of carbapenemase genes in 92.86% of studied Indian isolates. These included blaOXA-51 , blaOXA-23 , blaOXA-58 , and blaNDM genes, with over a third expressing dual carbapenemase genes. As per the MLST scheme, IC2Oxf /CC2Pas was the predominant clone, with 57.14% isolates belonging to this lineage. The presence of these carbapenemase genes resulted in sulbactam (SUL) resistance (MIC: 16-256 µg/ml) in all of the studied isolates. The efficacy of durlobactam (DUR), a novel β-lactamase inhibitor that also inhibits PBP2 was assessed through in silico intermolecular interaction analysis. Several nonsynonymous single nucleotide polymorphisms were identified in PBP2 (G264S, I108V, S259T) and PBP3 (A515V, T526S) sequences. Minimal variations were recorded in the protein backbone dynamics in active-site motifs of wild-type and mutants, which correlated with negligible binding energy fluctuations for the PBP3-SUL (-5.85 ± 0.04 kcal/mol) and PBP2-DUR (-5.16 ± 0.66 kcal/mol) complexes. Furthermore, higher binding affinities and low inhibition constants were noted in OXA23-DUR (-7.36 kcal/mol; 4.01 µM), OXA58-DUR (-6.44 kcal/mol; 19.07 µM), and NDM-DUR (-6.82 kcal/mol; 10.01 µM) complexes when compared with the conventional drugs avibactam and aztreonam. Stable interaction profiles of DUR with carbapenemases can possibly restore SUL activity against both PBP3WT and PBP3MTs . The study establishes the efficacy of the novel SUL-DUR combination as a successful treatment strategy in combating emerging XDR strains of A. baumannii.

Inherent backbone dynamics fine-tune the functional plasticity of Tudor domains

Tudor domains are crucial for mediating a diversity of protein-protein or protein-DNA interactions involved in nucleic acid metabolism. Using solution NMR spectroscopy, we assess the comprehensive understanding of the dynamical properties of the respective Tudor domains from four different bacterial (Escherichia coli) proteins UvrD, Mfd, RfaH, and NusG involved in different aspects of bacterial transcription regulation and associated processes. These proteins are benchmarked to the canonical Tudor domain fold from the human SMN protein. The detailed analysis of protein backbone dynamics and subsequent analysis by the Lipari-Szabo model-free approach revealed subtle differences in motions of the amide-bond vector on both pico- to nanosecond and micro- to millisecond timescales. On these timescales, our comparative approach reveals the usefulness of discrete amplitudes of dynamics to discern the different functionalities for Tudor domains exhibiting promiscuous binding, including the metamorphic Tudor domain included in the study.

Early Strides in NMR Dynamics Measurements.

The study of protein dynamics using the measurement of relaxation times by NMR was based on a set of studies in the mid-20th century that outlined theories and methods. However, the complexity of protein NMR was such that these simple experiments were not practical for application to proteins. The advent of techniques in the 1980s for isotopic labeling of proteins meant that pulse sequences could now be applied in multidimensional NMR experiments to enable per-residue information about the local relaxation times. One of the earliest advances was published in Biochemistry in 1989. The paper "Backbone dynamics of proteins as studied by 15N inverse detected heteronuclear NMR spectroscopy: application to staphylococcal nuclease" by Lewis Kay, Dennis Torchia, and Ad Bax delineated a set of pulse sequences that are used with minor modifications even today. This paper, with others from a limited number of other laboratories, forms the basis for the experimental determination of the backbone dynamics of proteins. The biological insights obtained from such measurements have only increased in the past 30 years. Sometimes, the best and perhaps only way to advance a field is an advancement in the technical capabilities that allows new perspectives to be reached.

Going the dHis-tance: Site-Directed Cu2+ Labeling of Proteins and Nucleic Acids.

In this Account, we showcase site-directed Cu2+ labeling in proteins and DNA, which has opened new avenues for the measurement of the structure and dynamics of biomolecules using electron paramagnetic resonance (EPR) spectroscopy. In proteins, the spin label is assembled in situ from natural amino acid residues and a metal complex and requires no post-expression synthetic modification or purification procedures. The labeling scheme exploits a double histidine (dHis) motif, which utilizes endogenous or site-specifically mutated histidine residues to coordinate a Cu2+ complex. Pulsed EPR measurements on such Cu2+-labeled proteins potentially yield distance distributions that are up to 5 times narrower than the common protein spin label-the approach, thus, overcomes the inherent limitation of the current technology, which relies on a spin label with a highly flexible side chain. This labeling scheme provides a straightforward method that elucidates biophysical information that is costly, complicated, or simply inaccessible by traditional EPR labels. Examples include the direct measurement of protein backbone dynamics at β-sheet sites, which are largely inaccessible through traditional spin labels, and rigid Cu2+-Cu2+ distance measurements that enable higher precision in the analysis of protein conformations, conformational changes, interactions with other biomolecules, and the relative orientations of two labeled protein subunits. Likewise, a Cu2+ label has been developed for use in DNA, which is small, is nucleotide independent, and is positioned within the DNA helix. The placement of the Cu2+ label directly reports on the biologically relevant backbone distance. Additionally, for both of these labeling techniques, we have developed models for interpretation of the EPR distance information, primarily utilizing molecular dynamics (MD) simulations. Initial results using force fields developed for both protein and DNA labels have agreed with experimental results, which has been a major bottleneck for traditional spin labels. Looking ahead, we anticipate new combinations of MD and EPR to further our understanding of protein and DNA conformational changes, as well as working synergistically to investigate protein-DNA interactions.

Backbone resonance assignment and dynamics of 110 kDa hexameric inorganic pyrophosphatase from Mycobacterium tuberculosis.

Family I soluble inorganic pyrophosphatases (PPases; EC 3.6.1.1) are enzymes essential for all organisms. They hydrolyze inorganic pyrophosphate, thus providing the driving force for numerous biosynthetic reactions. Soluble PPases retain enzymatic activity only in multimeric forms. PPases from various organisms are extensively studied by X-ray crystallography but until now there was no information on their structure and dynamics in solution. Hexameric 110 kDa (6×18.3 kDa) PPase from Mycobacterium tuberculosis (Mt-PPase) is a promising target for the rational design of potential anti-tuberculosis agents. In order to use NMR techniques in functional studies of Mt-PPase and rational design of the inhibitors for this enzyme, it is necessary to have information on the backbone 1H, 13C and 15N resonance assignments. Samples of Mt-PPase enriched with 99% of 13C and 15N isotopes, and 95% of 2H were obtained using recombinant protein expression in an isotopically-labeled medium and effective heat-shock protocol for the deuterium-to-hydrogen exchange of the amide groups. Backbone resonance assignment was achieved for more than 95% of the residues. It was found that the secondary structure of Mt-PPase in solution corresponds well to the crystal structure of this protein. Protein backbone dynamics were studied using 15N NMR relaxation experiments. Determined resonance assignments and dynamic properties provide the basis for the subsequent structure-based design of novel inhibitors of Mt-PPase-potential anti-tuberculosis drugs.

Alpha-Helical Membrane Protein Folding in “Mixed” and “Ideal” Bicelles

Lipid-detergent bicelles are often used in membrane protein studies as an alternative to micelles. The ideal bicelle model contains a segregated lipid core surrounded by a detergent rim, and their morphologies depend on the ratio of lipid to detergent (q-ratio). There are several reports that the most widely used bicelle, containing 1,2-dimyristoyl-sn-glycero-3-phosphocholine and 1,2-dihexanoyl-sn-glycero-3-phosphocholine, forms ideally mixed micelles at q-ratios below 0.5. To date, almost all membrane protein studies using isotropic bicelles, particularly solution NMR spectroscopy, have relied on q-ratios of ∼ 0.3 due to favorable tumbling properties. Designing bicelle mixtures with segregated lipid and detergent across a broad range of q-ratios is, therefore, of high interest for membrane protein applications. We investigated common bicelle mixtures with q-ratios of 0.0-1.0 and measured the short diameter, L, with SAXS. The diameter was previously shown to correspond with the mixing behavior reported by SANS and molecular dynamics simulations, with L equal to two lipid lengths in a lipid-segregated bicelle. All alkyl-chain lipid and detergent combinations with phosphocholine-like head groups form mixed micelles up to ∼ q-ratio = 0.5. Bicelles containing bulky moieties at detergent alkyl tails and head groups were additionally investigated, which we hypothesize will increase segregation between lipid and detergent components at q ratios less than 0.5. The mixtures most promising for ideal bicelle morphologies at low q-ratios were then used to investigate protein fold and backbone dynamics of a polytopic alpha-helical transmembrane protein using continuous-wave EPR and solution NMR spectroscopy. Our results suggest that lipid-protein interactions influence protein fold distinct from detergent-protein interactions.

Protein Backbone and Average Particle Dynamics in Reconstituted Discoidal and Spherical HDL Probed by Hydrogen Deuterium Exchange and Elastic Incoherent Neutron Scattering

Lipoproteins are supramolecular assemblies of proteins and lipids with dynamic characteristics critically linked to their biological functions as plasma lipid transporters and lipid exchangers. Among them, spherical high-density lipoproteins are the most abundant forms of high-density lipoprotein (HDL) in human plasma, active participants in reverse cholesterol transport, and associated with reduced development of atherosclerosis. Here, we employed elastic incoherent neutron scattering (EINS) and hydrogen-deuterium exchange mass spectrometry (HDX-MS) to determine the average particle dynamics and protein backbone local mobility of physiologically competent discoidal and spherical HDL particles reconstituted with human apolipoprotein A-I (apoA-I). Our EINS measurements indicated that discoidal HDL was more dynamic than spherical HDL at ambient temperatures, in agreement with their lipid-protein composition. Combining small-angle neutron scattering (SANS) with contrast variation and MS cross-linking, we showed earlier that the most likely organization of the three apolipoprotein A-I (apoA-I) chains in spherical HDL is a combination of a hairpin monomer and a helical antiparallel dimer. Here, we corroborated those findings with kinetic studies, employing hydrogen-deuterium exchange mass spectrometry (HDX-MS). Many overlapping apoA-I digested peptides exhibited bimodal HDX kinetics behavior, suggesting that apoA-I regions with the same amino acid composition located on different apoA-I chains had different conformations and/or interaction environments.

Catcher/Tag cyclization introduces electrostatic interaction mediated protein-protein interactions to enhance the thermostability of luciferase

Efficiently engineering proteins to enhance their thermostability and other features hinges on the correct understanding and prediction of the effects of protein backbone topology and dynamics. Here, we fused a Catcher module to the C-terminus of luciferase and a Tag module to the N-terminus, so that the spontaneous reaction between the Catcher and Tag covalently cyclized the luciferase. The thermal stability of cyclized luciferase improved, with T50, t1/2, and Tm values higher by 8.9 °C, 2.1-fold and 7.3 °C respectively compared with the linear enzyme. Furthermore, MD simulations suggested that the cyclizations resulted in a tight docking between Catcher/ Tag and luciferase. A series of interactions on the interface between Catcher/Tag and the C-domain of luciferase were identified. Compared with the uncyclized enzyme, these interactions appear to impose remarkable changes to the backbone dynamics of the enzyme. Our data showed that the cyclization increased the gyration radii, and decreased both RMSD values and Cα-fluctuations of the enzyme, which was consistent with experimental observations and indicates that the cyclized luciferase adopted a more stable conformation enssemble. Taken together, this study provides a theoretical understanding of how Catcher/Tag cyclization can enhance the thermal resilience of enzymes, and might be useful in guiding protein design.

Accurate calculation of side chain packing and free energy with applications to protein molecular dynamics.

To address the large gap between time scales that can be easily reached by molecular simulations and those required to understand protein dynamics, we present a rapid self-consistent approximation of the side chain free energy at every integration step. In analogy with the adiabatic Born-Oppenheimer approximation for electronic structure, the protein backbone dynamics are simulated as preceding according to the dictates of the free energy of an instantaneously-equilibrated side chain potential. The side chain free energy is computed on the fly, allowing the protein backbone dynamics to traverse a greatly smoothed energetic landscape. This computation results in extremely rapid equilibration and sampling of the Boltzmann distribution. Our method, termed Upside, employs a reduced model involving the three backbone atoms, along with the carbonyl oxygen and amide proton, and a single (oriented) side chain bead having multiple locations reflecting the conformational diversity of the side chain’s rotameric states. We also introduce a novel, maximum-likelihood method to parameterize the side chain interactions using protein structures. We demonstrate state-of-the-art accuracy for predicting χ1 rotamer states while consuming only milliseconds of CPU time. Our method enables rapidly equilibrating coarse-grained simulations that can nonetheless contain significant molecular detail. We also show that the resulting free energies of the side chains are sufficiently accurate for de novo folding of some proteins.

Lipid- and Cholesterol-Mediated Time-Scale-Specific Modulation of the Outer Membrane Protein X Dynamics in Lipid Bilayers.

Membrane protein function fundamentally depends on lipid-bilayer fluidity and the composition of the biological membrane. Although dynamic interdependencies between membrane proteins and the surrounding lipids are suspected, a detailed description is still missing. To uncover lipid-modulated membrane protein backbone dynamics, time-scale-specific NMR relaxation experiments with residue-resolution were recorded. The data revealed that lipid order, modified either biochemically or biophysically, changes the dynamics of the immersed membrane protein in a specific and time-scale-dependent manner. A temperature-dependent dynamics analysis furthermore suggests a direct coupling between lipid and protein dynamics in the picosecond-nanosecond, microsecond, and millisecond time scales, caused by the lipid's trans-gauche isomerization, the segmental and rotational motion of lipids, and the fluidity of the lipid phase, respectively. These observations provide evidence of a direct modulatory capability of the membrane to regulate protein function through lipid dynamics ranging from picoseconds to milliseconds.

NMR assignments of the WBSCR27 protein related to Williams-Beuren syndrome.

Williams-Beuren syndrome is a genetic disorder characterized by physiological and mental abnormalities, and is caused by hemizygous deletion of several genes in chromosome 7. One of the removed genes encodes the WBSCR27 protein. Bioinformatic analysis of the sequence of WBSCR27 indicates that it belongs to the family of SAM-dependent methyltransferases. However, exact cellular functions of this protein or phenotypic consequences of its deficiency are still unknown. Here we report nearly complete 1H, 15N, and 13C chemical shifts assignments of the 26kDa WBSCR27 protein from Mus musculus in complex with the cofactor S-adenosyl-L-methionine (SAM). Analysis of the assigned chemical shifts allowed us to characterize the protein's secondary structure and backbone dynamics. The topology of the protein's fold confirms the assumption that the WBSCR27 protein belongs to the family of class I methyltransferases.

Dynamics of dehaloperoxidase-hemoglobin A derived from NMR relaxation spectroscopy and molecular dynamics simulation

Dehaloperoxidase-hemoglobin is the first hemoglobin identified with biologically-relevant oxidative functions, which include peroxidase, peroxygenase and oxidase activities. Herein we report a study of the protein backbone dynamics of DHP using heteronuclear NMR relaxation methods and molecular dynamics (MD) simulations to address the role of protein dynamics in switching from one function to another. The results show that DHP's backbone helical regions and turns have average order parameters of S2 = 0.87 ± 0.03 and S2 = 0.76 ± 0.08, respectively. Furthermore, DHP is primarily a monomer in solution based on the overall tumbling correlation time τm is 9.49 ± 1.65 ns calculated using the prolate diffusion tensor model in the program relax. A number of amino acid residues have significant Rex using the Lipari-Szabo model-free formalism. These include Lys3, Ile6, Leu13, Gln18, Arg32, Ser48, Met49, Thr56, Phe60, Arg69, Thr71 Cys73, Ala77, Asn81, Gly95, Arg109, Phe115, Leu127 and Met136, which may experience slow conformational motions on the microseconds-milliseconds time scale according to the model. Caution should be used when the model contains >4 fitting parameters. The program caver3.0 was used to identify tunnels inside DHP obtained from MD simulation snapshots that are consistent with the importance of the Xe binding site, which is located at the central intersection of the tunnels. These tunnels provide diffusion pathways for small ligands such as O2, H2O and H2O2 to enter the distal pocket independently of the trajectory of substrates and inhibitors, both of which are aromatic molecules.

Binding and backbone dynamics of protein under topological constraint: calmodulin as a model system

Herein we present the effect of artificially imposed topological constraint on calmodulin (CaM) backbone dynamics and its molecular recognition behavior. While backbone dynamics of CaM remain largely unperturbed, the thermodynamic profile of CaM binding to the smooth-muscle myosin light-chain kinase (smMLCK) peptide is modulated significantly.

Characterization of Bifunctional Spin Labels for Investigating the Structural and Dynamic Properties of Membrane Proteins Using EPR Spectroscopy.

Site-directed spin labeling (SDSL) coupled with electron paramagnetic resonance (EPR) spectroscopy is a very powerful technique to study structural and dynamic properties of membrane proteins. The most widely used spin label is methanthiosulfonate (MTSL). However, the flexibility of this spin label introduces greater uncertainties in EPR measurements obtained for determining structures, side-chain dynamics, and backbone motion of membrane protein systems. Recently, a newer bifunctional spin label (BSL), 3,4-bis(methanethiosulfonylmethyl)-2,2,5,5-tetramethyl-2,5-dihydro-1H-pyrrol-1-yloxy, has been introduced to overcome the dynamic limitations associated with the MTSL spin label and has been invaluable in determining protein backbone dynamics and inter-residue distances due to its restricted internal motion and fewer size restrictions. While BSL has been successful in providing more accurate information about the structure and dynamics of several proteins, a detailed characterization of the spin label is still lacking. In this study, we characterized BSLs by performing CW-EPR spectral line shape analysis as a function of temperature on spin-labeled sites inside and outside of the membrane for the integral membrane protein KCNE1 in POPC/POPG lipid bilayers and POPC/POPG lipodisq nanoparticles. The experimental data revealed a powder pattern spectral line shape for all of the KCNE1-BSL samples at 296 K, suggesting the motion of BSLs approaches the rigid limit regime for these series of samples. BSLs were further utilized to report for the first time the distance measurement between two BSLs attached on an integral membrane protein KCNE1 in POPC/POPG lipid bilayers at room temperature using dipolar line broadening CW-EPR spectroscopy. The CW dipolar line broadening EPR data revealed a 15 ± 2 Å distance between doubly attached BSLs on KCNE1 (53/57-63/67) which is consistent with molecular dynamics modeling and the solution NMR structure of KCNE1 which yielded a distance of 17 Å. This study demonstrates the utility of investigating the structural and dynamic properties of membrane proteins in physiologically relevant membrane mimetics using BSLs.

Microsecond Dynamics in Ubiquitin Probed by Solid-State 15 N NMR Spectroscopy R1ρ Relaxation Experiments under Fast MAS (60-110 kHz).

15 N R1ρ relaxation experiments in solid-state NMR spectroscopy are sensitive to timescales and amplitudes of internal protein motions in the hundreds of nano- to microsecond time window, which is difficult to probe by solution-state NMR spectroscopy. By using 15 N R1ρ relaxation experiments, a simplified approach to detect low microsecond protein dynamics is described and residue-specific correlation times are determined from the ratio of 15 N R1ρ rate constants at different magic angle spinning frequencies. Microcrystalline ubiquitin exhibits small-amplitude dynamics on a timescale of about 1 μs across the entire protein, and larger amplitude motions, also on the 1 μs timescale, for several sites, including the β1 -β2 turn and the N terminus of the α helix. According to the analysis, the microsecond protein backbone dynamics are of lower amplitude than that concluded in previous solid-state NMR spectroscopy studies, but persist across the entire protein with a rather uniform timescale of 1 μs.

Evolution of P450 Monooxygenases toward Formation of Transient Channels and Exclusion of Nonproductive Gases

Gaseous molecules are essential for biological processes, yet mapping the migration of gas molecules into and out of proteins represents a significant challenge. Cytochrome P450 enzymes (P450s) contain numerous channels thought to be populated by their substrates, products, solvents, and gases, yet the principles underlying channel preference are unknown. We identified multiple putative ligand migration channels of two bacterial P450s, CYP102A1 (BM3) from Bacillus megaterium and CYP102A5 from Bacillus cereus, using implicit ligand sampling and free molecular dynamics simulations and furthermore characterized the energy of gas migration through each. We observed strong discrimination between preferred gas migration channels, previously identified substrate/product migration channels and water channels, and mapped putative O2 reservoirs in the enzyme core. The protein backbone dynamics (S2 order parameter) unexpectedly revealed that some channels are transient in nature, with subchannels forming and merging...

Direct Observation of the Intrinsic Backbone Torsional Mobility of Disordered Proteins

The fundamental backbone dynamics of unfolded proteins arising due to intrinsic ϕ-ψ dihedral angle fluctuations dictate the course of protein folding, binding, assembly, and function. These internal fluctuations are also critical for protein misfolding associated with a range of human diseases. However, direct observation and unambiguous assignment of this inherent dynamics in chemically denatured proteins is extremely challenging due to various experimental limitations. To directly map the backbone torsional mobility in the ϕ-ψ dihedral angle space, we used a model intrinsically disordered protein, namely, α-synuclein, that adopts an expanded state under native conditions. We took advantage of nonoccurrence of tryptophan in α-synuclein and created a number of single-tryptophan variants encompassing the entire polypeptide chain. We then utilized highly sensitive picosecond time-resolved fluorescence depolarization measurements that allowed us to discern the site-specific torsional relaxation at a low protein concentration under physiological conditions. For all the locations, the depolarization kinetics exhibited two well-separated rotational-correlation-time components. The shorter, subnanosecond component arises due to the local mobility of the indole side chain, whereas the longer rotational-correlation-time component (1.37 ± 0.15 ns), independent of global tumbling, represents a characteristic timescale for short-range conformational exchange in the ϕ-ψ dihedral space. This correlation time represents an intrinsic timescale for torsional relaxation and is independent of position, which is expected for an extended polypeptide chain having little or no propensity to form persistent structures. We were also able to capture this intrinsic timescale at the N-terminal unstructured domain of the prion protein. Our estimated timescale of the segmental mobility is similar to that of unfolded proteins studied by nuclear magnetic resonance in conjunction with molecular dynamics simulations. Our results have broader implications for a diverse range of functionally and pathologically important intrinsically disordered proteins and disordered regions.

NMR Backbone Assignment of Large Proteins by Using 13Cα‐Only Triple‐Resonance Experiments

Nuclear magnetic resonance (NMR) is a powerful tool to interrogate protein structure and dynamics residue by residue. However, the prerequisite chemical-shift assignment remains a bottleneck for large proteins due to the fast relaxation and the frequency degeneracy of the (13) Cα nuclei. Herein, we present a covariance NMR strategy to assign the backbone chemical shifts by using only HN(CO)CA and HNCA spectra that has a high sensitivity even for large proteins. By using the peak linear correlation coefficient (LCC), which is a sensitive probe even for tiny chemical-shift displacements, we correctly identify the fidelity of approximately 92 % cross-peaks in the covariance spectrum, which is thus a significant improvement on the approach developed by Snyder and Brüschweiler (66 %) and the use of spectral derivatives (50 %). Thus, we calculate the 4D covariance spectrum from HN(CO)CA and HNCA experiments, in which cross-peaks with LCCs above a universal threshold are considered as true correlations. This 4D covariance spectrum enables the sequential assignment of a 42 kDa maltose binding protein (MBP), in which about 95 % residues are successfully assigned with a high accuracy of 98 %. Our LCC approach, therefore, paves the way for a residue-by-residue study of the backbone structure and dynamics of large proteins.