Structural basis of hAT transposon end recognition by Hermes, an octameric DNA transposase from Musca domestica.

Hermes is a member of the hAT transposon superfamily that has active representatives, including McClintock's archetypal Ac mobile genetic element, in many eukaryotic species. The crystal structure of the Hermes transposase-DNA complex reveals that Hermes forms an octameric ring organized as a tetramer of dimers. Although isolated dimers are active in vitro for all the chemical steps of transposition, only octamers are active in vivo. The octamer can provide not only multiple specific DNA-binding domains to recognize repeated subterminal sequences within the transposon ends, which are important for activity, but also multiple nonspecific DNA binding surfaces for target capture. The unusual assembly explains the basis of bipartite DNA recognition at hAT transposon ends, provides a rationale for transposon end asymmetry, and suggests how the avidity provided by multiple sites of interaction could allow a transposase to locate its transposon ends amidst a sea of chromosomal DNA.

Conformational heterogeneity of UCAAUC RNA oligonucleotide from molecular dynamics simulations, SAXS, and NMR experiments.

We describe the conformational ensemble of the single-stranded r(UCAAUC) oligonucleotide obtained using extensive molecular dynamics (MD) simulations and Rosetta's FARFAR2 algorithm. The conformations observed in MD consist of A-form-like structures and variations thereof. These structures are not present in the pool generated using FARFAR2. By comparing with available nuclear magnetic resonance (NMR) measurements, we show that the presence of both A-form-like and other extended conformations is necessary to quantitatively explain experimental data. To further validate our results, we measure solution X-ray scattering (SAXS) data on the RNA hexamer and find that simulations result in more compact structures than observed from these experiments. The integration of simulations with NMR via a maximum entropy approach shows that small modifications to the MD ensemble lead to an improved description of the conformational ensemble. Nevertheless, we identify persisting discrepancies in matching experimental SAXS data.

Extended q-Range X-Ray Scattering Reveals High-Resolution Structural Details of Biomacromolecules in Aqueous Solutions.

We communicate a feasibility study for high-resolution structural characterization of biomacromolecules in aqueous solution from X-ray scattering experiments measured over a range of scattering vectors (q) that is approximately two orders of magnitude wider than used previously for such systems. Scattering data with such an extended q-range enables the recovery of the underlying real-space atomic pair distribution function, which facilitates structure determination. We demonstrate the potential of this method for biomacromolecules using several types of cyclodextrins (CD) as model systems. We successfully identified deviations of the tilting angles for the glycosidic units in CDs in aqueous solutions relative to their values in the crystalline forms of these molecules. Such level of structural detail is inaccessible from standard small angle scattering measurements. Our results call for further exploration of ultra-wide-angle X-ray scattering measurements for biomacromolecules.

Structure of the cell-binding component of the Clostridium difficile binary toxin reveals a di-heptamer macromolecular assembly.

Targeting Clostridium difficile infection is challenging because treatment options are limited, and high recurrence rates are common. One reason for this is that hypervirulent C. difficile strains often have a binary toxin termed the C. difficile toxin, in addition to the enterotoxins TsdA and TsdB. The C. difficile toxin has an enzymatic component, termed CDTa, and a pore-forming or delivery subunit termed CDTb. CDTb was characterized here using a combination of single-particle cryoelectron microscopy, X-ray crystallography, NMR, and other biophysical methods. In the absence of CDTa, 2 di-heptamer structures for activated CDTb (1.0 MDa) were solved at atomic resolution, including a symmetric (SymCDTb; 3.14 Å) and an asymmetric form (AsymCDTb; 2.84 Å). Roles played by 2 receptor-binding domains of activated CDTb were of particular interest since the receptor-binding domain 1 lacks sequence homology to any other known toxin, and the receptor-binding domain 2 is completely absent in other well-studied heptameric toxins (i.e., anthrax). For AsymCDTb, a Ca2+ binding site was discovered in the first receptor-binding domain that is important for its stability, and the second receptor-binding domain was found to be critical for host cell toxicity and the di-heptamer fold for both forms of activated CDTb. Together, these studies represent a starting point for developing structure-based drug-design strategies to target the most severe strains of C. difficile.

Structural and biophysical properties of farnesylated KRas interacting with the chaperone SmgGDS-558.

KRas is a small GTPase and membrane-bound signaling protein. Newly synthesized KRas is post-translationally modified with a membrane-anchoring prenyl group. KRas chaperones are therapeutic targets in cancer due to their participation in trafficking oncogenic KRas to membranes. SmgGDS splice variants are chaperones for small GTPases with basic residues in their hypervariable domain (HVR), including KRas. SmgGDS-607 escorts pre-prenylated small GTPases, while SmgGDS-558 escorts prenylated small GTPases. We provide a structural description of farnesylated and fully processed KRas (KRas-FMe) in complex with SmgGDS-558 and define biophysical properties of this interaction. Surface plasmon resonance measurements on biomimetic model membranes quantified the thermodynamics of the interaction of SmgGDS with KRas, and small-angle x-ray scattering was used to characterize complexes of SmgGDS-558 and KRas-FMe structurally. Structural models were refined using Monte Carlo and molecular dynamics simulations. Our results indicate that SmgGDS-558 interacts with the HVR and the farnesylated C-terminus of KRas-FMe, but not its G-domain. Therefore, SmgGDS-558 interacts differently with prenylated KRas than prenylated RhoA, whose G-domain was found in close contact with SmgGDS-558 in a recent crystal structure. Using immunoprecipitation assays, we show that SmgGDS-558 binds the GTP-bound, GDP-bound, and nucleotide-free forms of farnesylated and fully processed KRas in cells, consistent with SmgGDS-558 not engaging the G-domain of KRas. We found that the dissociation constant, Kd, for KRas-FMe binding to SmgGDS-558 is comparable with that for the KRas complex with PDEδ, a well-characterized KRas chaperone that also does not interact with the KRas G-domain. These results suggest that KRas interacts in similar ways with the two chaperones SmgGDS-558 and PDEδ. Therapeutic targeting of the SmgGDS-558/KRas complex might prove as useful as targeting the PDEδ/KRas complex in KRas-driven cancers.

Chemical shifts-based similarity restraints improve accuracy of RNA structures determined via NMR.

Determination of structure of RNA via NMR is complicated in large part by the lack of a precise parameterization linking the observed chemical shifts to the underlying geometric parameters. In contrast to proteins, where numerous high-resolution crystal structures serve as coordinate templates for this mapping, such models are rarely available for smaller oligonucleotides accessible via NMR, or they exhibit crystal packing and counter-ion binding artifacts that prevent their use for the chemical shifts analysis. On the other hand, NMR-determined structures of RNA often are not solved at the density of restraints required to precisely define the variable degrees of freedom. In this study we sidestep the problems of direct parameterization of the RNA chemical shifts/structure relationship and examine the effects of imposing local fragmental coordinate similarity restraints based on similarities of the experimental secondary ribose 13C/1H chemical shifts instead. The effect of such chemical shift similarity (CSS) restraints on the structural accuracy is assessed via residual dipolar coupling (RDC)-based cross-validation. Improvements in the coordinate accuracy are observed for all of the six RNA constructs considered here as test cases, which argues for routine inclusion of these terms during NMR-based oligonucleotide structure determination. Such accuracy improvements are expected to facilitate derivation of the chemical shift/structure relationships for RNA.

HIV-1 gp120-CD4-Induced Antibody Complex Elicits CD4 Binding Site-Specific Antibody Response in Mice.

Elicitation of broadly neutralizing Ab (bNAb) responses toward the conserved HIV-1 envelope (Env) CD4 binding site (CD4bs) by vaccination is an important goal for vaccine development and yet to be achieved. The outcome of previous immunogenicity studies suggests that the limited accessibility of the CD4bs and the presence of predominant nonneutralizing determinants (nND) on Env may impede the elicitation of bNAbs and their precursors by vaccination. In this study, we designed a panel of novel immunogens that 1) preferentially expose the CD4bs by selective elimination of glycosylation sites flanking the CD4bs, and 2) minimize the nND immune response by engineering fusion proteins consisting of gp120 Core and one or two CD4-induced (CD4i) mAbs for masking nND epitopes, referred to as gp120-CD4i fusion proteins. As expected, the fusion proteins possess improved antigenicity with retained affinity for VRC01-class, CD4bs-directed bNAbs and dampened affinity for nonneutralizing Abs. We immunized C57BL/6 mice with these fusion proteins and found that overall the fusion proteins elicit more focused CD4bs Ab response than prototypical gp120 Core by serological analysis. Consistently, we found that mice immunized with selected gp120-CD4i fusion proteins have higher frequencies of germinal center-activated B cells and CD4bs-directed memory B cells than those inoculated with parental immunogens. We isolated three mAbs from mice immunized with selected gp120-CD4i fusion proteins and found that their footprints on Env are similar to VRC01-class bNAbs. Thus, using gp120-CD4i fusion proteins with selective glycan deletion as immunogens could focus Ab response toward CD4bs epitope.

Structural Characterization and Modeling of a Respiratory Syncytial Virus Fusion Glycoprotein Nanoparticle Vaccine in Solution.

The respiratory syncytial virus (RSV) fusion (F) protein/polysorbate 80 (PS80) nanoparticle vaccine is the most clinically advanced vaccine for maternal immunization and protection of newborns against RSV infection. It is composed of a near-full-length RSV F glycoprotein, with an intact membrane domain, formulated into a stable nanoparticle with PS80 detergent. To understand the structural basis for the efficacy of the vaccine, a comprehensive study of its structure and hydrodynamic properties in solution was performed. Small-angle neutron scattering experiments indicate that the nanoparticle contains an average of 350 PS80 molecules, which form a cylindrical micellar core structure and five RSV F trimers that are arranged around the long axis of the PS80 core. All-atom models of full-length RSV F trimers were built from crystal structures of the soluble ectodomain and arranged around the long axis of the PS80 core, allowing for the generation of an ensemble of conformations that agree with small-angle neutron and X-ray scattering data as well as transmission electron microscopy (TEM) images. Furthermore, the hydrodynamic size of the RSV F nanoparticle was found to be modulated by the molar ratio of PS80 to protein, suggesting a mechanism for nanoparticle assembly involving addition of RSV F trimers to and growth along the long axis of the PS80 core. This study provides structural details of antigen presentation and conformation in the RSV F nanoparticle vaccine, helping to explain the induction of broad immunity and observed clinical efficacy. Small-angle scattering methods provide a general strategy to visualize surface glycoproteins from other pathogens and to structurally characterize nanoparticle vaccines.

Phosphorylation-induced conformational dynamics in an intrinsically disordered protein and potential role in phenotypic heterogeneity.

Intrinsically disordered proteins (IDPs) that lack a unique 3D structure and comprise a large fraction of the human proteome play important roles in numerous cellular functions. Prostate-Associated Gene 4 (PAGE4) is an IDP that acts as a potentiator of the Activator Protein-1 (AP-1) transcription factor. Homeodomain-Interacting Protein Kinase 1 (HIPK1) phosphorylates PAGE4 at S9 and T51, but only T51 is critical for its activity. Here, we identify a second kinase, CDC-Like Kinase 2 (CLK2), which acts on PAGE4 and hyperphosphorylates it at multiple S/T residues, including S9 and T51. We demonstrate that HIPK1 is expressed in both androgen-dependent and androgen-independent prostate cancer (PCa) cells, whereas CLK2 and PAGE4 are expressed only in androgen-dependent cells. Cell-based studies indicate that PAGE4 interaction with the two kinases leads to opposing functions. HIPK1-phosphorylated PAGE4 (HIPK1-PAGE4) potentiates c-Jun, whereas CLK2-phosphorylated PAGE4 (CLK2-PAGE4) attenuates c-Jun activity. Consistent with the cellular data, biophysical measurements (small-angle X-ray scattering, single-molecule fluorescence resonance energy transfer, and NMR) indicate that HIPK1-PAGE4 exhibits a relatively compact conformational ensemble that binds AP-1, whereas CLK2-PAGE4 is more expanded and resembles a random coil with diminished affinity for AP-1. Taken together, the results suggest that the phosphorylation-induced conformational dynamics of PAGE4 may play a role in modulating changes between PCa cell phenotypes. A mathematical model based on our experimental data demonstrates how differential phosphorylation of PAGE4 can lead to transitions between androgen-dependent and androgen-independent phenotypes by altering the AP-1/androgen receptor regulatory circuit in PCa cells.

Accuracy of MD solvent models in RNA structure refinement assessed via liquid-crystal NMR and spin relaxation data.

Molecular dynamics (MD) simulations play an important role in characterizing Ribonucleic Acid (RNA) structure, augmenting information from experimental techniques such as Nuclear Magnetic Resonance (NMR). In this work, we examine the accuracy of structural representation resulting from application of a number of explicit and implicit solvent models and refinement protocols against experimental data ranging from high density of residual dipolar coupling (RDC) restraints to completely unrestrained simulations. For a prototype A-form RNA helix, our results indicate that AMBER RNA force field with either implicit or explicit solvent can produce a realistic dynamic representation of RNA helical structure, accurately cross-validating with respect to a diverse array of NMR observables. In refinement against NMR distance restraints, modern MD force fields are found to be equally adequate, with high fidelity cross-validation to the residual dipolar couplings (RDCs) and residual chemical shift anisotropies (RCSAs), while slightly over-estimating structural order as monitored via NMR relaxation data. With restraints trimmed to encode only for base pairing information, cross-validation quality significantly deteriorates, now exhibiting a pronounced dependence on the choice of the solvent model. This deterioration is found to be partially reversible by increasing planarity restraints on the nucleobase geometry. For completely unrestrained MD simulations, the choice of water model becomes very important, with the best-performing TIP4P-Ew accurately reproducing both the RDC and RCSA data, while closely matching the NMR-derived order parameters. The information provided here will serve as a foundation for MD-based refinement of solution state NMR structures of RNA.

A trapped human PPM1A-phosphopeptide complex reveals structural features critical for regulation of PPM protein phosphatase activity.

Metal-dependent protein phosphatases (PPM) are evolutionarily unrelated to other serine/threonine protein phosphatases and are characterized by their requirement for supplementation with millimolar concentrations of Mg2+ or Mn2+ ions for activity in vitro The crystal structure of human PPM1A (also known as PP2Cα), the first PPM structure determined, displays two tightly bound Mn2+ ions in the active site and a small subdomain, termed the Flap, located adjacent to the active site. Some recent crystal structures of bacterial or plant PPM phosphatases have disclosed two tightly bound metal ions and an additional third metal ion in the active site. Here, the crystal structure of the catalytic domain of human PPM1A, PPM1Acat, complexed with a cyclic phosphopeptide, c(MpSIpYVA), a cyclized variant of the activation loop of p38 MAPK (a physiological substrate of PPM1A), revealed three metal ions in the active site. The PPM1Acat D146E-c(MpSIpYVA) complex confirmed the presence of the anticipated third metal ion in the active site of metazoan PPM phosphatases. Biophysical and computational methods suggested that complex formation results in a slightly more compact solution conformation through reduced conformational flexibility of the Flap subdomain. We also observed that the position of the substrate in the active site allows solvent access to the labile third metal-binding site. Enzyme kinetics of PPM1Acat toward a phosphopeptide substrate supported a random-order, bi-substrate mechanism, with substantial interaction between the bound substrate and the labile metal ion. This work illuminates the structural and thermodynamic basis of an innate mechanism regulating the activity of PPM phosphatases.

Maximizing accuracy of RNA structure in refinement against residual dipolar couplings.

Structural information about ribonucleic acid (RNA) is lagging behind that of proteins, in part due to its high charge and conformational variability. Molecular dynamics (MD) has played an important role in describing RNA structure, complementing information from both nuclear magnetic resonance (NMR), or X-ray crystallography. We examine the impact of the choice of the empirical force field for RNA structure refinement using cross-validation against residual dipolar couplings (RDCs) as structural accuracy reporter. Four force fields, representing both the state-of-the art in RNA simulation and the most popular selections in NMR structure determination, are compared for a prototypical A-RNA helix. RNA structural accuracy is also evaluated as a function of both density and nature of input NMR data including RDCs, anisotropic chemical shifts, and distance restraints. Our results show a complex interplay between the experimental restraints and the force fields indicating two best-performing choices: high-fidelity refinement in explicit solvent, and the conformational database-derived potentials. Accuracy of RNA models closely tracks the density of 1-bond C-H RDCs, with other data types having beneficial, but smaller effects. At lower RDC density, or when refining against NOEs only, the two selected force fields are capable of accurately describing RNA helices with little or no experimental RDC data, making them available for the higher order structure assembly or better quantification of the intramolecular dynamics. Unrestrained simulations of simple RNA motifs with state-of-the art MD force fields appear to capture the flexibility inherent in nucleic acids while also maintaining a good agreement with the experimental observables.

Dynamic equilibrium between closed and partially closed states of the bacterial Enzyme I unveiled by solution NMR and X-ray scattering.

Enzyme I (EI) is the first component in the bacterial phosphotransferase system, a signal transduction pathway in which phosphoryl transfer through a series of bimolecular protein-protein interactions is coupled to sugar transport across the membrane. EI is a multidomain, 128-kDa homodimer that has been shown to exist in two conformational states related to one another by two large (50-90°) rigid body domain reorientations. The open conformation of apo EI allows phosphoryl transfer from His189 located in the N-terminal domain α/ß (EIN(α/ß)) subdomain to the downstream protein partner bound to the EIN(α) subdomain. The closed conformation, observed in a trapped phosphoryl transfer intermediate, brings the EIN(α/ß) subdomain into close proximity to the C-terminal dimerization domain (EIC), thereby permitting in-line phosphoryl transfer from phosphoenolpyruvate (PEP) bound to EIC to His189. Here, we investigate the solution conformation of a complex of an active site mutant of EI (H189A) with PEP. Simulated annealing refinement driven simultaneously by solution small angle X-ray scattering and NMR residual dipolar coupling data demonstrates unambiguously that the EI(H189A)-PEP complex exists in a dynamic equilibrium between two approximately equally populated conformational states, one corresponding to the closed structure and the other to a partially closed species. The latter likely represents an intermediate in the open-to-closed transition.

Hybrid Applications of Solution Scattering to Aid Structural Biology.

Biomolecular applications of solution X-ray and neutron scattering (SAXS and SANS, respectively) started in late 1960s - early 1970s but were relatively limited in their ability to provide a detailed structural picture and lagged behind what became the two primary methods of experimental structural biology - X-ray crystallography and NMR. However, improvements in both data analysis and instrumentation led to an explosive growth in the number of studies that used small-angle scattering (SAS) for investigation of macromolecular structure, often in combination with other biophysical techniques. Such hybrid applications are nowadays quickly becoming a norm whenever scattering data are used for two reasons. First, it is generally accepted that SAS data on their own cannot lead to a uniquely defined high-resolution structural model, creating a need for supplementing them with information from complementary techniques. Second, solution scattering data are frequently applied in situations when a method such NMR or X-ray crystallography cannot provide a satisfactory structural picture, which makes these additional restraints highly desirable. Maturation of the hybrid bio-SAS approaches brings to light new questions including completeness of the conformational space sampling, model validation, and data compatibility.

Dissociation of the trimeric gp41 ectodomain at the lipid-water interface suggests an active role in HIV-1 Env-mediated membrane fusion.

The envelope glycoprotein gp41 mediates the process of membrane fusion that enables entry of the HIV-1 virus into the host cell. The actual fusion process involves a switch from a homotrimeric prehairpin intermediate conformation, consisting of parallel coiled-coil helices, to a postfusion state where the ectodomains are arranged as a trimer of helical hairpins, adopting a six-helix bundle (6HB) state. Here, we show by solution NMR spectroscopy that a water-soluble 6HB gp41 ectodomain binds to zwitterionic detergents that contain phosphocholine or phosphatidylcholine head groups and phospholipid vesicles that mimic T-cell membrane composition. Binding results in the dissociation of the 6HB and the formation of a monomeric state, where its two α-helices, N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR), become embedded in the lipid-water interface of the virus and host cell. The atomic structure of the gp41 ectodomain monomer, based on NOE distance restraints and residual dipolar couplings, shows that the NHR and CHR helices remain mostly intact, but they completely lose interhelical contacts. The high affinity of the ectodomain helices for phospholipid surfaces suggests that unzippering of the prehairpin intermediate leads to a state where the NHR and CHR helices become embedded in the host cell and viral membranes, respectively, thereby providing a physical force for bringing these membranes into close juxtaposition before actual fusion.

Probing the Action of Chemical Denaturant on an Intrinsically Disordered Protein by Simulation and Experiment.

Chemical denaturants are the most commonly used agents for unfolding proteins and are thought to act by better solvating the unfolded state. Improved solvation is expected to lead to an expansion of unfolded chains with increasing denaturant concentration, providing a sensitive probe of the denaturant action. However, experiments have so far yielded qualitatively different results concerning the effects of chemical denaturation. Studies using Förster resonance energy transfer (FRET) and other methods found an increase in radius of gyration with denaturant concentration, but most small-angle X-ray scattering (SAXS) studies found no change. This discrepancy therefore challenges our understanding of denaturation mechanism and more generally the accuracy of these experiments as applied to unfolded or disordered proteins. Here, we use all-atom molecular simulations to investigate the effect of urea and guanidinium chloride on the structure of the intrinsically disordered protein ACTR, which can be studied by experiment over a wide range of denaturant concentration. Using unbiased molecular simulations with a carefully calibrated denaturant model, we find that the protein chain indeed swells with increasing denaturant concentration. This is due to the favorable association of urea or guanidinium chloride with the backbone of all residues and with the side-chains of almost all residues, with denaturant-water transfer free energies inferred from this association in reasonable accord with experimental estimates. Interactions of the denaturants with the backbone are dominated by hydrogen bonding, while interactions with side-chains include other contributions. By computing FRET efficiencies and SAXS intensities at each denaturant concentration, we show that the simulation trajectories are in accord with both experiments on this protein, demonstrating that there is no fundamental inconsistency between the two types of experiment. Agreement with experiment also supports the picture of chemical denaturation described in our simulations, driven by weak association of denaturant with the protein. Our simulations support some assumptions needed for each experiment to accurately reflect changes in protein size, namely, that the commonly used FRET chromophores do not qualitatively alter the results and that possible effects such as preferential solvent partitioning into the interior of the chain do not interfere with the determination of radius of gyration from the SAXS experiments.

Consistent View of Polypeptide Chain Expansion in Chemical Denaturants from Multiple Experimental Methods.

There has been a long-standing controversy regarding the effect of chemical denaturants on the dimensions of unfolded and intrinsically disordered proteins: A wide range of experimental techniques suggest that polypeptide chains expand with increasing denaturant concentration, but several studies using small-angle X-ray scattering (SAXS) have reported no such increase of the radius of gyration (Rg). This inconsistency challenges our current understanding of the mechanism of chemical denaturants, which are widely employed to investigate protein folding and stability. Here, we use a combination of single-molecule Förster resonance energy transfer (FRET), SAXS, dynamic light scattering (DLS), and two-focus fluorescence correlation spectroscopy (2f-FCS) to characterize the denaturant dependence of the unfolded state of the spectrin domain R17 and the intrinsically disordered protein ACTR in two different denaturants. Standard analysis of the primary data clearly indicates an expansion of the unfolded state with increasing denaturant concentration irrespective of the protein, denaturant, or experimental method used. This is the first case in which SAXS and FRET have yielded even qualitatively consistent results regarding expansion in denaturant when applied to the same proteins. To more directly illustrate this self-consistency, we used both SAXS and FRET data in a Bayesian procedure to refine structural ensembles representative of the observed unfolded state. This analysis demonstrates that both of these experimental probes are compatible with a common ensemble of protein configurations for each denaturant concentration. Furthermore, the resulting ensembles reproduce the trend of increasing hydrodynamic radius with denaturant concentration obtained by 2f-FCS and DLS. We were thus able to reconcile the results from all four experimental techniques quantitatively, to obtain a comprehensive structural picture of denaturant-induced unfolded state expansion, and to identify the most likely sources of earlier discrepancies.

Quantitative Characterization of Configurational Space Sampled by HIV-1 Nucleocapsid Using Solution NMR, X-ray Scattering and Protein Engineering.

Nucleic-acid-related events in the HIV-1 replication cycle are mediated by nucleocapsid, a small protein comprising two zinc knuckles connected by a short flexible linker and flanked by disordered termini. Combining experimental NMR residual dipolar couplings, solution X-ray scattering and protein engineering with ensemble simulated annealing, we obtain a quantitative description of the configurational space sampled by the two zinc knuckles, the linker and disordered termini in the absence of nucleic acids. We first compute the conformational ensemble (with an optimal size of three members) of an engineered nucleocapsid construct lacking the N- and C-termini that satisfies the experimental restraints, and then validate this ensemble, as well as characterize the disordered termini, using the experimental data from the full-length nucleocapsid construct. The experimental and computational strategy is generally applicable to multidomain proteins. Differential flexibility within the linker results in asymmetric motion of the zinc knuckles which may explain their functionally distinct roles despite high sequence identity. One of the configurations (populated at a level of ≈40 %) closely resembles that observed in various ligand-bound forms, providing evidence for conformational selection and a mechanistic link between protein dynamics and function.

Side Chain Conformational Distributions of a Small Protein Derived from Model-Free Analysis of a Large Set of Residual Dipolar Couplings.

Accurate quantitative measurement of structural dispersion in proteins remains a prime challenge to both X-ray crystallography and NMR spectroscopy. Here we use a model-free approach based on measurement of many residual dipolar couplings (RDCs) in differentially orienting aqueous liquid crystalline solutions to obtain the side chain χ1 distribution sampled by each residue in solution. Applied to the small well-ordered model protein GB3, our approach reveals that the RDC data are compatible with a single narrow distribution of side chain χ1 angles for only about 40% of the residues. For more than half of the residues, populations greater than 10% for a second rotamer are observed, and four residues require sampling of three rotameric states to fit the RDC data. In virtually all cases, sampled χ1 values are found to center closely around ideal g(-), g(+) and t rotameric angles, even though no rotamer restraint is used when deriving the sampled angles. The root-mean-square difference between experimental (3)JHαHß couplings and those predicted by the Haasnoot-parametrized, motion-adjusted Karplus equation reduces from 2.05 to 0.75 Hz when using the new rotamer analysis instead of the 1.1-Å X-ray structure as input for the dihedral angles. A comparison between observed and predicted (3)JHαHß values suggests that the root-mean-square amplitude of χ1 angle fluctuations within a given rotamer well is ca. 20°. The quantitatively defined side chain rotamer equilibria obtained from our study set new benchmarks for evaluating improved molecular dynamics force fields, and also will enable further development of quantitative relations between side chain chemical shift and structure.

Quantitative residue-specific protein backbone torsion angle dynamics from concerted measurement of 3J couplings.

Three-bond (3)J(C'C') and (3)J(HNHα) couplings in peptides and proteins are functions of the intervening backbone torsion angle ϕ. In well-ordered regions, (3)J(HNHα) is tightly correlated with (3)J(C'C'), but the presence of large ϕ angle fluctuations differentially affects the two types of couplings. Assuming the ϕ angles follow a Gaussian distribution, the width of this distribution can be extracted from (3)J(C'C') and (3)J(HNHα), as demonstrated for the folded proteins ubiquitin and GB3. In intrinsically disordered proteins, slow transverse relaxation permits measurement of (3)J(C'C') and (3)J(HNH) couplings at very high precision, and impact of factors other than the intervening torsion angle on (3)J will be minimal, making these couplings exceptionally valuable structural reporters. Analysis of α-synuclein yields rather homogeneous widths of 69 ± 6° for the ϕ angle distributions and (3)J(C'C') values that agree well with those of a recent maximum entropy analysis of chemical shifts, J couplings, and (1)H-(1)H NOEs. Data are consistent with a modest (≤30%) population of the polyproline II region.