Structure and Disorder in an Unfolded State under Nondenaturing Conditions from Ensemble Models Consistent with a Large Number of Experimental Restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include ¹H-¹⁵N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O₂-induced ¹³C paramagnetic shifts, hydrogen-exchange protection factors, and ¹⁵N R₂ data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu²⁺-Ni²⁺ binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative α-helical structure and both native-like and nonnative tertiary structure.     

Retinal dynamics during light activation of rhodopsin revealed by solid-state NMR spectroscopy

Rhodopsin is a canonical member of class A of the G protein-coupled receptors (GPCRs) that are implicated in many of the drug interventions in humans and are of great pharmaceutical interest. The molecular mechanism of rhodopsin activation remains unknown as atomistic structural information for the active metarhodopsin II state is currently lacking. Solid-state ²H NMR constitutes a powerful approach to study atomic-level dynamics of membrane proteins. In the present application, we describe how information is obtained about interactions of the retinal cofactor with rhodopsin that change with light activation of the photoreceptor. The retinal methyl groups play an important role in rhodopsin function by directing conformational changes upon transition into the active state. Site-specific ²H labels have been introduced into the methyl groups of retinal and solid-state ²H NMR methods applied to obtain order parameters and correlation times that quantify the mobility of the cofactor in the inactive dark state, as well as the cryotrapped metarhodopsin I and metarhodopsin II states. Analysis of the angular-dependent ²H NMR line shapes for selectively deuterated methyl groups of rhodopsin in aligned membranes enables determination of the average ligand conformation within the binding pocket. The relaxation data suggest that the β-ionone ring is not expelled from its hydrophobic pocket in the transition from the pre-activated metarhodopsin I to the active metarhodopsin II state. Rather, the major structural changes of the retinal cofactor occur already at the metarhodopsin I state in the activation process. The metarhodopsin I to metarhodopsin II transition involves mainly conformational changes of the protein within the membrane lipid bilayer rather than the ligand. The dynamics of the retinylidene methyl groups upon isomerization are explained by an activation mechanism involving cooperative rearrangements of extracellular loop E2 together with transmembrane helices H5 and H6. These activating movements are triggered by steric clashes of the isomerized all-trans retinal with the β4 strand of the E2 loop and the side chains of Glu¹²² and Trp²⁶⁵ within the binding pocket. The solid-state ²H NMR data are discussed with regard to the pathway of the energy flow in the receptor activation mechanism.     

Point mutations in the HIV-1 matrix protein turn off the myristyl switch

During the late phase of human immunodeficiency virus type-1 (HIV-1) replication, newly synthesized retroviral Gag proteins are targeted to lipid raft regions of specific cellular membranes, where they assemble and bud to form new virus particles. Gag binds preferentially to the plasma membrane (PM) of most hematopoietic cell types, a process mediated by interactions between the cellular PM marker phosphatidylinositol-(4,5)-bisphosphate (PI(4,5)P(2)) and Gag's N-terminally myristoylated matrix (MA) domain. We recently demonstrated that PI(4,5)P(2) binds to a conserved cleft on MA and promotes myristate exposure, suggesting a role as both a direct membrane anchor and myristyl switch trigger. Here we show that PI(4,5)P(2) is also capable of binding to MA proteins containing point mutations that inhibit membrane binding in vitro, and in vivo, including V7R, L8A and L8I. However, these mutants do not exhibit PI(4,5)P(2) or concentration-dependent myristate exposure. NMR studies of V7R and L8A MA reveal minor structural changes that appear to be responsible for stabilizing the myristate-sequestered (myr(s)) species and inhibiting exposure. Unexpectedly, the myristyl group of a revertant mutant with normal PM targeting properties (V7R,L21K) is also tightly sequestered and insensitive to PI(4,5)P(2) binding. This mutant binds PI(4,5)P(2) with twofold higher affinity compared with the native protein, suggesting a potential compensatory mechanism for membrane binding.     

Accurate measurement of <Superscript>15</Superscript>N–<Superscript>13</Superscript>C residual dipolar couplings in nucleic acids

New 3D HCN quantitative J (QJ) pulse schemes are presented for the precise and accurate measurement of one-bond ¹⁵N_1/9–¹³C₁, ¹⁵N_1/9–¹³C_6/8, and ¹⁵N_1/9–¹³C_2/4 residual dipolar couplings (RDCs) in weakly aligned nucleic acids. The methods employ ¹H–¹³C multiple quantum (MQ) coherence or TROSY-type pulse sequences for optimal resolution and sensitivity. RDCs are obtained from the intensity ratio of H₁–C₁–N_1/9 (MQ-HCN-QJ) or H_6/8–C_6/8–N_1/9 (TROSY-HCN-QJ) correlations in two interleaved 3D NMR spectra, with dephasing intervals of zero (reference spectrum) and 1/(2J_NC) (attenuated spectrum). The different types of ¹⁵N–¹³C couplings can be obtained by using either the 3D MQ-HCN-QJ or TROSY-HCN-QJ pulse scheme, with the appropriate setting of the duration of the constant-time ¹⁵N evolution period and the offset of two frequency-selective ¹³C pulses. The methods are demonstrated for a uniformly ¹³C, ¹⁵N-enriched 24-nucleotide stem-loop RNA sequence, helix-35, aligned in the magnetic field using phage Pf1. For measurements of RDCs systematic errors are found to be negligible, and experiments performed on a 1.5 mM helix-35 sample result in an estimated precision of ca. 0.07 Hz for ¹D_NC, indicating the utility of the measured RDCs in structure validation and refinement. Indeed, for a complete set of ¹⁵N_1/9–¹³C₁, ¹⁵N_1/9–¹³C_6/8, and ¹⁵N_1/9–¹³C_2/4 dipolar couplings obtained for the stem nucleotides, the measured RDCs are in excellent agreement with those predicted for an NMR structure of helix-35, refined using independently measured observables, including ¹³C–¹H, ¹³C–¹³C and ¹H–¹H dipolar couplings.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-005-0646-2.     

Structural Characterization of Protein-Protein Complexes by Integrating Computational Docking with Small-angle Scattering Data

X-ray crystallography and NMR can provide detailed structural information of protein-protein complexes, but technical problems make their application challenging in the high-throughput regime. Other methods such as small-angle X-ray scattering (SAXS) are more promising for large-scale application, but at the cost of lower resolution, which is a problem that can be solved by complementing SAXS data with theoretical simulations. Here, we propose a novel strategy that combines SAXS data and accurate protein-protein docking simulations. The approach has been benchmarked on a large pool of known structures with synthetic SAXS data, and on three experimental examples. The combined approach (pyDockSAXS) provided a significantly better success rate (43% for the top 10 predictions) than either of the two methods alone. Further analysis of the influence of different docking parameters made it possible to increase the success rates for specific cases, and to define guidelines for improving the data-driven protein-protein docking protocols.     

Solution NMR Studies of Chlorella Virus DNA Ligase-adenylate

DNA ligases are essential guardians of genome integrity by virtue of their ability to recognize and seal 3′-OH/5′-phosphate nicks in duplex DNA. The substrate binding and three chemical steps of the ligation pathway are coupled to global and local changes in ligase structure, involving both massive protein domain movements and subtle remodeling of atomic contacts in the active site. Here we applied solution NMR spectroscopy to study the conformational dynamics of the Chlorella virus DNA ligase (ChVLig), a minimized eukaryal ATP-dependent ligase consisting of nucleotidyltransferase, OB, and latch domains. Our analysis of backbone ¹⁵N spin relaxation and ¹⁵N,¹H residual dipolar couplings of the covalent ChVLig-AMP intermediate revealed conformational sampling on fast (picosecond to nanosecond) and slow timescales (microsecond to millisecond), indicative of interdomain and intradomain flexibility. We identified local and global changes in ChVLig-AMP structure and dynamics induced by phosphate. In particular, the chemical shift perturbations elicited by phosphate were clustered in the peptide motifs that comprise the active site. We hypothesize that phosphate anion mimics some of the conformational transitions that occur when ligase-adenylate interacts with the nick 5′-phosphate.     

Novel Protein-Protein Contacts Facilitate mRNA 3′-Processing Signal Recognition by Rna15 and Hrp1

Precise 3′-end processing of mRNA is essential for correct gene expression, yet in yeast, 3′-processing signals consist of multiple ambiguous sequence elements. Two neighboring elements upstream of the cleavage site are particularly important for the accuracy (positioning element) and efficiency (efficiency element) of 3′-processing and are recognized by the RNA-binding proteins Rna15 and Hrp1, respectively. In vivo, these interactions are strengthened by the scaffolding protein Rna14 that stabilizes their association. The NMR structure of the 34 -kDa ternary complex of the RNA recognition motif (RRM) domains of Hrp1 and Rna15 bound to this pair of RNA elements was determined by residual dipolar coupling and paramagnetic relaxation experiments. It reveals how each of the proteins binds to RNA and introduces a novel class of protein-protein contact in regions of previously unknown function. These interdomain contacts had previously been overlooked in other multi-RRM structures, although a careful analysis suggests that they may be frequently present. Mutations in the regions of these contacts disrupt 3′-end processing, suggesting that they may structurally organize the ribonucleoprotein complexes responsible for RNA processing.     

Solution Structure of Histone Chaperone ANP32B: Interaction with Core Histones H3-H4 through Its Acidic Concave Domain

Eukaryotic gene expression is regulated by histone deposition onto and eviction from nucleosomes, which are mediated by several chromatin-modulating factors. Among them, histone chaperones are key factors that facilitate nucleosome assembly. Acidic nuclear phosphoprotein 32B (ANP32B) belongs to the ANP32 family, which shares N-terminal leucine-rich repeats (LRRs) and a C-terminal variable anionic region. The C-terminal region functions as an inhibitor of histone acetylation, but the functional roles of the LRR domain in chromatin regulation have remained elusive. Here, we report that the LRR domain of ANP32B possesses histone chaperone activity and forms a curved structure with a parallel β-sheet on the concave side and mostly helical elements on the convex side. Our analyses revealed that the interaction of ANP32B with the core histones H3-H4 occurs on its concave side, and both the acidic and hydrophobic residues that compose the concave surface are critical for histone binding. These results provide a structural framework for understanding the functional mechanisms of acidic histone chaperones.     

Conformational Dynamics and Structural Plasticity Play Critical Roles in the Ubiquitin Recognition of a UIM Domain

Ubiquitin-interacting motifs (UIMs) are an important class of protein domains that interact with ubiquitin or ubiquitin-like proteins. These approximately 20-residue-long domains are found in a variety of ubiquitin receptor proteins and serve as recognition modules towards intracellular targets, which may be individual ubiquitin subunits or polyubiquitin chains attached to a variety of proteins. Previous structural studies of interactions between UIMs and ubiquitin have shown that UIMs adopt an extended structure of a single α-helix, containing a hydrophobic surface with a conserved sequence pattern that interacts with key hydrophobic residues on ubiquitin. In light of this large body of structural studies, details regarding the presence and the roles of structural dynamics and plasticity are surprisingly lacking. In order to better understand the structural basis of ubiquitin-UIM recognition, we have characterized changes in the structure and dynamics of ubiquitin upon binding of a UIM domain from the yeast Vps27 protein. The solution structure of a ubiquitin-UIM fusion protein designed to study these interactions is reported here and found to consist of a well-defined ubiquitin core and a bipartite UIM helix. Moreover, we have studied the plasticity of the docking interface, as well as global changes in ubiquitin due to UIM binding at the picoseconds-to-nanoseconds and microseconds-to-milliseconds protein motions by nuclear magnetic resonance relaxation. Changes in generalized-order parameters of amide groups show a distinct trend towards increased structural rigidity at the UIM-ubiquitin interface relative to values determined in unbound ubiquitin. Analysis of ¹⁵N Carr-Purcell-Meiboom-Gill relaxation dispersion measurements suggests the presence of two types of motions: one directly related to the UIM-binding interface and the other induced to distal parts of the protein. This study demonstrates a case where localized interactions among protein domains have global effects on protein motions at timescales ranging from picoseconds to milliseconds.     

Molecular crowding inhibits intramolecular breathing motions in proteins

In aqueous solution some proteins undergo large-scale movements of secondary structures, subunits or domains, referred to as protein “breathing”, that define a native-state ensemble of structures. These fluctuations are sensitive to the nature and concentration of solutes and other proteins and are thereby expected to be different in the crowded interior of a cell than in dilute solution. Here we use a combination of wide angle X-ray scattering (WAXS) and computational modeling to derive a quantitative measure of the spatial scale of conformational fluctuations in a protein solution. Concentration-dependent changes in the observed scattering intensities are consistent with a model of structural fluctuations in which secondary structures undergo rigid-body motions relative to one another. This motion increases with decreasing protein concentration or increasing temperature. Analysis of a set of five structurally and functionally diverse proteins reveals a diversity of kinetic behaviors. Proteins with multiple disulfide bonds exhibit little or no increase in breathing in dilute solutions. The spatial extent of structural fluctuations appears highly dependent on both protein structure and concentration and is universally suppressed at very high protein concentrations.