Structure and disorder in the ribonuclease S-peptide probed by NMR residual dipolar couplings

Monitoring protein folding through high pressure NMR spectroscopy

High-pressure is a well-known perturbation method used to destabilize globular proteins. It is perfectly reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. In contrast to other perturbation methods such as heat or chemical denaturant that destabilize protein structures uniformly, pressure exerts local effects on regions or domains of a protein containing internal cavities. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility to monitor at a residue level the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. High-pressure NMR experiments can now be routinely performed, owing to the recent development of commercially available high-pressure sample cells. This review summarizes recent advances and some future directions of high-pressure NMR techniques for the characterization at atomic resolution of the energy landscape of protein folding.

RDC-enhanced structure calculation of a β-heptapeptide in methanol

Residual dipolar couplings (RDCs) are a rich source of structural information that goes beyond the range covered by the nuclear Overhauser effect or scalar coupling constants. They can only be measured in partially oriented samples. RDC studies of peptides in organic solvents have so far been focused on samples in chloroform or DMSO. Here, we show that stretched poly(vinyl acetate) can be used for the partial alignment of a linear β-peptide with proteinogenic side chains in methanol. ¹ D_CH , ¹ D_NH , and ² D_HH RDCs were collected with this sample and included as restraints in a simulated annealing calculation. Incorporation of RDCs in the structure calculation process improves the long-range definition in the backbone of the resulting 3₁₄ -helix and uncovers side-chain mobility. Experimental side-chain RDCs of the central leucine and valine residues are in good agreement with predicted values from a local three-state model. Copyright © 2016 John Wiley & Sons, Ltd.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

The illusive search for the lowest free energy state of globular proteins and RNAs

As a consequence of the one-dimensional storage and transfer of genetic information, DNA→RNA→protein, the process by which globular proteins and RNAs achieve their three-dimensional structure involves folding of a linear chain. The folding process itself could create massive activation barriers that prevent the attainment of many stable protein and RNA structures. We consider several kinds of energy barriers inherent in folding that might serve as kinetic constraints to achieving the lowest energy state. Alternative approaches to forming 3D structure, where a substantial number of weak interactions would be created prior to the formation of all the peptide (or phosphodiester) bonds, might not be subjected to such high barriers. This could lead to unique 3D conformational states, potentially more stable than "native" proteins and RNAs, with new functionalities.

Inhibition of semen-derived enhancer of virus infection (SEVI) fibrillogenesis by zinc and copper

Semen-derived enhancer of virus infection (SEVI), a naturally occurring peptide fragment of prostatic acid phosphatase, enhances HIV infectivity by forming cationic amyloid fibrils that aid the fusion of negatively charged virion and target cell membranes. Cu(II) and Zn(II) inhibit fibrillization of SEVI in a kinetic assay using the fibril-specific dye ThT. TEM suggests that the metals do not affect fibril morphology. NMR shows that the metals bind to histidines 3 and 23 in the SEVI sequence. ITC experiments indicate that SEVI forms oligomeric complexes with the metals. Dissociation constants are micromolar for Cu(II) and millimolar for Zn(II). Because the Cu(II) and Zn(II) concentrations that inhibit fibrillization are comparable with those found in seminal fluid the metals may modulate SEVI fibrillization under physiological conditions.

Domain cooperativity in multidomain proteins: what can we learn from molecular alignment in anisotropic media?

Many proteins have modular design with multiple globular domains connected via flexible linkers. As a simple model of such system, we study a tandem construct consisting of two identical SH3 domains and a variable-length Gly/Ser linker. When the linker is short, this construct represents a dumbbell-shaped molecule with limited amount of domain-domain mobility. Due to its elongated shape, this molecule efficiently aligns in steric alignment media. As the length of the linker increases, the two domains become effectively uncoupled and begin to behave as independent entities. Consequently, their degree of alignment drops, approaching that found in the (near-spherical) isolated SH3 domains. To model the dependence of alignment parameters on the length of the interdomain linker, we have generated in silico a series of conformational ensembles representing SH3 tandems with different linker length. These ensembles were subsequently used as input for alignment prediction software PALES. The predicted alignment tensors were compared with the results of experimental measurements using a series of tandem-SH3 samples in PEG/hexanol alignment media. This comparison broadly confirmed the expected trends. At the same time, it has been found that the isolated SH3 domain aligns much stronger than expected. This finding can be attributed to complex morphology of the PEG/hexanol media and/or to weak site-specific interactions between the protein and the media. In the latter case, there are strong indications that electrostatic interactions may play a role. The fact that PEG/hexanol does not behave as a simple steric media should serve as a caution for studies that use PALES as a quantitative prediction tool (especially for disordered proteins). Further progress in this area depends on our ability to accurately model the anisotropic media and its site-specific interactions with protein molecules. Once this ability is improved, it should be possible to use the alignment parameters as a measure of domain-domain cooperativity, thus identifying the situations where two domains transiently interact with each other or become coupled through a partially structured linker.

NMR characterization of long-range order in intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are predicted to represent a significant fraction of the human genome, and the development of meaningful molecular descriptions of these proteins remains a key challenge for contemporary structural biology. In order to describe the conformational behavior of IDPs, a molecular representation of the disordered state based on diverse sources of structural data that often exhibit complex and very different averaging behavior is required. In this study, we propose a combination of paramagnetic relaxation enhancements (PREs) and residual dipolar couplings (RDCs) to define both long-range and local structural features of IDPs in solution. We demonstrate that ASTEROIDS, an ensemble selection algorithm, faithfully reproduces intramolecular contacts, even in the presence of highly diffuse, ill-defined target interactions. We also show that explicit modeling of spin-label mobility significantly improves the reproduction of experimental PRE data, even in the case of highly disordered proteins. Prediction of the effects of transient long-range contacts on RDC profiles reveals that weak intramolecular interactions can induce a severe distortion of the profiles that compromises the description of local conformational sampling if it is not correctly taken into account. We have developed a solution to this problem that involves efficiently combining RDC and PRE data to simultaneously determine long-range and local structure in highly flexible proteins. This combined analysis is shown to be essential for the accurate interpretation of experimental data from alpha-synuclein, an important IDP involved in human neurodegenerative disease, confirming the presence of long-range order between distant regions in the protein.

Quantitative description of backbone conformational sampling of unfolded proteins at amino acid resolution from NMR residual dipolar couplings

An atomic resolution characterization of the structural properties of unfolded proteins that explicitly invokes the highly dynamic nature of the unfolded state will be extremely important for the development of a quantitative understanding of the thermodynamic basis of protein folding and stability. Here we develop a novel approach using residual dipolar couplings (RDCs) from unfolded proteins to determine conformational behavior on an amino acid specific basis. Conformational sampling is described in terms of ensembles of structures selected from a large pool of conformers. We test this approach, using extensive simulation, to determine how well the fitting of RDCs to reduced conformational ensembles containing few copies of the molecule can correctly reproduce the backbone conformational behavior of the protein. Having established approaches that allow accurate mapping of backbone dihedral angle conformational space from RDCs, we apply these methods to obtain an amino acid specific description of ubiquitin denatured in 8 M urea at pH 2.5. Cross-validation of data not employed in the fit verifies that an ensemble size of 200 structures is appropriate to characterize the highly fluctuating backbone. This approach allows us to identify local conformational sampling properties of urea-unfolded ubiquitin, which shows that the backbone sampling of certain types of charged or polar amino acids, in particular threonine, glutamic acid, and arginine, is affected more strongly by urea binding than amino acids with hydrophobic side chains. In general, the approach presented here establishes robust procedures for the study of all denatured and intrinsically disordered states.

Probing the urea dependence of residual structure in denatured human alpha-lactalbumin

Backbone (15)N relaxation parameters and (15)N-(1)H(N) residual dipolar couplings (RDCs) have been measured for a variant of human alpha-lactalbumin (alpha-LA) in 4, 6, 8 and 10 M urea. In the alpha-LA variant, the eight cysteine residues in the protein have been replaced by alanines (all-Ala alpha-LA). This protein is a partially folded molten globule at pH 2 and has been shown previously to unfold in a stepwise non-cooperative manner on the addition of urea. (15)N R(2) values in some regions of all-Ala alpha-LA show significant exchange broadening which is reduced as the urea concentration is increased. Experimental RDC data are compared with RDCs predicted from a statistical coil model and with bulkiness, average area buried upon folding and hydrophobicity profiles in order to identify regions of non-random structure. Residues in the regions corresponding to the B, D and C-terminal 3(10) helices in native alpha-LA show R(2) values and RDC data consistent with some non-random structural propensities even at high urea concentrations. Indeed, for residues 101-106 the residual structure persists in 10 M urea and the RDC data suggest that this might include the formation of a turn-like structure. The data presented here allow a detailed characterization of the non-cooperative unfolding of all-Ala alpha-LA at higher concentrations of denaturant and complement previous studies which focused on structural features of the molten globule which is populated at lower concentrations of denaturant.

Insights into the dynamics of specific telomeric single-stranded DNA recognition by Pot1pN

The N-terminal oligonucleotide/oligosaccharide-binding fold domain of the Schizosaccharomyces pombe protection of telomeres 1 (Pot1) protein, Pot1pN (residues 1-187 of full-length Pot1), specifically recognizes telomeric single-stranded DNA (ssDNA) via a complex series of molecular interactions that are punctuated by unusual internucleotide hydrogen bonds. While the structure of ssDNA-bound Pot1pN provides an initial model for understanding how the Pot1pN-ssDNA complex is assembled and how specific nucleotide recognition occurs, further refinement requires knowledge of the ssDNA-free state of Pot1pN and the dynamic changes that accompany the binding of ssDNA. Using NMR strategies, we found that ssDNA-free Pot1pN adopts a similar overall protein backbone topology as ssDNA-bound Pot1pN does. Although the backbone structure remained relatively unchanged, we observed unexpected differential dynamic changes within the ssDNA-binding pockets of Pot1pN upon binding of cognate ssDNA. These studies support a model in which conformational selection and induced fit play important roles in the recognition of ssDNA by Pot1pN. Furthermore, the studies presented here provide a more comprehensive understanding of how specific nucleotide recognition is achieved by the telomere-end protection family of essential proteins.

Model-independent interpretation of NMR relaxation data for unfolded proteins: the acid-denatured state of ACBP

We have investigated the acid-unfolded state of acyl-coenzyme A binding protein (ACBP) using 15N laboratory frame nuclear magnetic resonance (NMR) relaxation experiments at three magnetic field strengths. The data have been analyzed using standard model-free fitting and models involving distribution of correlation times. In particular, a model-independent method of analysis that does not assume any analytical form for the correlation time distribution is proposed. This method explains correlations between model-free parameters and the analytical distribution parameters found by other authors. The analysis also shows that the relaxation data are consistent with and complementary to information obtained from other parameters, especially secondary chemical shifts and residual dipolar couplings, and strengthens the conclusions of previous observations that three out of the four regions that form helices in the native structure appear to contain residual secondary structure also in the acid-denatured state.

Back to the future: Ribonuclease A

Pancreatic ribonuclease A (EC 3.1.27.5, RNase) is, perhaps, the best-studied enzyme of the 20th century. It was isolated by René Dubos, crystallized by Moses Kunitz, sequenced by Stanford Moore and William Stein, and synthesized in the laboratory of Bruce Merrifield, all at the Rockefeller Institute/University. It has proven to be an excellent model system for many different types of experiments, both as an enzyme and as a well-characterized protein for biophysical studies. Of major significance was the demonstration by Chris Anfinsen at NIH that the primary sequence of RNase encoded the three-dimensional structure of the enzyme. Many other prominent protein chemists/enzymologists have utilized RNase as a dominant theme in their research. In this review, the history of RNase and its offspring, RNase S (S-protein/S-peptide), will be considered, especially the work in the Merrifield group, as a preface to preliminary data and proposed experiments addressing topics of current interest. These include entropy-enthalpy compensation, entropy of ligand binding, the impact of protein modification on thermal stability, and the role of protein dynamics in enzyme action. In continuing to use RNase as a prototypical enzyme, we stand on the shoulders of the giants of protein chemistry to survey the future.

Theoretical framework for NMR residual dipolar couplings in unfolded proteins

A theoretical framework for the prediction of nuclear magnetic resonance (NMR) residual dipolar couplings (RDCs) in unfolded proteins under weakly aligning conditions is presented. The unfolded polypeptide chain is modeled as a random flight chain while the alignment medium is represented by a set of regularly arranged obstacles. For the case of bicelles oriented perpendicular to the magnetic field, a closed-form analytical result is derived. With the obtained analytical expression the RDCs are readily accessible for any locus along the chain, for chains of differing length, and for varying bicelle concentrations. The two general features predicted by the model are (i) RDCs in the center segments of a polypeptide chain are larger than RDCs in the end segments, resulting in a bell-shaped sequential distribution of RDCs, and (ii) couplings are larger for shorter chains than for longer chains at a given bicelle concentration. Experimental data available from the literature confirm the first prediction of the model, providing a tool for recognizing fully unfolded polypeptide chains. With less certainty experimental data appear to support the second prediction as well. However, more systematic experimental studies are needed in order to validate or disprove the predictions of the model. The presented framework is an important step towards a solid theoretical foundation for the analysis of experimentally measured RDCs in unfolded proteins in the case of alignment media such as polyacrylamide gels and neutral bicelle systems which align biomacromolecules by a steric mechanism. Various improvements and generalizations are possible within the suggested approach.

The Cytoplasmic Domain of the Chloride Channel ClC-0: Structural and Dynamic Characterization of Flexible Regions

Eukaryotic members of the ClC family of chloride channels and transporters are composed of a transmembrane ion transport domain followed by a cytoplasmic domain, which is believed to be involved in the modulation of ClC function. In some family members this putative regulatory domain contains next to a well-folded structured part, long sequence stretches with low sequence complexity. These regions, a 96 residue long linker connecting two structured sub-domains, and 35 residues on the C teminus of the domain were found disordered in a recent crystal structure of this domain in ClC-0. Both regions have a large influence on the modulation of channel function in closely related family members. Here we describe a NMR study to characterize the structural and dynamic properties of these putatively unstructured stretches. Our study reveals that the two regions indeed show large conformational flexibility with dynamics on the nanosecond timescale. However, small islands of secondary structure are found interdispersed between the unfolded regions. This study characterizes for the first time the biophysical properties of these protein segments, which may become important for the understanding of novel regulatory mechanisms within the ClC family.

Conformational Preference of Polyglycine in Solution to Elongated Structure

Do polypeptide chains ever behave like a random coil? In this report we demonstrate that glycine, the residue with the fewest backbone restrictions, exhibits a strong preference for an extended conformation in solution when polymerized in short segments of polyglycine. A model peptide system comprised of two unique tripeptide units, between which 1 to 18 glycine residues are inserted, is characterized by NMR and by small-angle X-ray scattering (SAXS). The residual dipolar coupling (RDC) values of the two tripeptide units are insensitive to changes in number of intervening glycines, suggesting that extension of the linker does not alter the average angular relationship between the tripeptides. Polyglycine segments longer than nine residues form insoluble aggregates. SAXS measurements using synchrotron radiation provide direct evidence that polyglycine peptides adopt elongated conformations. In particular, the construct with a linker with six glycines showed a scattering profile indicative of a monomeric state with a radius of gyration and the maximum dimension of 9.1 A and approximately 34 A, respectively. The ensemble averaged global structure of this 12-mer peptide can best be approximated by a cylinder with a radius of 4 A and a length of approximately 33 A, making it intermediate in extension between a beta strand and an alpha helix.

Foldon, the natural trimerization domain of T4 fibritin, dissociates into a monomeric A-state form containing a stable beta-hairpin: atomic details of trimer dissociation and local beta-hairpin stability from residual dipolar couplings

The C-terminal domain of T4 fibritin (foldon) is obligatory for the formation of the fibritin trimer structure and can be used as an artificial trimerization domain. Its native structure consists of a trimeric beta-hairpin propeller. At low pH, the foldon trimer disintegrates into a monomeric (A-state) form that has similar properties as that of an early intermediate of the trimer folding pathway. The formation of this A-state monomer from the trimer, its structure, thermodynamic stability, equilibrium association and folding dynamics have been characterized to atomic detail by modern high-resolution NMR techniques. The foldon A-state monomer forms a beta-hairpin with intact and stable H-bonds that is similar to the monomer in the foldon trimer, but lacks a defined structure in its N and C-terminal parts. Its thermodynamic stability in pure water is comparable to designed hairpins stabilized in alcohol/water mixtures. Details of the thermal unfolding of the foldon A-state have been characterized by chemical shifts and residual dipolar couplings (RDCs) detected in inert, mechanically stretched polyacrylamide gels. At the onset of the thermal transition, uniform relative changes in RDC values indicate a uniform decrease of local N-HN and Calpha-Halpha order parameters for the hairpin strand residues. In contrast, near-turn residues show particular thermal stability in RDC values and hence in local order parameters. This coincides with increased transition temperatures of the beta-turn residues observed by chemical shifts. At high temperatures, the RDCs converge to non-zero average values consistent with predictions from random chain polymer models. Residue-specific deviations above the unfolding transition reveal the persistence of residual order around proline residues, large hydrophobic residues and at the beta-turn.

Determination of protein backbone structures from residual dipolar couplings

There are a number of circumstances in which a focus on determination of the backbone structure of a protein, as opposed to a complete all-atom structure, may be appropriate. This is particularly the case for structures determined as a part of a structural genomics initiative in which computational modeling of many sequentially related structures from the backbone of a single family representative is anticipated. It is, however, also the case when the backbone may be a stepping-stone to more targeted studies of ligand interaction or protein-protein interaction. Here an NMR protocol is described that can produce a backbone structure of a protein without the need for extensive experiments directed at side chain resonance assignment or the collection of structural information on side chains. The procedure relies primarily on orientational constraints from residual dipolar couplings as opposed to distance constraints from NOEs. Procedures for sample preparation, data acquisition, and data analysis are described, along with examples from application to small target proteins of a structural genomics project.

Structural Characterization of Unfolded States of Apomyoglobin using Residual Dipolar Couplings

The conformational propensities of unfolded states of apomyoglobin have been investigated by measurement of residual dipolar couplings between (15)N and (1)H in backbone amide groups. Weak alignment of apomyoglobin in acid and urea-unfolded states was induced with both stretched and compressed polyacrylamide gels. In 8 M urea solution at pH 2.3, conditions under which apomyoglobin contains no detectable secondary or tertiary structure, significant residual dipolar couplings of uniform sign were observed for all residues. At pH 2.3 in the absence of urea, a change in the magnitude and/or sign of the residual dipolar couplings occurs in local regions of the polypeptide where there is a high propensity for helical secondary structure. These results are interpreted on the basis of the statistical properties of the unfolded polypeptide chain, viewed as a polymer of statistical segments. For a folded protein, the magnitude and sign of the residual dipolar couplings depend on the orientation of each bond vector relative to the alignment tensor of the entire molecule, which reorients as a single entity. For unfolded proteins, there is no global alignment tensor; instead, residual dipolar couplings are attributed to alignment of the statistical segments or of transient elements of secondary structure. For apomyoglobin in 8 M urea, the backbone is highly extended, with phi and psi dihedral angles favoring the beta or P(II) regions. Each statistical segment has a highly anisotropic shape, with the N-H bond vectors approximately perpendicular to the long axis, and becomes weakly aligned in the anisotropic environment of the strained acrylamide gels. Local regions of enhanced flexibility or chain compaction are characterized by a decrease in the magnitude of the residual dipolar couplings. The formation of a small population of helical structure in the acid-denatured state of apomyoglobin leads to a change in sign of the residual dipolar couplings in local regions of the polypeptide; the population of helix estimated from the residual dipolar couplings is in excellent agreement with that determined from chemical shifts. The alignment model described here for apomyoglobin can also explain the pattern of residual dipolar couplings reported previously for denatured states of staphylococcal nuclease and other proteins. In conjunction with other NMR experiments, residual dipolar couplings can provide valuable insights into the dynamic conformational propensities of unfolded and partly folded states of proteins and thereby help to chart the upper reaches of the folding landscape.

Short-range, Long-range and Transition State Interactions in the Denatured State of ACBP from Residual Dipolar Couplings

Residual dipolar couplings in the denatured state of bovine acyl-coenzyme A binding protein (ACBP) oriented in strained polyacrylamide gels have been shown to be a sensitive, sequence-specific probe for residual secondary structure. Results supporting this were obtained by comparing residual dipolar couplings under different denaturing conditions. The data were analyzed using the program molecular fragment replacement (MFR), which demonstrated alpha-helix propensity in four isolated stretches along the protein backbone, and these coincide with the location of native helices. This is in full agreement with earlier findings based on secondary chemical shift values. Furthermore, N-H residual dipolar couplings provided direct evidence for the existence of native-like hydrophobic interactions in the acid-denatured state of ACBP at pH 2.3. It was shown that replacement of the hydrophobic side-chain of residue Ile27 with alanine in helix A2 leads to large decreases of residual dipolar couplings in residues that form helix A4 in the native state. It is suggested that the Ile to Ala mutation changes the probability for the formation of long-range interactions, which are present in the acid-denatured state of the wild-type protein. These long-range interactions are similar to those proposed to form in the transition state of folding of ACBP. Therefore, the application of residual dipolar couplings in combination with a comparative mutation study has demonstrated the presence of precursors to the folding transition state under acid-unfolding conditions.