Empirical correlation between protein backbone conformation and C.alpha. and C.beta. 13C nuclear magnetic resonance chemical shifts

NMR Methods for the Study of Instrinsically Disordered Proteins Structure, Dynamics, and Interactions: General Overview and Practical Guidelines

Thanks to recent improvements in NMR instrumentation, pulse sequence design, and sample preparation, a panoply of new NMR tools has become available for atomic resolution characterization of intrinsically disordered proteins (IDPs) that are optimized for the particular chemical and spectroscopic properties of these molecules. A wide range of NMR observables can now be measured on increasingly complex IDPs that report on their structural and dynamic properties in isolation, as part of a larger complex, or even inside an entire living cell. Herein we present basic NMR concepts, as well as optimised tools available for the study of IDPs in solution. In particular, the following sections are discussed hereafter: a short introduction to NMR spectroscopy and instrumentation (Sect. 3.1), the effect of order and disorder on NMR observables (Sect. 3.2), particular challenges and bottlenecks for NMR studies of IDPs (Sect. 3.3), 2D HN and CON NMR experiments: the fingerprint of an IDP (Sect. 3.4), tools for overcoming major bottlenecks of IDP NMR studies (Sect. 3.5), 13C detected experiments (Sect. 3.6), from 2D to 3D: from simple snapshots to site-resolved characterization of IDPs (Sect. 3.7), sequential NMR assignment: 3D experiments (Sect. 3.8), high-dimensional NMR experiments (nD, with n>3) (Sect. 3.9) and conclusions and perspectives (Sect. 3.10).

Protein structural information derived from NMR chemical shift with the neural network program TALOS-N

Chemical shifts are obtained at the first stage of any protein structural study by NMR spectroscopy. Chemical shifts are known to be impacted by a wide range of structural factors, and the artificial neural network based TALOS-N program has been trained to extract backbone and side-chain torsion angles from (1)H, (15)N, and (13)C shifts. The program is quite robust and typically yields backbone torsion angles for more than 90 % of the residues and side-chain χ 1 rotamer information for about half of these, in addition to reliably predicting secondary structure. The use of TALOS-N is illustrated for the protein DinI, and torsion angles obtained by TALOS-N analysis from the measured chemical shifts of its backbone and (13)C(β) nuclei are compared to those seen in a prior, experimentally determined structure. The program is also particularly useful for generating torsion angle restraints, which then can be used during standard NMR protein structure calculations.

Structure calculation, refinement and validation using CcpNmr Analysis

CcpNmr Analysis provides a streamlined pipeline for both NMR chemical shift assignment and structure determination of biological macromolecules. In addition, it encompasses tools to analyse the many additional experiments that make NMR such a pivotal technique for research into complex biological questions. This report describes how CcpNmr Analysis can seamlessly link together all of the tasks in the NMR structure-determination process. It details each of the stages from generating NMR restraints [distance, dihedral, hydrogen bonds and residual dipolar couplings (RDCs)], exporting these to and subsequently re-importing them from structure-calculation software (such as the programs CYANA or ARIA) and analysing and validating the results obtained from the structure calculation to, ultimately, the streamlined deposition of the completed assignments and the refined ensemble of structures into the PDBe repository. Until recently, such solution-structure determination by NMR has been quite a laborious task, requiring multiple stages and programs. However, with the new enhancements to CcpNmr Analysis described here, this process is now much more intuitive and efficient and less error-prone.

Correlation of chemical shifts predicted by molecular dynamics simulations for partially disordered proteins

There has been a longstanding interest in being able to accurately predict NMR chemical shifts from structural data. Recent studies have focused on using molecular dynamics (MD) simulation data as input for improved prediction. Here we examine the accuracy of chemical shift prediction for intein systems, which have regions of intrinsic disorder. We find that using MD simulation data as input for chemical shift prediction does not consistently improve prediction accuracy over use of a static X-ray crystal structure. This appears to result from the complex conformational ensemble of the disordered protein segments. We show that using accelerated molecular dynamics (aMD) simulations improves chemical shift prediction, suggesting that methods which better sample the conformational ensemble like aMD are more appropriate tools for use in chemical shift prediction for proteins with disordered regions. Moreover, our study suggests that data accurately reflecting protein dynamics must be used as input for chemical shift prediction in order to correctly predict chemical shifts in systems with disorder.

Case study of hydrogen bonding in a hydrophobic cavity

Protein internal hydrogen bonds and hydrophobicity determine protein folding and structure stabilization, and the introduction of a hydrogen bond has been believed to represent a better interaction for consolidating protein structure. We observed an alternative example for chicken IL-1β. The native IL-1β contains a hydrogen bond between the Y157 side-chain OηH and I133 backbone CO, whereby the substitution from Tyr to Phe abolishes the connection and the mutant without the hydrogen bond is more stable. An attempt to explain the energetic view of the presence of the hydrogen bond fails when only considering the nearly identical X-ray structures. Here, we resolve the mechanism by monitoring the protein backbone dynamics and interior hydrogen bond network. IL-1β contains a hydrophobic cavity in the protein interior, and Y157 is one of the surrounding residues. The Y157 OηH group introduces an unfavorable energy in the hydrophobic cavity, therefore sequestering itself by forming a hydrogen bond with the proximate residue I133. The hydrogen bonding confines Y157 orientation but exerts a force to disrupt the hydrogen bond network surrounding the cavity. The effect propagates over the entire protein and reduces the stability, as reflected in the protein backbone dynamics observed by an NMR hydrogen-deuterium (H/D) exchange experiment. We describe the particular case in which a hydrogen bond does not necessarily confer enhanced protein stability while the disruption of hydrophobicity must be integrally considered.

Triple resonance-based ¹³C(α) and ¹³C(β) CEST experiments for studies of ms timescale dynamics in proteins

A pair of triple resonance based CEST pulse schemes are presented for measuring ¹³C(α) and ¹³C(β) chemical shifts of sparsely populated and transiently formed conformers that are invisible to traditional NMR experiments. CEST profiles containing dips at resonance positions of ¹³C(α) or ¹³C(β) spins of major (ground) and minor (excited) conformers are obtained in a pseudo 3rd dimension that is generated by quantifying modulations of cross peaks in ¹⁵N, ¹H(N) correlation spectra. An application to the folding reaction of a G48A mutant of the Fyn SH3 domain is presented, illustrating and validating the methodology.

Molecular insights into the interaction between Plasmodium falciparum apical membrane antigen 1 and an invasion-inhibitory peptide

Apical membrane antigen 1 (AMA1) of the human malaria parasite Plasmodium falciparum has been implicated in invasion of the host erythrocyte. It interacts with malarial rhoptry neck (RON) proteins in the moving junction that forms between the host cell and the invading parasite. Agents that block this interaction inhibit invasion and may serve as promising leads for anti-malarial drug development. The invasion-inhibitory peptide R1 binds to a hydrophobic cleft on AMA1, which is an attractive target site for small molecules that block parasite invasion. In this work, truncation and mutational analyses show that Phe5-Phe9, Phe12 and Arg15 in R1 are the most important residues for high affinity binding to AMA1. These residues interact with two well-defined binding hot spots on AMA1. Computational solvent mapping reveals that one of these hot spots is suitable for small molecule targeting. We also confirm that R1 in solution binds to AMA1 with 1∶1 stoichiometry and adopts a secondary structure consistent with the major form of R1 observed in the crystal structure of the complex. Our results provide a basis for designing high affinity inhibitors of the AMA1-RON2 interaction.

Backbone ¹H, ¹⁵N, and ¹³C and side chain ¹³Cβ NMR chemical shift assignment of murine interleukin-10

Almost complete assignment of backbone (1)H, (13)C, (15)N and side chain (13)Cβ resonances for the immune-regulatory cytokine IL-10 is reported. The protein was overexpressed in Escherichia coli and was refolded from inclusion bodies. The point mutation C149Y was introduced to suppress incorrect disulfide bond formation and to improve protein refolding.

¹H, ¹³C and ¹⁵N NMR assignments of a mutant of UV inducible transcript (S55A-UVI31+) from Chlamydomonas reinhardtii

Almost complete sequence specific (1)H, (13)C and (15)N resonance assignments of a mutant of UV inducible transcript (S55A-UVI31+) from Chlamydomonas reinhardtii are reported, as a prelude to its structural and functional characterization. Site-directed mutagenesis of uvi31+ was carried out using complementary mutation harbouring oligonucleotides for S55A mutant. The resulting mutant was sequenced and then S55A mutant of uvi31+ cDNA were expressed in E. coli, and purification were carried out from where the protein was purified to high homogeneity. The point mutation S55A in UVI31+ results in a significant enhancement in its DNA endonuclease activity as compared to its wild type protein.

Discriminating protein structure classes by incorporating Pseudo Average Chemical Shift to Chou's general PseAAC and Support Vector Machine

Proteins control all biological functions in living species. Protein structure is comprised of four major classes including all-α class, all-β class, α+β, and α/β. Each class performs different function according to their nature. Owing to the large exploration of protein sequences in the databanks, the identification of protein structure classes is difficult through conventional methods with respect to cost and time. Looking at the importance of protein structure classes, it is thus highly desirable to develop a computational model for discriminating protein structure classes with high accuracy. For this purpose, we propose a silco method by incorporating Pseudo Average Chemical Shift and Support Vector Machine. Two feature extraction schemes namely Pseudo Amino Acid Composition and Pseudo Average Chemical Shift are used to explore valuable information from protein sequences. The performance of the proposed model is assessed using four benchmark datasets 25PDB, 1189, 640 and 399 employing jackknife test. The success rates of the proposed model are 84.2%, 85.0%, 86.4%, and 89.2%, respectively on the four datasets. The empirical results reveal that the performance of our proposed model compared to existing models is promising in the literature so far and might be useful for future research.

¹H, ¹⁵N, and ¹³C chemical shift assignments for the winged helix domains of two archeal MCM C-termini

High-fidelity replication guarantees the stable inheritance of genetic information stored in the DNA of living organisms. The minichromosome maintenance (MCM) complex functions as replicative DNA-unwinding helicase and has been identified as one key player in the replication process of archea and eukarya. Despite the availability of considerable structural information on archeal MCMs, such information was missing for their C-terminal domain. In order to obtain more detailed structural information, we assigned the NMR chemical shifts for backbone and side chain nuclei for the MCM C-terminal winged helix domains of the archeal species Methanothermobacter thermautrophicus and Sulfolobus solfataricus.

Inhibition of the human respiratory syncytial virus small hydrophobic protein and structural variations in a bicelle environment

The small hydrophobic (SH) protein is a 64-amino-acid polypeptide encoded by the human respiratory syncytial virus (hRSV). SH protein has a single α-helical transmembrane (TM) domain that forms pentameric ion channels. Herein, we report the first inhibitor of the SH protein channel, pyronin B, and we have mapped its binding site to a conserved surface of the RSV SH pentamer, at the C-terminal end of the transmembrane domain. The validity of the SH protein structural model used has been confirmed by using a bicellar membrane-mimicking environment. However, in bicelles the α-helical stretch of the TM domain extends up to His-51, and by comparison with previous models both His-22 and His-51 adopt an interhelical/lumenal orientation relative to the channel pore. Neither His residue was found to be essential for channel activity although His-51 protonation reduced channel activity at low pH, with His-22 adopting a more structural role. The latter results are in contrast with previous patch clamp data showing channel activation at low pH, which could not be reproduced in the present work. Overall, these results establish a solid ground for future drug development targeting this important viroporin. Importance: The human respiratory syncytial virus (hRSV) is responsible for 64 million reported cases of infection and 160,000 deaths each year. Lack of adequate antivirals fuels the search for new targets for treatment. The small hydrophobic (SH) protein is a 64-amino-acid polypeptide encoded by hRSV and other paramyxoviruses, and its absence leads to viral attenuation in vivo and early apoptosis in infected cells. SH protein forms pentameric ion channels that may constitute novel drug targets, but no inhibitor for this channel activity has been reported so far. A small-molecule inhibitor, pyronin B, can reduce SH channel activity, and its likely binding site on the SH protein channel has been identified. Black lipid membrane (BLM) experiments confirm that protonation of both histidine residues reduces stability and channel activity. These results contrast with previous patch clamp data that showed low-pH activation, which we have not been able to reproduce.

acACS: improving the prediction accuracy of protein subcellular locations and protein classification by incorporating the average chemical shifts composition

The chemical shift is sensitive to changes in the local environments and can report the structural changes. The structure information of a protein can be represented by the average chemical shifts (ACS) composition, which has been broadly applied for enhancing the prediction accuracy in protein subcellular locations and protein classification. However, different kinds of ACS composition can solve different problems. We established an online web server named acACS, which can convert secondary structure into average chemical shift and then compose the vector for representing a protein by using the algorithm of auto covariance. Our solution is easy to use and can meet the needs of users.

Probing the non-native H helix translocation in apomyoglobin folding intermediates

Apomyoglobin folds via sequential helical intermediates that are formed by rapid collapse of the A, B, G, and H helix regions. An equilibrium molten globule with a similar structure is formed near pH 4. Previous studies suggested that the folding intermediates are kinetically trapped states in which folding is impeded by non-native packing of the G and H helices. Fluorescence spectra of mutant proteins in which cysteine residues were introduced at several positions in the G and H helices show differential quenching of W14 fluorescence, providing direct evidence of translocation of the H helix relative to helices A and G in both the kinetic and equilibrium intermediates. Förster resonance energy transfer measurements show that a 5-({2-[(acetyl)amino]ethyl}amino)naphthalene-1-sulfonic acid acceptor coupled to K140C (helix H) is closer to Trp14 (helix A) in the equilibrium molten globule than in the native state, by a distance that is consistent with sliding of the H helix in an N-terminal direction by approximately one helical turn. Formation of an S108C–L135C disulfide prevents H helix translocation in the equilibrium molten globule by locking the G and H helices into their native register. By enforcing nativelike packing of the A, G, and H helices, the disulfide resolves local energetic frustration and facilitates transient docking of the E helix region onto the hydrophobic core but has only a small effect on the refolding rate. The apomyoglobin folding landscape is highly rugged, with several energetic bottlenecks that frustrate folding; relief of any one of the major identified bottlenecks is insufficient to speed progression to the transition state.

NMR spectroscopy of native and in vitro tissues implicates polyADP ribose in biomineralization

Nuclear magnetic resonance (NMR) spectroscopy is useful to determine molecular structure in tissues grown in vitro only if their fidelity, relative to native tissue, can be established. Here, we use multidimensional NMR spectra of animal and in vitro model tissues as fingerprints of their respective molecular structures, allowing us to compare the intact tissues at atomic length scales. To obtain spectra from animal tissues, we developed a heavy mouse enriched by about 20% in the NMR-active isotopes carbon-13 and nitrogen-15. The resulting spectra allowed us to refine an in vitro model of developing bone and to probe its detailed structure. The identification of an unexpected molecule, poly(adenosine diphosphate ribose), that may be implicated in calcification of the bone matrix, illustrates the analytical power of this approach.

The arginine-rich RNA-binding motif of HIV-1 Rev is intrinsically disordered and folds upon RRE binding

Arginine-rich motifs (ARMs) capable of binding diverse RNA structures play critical roles in transcription, translation, RNA trafficking, and RNA packaging. The regulatory HIV-1 protein Rev is essential for viral replication and belongs to the ARM family of RNA-binding proteins. During the early stages of the HIV-1 life cycle, incompletely spliced and full-length viral mRNAs are very inefficiently recognized by the splicing machinery of the host cell and are subject to degradation in the cell nucleus. These transcripts harbor the Rev Response Element (RRE), which orchestrates the interaction with the Rev ARM and the successive Rev-dependent mRNA export pathway. Based on established criteria for predicting intrinsic disorder, such as hydropathy, combined with significant net charge, the very basic primary sequences of ARMs are expected to adopt coil-like structures. Thus, we initiated this study to investigate the conformational changes of the Rev ARM associated with RNA binding. We used multidimensional NMR and circular dichroism spectroscopy to monitor the observed structural transitions, and described the conformational landscapes using statistical ensemble and molecular-dynamics simulations. The combined spectroscopic and simulated results imply that the Rev ARM is intrinsically disordered not only as an isolated peptide but also when it is embedded into an oligomerization-deficient Rev mutant. RRE recognition triggers a crucial coil-to-helix transition employing an induced-fit mechanism.

¹H, ¹³C and ¹⁵N resonance assignments of S114A mutant of UVI31+ from Chlamydomonas reinhardtii

Almost complete sequence specific (1)H, (13)C and (15)N resonance assignments of S114A mutant of UVI31+ from Chlamydomonas reinhardtii are reported. The cDNA of S114A mutant of UVI31+ was cloned from a eukaryotic green algae (C. reinhardtii) and overexpressed in E.coli from where the protein was purified to homogeneity. The point mutation S114A in UVI31+ reduces its DNA endonuclease activity substantially as compared with its wild type. As a prelude to the structural characterization of S114A mutant of UVI31+, we report here complete sequence-specific (1)H, (13)C and (15)N NMR assignments.

Structure of a conserved Golgi complex-targeting signal in coronavirus envelope proteins

Coronavirus envelope (CoV E) proteins are ∼100-residue polypeptides with at least one channel-forming α-helical transmembrane (TM) domain. The extramembrane C-terminal tail contains a completely conserved proline, at the center of a predicted β-coil-β motif. This hydrophobic motif has been reported to constitute a Golgi-targeting signal or a second TM domain. However, no structural data for this or other extramembrane domains in CoV E proteins is available. Herein, we show that the E protein in the severe acute respiratory syndrome virus has only one TM domain in micelles, whereas the predicted β-coil-β motif forms a short membrane-bound α-helix connected by a disordered loop to the TM domain. However, complementary results suggest that this motif is potentially poised for conformational change or in dynamic exchange with other conformations.

NMR characterization of the near native and unfolded states of the PTB domain of Dok1: alternate conformations and residual clusters

BACKGROUND

Phosphotyrosine binding (PTB) domains are critically involved in cellular signaling and diseases. PTB domains are categorized into three distinct structural classes namely IRS-like, Shc-like and Dab-like. All PTB domains consist of a core pleckstrin homology (PH) domain with additional structural elements in Shc and Dab groups. The core PH fold of the PTB domain contains a seven stranded β-sheet and a long C-terminal helix.

PRINCIPAL FINDINGS

In this work, the PTB domain of Dok1 protein has been characterized, by use of NMR spectroscopy, in solutions containing sub-denaturing and denaturing concentrations of urea. We find that the Dok1 PTB domain displays, at sub-denaturing concentrations of urea, alternate conformational states for residues located in the C-terminal helix and in the β5 strand of the β-sheet region. The β5 strand of PTB domain has been found to be experiencing significant chemical shift perturbations in the presence of urea. Notably, many of these residues in the helix and in the β5 strand are also involved in ligand binding. Structural and dynamical analyses at 7 M urea showed that the PTB domain is unfolded with islands of motionally restricted regions in the polypeptide chain. Further, the C-terminal helix appears to be persisted in the unfolded state of the PTB domain. By contrast, residues encompassing β-sheets, loops, and the short N-terminal helix lack any preferred secondary structures. Moreover, these residues demonstrated an intimate contact with the denaturant.

SIGNIFICANCE

This study implicates existence of alternate conformational states around the ligand binding pocket of the PTB domain either in the native or in the near native conditions. Further, the current study demonstrates that the C-terminal helical region of PTB domain may be considered as a potential site for the initiation of folding.

Capping parallel β-sheets of acetyl(Ala)6NH2 with an acetyl(Ala)5ProNH2 can arrest the growth of the sheet, suggesting a potential for curtailing amyloid growth. An ONIOM and density functional theory study

We present ONIOM calculations using B3LYP/d95(d,p) as the high level and AM1 as the medium level on parallel β-sheets containing four strands of Ac-AAAAAA-NH2 capped with either Ac-AAPAAA-NH2 or Ac-AAAPAA-NH2. Because Pro can form H-bonds from only one side of the peptide linkage (that containing the C═O H-bond acceptor), only one of the two Pro-containing strands can favorably add to the sheet on each side. Surprisingly, when the sheet is capped with AAPAAA-NH2 at one edge, the interaction between the cap and sheet is slightly more stabilizing than that of another all Ala strand. Breaking down the interaction enthalpies into H-bonding and distortion energies shows the favorable interaction to be due to lower distortion energies in both the strand and the four-stranded sheet. Because another strand would be inhibited for attachment to the other side of the capping (Pro-containing) strand, we suggest the possible use of Pro residues in peptides designed to arrest the growth of many amyloids.

"Structural characterization of the minimal segment of TDP-43 competent for aggregation"

TDP-43 is a nuclear protein whose abnormal aggregates are implicated in ALS and FTLD. Recently, an Asn/Gln rich C-terminal segment of TDP-43 has been shown to produce aggregation in vitro and reproduce most of the protein's pathological hallmarks in cells, but little is known about this segment's structure. Here, CD and 2D heteronuclear NMR spectroscopies provide evidence that peptides corresponding to the wild type and mutated sequences of this segment adopt chiefly disordered conformations that, in the case of the wild type sequence, spontaneously forms a β-sheet rich oligomer. Moreover, MD simulation provides evidence for a structure consisting of two β-strands and a well-defined, yet non-canonical structural element. Furthermore, MD simulations of four pathological mutations (Q343R, N345K, G348V and N352S) occurring in this segment predict that all of them could affect this region's structure. In particular, the Q343R variant tends to stabilize disordered conformers, N345K permits the formation of longer, more stable β-strands, and G348V tends to shorten and destabilize them. Finally, N352S acts to alter the β-stand register and when S352 is phosphorylated, it induces partial unfolding. Our results provide a better understanding of TDP-43 aggregation process and will be useful to design effectors capable to modulate its progression.

13C Chemical Shifts in Proteins: A Rich Source of Encoded Structural Information

Despite the formidable progress in Nuclear Magnetic Resonance (NMR) spectroscopy, quality assessment of NMR-derived structures remains as an important problem. Thus, validation of protein structures is essential for the spectroscopists, since it could enable them to detect structural flaws and potentially guide their efforts in further refinement. Moreover, availability of accurate and efficient validation tools would help molecular biologists and computational chemists to evaluate quality of available experimental structures and to select a protein model which is the most suitable for a given scientific problem. The ¹³C^α nuclei are ubiquitous in proteins, moreover, their shieldings are easily obtainable from NMR experiments and represent a rich source of encoded structural information that makes ¹³C^α chemical shifts an attractive candidate for use in computational methods aimed at determination and validation of protein structures. In this chapter, the basis of a novel methodology of computing, at the quantum chemical level of theory, the ¹³C^α shielding for the amino acid residues in proteins is described. We also identify and examine the main factors affecting the ¹³C^α-shielding computation. Finally, we illustrate how the information encoded in the ¹³C chemical shifts can be used for a number of applications, viz., from protein structure prediction of both α-helical and β-sheet conformations, to determination of the fraction of the tautomeric forms of the imidazole ring of histidine in proteins as a function of pH or to accurate detection of structural flaws, at a residue-level, in NMR-determined protein models.

Local interactions influence the fibrillation kinetics, structure and dynamics of Aβ(1–40) but leave the general fibril structure unchanged

Selective point mutations introducing local fields do not alter the overall structure and morphology of Aβ(1–40) fibrils.

Conformational dynamics is more important than helical propensity for the folding of the all α-helical protein Im7

Im7 folds via an on-pathway intermediate that contains three of the four native α-helices. The missing helix, helix III, is the shortest and its failure to be formed until late in the pathway is related to frustration in the structure. Im7H3M3, a 94-residue variant of the 87-residue Im7 in which helix III is the longest of the four native helices, also folds via an intermediate. To investigate the structural basis for this we calculated the frustration in the structure of Im7H3M3 and used NMR to investigate its dynamics. We found that the native state of Im7H3M3 is highly frustrated and in equilibrium with an intermediate state that lacks helix III, similar to Im7. Model-free analysis identified residues with chemical exchange contributions to their relaxation that aligned with the residues predicted to have highly frustrated interactions, also like Im7. Finally, we determined properties of urea-denatured Im7H3M3 and identified four clusters of interacting residues that corresponded to the α-helices of the native protein. In Im7 the cluster sizes were related to the lengths of the α-helices with cluster III being the smallest but in Im7H3M3 cluster III was also the smallest, despite this region forming the longest helix in the native state. These results suggest that the conformational properties of the urea-denatured states promote formation of a three-helix intermediate in which the residues that form helix III remain non-helical. Thus it appears that features of the native structure are formed early in folding linked to collapse of the unfolded state.

Global fold of human cannabinoid type 2 receptor probed by solid-state 13C-, 15N-MAS NMR and molecular dynamics simulations

The global fold of human cannabinoid type 2 (CB2 ) receptor in the agonist-bound active state in lipid bilayers was investigated by solid-state (13)C- and (15)N magic-angle spinning (MAS) NMR, in combination with chemical-shift prediction from a structural model of the receptor obtained by microsecond-long molecular dynamics (MD) simulations. Uniformly (13)C- and (15)N-labeled CB2 receptor was expressed in milligram quantities by bacterial fermentation, purified, and functionally reconstituted into liposomes. (13)C MAS NMR spectra were recorded without sensitivity enhancement for direct comparison of Cα, Cβ, and C=O bands of superimposed resonances with predictions from protein structures generated by MD. The experimental NMR spectra matched the calculated spectra reasonably well indicating agreement of the global fold of the protein between experiment and simulations. In particular, the (13) C chemical shift distribution of Cα resonances was shown to be very sensitive to both the primary amino acid sequence and the secondary structure of CB2. Thus the shape of the Cα band can be used as an indicator of CB2 global fold. The prediction from MD simulations indicated that upon receptor activation a rather limited number of amino acid residues, mainly located in the extracellular Loop 2 and the second half of intracellular Loop 3, change their chemical shifts significantly (≥ 1.5 ppm for carbons and ≥ 5.0 ppm for nitrogens). Simulated two-dimensional (13) Cα(i)-(13)C=O(i) and (13)C=O(i)-(15)NH(i + 1) dipolar-interaction correlation spectra provide guidance for selective amino acid labeling and signal assignment schemes to study the molecular mechanism of activation of CB2 by solid-state MAS NMR.

Physics-based method to validate and repair flaws in protein structures

A method that makes use of information provided by the combination of (13)C(α) and (13)C(β) chemical shifts, computed at the density functional level of theory, enables one to (i) validate, at the residue level, conformations of proteins and detect backbone or side-chain flaws by taking into account an ensemble average of chemical shifts over all of the conformations used to represent a protein, with a sensitivity of ∼90%; and (ii) provide a set of (χ1/χ2) torsional angles that leads to optimal agreement between the observed and computed (13)C(α) and (13)C(β) chemical shifts. The method has been incorporated into the CheShift-2 protein validation Web server. To test the reliability of the provided set of (χ1/χ2) torsional angles, the side chains of all reported conformations of five NMR-determined protein models were refined by a simple routine, without using NOE-based distance restraints. The refinement of each of these five proteins leads to optimal agreement between the observed and computed (13)C(α) and (13)C(β) chemical shifts for ∼94% of the flaws, on average, without introducing a significantly large number of violations of the NOE-based distance restraints for a distance range ≤ 0.5 , in which the largest number of distance violations occurs. The results of this work suggest that use of the provided set of (χ1/χ2) torsional angles together with other observables, such as NOEs, should lead to a fast and accurate refinement of the side-chain conformations of protein models.

Probing the residual structure in avian prion hexarepeats by CD, NMR and MD techniques

Many proteins perform essential biological functions by means of regions that lacking specific organized structure exist as an ensemble of interconverting transient conformers. The characterization of such regions, including the description of their structural propensities, number of conformations and relative populations can provide useful insights. Prion diseases result from the conversion of a normal glycoprotein into a misfolded pathogenic isoform. The structures of mammal and chicken prion proteins show a similar fold with a globular domain and a flexible N-terminal portion that contains different repeated regions: octarepeats (PHGGGWGQ) in mammals and hexarepeats (PHNPGY) in chickens. The higher number of prolines in the hexarepeat region suggests that this region may retain a significant amount of residual secondary structure. Here, we report the CD, NMR and MD characterization of a peptide (2-HexaPY) composed of two hexarepeats. We combine experimental NMR data and MD to investigate at atomic level its ensemble-averaged structural properties, demonstrating how each residue of both repeats has a different quantified PPII propensity that shows a periodicity along the sequence. This feature explains the absence of cooperativity to stabilize a PPII conformation. Nonetheless, such residual structure can play a role in nucleating local structural transitions as well as modulating intra-molecular or inter-molecular interactions.

13C structuring shifts for the analysis of model β-hairpins and β-sheets in proteins: diagnostic shifts appear only at the cross-strand H-bonded residues

The present studies have shown that (13)C=O, (13)C(α) and (13)C(β) of H-bonded strand residues in β-hairpins provide additional probes for quantitating the extent of folding in β-hairpins and other β-sheet models. Large differences in the structuring shifts (CSDs) of these (13)C sites in H-bonded versus non-H-bonded sites are observed: the differences between H-bonded and non-H-bonded sites are greater than 1.2 ppm for all three (13)C probes. This prompts us to suggest that efforts to determine the extent of hairpin folding from (13)C shifts should be based exclusively on the observation at the cross-strand H-bonded sites. Furthermore, the statistics suggest the (13)C' and (13)C(β) CSDs will provide the best differentiation with 100%-folded CSD values approaching -2.6 and +3 ppm, respectively, for the H-bonded sites. These conclusions can be extended to edge-strands of protein β-sheets. Our survey of reported (13)C shifts in β-proteins indicates that some of the currently employed random coil values need to be adjusted, particularly for ionization-induced effects.

Protein backbone and sidechain torsion angles predicted from NMR chemical shifts using artificial neural networks

A new program, TALOS-N, is introduced for predicting protein backbone torsion angles from NMR chemical shifts. The program relies far more extensively on the use of trained artificial neural networks than its predecessor, TALOS+. Validation on an independent set of proteins indicates that backbone torsion angles can be predicted for a larger, ≥90 % fraction of the residues, with an error rate smaller than ca 3.5 %, using an acceptance criterion that is nearly two-fold tighter than that used previously, and a root mean square difference between predicted and crystallographically observed (ϕ, ψ) torsion angles of ca 12º. TALOS-N also reports sidechain χ(1) rotameric states for about 50 % of the residues, and a consistency with reference structures of 89 %. The program includes a neural network trained to identify secondary structure from residue sequence and chemical shifts.

Structural transitions in tau k18 on micelle binding suggest a hierarchy in the efficacy of individual microtubule-binding repeats in filament nucleation

The protein tau is found in an aggregated filamentous state in the intraneuronal paired helical filament deposits characteristic of Alzheimer's disease and other related dementias and mutations in tau protein and mRNA cause frontotemproal dementia. Tau isoforms include a microtubule-binding domain containing either three or four imperfect tandem microtubule binding repeats that also form the core of tau filaments and contain hexapaptide motifs that are critical for tau aggregation. The tau microtubule-binding domain can also engage in direct interactions with detergents, fatty acids, or membranes, which can greatly facilitate tau aggregation and may also mediate some tau functions. Here, we show that the alternatively spliced second microtubule-binding repeat exhibits significantly different structural characteristics compared with the other three repeats in the context of the intact repeat domain. Most notably, the PHF6* hexapeptide motif located at the N-terminus of repeat 2 has a lower propensity to form strand-like structure than the corresponding PHF6 motif in repeat 3, and unlike PHF6 converts to partially helical structure in the micelle-bound state. Interestingly, the behavior of the Module-B motif, located at the beginning of repeat 4, resembles that of PHF6* rather than PHF6. Our observations, combined with previous results showing that PHF6* and Module-B are both less effective than PHF6 in nucleating tau aggregation, suggest a hierarchy in the efficacy of these motifs in nucleating tau aggregation that originates in differences in their intrinsic propensities for extended strand-like structure and the resistance of these propensities to changes in tau's environment.

Functional manipulation of a calcium-binding protein from Entamoeba histolytica guided by paramagnetic NMR

EhCaBP1, one of the calcium-binding proteins from Entamoeba histolytica, is a two-domain EF-hand protein. The two domains of EhCaBP1 are structurally and functionally different from each other. However, both domains are required for structural stability and a full range of functional diversity. Analysis of sequence and structure of EhCaBP1 and other CaBPs indicates that the C-terminal domain of EhCaBP1 possesses a unique structure compared with other family members. This had been attributed to the absence of a Phe-Phe interaction between highly conserved Phe residues at the -4 position in EF-hand III (F[-4]; Tyr(81)) and at the 13th position in EF-hand IV (F[+13]; Phe(129)) of the C-terminal domain. Against this backdrop, we mutated the Tyr residue at the -4th position of EF III to the Phe residue (Y81F), to bring in the Phe-Phe interaction and understand the nature of structural and functional changes in the protein by NMR spectroscopy, molecular dynamics (MD) simulation, isothermal titration calorimetry (ITC), and biological assays, such as imaging and actin binding. The Y81F mutation in EhCaBP1 resulted in a more compact structure for the C-terminal domain of the mutant as in the case of calmodulin and troponin C. The compact structure is favored by the presence of a π-π interaction between Phe(81) and Phe(129) along with several hydrophobic interactions of Phe(81), which are not seen in the wild-type protein. Furthermore, the biological assays reveal preferential membrane localization of the mutant, loss of its colocalization with actin in the phagocytic cups, whereas retaining its ability to bind G- and F-actin.

Structural characterization by NMR of a double phosphorylated chimeric peptide vaccine for treatment of Alzheimer's disease

Rational design of peptide vaccines becomes important for the treatment of some diseases such as Alzheimer's disease (AD) and related disorders. In this study, as part of a larger effort to explore correlations of structure and activity, we attempt to characterize the doubly phosphorylated chimeric peptide vaccine targeting a hyperphosphorylated epitope of the Tau protein. The 28-mer linear chimeric peptide consists of the double phosphorylated B cell epitope Tau₂₂₉₋₂₃₇[pThr231/pSer235] and the immunomodulatory T cell epitope Ag85B₂₄₁₋₂₅₅ originating from the well-known antigen Ag85B of the Mycobacterium tuberculosis, linked by a four amino acid sequence -GPSL-. NMR chemical shift analysis of our construct demonstrated that the synthesized peptide is essentially unfolded with a tendency to form a β-turn due to the linker. In conclusion, the -GPSL- unit presumably connects the two parts of the vaccine without transferring any structural information from one part to the other. Therefore, the double phosphorylated epitope of the Tau peptide is flexible and accessible.

Solution nuclear magnetic resonance spectroscopy

Solution nuclear magnetic resonance (NMR) spectroscopy has come a long way in characterizing the structure and function of biological molecules since the first one-dimensional spectrum of protein was recorded about 30 years ago. To date (September 1, 2012), there are 9,521 solution NMR structures in the Protein Data Bank, compared to 74,009 determined by crystallographic methods. Unlike X-ray and electron microscopy (EM) methods, which are based on the concepts of Fourier optics and image reconstruction, structure determination by NMR involves measuring structural restraints and finding structural solutions that satisfy the restraints. Although the NMR approach is much less direct in a physical sense, it has proven itself over the years to be capable of de novo structure determination at high precision. Moreover, the method is highly versatile and can be used in a variety of ways for addressing mechanistic questions. NMR measurements of protein internal dynamics and protein-protein or protein-ligand interaction are directly relevant to function in vivo because the molecules are often in physiological buffer conditions. The method can also be applied to investigate protein-folding intermediates, conformational changes, as well as intrinsically unfolded proteins. Recently, along with X-ray and EM, solution NMR has entered a state of rapid growth for structural studies of membrane proteins, already demonstrating its feasibility in de novo structure determination of membrane-embedded ion channels and receptors. As the hardware advances rapidly, especially in cryogenic probes that have much higher sensitivity, the sample concentration required for solution NMR investigation is decreasing, hopefully soon to a concentration level at which nonspecific protein aggregation is no longer an issue. After three decades of improvement in spectrometer technology, NMR pulse experiments, isotope labeling schemes, and structure determination software, we believe that solution NMR will truly enter the production phase in the next decade to answer biological questions of high impact, and to become more versatile than ever in complementing X-ray and EM in investigating protein structure and function.

Conformational coupling across the plasma membrane in activation of the EGF receptor

How the epidermal growth factor receptor (EGFR) activates is incompletely understood. The intracellular portion of the receptor is intrinsically active in solution, and to study its regulation, we measured autophosphorylation as a function of EGFR surface density in cells. Without EGF, intact EGFR escapes inhibition only at high surface densities. Although the transmembrane helix and the intracellular module together suffice for constitutive activity even at low densities, the intracellular module is inactivated when tethered on its own to the plasma membrane, and fluorescence cross-correlation shows that it fails to dimerize. NMR and functional data indicate that activation requires an N-terminal interaction between the transmembrane helices, which promotes an antiparallel interaction between juxtamembrane segments and release of inhibition by the membrane. We conclude that EGF binding removes steric constraints in the extracellular module, promoting activation through N-terminal association of the transmembrane helices.

Probing slowly exchanging protein systems via ¹³Cα-CEST: monitoring folding of the Im7 protein

A¹³C(α) chemical exchange saturation transfer based experiment is presented for the study of protein systems undergoing slow interconversion between an 'observable' ground state and one or more 'invisible' excited states. Here a labeling strategy whereby [2-(13)C]-glucose is the sole carbon source is exploited, producing proteins with ¹³C at the C(α) position, while the majority of residues remain unlabeled at CO or C(β). The new experiment is demonstrated with an application to the folding reaction of the Im7 protein that involves an on-pathway excited state. The obtained excited state (13)C(α) chemical shifts are cross validated by comparison to values extracted from analysis of CPMG relaxation dispersion profiles, establishing the utility of the methodology.

(13)C-(13)c homonuclear recoupling in solid-state nuclear magnetic resonance at a moderately high magic-angle-spinning frequency

Two-dimensional (13)C-(13)C correlation experiments are widely employed in structure determination of protein assemblies using solid-state nuclear magnetic resonance. Here, we investigate the process of (13)C-(13)C magnetisation transfer at a moderate magic-angle-spinning frequency of 30 kHz using some of the prominent second-order dipolar recoupling schemes. The effect of isotropic chemical-shift difference and spatial distance between two carbons and amplitude of radio frequency on (1)H channel on the magnetisation transfer efficiency of these schemes is discussed in detail.

High-pressure NMR reveals close similarity between cold and alcohol protein denaturation in ubiquitin

Proteins denature not only at high, but also at low temperature as well as high pressure. These denatured states are not easily accessible for experiment, because usually heat denaturation causes aggregation, whereas cold or pressure denaturation occurs at temperatures well below the freezing point of water or pressures above 5 kbar, respectively. Here we have obtained atomic details of the pressure-assisted, cold-denatured state of ubiquitin at 2,500 bar and 258 K by high-resolution NMR techniques. Under these conditions, a folded, native-like and a disordered state exist in slow exchange. Secondary chemical shifts show that the disordered state has structural propensities for a native-like N-terminal β-hairpin and α-helix and a nonnative C-terminal α-helix. These propensities are very similar to the previously described alcohol-denatured (A-)state. Similar to the A-state, (15)N relaxation data indicate that the secondary structure elements move as independent segments. The close similarity of pressure-assisted, cold-denatured, and alcohol-denatured states with native and nonnative secondary elements supports a hierarchical mechanism of folding and supports the notion that similar to alcohol, pressure and cold reduce the hydrophobic effect. Indeed, at nondenaturing concentrations of methanol, a complete transition from the native to the A-state can be achieved at ambient temperature by varying the pressure from 1 to 2,500 bar. The methanol-assisted pressure transition is completely reversible and can also be induced in protein G. This method should allow highly detailed studies of protein-folding transitions in a continuous and reversible manner.

Structural ensemble of an intrinsically disordered polypeptide

Intrinsically disordered proteins (IDPs), which play key roles in cell signaling and regulation, do not display specific tertiary structure when isolated in solution. Instead, they dynamically explore an ensemble of unfolded configurations, adopting more stable, ordered structures only after binding to their ligands. Whether ligands induce IDP structural changes upon binding or simply bind to pre-existing conformers that are populated within the IDP's structural ensemble is not well understood. Molecular simulations can provide information with the spatiotemporal resolution necessary to resolve these issues. Here, we report on the conformational ensemble of a 15-residue wild-type p53 fragment from the TAD domain and its mutant (TAD-P27L) obtained by replica exchange molecular dynamics simulation using an optimized (fully atomistic, explicit solvent) protein model and the experimental validation of the simulation results. We use a clustering method based on structural similarity to identify conformer states populated by the peptides in solution from the simulated ensemble. We show that p53 populates solution structures that strongly resemble the ligand (MDM2)-bound structure, but at the same time, the conformational free-energy landscape is relatively flat in the absence of the ligand.

Comparison of some dispersion-corrected and traditional functionals as applied to peptides and conformations of cyclohexane derivatives

We compare the energetic and structural properties of fully optimized α-helical and antiparallel β-sheet polyalanines and the energetic differences between axial and equatorial conformations of three cyclohexane derivatives (methyl, fluoro, and chloro) as calculated using several functionals designed to treat dispersion (B97-D, ωB97x-D, M06, M06L, and M06-2X) with other traditional functionals not specifically parametrized to treat dispersion (B3LYP, X3LYP, and PBE1PBE) and with experimental results. Those functionals developed to treat dispersion significantly overestimate interaction enthalpies of folding for the α-helix and predict unreasonable structures that contain Ramachandran φ and ψ and C = O...N H-bonding angles that are out of the bounds of databases compiled the β-sheets. These structures are consistent with overestimation of the interaction energies. For the cyclohexanes, these functionals overestimate the stabilities of the axial conformation, especially when used with smaller basis sets. Their performance improves when the basis set is improved from D95∗∗ to aug-cc-pVTZ (which would not be possible with systems as large as the peptides).

The maturational refolding of the β-hairpin motif of equine infectious anemia virus capsid protein extends its helix α1 at capsid assembly locus

A retroviral capsid (CA) protein consists of two helical domains, CA^N and CA^C, which drive hexamer and dimer formations, respectively, to form a capsid lattice. The N-terminal 13 residues of CA refold to a β-hairpin motif upon processing from its precursor polyprotein Gag. The β-hairpin is essential for correct CA assembly but unexpectedly it is not within any CA oligomeric interfaces. To understand the β-hairpin function we studied the full-length CA protein from equine infectious anemia virus (EIAV), a lentivirus sharing the same cone-shaped capsid core as HIV-1. Solution NMR spectroscopy is perfectly suited to study EIAV-CA that dimerizes weaker than HIV-1-CA. Comparison between the wild-type (wt) EIAV-CA and a variant lacking the β-hairpin structure demonstrated that folding of the β-hairpin specifically extended the N terminus of helix α1 from Tyr²⁰ to Pro¹⁷. This coil to helix transition involves the conserved sequence of Thr¹⁶-Pro¹⁷-Arg¹⁸ (Ser¹⁶-Pro¹⁷-Arg¹⁸ in HIV-1-CA). The extended region of helix α1 constituted an expanded EIAV-CA^N oligomeric interface and overlapped with the HIV-1-CA hexamer-core residue Arg¹⁸, helical in structure and pivotal in assembly. Therefore we propose the function of the maturational refolding of the β-hairpin in CA assembly is to extend helix α1 at the N terminus to enhance the CA^N oligomerization along the capsid assembly core interface. In addition, NMR resonance line broadening indicated the presence of micro-millisecond exchange kinetics due to the EIAV-CA^N domain oligomerization, independent to the faster EIAV-CA^C domain dimerization.

The pore size of PLGA bone implants determines the de novo formation of bone tissue in tibial head defects in rats

PURPOSE

The influence of the pore size of biodegradable poly(lactic-co-glycolic acid) scaffolds on bone regeneration was investigated.

METHODS

Cylindrical poly(lactic-co-glycolic acid) scaffolds were implanted into a defect in the tibial head of rats. Pore sizes of 100-300, 300-500, and 500-710 μm were tested and compared to untreated defects as control. Two and four weeks after implantation, the specimens were explanted and defect regeneration and de novo extracellular matrix generation were investigated by MRI, quantitative solid-state NMR, and mass spectrometry.

RESULTS

The pore size of the scaffolds had a pronounced influence on the quantity of the extracellular matrix synthesized in the graft; most collagen was synthesized within the first 2 weeks of implantation, while the amount of hydroxyapatite increased in the second 2 weeks. After 4 weeks, the scaffolds contained large quantities of newly formed lamellar bone while the control defects were filled by inhomogenous woven bone. Best results were obtained for scaffolds of a pore size of 300-500 μm.

CONCLUSION

Our analysis showed that the structure and dynamics of the regenerated extracellular matrix was very similar to that of the native bone, suggesting that biomineralization was significantly enhanced by the choice of the most appropriate implant material.

Improved chemical shift prediction by Rosetta conformational sampling

Chemical shift frequencies represent a time-average of all the conformational states populated by a protein. Thus, chemical shift prediction programs based on sequence and database analysis yield higher accuracy for rigid rather than flexible protein segments. Here we show that the prediction accuracy can be significantly improved by averaging over an ensemble of structures, predicted solely from amino acid sequence with the Rosetta program. This approach to chemical shift and structure prediction has the potential to be useful for guiding resonance assignments, especially in solid-state NMR structural studies of membrane proteins in proteoliposomes.

3D-TROSY-based backbone and ILV-methyl resonance assignments of a 319-residue homodimer from a single protein sample

The feasibility of practically complete backbone and ILV methyl chemical shift assignments from a single [U-(2)H,(15)N,(13)C; Ileδ1-{(13)CH(3)}; Leu,Val-{(13)CH(3)/(12)CD(3)}]-labeled protein sample of the truncated form of ligand-free Bst-Tyrosyl tRNA Synthetase (Bst-ΔYRS), a 319-residue predominantly helical homodimer, is established. Protonation of ILV residues at methyl positions does not appreciably detract from the quality of TROSY triple resonance data. The assignments are performed at 40 °C to improve the sensitivity of the measurements and alleviate the overlap of (1)H-(15)N correlations in the abundant α-helical segments of the protein. A number of auxiliary approaches are used to assist in the assignment process: (1) selection of (1)H-(15)N amide correlations of certain residue types (Ala, Thr/Ser) that simplifies 2D (1)H-(15)N TROSY spectra, (2) straightforward identification of ILV residue types from the methyl-detected 'out-and-back' HMCM(CG)CBCA experiment, and (3) strong sequential HN-HN NOE connectivities in the helical regions. The two subunits of Bst-YRS were predicted earlier to exist in two different conformations in the absence of ligands. In agreement with our earlier findings (Godoy-Ruiz in J Am Chem Soc 133:19578-195781, 2011), no evidence of dimer asymmetry has been observed in either amide- or methyl-detected experiments.

Assessing and refining molecular dynamics simulations of proteins with nuclear magnetic resonance data

The sophistication of the force fields, algorithms and hardware used for molecular dynamics (MD) simulations of proteins is continuously increasing. No matter how advanced the methodology, however, it is essential to evaluate the appropriateness of the structures sampled in a simulation by comparison with quantitative experimental data. Solution nuclear magnetic resonance (NMR) data are particularly useful for checking the quality of protein simulations, as they provide both structural and dynamic information on a variety of temporal and spatial scales. Here, various features and implications of using NMR data to validate and bias MD simulations are outlined, including an overview of the different types of NMR data that report directly on structural properties and of relevant simulation techniques. The focus throughout is on how to properly account for conformational averaging, particularly within the context of the assumptions inherent in the relationships that link NMR data to structural properties.