Validation of Protein Structure from Anisotropic Carbonyl Chemical Shifts in a Dilute Liquid Crystalline Phase

On the structural denaturation of biological analytes in trapped ion mobility spectrometry – mass spectrometry

Trapped ion mobility spectra recorded for ubiquitin are consistent with structures reported for the native state by NMR.

Kinetics from Implicit Solvent Simulations of Biomolecules as a Function of Viscosity

Kinetic properties of alanine dipeptide, the B1 domain of streptococcal protein G, and ubiquitin are compared between explicit solvent and implicit solvent simulations with the generalized Born molecular volume (GBMV) method. The results indicate that kinetics from explicit solvent simulations and experiments can be matched closely when the implicit solvent simulations are combined with Langevin dynamics and a friction coefficient near 10 ps(-1). Smaller and larger friction coefficients accelerate and slow down conformational sampling. It is found that local conformational exploration without the crossing of significant barriers can be accelerated by a factor of 4-5; however, the acceleration of barrier crossings is limited to about a factor of 2. The use of a Nosé-Hoover thermostat instead of Langevin dynamics greatly enhances local conformational sampling but slows down the crossing of barriers by at least an order of magnitude because of the lack of solute-solvent stochastic collisions.

Direct Comparison of Experimental and Calculated NMR Scalar Coupling Constants for Force Field Validation and Adaptation

The ability to measure scalar coupling constants across hydrogen bonds ((3h)JNC') from high-resolution NMR experiments allows the characterization of detailed structural properties of biomolecules. To analyze those, a parametrized model based on the linear combination of atomic orbitals relates H-bond geometries with the measured (3h)JNC' coupling magnitude. In the present study the dependence of calculated (3h)JNC' coupling constants on force field parameters is assessed. It is shown that increased polarity of the hydrogen bond improves the calculated (3h)JNC' coupling constants and shifts the conformational ensemble sampled from the molecular dynamics (MD) simulations toward the experimentally measured one. Increased charges lead to more narrow distance and angle distributions and improve the agreement between calculated and measured (3h)JNC' couplings. However, different secondary structures are better represented by different magnitudes of electrostatic interactions-different atomic partial charges in the present work-as indicated by root-mean square deviations (rsmds) between observed and calculated coupling constants (3h)JNC'. The parametrization of the empirical formula is found to be meaningful and robust, but the parameter values are not universal across different proteins and different secondary structural elements (α-helices, β-sheets and loops). Using standard and slightly increased CHARMM charges, predictions for the as-yet unknown scalar coupling constants for the V54A and I6A mutants of protein G are made.

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.

Conformational analysis of small organic molecules using NOE and RDC data: A discussion of strychnine and α-methylene-γ-butyrolactone

To understand the properties and/or reactivity of an organic molecule, an understanding of its three-dimensional structure is necessary. Simultaneous determination of configuration and conformation often poses a daunting challenge. Thus, the more information accessible for a given molecule, the better. Additionally to (3)J-couplings, two sources of information, quantitative NOE and more recently also RDCs, are used for conformational analysis by NMR spectroscopy. In this paper, we compare these sources of conformational information in two molecules: the configurationally well-characterized strychnine 1, and the only recently configurationally and conformationally characterized α-methylene-γ-butyrolactone 2. We discuss possible sources of error in the measurement and analysis process, and how to exclude them. By this means, we are able to bolster the previously proposed flexibility for these two molecules.

Ensembles of a small number of conformations with relative populations

In our previous work, we proposed a new way to represent protein native states, using ensembles of a small number of conformations with relative Populations, or ESP in short. Using Ubiquitin as an example, we showed that using a small number of conformations could greatly reduce the potential of overfitting and assigning relative populations to protein ensembles could significantly improve their quality. To demonstrate that ESP indeed is an excellent alternative to represent protein native states, in this work we compare the quality of two ESP ensembles of Ubiquitin with several well-known regular ensembles or average structure representations. Extensive amount of significant experimental data are employed to achieve a thorough assessment. Our results demonstrate that ESP ensembles, though much smaller in size comparing to regular ensembles, perform equally or even better sometimes in all four different types of experimental data used in the assessment, namely, the residual dipolar couplings, residual chemical shift anisotropy, hydrogen exchange rates, and solution scattering profiles. This work further underlines the significance of having relative populations in describing the native states.

Spatial reorientation experiments for NMR of solids and partially oriented liquids

Motional reorientation experiments are extensions of Magic Angle Spinning (MAS) where the rotor axis is changed in order to average out, reintroduce, or scale anisotropic interactions (e.g. dipolar couplings, quadrupolar interactions or chemical shift anisotropies). This review focuses on Variable Angle Spinning (VAS), Switched Angle Spinning (SAS), and Dynamic Angle Spinning (DAS), all of which involve spinning at two or more different angles sequentially, either in successive experiments or during a multidimensional experiment. In all of these experiments, anisotropic terms in the Hamiltonian are scaled by changing the orientation of the spinning sample relative to the static magnetic field. These experiments vary in experimental complexity and instrumentation requirements. In VAS, many one-dimensional spectra are collected as a function of spinning angle. In SAS, dipolar couplings and/or chemical shift anisotropies are reintroduced by switching the sample between two different angles, often 0° or 90° and the magic angle, yielding a two-dimensional isotropic-anisotropic correlation spectrum. Dynamic Angle Spinning (DAS) is a related experiment that is used to simultaneously average out the first- and second-order quadrupolar interactions, which cannot be accomplished by spinning at any unique rotor angle in physical space. Although motional reorientation experiments generally require specialized instrumentation and data analysis schemes, some are accessible with only minor modification of standard MAS probes. In this review, the mechanics of each type of experiment are described, with representative examples. Current and historical probe and coil designs are discussed from the standpoint of how each one accomplishes the particular objectives of the experiment(s) it was designed to perform. Finally, applications to inorganic materials and liquid crystals, which present very different experimental challenges, are discussed. The review concludes with perspectives on how motional reorientation experiments can be applied to current problems in chemistry, molecular biology, and materials science, given the many advances in high-field NMR magnets, fast spinning, and sample preparation realized in recent years.

ProCS15: a DFT-based chemical shift predictor for backbone and Cβ atoms in proteins

We present ProCS15: a program that computes the isotropic chemical shielding values of backbone and Cβ atoms given a protein structure in less than a second. ProCS15 is based on around 2.35 million OPBE/6-31G(d,p)//PM6 calculations on tripeptides and small structural models of hydrogen-bonding. The ProCS15-predicted chemical shielding values are compared to experimentally measured chemical shifts for Ubiquitin and the third IgG-binding domain of Protein G through linear regression and yield RMSD values of up to 2.2, 0.7, and 4.8 ppm for carbon, hydrogen, and nitrogen atoms. These RMSD values are very similar to corresponding RMSD values computed using OPBE/6-31G(d,p) for the entire structure for each proteins. These maximum RMSD values can be reduced by using NMR-derived structural ensembles of Ubiquitin. For example, for the largest ensemble the largest RMSD values are 1.7, 0.5, and 3.5 ppm for carbon, hydrogen, and nitrogen. The corresponding RMSD values predicted by several empirical chemical shift predictors range between 0.7–1.1, 0.2–0.4, and 1.8–2.8 ppm for carbon, hydrogen, and nitrogen atoms, respectively.

The Structural Role of Antibody N-Glycosylation in Receptor Interactions

Asparagine(N)297-linked glycosylation of immunoglobulin G (IgG) Fc is required for binding to FcγRIIa, IIb, and IIIa, although it is unclear how it contributes. We found the quaternary structure of glycosylated Fc was indistinguishable from aglycosylated Fc, indicating that N-glycosylation does not maintain relative Fc Cγ2/Cγ3 domain orientation. However, the conformation of the C'E loop, which contains N297, was significantly perturbed in the aglycosylated Fc variant. The conformation of the C'E loop as measured with a range of Fc variants shows a strong correlation with FcγRIIIa affinity. These results indicate that the primary role of the IgG1 Fc N-glycan is to stabilize the C'E loop through intramolecular interactions between carbohydrate and amino acid residues, and preorganize the FcγRIIIa interface for optimal binding affinity. The features that contribute to the capacity of the IgG1 Fc N-glycan to restrict protein conformation and tune binding affinity are conserved in other antibodies including IgG2-IgG4, IgD, IgE, and IgM.

Assessment of the utility of contact-based restraints in accelerating the prediction of protein structure using molecular dynamics simulations

Molecular dynamics (MD) simulation is a well-established tool for the computational study of protein structure and dynamics, but its application to the important problem of protein structure prediction remains challenging, in part because extremely long timescales can be required to reach the native structure. Here, we examine the extent to which the use of low-resolution information in the form of residue-residue contacts, which can often be inferred from bioinformatics or experimental studies, can accelerate the determination of protein structure in simulation. We incorporated sets of 62, 31, or 15 contact-based restraints in MD simulations of ubiquitin, a benchmark system known to fold to the native state on the millisecond timescale in unrestrained simulations. One-third of the restrained simulations folded to the native state within a few tens of microseconds-a speedup of over an order of magnitude compared with unrestrained simulations and a demonstration of the potential for limited amounts of structural information to accelerate structure determination. Almost all of the remaining ubiquitin simulations reached near-native conformations within a few tens of microseconds, but remained trapped there, apparently due to the restraints. We discuss potential methodological improvements that would facilitate escape from these near-native traps and allow more simulations to quickly reach the native state. Finally, using a target from the Critical Assessment of protein Structure Prediction (CASP) experiment, we show that distance restraints can improve simulation accuracy: In our simulations, restraints stabilized the native state of the protein, enabling a reasonable structural model to be inferred.

The Quest for Simplicity: Remarks on the Free-Approach Models

Nuclear magnetic relaxation provides a powerful method giving insight into molecular motions at atomic resolution on a broad time scale. Dynamics of biological macromolecules has been widely exploited by measuring (15)N and (13)C relaxation data. Interpretation of these data relies almost exclusively on the use of the model-free approach (MFA) and its extended version (EMFA) which requires no particular physical model of motion and a small number of parameters. It is shown that EMFA is often unable to cope with three different time scales and fails to describe slow internal motions properly. In contrast to EMFA, genuine MFA with two time scales can reproduce internal motions slower than the overall tumbling. It is also shown that MFA and simplified EMFA are equivalent with respect to the values of the N-H bond length and chemical shift anisotropy. Therefore, the vast majority of (15)N relaxation data for proteins can be satisfactorily interpreted solely with MFA.

Validating solution ensembles from molecular dynamics simulation by wide-angle X-ray scattering data

Wide-angle x-ray scattering (WAXS) experiments of biomolecules in solution have become increasingly popular because of technical advances in light sources and detectors. However, the structural interpretation of WAXS profiles is problematic, partly because accurate calculations of WAXS profiles from structural models have remained challenging. In this work, we present the calculation of WAXS profiles from explicit-solvent molecular dynamics (MD) simulations of five different proteins. Using only a single fitting parameter that accounts for experimental uncertainties because of the buffer subtraction and dark currents, we find excellent agreement to experimental profiles both at small and wide angles. Because explicit solvation eliminates free parameters associated with the solvation layer or the excluded solvent, which would require fitting to experimental data, we minimize the risk of overfitting. We further find that the influence from water models and protein force fields on calculated profiles are insignificant up to q≈15nm(-1). Using a series of simulations that allow increasing flexibility of the proteins, we show that incorporating thermal fluctuations into the calculations significantly improves agreement with experimental data, demonstrating the importance of protein dynamics in the interpretation of WAXS profiles. In addition, free MD simulations up to one microsecond suggest that the calculated profiles are highly sensitive with respect to minor conformational rearrangements of proteins, such as an increased flexibility of a loop or an increase of the radius of gyration by < 1%. The present study suggests that quantitative comparison between MD simulations and experimental WAXS profiles emerges as an accurate tool to validate solution ensembles of biomolecules.

Less is more: structures of difficult targets with minimal constraints

By merging recent experimental and computational methodology advances, resolution-adapted structural recombination Rosetta has emerged as a powerful strategy for solving the structure of traditionally challenging targets. In this issue of Structure, Sgourakis and colleagues solve the structure of one such target, the immunoevasin protein m04, using this approach.

Capturing conformational States in proteins using sparse paramagnetic NMR data

Capturing conformational changes in proteins or protein-protein complexes is a challenge for both experimentalists and computational biologists. Solution nuclear magnetic resonance (NMR) is unique in that it permits structural studies of proteins under greatly varying conditions, and thus allows us to monitor induced structural changes. Paramagnetic effects are increasingly used to study protein structures as they give ready access to rich structural information of orientation and long-range distance restraints from the NMR signals of backbone amides, and reliable methods have become available to tag proteins with paramagnetic metal ions site-specifically and at multiple sites. In this study, we show how sparse pseudocontact shift (PCS) data can be used to computationally model conformational states in a protein system, by first identifying core structural elements that are not affected by the environmental change, and then computationally completing the remaining structure based on experimental restraints from PCS. The approach is demonstrated on a 27 kDa two-domain NS2B-NS3 protease system of the dengue virus serotype 2, for which distinct closed and open conformational states have been observed in crystal structures. By changing the input PCS data, the observed conformational states in the dengue virus protease are reproduced without modifying the computational procedure. This data driven Rosetta protocol enables identification of conformational states of a protein system, which are otherwise difficult to obtain either experimentally or computationally.

Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering, Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation

The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.

A novel RNA binding surface of the TAM domain of TIP5/BAZ2A mediates epigenetic regulation of rRNA genes

The chromatin remodeling complex NoRC, comprising the subunits SNF2h and TIP5/BAZ2A, mediates heterochromatin formation at major clusters of repetitive elements, including rRNA genes, centromeres and telomeres. Association with chromatin requires the interaction of the TAM (TIP5/ARBP/MBD) domain of TIP5 with noncoding RNA, which targets NoRC to specific genomic loci. Here, we show that the NMR structure of the TAM domain of TIP5 resembles the fold of the MBD domain, found in methyl-CpG binding proteins. However, the TAM domain exhibits an extended MBD fold with unique C-terminal extensions that constitute a novel surface for RNA binding. Mutation of critical amino acids within this surface abolishes RNA binding in vitro and in vivo. Our results explain the distinct binding specificities of TAM and MBD domains to RNA and methylated DNA, respectively, and reveal structural features for the interaction of NoRC with non-coding RNA.

Structural heterogeneity in microcrystalline ubiquitin studied by solid-state NMR

By applying [1-(13) C]- and [2-(13) C]-glucose labeling schemes to the folded globular protein ubiquitin, a strong reduction of spectral crowding and increase in resolution in solid-state NMR (ssNMR) spectra could be achieved. This allowed spectral resonance assignment in a straightforward manner and the collection of a wealth of long-range distance information. A high precision solid-state NMR structure of microcrystalline ubiquitin was calculated with a backbone rmsd of 1.57 to the X-ray structure and 1.32 Å to the solution NMR structure. Interestingly, we can resolve structural heterogeneity as the presence of three slightly different conformations. Structural heterogeneity is most significant for the loop region β1-β2 but also for β-strands β1, β2, β3, and β5 as well as for the loop connecting α1 and β3. This structural polymorphism observed in the solid-state NMR spectra coincides with regions that showed dynamics in solution NMR experiments on different timescales.

Generation of pseudocontact shifts in proteins with lanthanides using small "clickable" nitrilotriacetic acid and iminodiacetic acid tags

Pseudocontact shifts (PCS) induced by paramagnetic lanthanide ions provide unique long-range structural information in nuclear magnetic resonance (NMR) spectra, but the site-specific attachment of lanthanide tags to proteins remains a challenge. Here we incorporated p-azido-phenylalanine (AzF) site-specifically into the proteins ubiquitin and GB1, and ligated the AzF residue with alkyne derivatives of small nitrilotriacetic acid and iminodiacetic acid tags using the Cu(I) -catalysed "click" reaction. These tags form lanthanide complexes with no or only a small net charge and produced sizeable PCSs with paramagnetic lanthanide ions in all mutants tested. The PCSs were readily fitted by single magnetic susceptibility anisotropy tensors. Protein precipitation during the click reaction was greatly alleviated by the presence of 150 mM NaCl.

MESMER: minimal ensemble solutions to multiple experimental restraints

MOTIVATION

Macromolecular structures and interactions are intrinsically heterogeneous, temporally adopting a range of configurations that can confound the analysis of data from bulk experiments. To obtain quantitative insights into heterogeneous systems, an ensemble-based approach can be employed, in which predicted data computed from a collection of models is compared to the observed experimental results. By simultaneously fitting orthogonal structural data (e.g. small-angle X-ray scattering, nuclear magnetic resonance residual dipolar couplings, dipolar electron-electron resonance spectra), the range and population of accessible macromolecule structures can be probed.

RESULTS

We have developed MESMER, software that enables the user to identify ensembles that can recapitulate experimental data by refining thousands of component collections selected from an input pool of potential structures. The MESMER suite includes a powerful graphical user interface (GUI) to streamline usage of the command-line tools, calculate data from structure libraries and perform analyses of conformational and structural heterogeneity. To allow for incorporation of other data types, modular Python plugins enable users to compute and fit data from nearly any type of quantitative experimental data.

RESULTS

Conformational heterogeneity in three macromolecular systems was analyzed with MESMER, demonstrating the utility of the streamlined, user-friendly software.

AVAILABILITY AND IMPLEMENTATION

https://code.google.com/p/mesmer/

A minor conformation of a lanthanide tag on adenylate kinase characterized by paramagnetic relaxation dispersion NMR spectroscopy

NMR relaxation dispersion techniques provide a powerful method to study protein dynamics by characterizing lowly populated conformations that are in dynamic exchange with the major state. Paramagnetic NMR is a versatile tool for investigating the structures and dynamics of proteins. These two techniques were combined here to measure accurate and precise pseudocontact shifts of a lowly populated conformation. This method delivers valuable long-range structural restraints for higher energy conformations of macromolecules in solution. Another advantage of combining pseudocontact shifts with relaxation dispersion is the increase in the amplitude of dispersion profiles. Lowly populated states are often involved in functional processes, such as enzyme catalysis, signaling, and protein/protein interactions. The presented results also unveil a critical problem with the lanthanide tag used to generate paramagnetic relaxation dispersion effects in proteins, namely that the motions of the tag can interfere severely with the observation of protein dynamics. The two-point attached CLaNP-5 lanthanide tag was linked to adenylate kinase. From the paramagnetic relaxation dispersion only motion of the tag is observed. The data can be described accurately by a two-state model in which the protein-attached tag undergoes a 23° tilting motion on a timescale of milliseconds. The work demonstrates the large potential of paramagnetic relaxation dispersion and the challenge to improve current tags to minimize relaxation dispersion from tag movements.

FANTEN: a new web-based interface for the analysis of magnetic anisotropy-induced NMR data

Pseudocontact shifts (PCSs) and residual dipolar couplings (RDCs) arising from the presence of paramagnetic metal ions in proteins as well as RDCs due to partial orientation induced by external orienting media are nowadays routinely measured as a part of the NMR characterization of biologically relevant systems. PCSs and RDCs are becoming more and more popular as restraints (1) to determine and/or refine protein structures in solution, (2) to monitor the extent of conformational heterogeneity in systems composed of rigid domains which can reorient with respect to one another, and (3) to obtain structural information in protein-protein complexes. The use of both PCSs and RDCs proceeds through the determination of the anisotropy tensors which are at the origin of these NMR observables. A new user-friendly web tool, called FANTEN (Finding ANisotropy TENsors), has been developed for the determination of the anisotropy tensors related to PCSs and RDCs and has been made freely available through the WeNMR ( http://fanten-enmr.cerm.unifi.it:8080 ) gateway. The program has many new features not available in other existing programs, among which the possibility of a joint analysis of several sets of PCS and RDC data and the possibility to perform rigid body minimizations.

Conformation of inhibitor-free HIV-1 protease derived from NMR spectroscopy in a weakly oriented solution

Flexibility of the glycine-rich flaps is known to be essential for catalytic activity of the HIV-1 protease, but their exact conformations at the different stages of the enzymatic pathway remain subject to much debate. Although hundreds of crystal structures of protease-inhibitor complexes have been solved, only about a dozen inhibitor-free protease structures have been reported. These latter structures reveal a large diversity of flap conformations, ranging from closed to semi-open to wide open. To evaluate the average structure in solution, we measured residual dipolar couplings (RDCs) and compared these to values calculated for crystal structures representative of the closed, semi-open, and wide-open states. The RDC data clearly indicate that the inhibitor-free protease, on average, adopts a closed conformation in solution that is very similar to the inhibitor-bound state. By contrast, a highly drug-resistant protease mutant, PR20, adopts the wide-open flap conformation.

Improvements to REDCRAFT: a software tool for simultaneous characterization of protein backbone structure and dynamics from residual dipolar couplings

Within the past two decades, there has been an increase in the acquisition of residual dipolar couplings (RDC) for investigations of biomolecular structures. Their use however is still not as widely adopted as the traditional methods of structure determination by NMR, despite their potential for extending the limits in studies that examine both the structure and dynamics of biomolecules. This is in part due to the difficulties associated with the analysis of this information-rich data type. The software analysis tool REDCRAFT was previously introduced to address some of these challenges. Here we describe and evaluate a number of additional features that have been incorporated in order to extend its computational and analytical capabilities. REDCRAFT's more traditional enhancements integrate a modified steric collision term, as well as structural refinement in the rotamer space. Other, non-traditional improvements include: the filtering of viable structures based on relative order tensor estimates, decimation of the conformational space based on structural similarity, and forward/reverse folding of proteins. Utilizing REDCRAFT's newest features we demonstrate de-novo folding of proteins 1D3Z and 1P7E to within less than 1.6 Å of the corresponding X-ray structures, using as many as four RDCs per residue and as little as two RDCs per residue, in two alignment media. We also show the successful folding of a structure to less than 1.6 Å of the X-ray structure using {C(i-1)-N(i), N(i)-H(i), and C(i-1)-H(i)} RDCs in one alignment medium, and only {N(i)-H(i)} in the second alignment medium (a set of data which can be collected on deuterated samples). The program is available for download from our website at http://ifestos.cse.sc.edu .

Dss1 is a 26S proteasome ubiquitin receptor

The ubiquitin-proteasome system is the major pathway for protein degradation in eukaryotic cells. Proteins to be degraded are conjugated to ubiquitin chains that act as recognition signals for the 26S proteasome. The proteasome subunits Rpn10 and Rpn13 are known to bind ubiquitin, but genetic and biochemical data suggest the existence of at least one other substrate receptor. Here, we show that the phylogenetically conserved proteasome subunit Dss1 (Sem1) binds ubiquitin chains linked by K63 and K48. Atomic resolution data show that Dss1 is disordered and binds ubiquitin by binding sites characterized by acidic and hydrophobic residues. The complementary binding region in ubiquitin is composed of a hydrophobic patch formed by I13, I44, and L69 flanked by two basic regions. Mutations in the ubiquitin-binding site of Dss1 cause growth defects and accumulation of ubiquitylated proteins.

•

Dss1 is a ubiquitin-binding protein

•

Dss1 binds ubiquitin via an intrinsically disordered region

•

The ubiquitin-binding activity of Dss1 is required for function

NMR structure validation in relation to dynamics and structure determination

NMR spectroscopy is a key technique for understanding the behaviour of proteins, especially highly dynamic proteins that adopt multiple conformations in solution. Overall, protein structures determined from NMR spectroscopy data constitute just over 10% of the Protein Data Bank archive. This review covers the validation of these NMR protein structures, but rather than describing currently available methodology, it focuses on concepts that are important for understanding where and how validation is most relevant. First, the inherent characteristics of the protein under study have an influence on quality and quantity of the distinct types of data that can be acquired from NMR experiments. Second, these NMR data are necessarily transformed into a model for use in a structure calculation protocol, and the protein structures that result from this reflect the types of NMR data used as well as the protein characteristics. The validation of NMR protein structures should therefore take account, wherever possible, of the inherent behavioural characteristics of the protein, the types of available NMR data, and the calculation protocol. These concepts are discussed in the context of 'knowledge based' and 'model versus data' validation, with suggestions for questions to ask and different validation categories to consider. The principal aim of this review is to stimulate discussion and to help the reader understand the relationships between the above elements in order to make informed decisions on which validation approaches are the most relevant in particular cases.

Towards understanding the molecular recognition process in prokaryotic zinc-finger domain

Eukaryotic Cys2His2 zinc finger domain is one of the most common and important structural motifs involved in protein-DNA interaction. The recognition motif is characterized by the tetrahedral coordination of a zinc ion by conserved cysteine and histidine residues. We have characterized the prokaryotic Cys2His2 zinc finger motif, included in the DNA binding region (Ros87) of Ros protein from Agrobacterium tumefaciens, demonstrating that, although possessing a similar zinc coordination sphere, this domain presents significant differences from its eukaryotic counterpart. Furthermore, basic residues flanking the zinc binding region on either side have been demonstrated, by Electrophoretic Mobility Shift Assay (EMSA) experiments, to be essential for Ros DNA binding. In spite of this wealth of knowledge, the structural details of the mechanism through which the prokaryotic zinc fingers recognize their target genes are still unclear. Here, to gain insights into the molecular DNA recognition process of prokaryotic zinc finger domains we applied a strategy in which we performed molecular docking studies using a combination of Nuclear Magnetic Resonance (NMR) and Molecular Dynamics (MD) simulations data. The results demonstrate that the MD ensemble provides a reasonable picture of Ros87 backbone dynamics in solution. The Ros87-DNA model indicates that the interaction involves the first two residue of the first α-helix, and several residues located in the basic regions flanking the zinc finger domain. Interestingly, the prokaryotic zinc finger domain, mainly with the C-terminal tail that is wrapped around the DNA, binds a more extended recognition site than the eukaryotic counterpart. Our analysis demonstrates that the introduction of the protein flexibility in docking studies can improve, in terms of accuracy, the quality of the obtained models and could be particularly useful for protein showing high conformational heterogeneity as well as for computational drug design applications.

Exploring regions of conformational space occupied by two-domain proteins

The presence of heterogeneity in the interdomain arrangement of several biomolecules is required for their function. Here we present a method to obtain crucial clues to distinguish between different kinds of protein conformational distributions based on experimental NMR data. The method explores subregions of the conformational space and provides both upper and lower bounds of probability for the system to be in each subregion.

The structure of mouse cytomegalovirus m04 protein obtained from sparse NMR data reveals a conserved fold of the m02-m06 viral immune modulator family

Immunoevasins are key proteins used by viruses to subvert host immune responses. Determining their high-resolution structures is key to understanding virus-host interactions toward the design of vaccines and other antiviral therapies. Mouse cytomegalovirus encodes a unique set of immunoevasins, the m02-m06 family, that modulates major histocompatibility complex class I (MHC-I) antigen presentation to CD8+ T cells and natural killer cells. Notwithstanding the large number of genetic and functional studies, the structural biology of immunoevasins remains incompletely understood, largely because of crystallization bottlenecks. Here we implement a technology using sparse nuclear magnetic resonance data and integrative Rosetta modeling to determine the structure of the m04/gp34 immunoevasin extracellular domain. The structure reveals a β fold that is representative of the m02-m06 family of viral proteins, several of which are known to bind MHC-I molecules and interfere with antigen presentation, suggesting its role as a diversified immune regulation module.

Correlated inter-domain motions in adenylate kinase

Correlated inter-domain motions in proteins can mediate fundamental biochemical processes such as signal transduction and allostery. Here we characterize at structural level the inter-domain coupling in a multidomain enzyme, Adenylate Kinase (AK), using computational methods that exploit the shape information encoded in residual dipolar couplings (RDCs) measured under steric alignment by nuclear magnetic resonance (NMR). We find experimental evidence for a multi-state equilibrium distribution along the opening/closing pathway of Adenylate Kinase, previously proposed from computational work, in which inter-domain interactions disfavour states where only the AMP binding domain is closed. In summary, we provide a robust experimental technique for study of allosteric regulation in AK and other enzymes.

Structural and functional implications of the QUA2 domain on RNA recognition by GLD-1

The STAR family comprises ribonucleic acid (RNA)-binding proteins that play key roles in RNA-regulatory processes. RNA recognition is achieved by a KH domain with an additional α-helix (QUA2) that seems to extend the RNA-binding surface to six nucleotides for SF1 (Homo sapiens) and seven nucleotides for GLD-1 (Caenorhabditis elegans). To understand the structural basis of this probable difference in specificity, we determined the solution structure of GLD-1 KH-QUA2 with the complete consensus sequence identified in the tra-2 gene. Compared to SF1, the GLD-1 KH-QUA2 interface adopts a different conformation resulting indeed in an additional sequence-specific binding pocket for a uracil at the 5'end. The functional relevance of this binding pocket is emphasized by our bioinformatics analysis showing that GLD-1 binding sites with this 5'end uracil are more predictive for the functional response of the messenger RNAs to gld-1 knockout. We further reveal the importance of the KH-QUA2 interface in vitro and that its alteration in vivo affects the level of translational repression dependent on the sequence of the GLD-1 binding motif. In conclusion, we demonstrate that the QUA2 domain distinguishes GLD-1 from other members of the STAR family and contributes more generally to the modulation of RNA-binding affinity and specificity of KH domain containing proteins.

The role of ubiquitin-binding domains in human pathophysiology

Ubiquitination, a fundamental post-translational modification (PTM) resulting in the covalent attachment of ubiquitin (Ub) to a target protein, is currently implicated in several key cellular processes. Although ubiquitination was initially associated with protein degradation, it is becoming increasingly evident that proteins labeled with polyUb chains of specific topology and length are activated in an ever-expanding repertoire of specific cellular processes. In addition to their involvement in the classical protein degradation pathways they are involved in DNA repair, kinase regulation and nuclear factor-κB (NF-κB) signaling. The sorting and processing of distinct Ub signals is mediated by small protein motifs, known as Ub-binding domains (UBDs), which are found in proteins that execute disparate biological functions. The involvement of UBDs in several biological pathways has been revealed by several studies which have highlighted the vital role of UBDs in cellular homeostasis. Importantly, functional impairment of UBDs in key regulatory pathways has been related to the development of pathophysiological conditions, including immune disorders and cancer. In this review, we present an up-to-date account of the crucial role of UBDs and their functions, with a special emphasis on their functional impairment in key biological pathways and the pathogenesis of several human diseases. The still under-investigated topic of Ub-UBD interactions as a target for developing novel therapeutic strategies against many diseases is also discussed.

A tensor-free method for the structural and dynamical refinement of proteins using residual dipolar couplings

Residual dipolar couplings (RDCs) are parameters measured in nuclear magnetic resonance spectroscopy that can provide exquisitely detailed information about the structure and dynamics of biological macromolecules. We describe here a method of using RDCs for the structural and dynamical refinement of proteins that is based on the observation that the RDC between two atomic nuclei depends directly on the angle ϑ between the internuclear vector and the external magnetic field. For every pair of nuclei for which an RDC is available experimentally, we introduce a structural restraint to minimize the deviation from the value of the angle ϑ derived from the measured RDC and that calculated in the refinement protocol. As each restraint involves only the calculation of the angle ϑ of the corresponding internuclear vector, the method does not require the definition of an overall alignment tensor to describe the preferred orientation of the protein with respect to the alignment medium. Application to the case of ubiquitin demonstrates that this method enables an accurate refinement of the structure and dynamics of this protein to be obtained.

Structural dynamics of a single-stranded RNA-helix junction using NMR

Many regulatory RNAs contain long single strands (ssRNA) that adjoin secondary structural elements. Here, we use NMR spectroscopy to study the dynamic properties of a 12-nucleotide (nt) ssRNA tail derived from the prequeuosine riboswitch linked to the 3' end of a 48-nt hairpin. Analysis of chemical shifts, NOE connectivity, (13)C spin relaxation, and residual dipolar coupling data suggests that the first two residues (A25 and U26) in the ssRNA tail stack onto the adjacent helix and assume an ordered conformation. The following U26-A27 step marks the beginning of an A6-tract and forms an acute pivot point for substantial motions within the tail, which increase toward the terminal end. Despite substantial internal motions, the ssRNA tail adopts, on average, an A-form helical conformation that is coaxial with the helix. Our results reveal a surprising degree of structural and dynamic complexity at the ssRNA-helix junction, which involves a fine balance between order and disorder that may facilitate efficient pseudoknot formation on ligand recognition.

A practical implicit solvent potential for NMR structure calculation

The benefits of protein structure refinement in water are well documented. However, performing structure refinement with explicit atomic representation of the solvent molecules is computationally expensive and impractical for NMR-restrained structure calculations that start from completely extended polypeptide templates. Here we describe a new implicit solvation potential, EEFx (Effective Energy Function for XPLOR-NIH), for NMR-restrained structure calculations of proteins in XPLOR-NIH. The key components of EEFx are an energy term for solvation energy that works together with other nonbonded energy functions, and a dedicated force field for conformational and nonbonded protein interaction parameters. The initial results obtained with EEFx show that significant improvements in structural quality can be obtained. EEFx is computationally efficient and can be used both to fold and refine structures. Overall, EEFx improves the quality of protein conformation and nonbonded atomic interactions. Moreover, such benefits are accompanied by enhanced structural precision and enhanced structural accuracy, reflected in improved agreement with the cross-validated dipolar coupling data. Finally, implementation of EEFx calculations is straightforward and computationally efficient. Overall, EEFx provides a useful method for the practical calculation of experimental protein structures in a physically realistic environment.

Dissecting the effects of concentrated carbohydrate solutions on protein diffusion, hydration, and internal dynamics

We present herein a thorough description of the effects of high glucose concentrations on the diffusion, hydration and internal dynamics of ubiquitin, as predicted from extensive molecular dynamics simulations on several systems described at fully atomistic level. We observe that the protein acts as a seed that speeds up the natural propensity of glucose to cluster at high concentration; the sugar molecules thus aggregate around the protein trapping it inside a dynamic cage. This process extensively dehydrates the protein surface, restricts the motions of the remaining water molecules, and drags the large-scale, collective motions of protein atoms slowing down the rate of exploration of the conformational space despite only a slight dampening of fast, local dynamics. We discuss how these effects could be relevant to the function of sugars as preservation agents in biological materials, and how crowding by small sticky molecules could modulate proteins across different reaction coordinates inside the cellular cytosol.

Simultaneous use of solution NMR and X-ray data in REFMAC5 for joint refinement/detection of structural differences

The program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic data and paramagnetic NMR data (pseudocontact shifts and self-orientation residual dipolar couplings) and/or diamagnetic residual dipolar couplings. Incorporation of these long-range NMR restraints in REFMAC5 can reveal differences between solid-state and solution conformations of molecules or, in their absence, can be used together with X-ray crystallographic data for structural refinement. Since NMR and X-ray data are complementary, when a single structure is consistent with both sets of data and still maintains reasonably `ideal' geometries, the reliability of the derived atomic model is expected to increase. The program was tested on five different proteins: the catalytic domain of matrix metalloproteinase 1, GB3, ubiquitin, free calmodulin and calmodulin complexed with a peptide. In some cases the joint refinement produced a single model consistent with both sets of observations, while in other cases it indicated, outside the experimental uncertainty, the presence of different protein conformations in solution and in the solid state.

High-resolution NMR spectroscopy of encapsulated proteins dissolved in low-viscosity fluids

High-resolution multi-dimensional solution NMR is unique as a biophysical and biochemical tool in its ability to examine both the structure and dynamics of macromolecules at atomic resolution. Conventional solution NMR approaches, however, are largely limited to examinations of relatively small (<25kDa) molecules, mostly due to the spectroscopic consequences of slow rotational diffusion. Encapsulation of macromolecules within the protective nanoscale aqueous interior of reverse micelles dissolved in low viscosity fluids has been developed as a means through which the 'slow tumbling problem' can be overcome. This approach has been successfully applied to diverse proteins and nucleic acids ranging up to 100kDa, considerably widening the range of biological macromolecules to which conventional solution NMR methodologies may be applied. Recent advances in methodology have significantly broadened the utility of this approach in structural biology and molecular biophysics.

Homonuclear decoupling for enhancing resolution and sensitivity in NOE and RDC measurements of peptides and proteins

Application of band-selective homonuclear (BASH) (1)H decoupling pulses during acquisition of the (1)H free induction decay is shown to be an efficient procedure for removal of scalar and residual dipolar couplings between amide and aliphatic protons. BASH decoupling can be applied in both dimensions of a homonuclear 2D NMR experiment and is particularly useful for enhancing spectral resolution in the H(N)-H(α) region of NOESY spectra of peptides and proteins, which contain important information on the backbone torsion angles. The method then also prevents generation of zero quantum and Hz(N)-Hz(α) terms, thereby facilitating analysis of intraresidue interactions. Application to the NOESY spectrum of a hexapeptide fragment of the intrinsically disordered protein α-synuclein highlights the considerable diffusion anisotropy present in linear peptides. Removal of residual dipolar couplings between H(N) and aliphatic protons in weakly aligned proteins increases resolution in the (1)H-(15)N HSQC region of the spectrum and allows measurement of RDCs in samples that are relatively strongly aligned. The approach is demonstrated for measurement of RDCs in protonated (15)N/(13)C-enriched ubiquitin, aligned in Pf1, yielding improved fitting to the ubiquitin structure.

Solution structure of lysine-free (K0) ubiquitin

Lysine-free ubiquitin (K0-Ub) is commonly used to study the ubiquitin-signaling pathway, where it is assumed to have the same structure and function as wild-type ubiquitin (wt-Ub). However, the K0-Ub (15) N heteronuclear single quantum correlation NMR spectrum differs significantly from wt-Ub and the melting temperature is depressed by 19°C, raising the question of the structural integrity and equivalence to wt-Ub. The three-dimensional structure of K0-Ub was determined by solution NMR, using chemical shift and residual dipolar coupling data. K0-Ub adopts the same backbone structure as wt-Ub, and all significant chemical shifts can be related to interactions impacted by the K to R mutations.

Improved cross validation of a static ubiquitin structure derived from high precision residual dipolar couplings measured in a drug-based liquid crystalline phase

The antibiotic squalamine forms a lyotropic liquid crystal at very low concentrations in water (0.3-3.5% w/v), which remains stable over a wide range of temperature (1-40 °C) and pH (4-8). Squalamine is positively charged, and comparison of the alignment of ubiquitin relative to 36 previously reported alignment conditions shows that it differs substantially from most of these, but is closest to liquid crystalline cetyl pyridinium bromide. High precision residual dipolar couplings (RDCs) measured for the backbone (1)H-(15)N, (15)N-(13)C', (1)H(α)-(13)C(α), and (13)C'-(13)C(α) one-bond interactions in the squalamine medium fit well to the static structural model previously derived from NMR data. Inclusion into the structure refinement procedure of these RDCs, together with (1)H-(15)N and (1)H(α)-(13)C(α) RDCs newly measured in Pf1, results in improved agreement between alignment-induced changes in (13)C' chemical shift, (3)JHNHα values, and (13)C(α)-(13)C(β) RDCs and corresponding values predicted by the structure, thereby validating the high quality of the single-conformer structural model. This result indicates that fitting of a single model to experimental data provides a better description of the average conformation than does averaging over previously reported NMR-derived ensemble representations. The latter can capture dynamic aspects of a protein, thus making the two representations valuable complements to one another.