Long-range motional restrictions in a multidomain zinc-finger protein from anisotropic tumbling

Allosteric autoregulation of DNA binding via a DNA-mimicking protein domain: a biophysical study of ZNF410-DNA interaction using small angle X-ray scattering

ZNF410 is a highly-conserved transcription factor, remarkable in that it recognizes a 15-base pair DNA element but has just a single responsive target gene in mammalian erythroid cells. ZNF410 includes a tandem array of five zinc-fingers (ZFs), surrounded by uncharacterized N- and C-terminal regions. Unexpectedly, full-length ZNF410 has reduced DNA binding affinity, compared to that of the isolated DNA binding ZF array, both in vitro and in cells. AlphaFold predicts a partially-folded N-terminal subdomain that includes a 30-residue long helix, preceded by a hairpin loop rich in acidic (aspartate/glutamate) and serine/threonine residues. This hairpin loop is predicted by AlphaFold to lie against the DNA binding interface of the ZF array. In solution, ZNF410 is a monomer and binds to DNA with 1:1 stoichiometry. Surprisingly, the single best-fit model for the experimental small angle X-ray scattering profile, in the absence of DNA, is the original AlphaFold model with the N-terminal long-helix and the hairpin loop occupying the ZF DNA binding surface. For DNA binding, the hairpin loop presumably must be displaced. After combining biophysical, biochemical, bioinformatic and artificial intelligence-based AlphaFold analyses, we suggest that the hairpin loop mimics the structure and electrostatics of DNA, and provides an additional mechanism, supplementary to sequence specificity, of regulating ZNF410 DNA binding.

Observing reorientation dynamics with Time-Resolved fluorescence and molecular dynamics in varying periodic boundary conditions

This work presents a combined study of time-resolved fluorescence spectroscopy and all-atom molecular dynamics simulation to investigate periodic boundary conditions' and water models' influence on the orientation dynamics and translational and rotational diffusion of peptides in solution. We have characterized the effects of solvent box size and water model choice on the dynamics of two peptide systems, NATA and WK5. Computationally, translational, and rotational diffusion and internal fluctuations are investigated through all-atom molecular dynamics simulation with two water models and different box sizes. These results are compared with time-resolved fluorescence anisotropy decay (FAD) measurements. The associated time constant and orientation dynamics from FAD measurement along the ¹L_b axis provided baseline data to validate molecular dynamics simulation. The modeling results show that diffusion rates vary roughly in inverse proportion to water model viscosity, as one would expect. Corrections for finite box size are significant for translational diffusion and insignificant for rotational diffusion. This study also finds that internal dynamics described by autocorrelation functions and kinetic network models are relatively insensitive to both box size and water model properties. Our observation suggests that different peptide properties respond differently to a change in simulation conditions.Communicated by Ramaswamy H. Sarma.

Explicit models of motions to analyze NMR relaxation data in proteins

Nuclear Magnetic Resonance (NMR) is a tool of choice to characterize molecular motions. In biological macromolecules, pico- to nano-second motions, in particular, can be probed by nuclear spin relaxation rates which depend on the time fluctuations of the orientations of spin interaction frames. For the past 40 years, relaxation rates have been successfully analyzed using the Model Free (MF) approach which makes no assumption on the nature of motions and reports on the effective amplitude and time-scale of the motions. However, obtaining a mechanistic picture of motions from this type of analysis is difficult at best, unless complemented with molecular dynamics (MD) simulations. In spite of their limited accuracy, such simulations can be used to obtain the information necessary to build explicit models of motions designed to analyze NMR relaxation data. Here, we present how to build such models, suited in particular to describe motions of methyl-bearing protein side-chains and compare them with the MF approach. We show on synthetic data that explicit models of motions are more robust in the presence of rotamer jumps which dominate the relaxation in methyl groups of protein side-chains. We expect this work to motivate the use of explicit models of motion to analyze MD and NMR data.

Zinc finger structure determination by NMR: Why zinc fingers can be a handful

Zinc fingers can be loosely defined as protein domains containing one or more tetrahedrally-co-ordinated zinc ions whose role is to stabilise the structure rather than to be involved in enzymatic chemistry; such zinc ions are often referred to as "structural zincs". Although structural zincs can occur in proteins of any size, they assume particular significance for very small protein domains, where they are often essential for maintaining a folded state. Such small structures, that sometimes have only marginal stability, can present particular difficulties in terms of sample preparation, handling and structure determination, and early on they gained a reputation for being resistant to crystallisation. As a result, NMR has played a more prominent role in structural studies of zinc finger proteins than it has for many other types of proteins. This review will present an overview of the particular issues that arise for structure determination of zinc fingers by NMR, and ways in which these may be addressed.

Molecular Dynamics-Assisted Optimization of Protein NMR Relaxation Analysis

NMR relaxation analysis of the mobile residues in globular proteins is sensitive to the form of the experimentally fitted internal autocorrelation function, which is used to represent that motion. Different order parameter representations can precisely fit the same set of ¹⁵N R₁, R₂, and heteronuclear NOE measurements while yielding significantly divergent predictions of the underlying autocorrelation functions, indicating the insufficiency of these experimental relaxation data for assessing which order parameter representation provides the most physically realistic predictions. Molecular dynamics simulations offer an unparalleled capability for discriminating among different order parameter representations to assess which representation can most accurately model a wide range of physically realistic autocorrelation functions. Six currently utilized AMBER and CHARMM force fields were applied to calculate autocorrelation functions for the backbone H–N bond vectors of ubiquitin as an operational test set. An optimized time constant-constrained triexponential (TCCT) representation was shown to markedly outperform the widely used (S_f²,τ_s,S²) extended Lipari–Szabo representation and the more closely related (S_f²,S_H², S_N²) Larmor frequency-selective representation. Optimization of the TCCT representation at both 600 and 900 MHz ¹H converged to the same parameterization. The higher magnetic field yielded systematically larger deviations in the back-prediction of the autocorrelation functions for the mobile amides, indicating little added benefit from multiple field measurements in analyzing amides that lack slower (∼ms) exchange line-broadening effects. Experimental ¹⁵N relaxation data efficiently distinguished among the different force fields with regard to their prediction of ubiquitin backbone conformational dynamics in the ps–ns time frame. While the earlier AMBER 99SB and CHARMM27 force fields underestimate the scale of backbone dynamics, which occur in this time frame, AMBER 14SB provided the most consistent predictions for the well-averaged highly mobile C-terminal residues of ubiquitin.

Zinc Fingers 10 and 11 of Miz-1 undergo conformational exchange to achieve specific DNA binding

Miz-1 (ZBTB17) is a poly-zinc finger BTB/POZ transcription factor with 12 consecutive C2H2 zinc fingers (ZFs) that binds transcriptional start sites (TSSs) to regulate the expression of genes involved in cell development and proliferation. As of now, it is not known which of the 12 consecutive ZFs are responsible for the recognition of the 24 base pair consensus sequence found at these TSSs. Evidence suggests ZFs 7-12 plays this role. We provide validation for this and describe the structural and dynamical characterization of unprecedented conformational exchange in the linker between ZFs 10 and 11. This conformational exchange uncouples ZFs 7-10 from 11 and 12 and promotes a scanning-recognition mechanism through which the two segments cooperate to bind two sub-sites at both ends of the consensus. We further show that this can result in the coiling of TSSs as part of Miz-1's mechanism of transcriptional transactivation.

Quantification of natural abundance NMR data differentiates the solution behavior of monoclonal antibodies and their fragments

Biotherapeutics are an important class of molecules for the treatment of a wide range of diseases. They include low molecular weight peptides, highly engineered protein scaffolds and monoclonal antibodies. During their discovery and development, assessments of the biophysical attributes is critical to understanding the solution behavior of therapeutic proteins and for de-risking liabilities. Thus, methods that can quantify, characterize, and provide a basis to inform risks and drive the selection of more optimal antibody and alternative scaffolds are needed. Nuclear Magnetic Resonance (NMR) spectroscopy is a technique that provides a means to probe antibody and antibody-like molecules in solution, at atomic resolution, under any formulated conditions. Here, all samples were profiled at natural abundance requiring no isotope enrichment. We present a numerical approach that quantitates two-dimensional methyl spectra. The approach was tested with a reference dataset that contained different types of antibody and antibody-like molecules. This dataset was processed through a procedure we call a Random Sampling of NMR Peaks for Covariance Analysis. This analysis revealed that the first two components were well correlated with the hydrodynamic radius of the molecules included in the reference set. Higher-order principal components were also linked to dynamic features between different tethered antibody-like molecules and contributed to decisions around candidate selection. The reference set provides a basis to characterize molecules with unknown solution behavior and is sensitive to the behavior of a molecule formulated under different conditions. The approach is independent of protein design, scaffold, formulation and provides a facile method to quantify solution behavior.

RNA Binding by the KTS Splice Variants of Wilms' Tumor Suppressor Protein WT1

Wilms' tumor suppressor protein WT1 regulates the expression of multiple genes through binding of the Cys₂-His₂ zinc finger domain to promoter sites. WT1 has also been proposed to be involved in post-transcriptional regulation, by binding to RNA using the same set of zinc fingers. WT1 has two major splice variants, where the Lys-Thr-Ser (KTS) tripeptide is inserted into the linker between the third and fourth zinc fingers. To obtain insights into the mechanism by which the different WT1 splice variants recognize both DNA and RNA, we have determined the solution structure of the WT1 (-KTS) zinc finger domain in complex with a 29mer stem-loop RNA. Zinc fingers 1-3 bind in a widened major groove favored by the presence of a bulge nucleotide in the double-stranded helical stem. Fingers 2 and 3 make specific contacts with the nucleobases in a conserved AUGG sequence in the helical stem. Nuclear magnetic resonance chemical shift mapping and relaxation analysis show that fingers 1-3 of the two splice variants (-KTS and +KTS) of WT1 form similar complexes with RNA. Finger 4 of the -KTS isoform interacts weakly with the RNA loop, an interaction that is abrogated in the +KTS isoform, and both isoforms bind with similar affinity to the RNA. In contrast, finger 4 is required for high-affinity binding to DNA and insertion of KTS into the linker of fingers 3 and 4 abrogates DNA binding. While finger 1 is required for RNA binding, it is dispensable for binding to consensus DNA sites.

13C NMR Relaxation Analysis of Protein GB3 for the Assessment of Side Chain Dynamics Predictions by Current AMBER and CHARMM Force Fields

Molecular simulations with seven current AMBER- and CHARMM-based force fields yield markedly differing internal bond vector autocorrelation function predictions for many of the 223 methine and methylene H-C bonds of the 56-residue protein GB3. To enable quantification of accuracy, ¹³C R₁, R₂, and heteronuclear NOE relaxation rates have been determined for the methine and stereochemically assigned methylene C^α and C^β positions. With only three experimental relaxation values for each bond vector, central to this analysis is the accuracy with which MD-derived autocorrelation curves can be represented by a 3-parameter equation which, in turn, maps onto the NMR relaxation values. In contrast to the more widely used extended Lipari-Szabo order parameter representation, 95% of these MD-derived internal autocorrelation curves for GB3 can be fitted to within 1.0% rmsd over the time frame from 30 ps to 4 ns by a biexponential Larmor frequency-selective representation (LF-S²). Applying the LF-S² representation to the experimental relaxation rates and uncertainties serves to determine the boundary range for the autocorrelation function of each bond vector consistent with the experimental data. Not surprisingly, all seven force fields predict the autocorrelation functions for the more motionally restricted ¹H^α-¹³C^α and ¹H^β-¹³C^β bond vectors with reasonable accuracy. However, for the ¹H^β-¹³C^β bond vectors exhibiting aggregate order parameter S² values less than 0.85, only 1% of the MD-derived predictions lie with 1 σ of the experimentally determined autocorrelation functions and only 7% within 2 σ. On the other hand, substantial residue type-specific improvements in predictive performance were observed among the recent AMBER force fields. This analysis indicates considerable potential for the use of ¹³C relaxation measurements in guiding the optimization of the side chain dynamics characteristics of protein molecular simulations.

Conformational Priming of RepA-WH1 for Functional Amyloid Conversion Detected by NMR Spectroscopy

How proteins with a stable globular fold acquire the amyloid state is still largely unknown. RepA, a versatile plasmidic DNA binding protein from Pseudomonas savastanoi, is functional as a transcriptional repressor or as an initiator or inhibitor of DNA replication, the latter via assembly of an amyloidogenic oligomer. Its N-terminal domain (WH1) is responsible for discrimination between these functional abilities by undergoing insufficiently understood structural changes. RepA-WH1 is a stable dimer whose conformational dynamics had not been explored. Here, we have studied it through NMR {¹H}-¹⁵N relaxation and H/D exchange kinetics measurements. The N- and the C-terminal α-helices, and the internal amyloidogenic loop, are partially unfolded in solution. S4-indigo, a small inhibitor of RepA-WH1 amyloidogenesis, binds to and tethers the N-terminal α-helix to a β-hairpin that is involved in dimerization, thus providing evidence for a priming role of fraying ends and dimerization switches in the amyloidogenesis of folded proteins.

Functional protein dynamics on uncharted time scales detected by nanoparticle-assisted NMR spin relaxation

Protein function depends critically on intrinsic internal dynamics, which is manifested in distinct ways, such as loop motions that regulate protein recognition and catalysis. Under physiological conditions, dynamic processes occur on a wide range of time scales from subpicoseconds to seconds. Commonly used NMR spin relaxation in solution provides valuable information on very fast and slow motions but is insensitive to the intermediate nanosecond to microsecond range that exceeds the protein tumbling correlation time. Presently, very little is known about the nature and functional role of these motions. It is demonstrated here how transverse spin relaxation becomes exquisitely sensitive to these motions at atomic resolution when studying proteins in the presence of nanoparticles. Application of this novel cross-disciplinary approach reveals large-scale dynamics of loops involved in functionally critical protein-protein interactions and protein-calcium ion recognition that were previously unobservable.

Characterization of Internal Protein Dynamics and Conformational Entropy by NMR Relaxation

Recent studies suggest that the fast timescale motion of methyl-bearing side chains may play an important role in mediating protein activity. These motions have been shown to encapsulate the residual conformational entropy of the folded state that can potentially contribute to the energetics of protein function. Here, we provide an overview of how to characterize these motions using nuclear magnetic resonance (NMR) spin relaxation methods. The strengths and limitations of several techniques are highlighted in order to assist with experimental design. Particular emphasis is placed on the practical aspects of sample preparation, data collection, data fitting, and statistical analysis. Additionally, discussion of the recently refined "entropy meter" is presented and its use in converting NMR observables to conformational entropy is illustrated. Taken together, these methods should yield new insights into the complex interplay between structure and dynamics in protein function.

Narrowing the gap between experimental and computational determination of methyl group dynamics in proteins

A synergistic analysis of methyl NMR relaxation data and MD simulations identifies ways to improve studies of protein dynamics.

Atomic resolution conformational dynamics of intrinsically disordered proteins from NMR spin relaxation

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental approaches for investigating the conformational behaviour of intrinsically disordered proteins (IDPs). IDPs represent a significant fraction of all proteomes, and, despite their importance for understanding fundamental biological processes, the molecular basis of their activity still remains largely unknown. The functional mechanisms exploited by IDPs in their interactions with other biomolecules are defined by their intrinsic dynamic modes and associated timescales, justifying the considerable interest over recent years in the development of technologies adapted to measure and describe this behaviour. NMR spin relaxation delivers information-rich, site-specific data reporting on conformational fluctuations occurring throughout the molecule. Here we review recent progress in the use of ¹⁵N relaxation to identify local backbone dynamics and long-range chain-like motions in unfolded proteins.

Comprehensive analysis of NMR data using advanced line shape fitting

NMR spectroscopy is uniquely suited for atomic resolution studies of biomolecules such as proteins, nucleic acids and metabolites, since detailed information on structure and dynamics are encoded in positions and line shapes of peaks in NMR spectra. Unfortunately, accurate determination of these parameters is often complicated and time consuming, in part due to the need for different software at the various analysis steps and for validating the results. Here, we present an integrated, cross-platform and open-source software that is significantly more versatile than the typical line shape fitting application. The software is a completely redesigned version of PINT ( https://pint-nmr.github.io/PINT/ ). It features a graphical user interface and includes functionality for peak picking, editing of peak lists and line shape fitting. In addition, the obtained peak intensities can be used directly to extract, for instance, relaxation rates, heteronuclear NOE values and exchange parameters. In contrast to most available software the entire process from spectral visualization to preparation of publication-ready figures is done solely using PINT and often within minutes, thereby, increasing productivity for users of all experience levels. Unique to the software are also the outstanding tools for evaluating the quality of the fitting results and extensive, but easy-to-use, customization of the fitting protocol and graphical output. In this communication, we describe the features of the new version of PINT and benchmark its performance.

Solution NMR views of dynamical ordering of biomacromolecules

BACKGROUND

To understand the mechanisms related to the 'dynamical ordering' of macromolecules and biological systems, it is crucial to monitor, in detail, molecular interactions and their dynamics across multiple timescales. Solution nuclear magnetic resonance (NMR) spectroscopy is an ideal tool that can investigate biophysical events at the atomic level, in near-physiological buffer solutions, or even inside cells.

SCOPE OF REVIEW

In the past several decades, progress in solution NMR has significantly contributed to the elucidation of three-dimensional structures, the understanding of conformational motions, and the underlying thermodynamic and kinetic properties of biomacromolecules. This review discusses recent methodological development of NMR, their applications and some of the remaining challenges.

MAJOR CONCLUSIONS

Although a major drawback of NMR is its difficulty in studying the dynamical ordering of larger biomolecular systems, current technologies have achieved considerable success in the structural analysis of substantially large proteins and biomolecular complexes over 1MDa and have characterised a wide range of timescales across which biomolecular motion exists. While NMR is well suited to obtain local structure information in detail, it contributes valuable and unique information within hybrid approaches that combine complementary methodologies, including solution scattering and microscopic techniques.

GENERAL SIGNIFICANCE

For living systems, the dynamic assembly and disassembly of macromolecular complexes is of utmost importance for cellular homeostasis and, if dysregulated, implied in human disease. It is thus instructive for the advancement of the study of the dynamical ordering to discuss the potential possibilities of solution NMR spectroscopy and its applications. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

Engineering double-stranded RNA binding activity into the Drosha double-stranded RNA binding domain results in a loss of microRNA processing function

Canonical processing of miRNA begins in the nucleus with the Microprocessor complex, which is minimally composed of the RNase III enzyme Drosha and two copies of its cofactor protein DGCR8. In structural analogy to most RNase III enzymes, Drosha possesses a modular domain with the double-stranded RNA binding domain (dsRBD) fold. Unlike the dsRBDs found in most members of the RNase III family, the Drosha-dsRBD does not display double-stranded RNA binding activity; perhaps related to this, the Drosha-dsRBD amino acid sequence does not conform well to the canonical patterns expected for a dsRBD. In this article, we investigate the impact on miRNA processing of engineering double-stranded RNA binding activity into Drosha's non-canonical dsRBD. Our findings corroborate previous studies that have demonstrated the Drosha-dsRBD is necessary for miRNA processing and suggest that the amino acid composition in the second α-helix of the domain is critical to support its evolved function.

Structural Insights into c-Myc-interacting Zinc Finger Protein-1 (Miz-1) Delineate Domains Required for DNA Scanning and Sequence-specific Binding

c-Myc-interacting zinc finger protein-1 (Miz-1) is a poly-Cys2His2 zinc finger (ZF) transcriptional regulator of many cell cycle genes. A Miz-1 DNA sequence consensus has recently been identified and has also unveiled Miz-1 functions in other cellular processes, underscoring its importance in the cell. Miz-1 contains 13 ZFs, but it is unknown why Miz-1 has so many ZFs and whether they recognize and bind DNA sequences in a typical fashion. Here, we used NMR to deduce the role of Miz-1 ZFs 1-4 in detecting the Miz-1 consensus sequence and preventing nonspecific DNA binding. In the construct containing the first 4 ZFs, we observed that ZFs 3 and 4 form an unusual compact and stable structure that restricts their motions. Disruption of this compact structure by an electrostatically mismatched A86K mutation profoundly affected the DNA binding properties of the WT construct. On the one hand, Miz1-4^WT was found to bind the Miz-1 DNA consensus sequence weakly and through ZFs 1-3 only. On the other hand, the four ZFs in the structurally destabilized Miz1-4^A86K mutant bound to the DNA consensus with a 30-fold increase in affinity (100 nm). The formation of such a thermodynamically stable but nonspecific complex is expected to slow down the rate of DNA scanning by Miz-1 during the search for its consensus sequence. Interestingly, we found that the motif stabilizing the compact structure between ZFs 3 and 4 is conserved and enriched in other long poly-ZF proteins. As discussed in detail, our findings support a general role of compact inter-ZF structures in minimizing the formation of off-target DNA complexes.

High-Resolution NMR Determination of the Dynamic Structure of Membrane Proteins

¹⁵N spin‐relaxation rates are demonstrated to provide critical information about the long‐range structure and internal motions of membrane proteins. Combined with an improved calculation method, the relaxation‐rate‐derived structure of the 283‐residue human voltage‐dependent anion channel revealed an anisotropically shaped barrel with a rigidly attached N‐terminal helix. Our study thus establishes an NMR spectroscopic approach to determine the structure and dynamics of mammalian membrane proteins at high accuracy and resolution.

Changes in conformational dynamics of basic side chains upon protein-DNA association

Basic side chains play major roles in recognition of nucleic acids by proteins. However, dynamic properties of these positively charged side chains are not well understood. In this work, we studied changes in conformational dynamics of basic side chains upon protein–DNA association for the zinc-finger protein Egr-1. By nuclear magnetic resonance (NMR) spectroscopy, we characterized the dynamics of all side-chain cationic groups in the free protein and in the complex with target DNA. Our NMR order parameters indicate that the arginine guanidino groups interacting with DNA bases are strongly immobilized, forming rigid interfaces. Despite the strong short-range electrostatic interactions, the majority of the basic side chains interacting with the DNA phosphates exhibited high mobility, forming dynamic interfaces. In particular, the lysine side-chain amino groups exhibited only small changes in the order parameters upon DNA-binding. We found a similar trend in the molecular dynamics (MD) simulations for the free Egr-1 and the Egr-1–DNA complex. Using the MD trajectories, we also analyzed side-chain conformational entropy. The interfacial arginine side chains exhibited substantial entropic loss upon binding to DNA, whereas the interfacial lysine side chains showed relatively small changes in conformational entropy. These data illustrate different dynamic characteristics of the interfacial arginine and lysine side chains.

A dynamic look backward and forward

The 2015 Gunther Laukien Prize recognized solution NMR studies of protein dynamics and thermodynamics. This Perspective surveys aspects of the development and application of NMR spin relaxation for investigations of protein flexibility and function over multiple time scales in solution. Methods highlighted include analysis of overall rotational diffusion, theoretical descriptions of R1ρ relaxation, and molecular dynamics simulations to interpret NMR spin relaxation. Applications are illustrated for the zinc-finger domain Xfin-31, the calcium-binding proteins calbindin D9k and calmodulin, and the bZip transcription factor of GCN4.

Structure and dynamics of the tetrameric mnt repressor and a model for its DNA complex

Abstract The tetrameric Mnt repressor of bacteriophage P22 consists of two dimeric DNA-binding domains and a tetramerization domain. The NOE and chemical shift data demonstrate that the structures of the domains in the wild-type repressor protein are similar to those of the separate domains, the three-dimensional structures of which have been determined previously. (15)N relaxation measurements show that the linker that connects the anti-parallel four-helix bundle with the two β-sheet DNA-binding dimers is highly flexible. No evidence was found for interactions between the distinct modules. The (15)N relaxation properties of the two domains differ substantially, confirming their structural independence. A model in which one two-stranded coiled coil of the four-helix bundle is attached to one N-terminal dimer is most consistent with the biochemical data and (15)N relaxation data. For the Mnt-DNA complex this geometry fits with a model in which the two β-sheet DNA-binding domains are bound at two successive major grooves of the Mnt operator and the tetramerization domain is packed between these two DNA-bound dimers. In such a model the two-fold symmetry axis of the four-helix bundle coincides with that of the operator sequence and the two bound dimers. Bending of the Mnt operator of approximately 30° upon binding of the tetramer, as measured by gel-shift assays, is in agreement with this model of the Mnt-DNA complex.

Movement and specificity in a modular DNA binding protein

The single-stranded DNA (ssDNA) binding protein RPA binds to and protects ssDNA while simultaneously recruiting numerous replication and repair proteins essential for genome integrity. In this issue of Structure, Brosey et al. (2015) show that the flexibility and interactions of the modular domains of RPA are altered by ssDNA binding and suggest that these changes in configurational freedom are important for the many functions of RPA.

Dynamics of GCN4 facilitate DNA interaction: a model-free analysis of an intrinsically disordered region

Intrinsically disordered proteins (IDPs) and proteins with intrinsically disordered regions (IDRs) are known to play important roles in regulatory and signaling pathways. A critical aspect of these functions is the ability of IDP/IDRs to form highly specific complexes with target molecules. However, elucidation of the contributions of conformational dynamics to function has been limited by challenges associated with structural heterogeneity of IDP/IDRs. Using NMR spin relaxation parameters ((15)N R1, (15)N R2, and {(1)H}-(15)N heteronuclear NOE) collected at four static magnetic fields ranging from 14.1 to 21.1 T, we have analyzed the backbone dynamics of the basic leucine-zipper (bZip) domain of the Saccharomyces cerevisiae transcription factor GCN4, whose DNA binding domain is intrinsically disordered in the absence of DNA substrate. We demonstrate that the extended model-free analysis can be applied to proteins with IDRs such as apo GCN4 and that these results significantly extend previous NMR studies of GCN4 dynamics performed using a single static magnetic field of 11.74 T [Bracken, et al., J. Mol. Biol., 1999, 285, 2133-2146] and correlate well with molecular dynamics simulations [Robustelli, et al., J. Chem. Theory Comput., 2013, 9, 5190-5200]. In contrast to the earlier work, data at multiple static fields allows the time scales of internal dynamics of GCN4 to be reliably quantified. Large amplitude dynamic fluctuations in the DNA-binding region have correlation times (τs ≈ 1.4-2.5 ns) consistent with a two-step mechanism in which partially ordered bZip conformations of GCN4 form initial encounter complexes with DNA and then rapidly rearrange to the high affinity state with fully formed basic region recognition helices.

Linkage-specific conformational ensembles of non-canonical polyubiquitin chains

Polyubiquitination is a critical protein post-translational modification involved in a variety of processes in eukaryotic cells. The molecular basis for selective recognition of the polyubiquitin signals by cellular receptors is determined by the conformations polyubiquitin chains adopt; this has been demonstrated for K48- and K63-linked chains. Recent studies of the so-called non-canonical chains (linked via K6, K11, K27, K29, or K33) suggest they play important regulatory roles in growth, development, and immune system pathways, but biophysical studies are needed to elucidate the physical/structural basis of their interactions with receptors. A first step towards this goal is characterization of the conformations these chains adopt in solution. We assembled diubiquitins (Ub2) comprised of every lysine linkage. Using solution NMR measurements, small-angle neutron scattering (SANS), and in silico ensemble generation, we determined population-weighted conformational ensembles that shed light on the structure and dynamics of the non-canonical polyubiquitin chains. We found that polyubiquitin is conformationally heterogeneous, and each chain type exhibits unique conformational ensembles. For example, K6-Ub2 and K11-Ub2 (at physiological salt concentration) are in dynamic equilibrium between at least two conformers, where one exhibits a unique Ub/Ub interface, distinct from that observed in K48-Ub2 but similar to crystal structures of these chains. Conformers for K29-Ub2 and K33-Ub2 resemble recent crystal structures in the ligand-bound state. Remarkably, a number of diubiquitins adopt conformers similar to K48-Ub2 or K63-Ub2, suggesting potential overlap of biological function among different lysine linkages. These studies highlight the potential power of determining function from elucidation of conformational states.

Halophilic Protein Adaptation Results from Synergistic Residue-Ion Interactions in the Folded and Unfolded States

Halophilic organisms thrive in environments with extreme salt concentrations and have adapted by allowing molar quantities of cosolutes, mainly KCl, to accumulate in their cytoplasm. To cope with this high intracellular salinity, halophilic organisms modified the chemical composition of their proteins to enrich their surface with acidic and short polar side chains, while lysines and bulky hydrophobic residues got depleted. We have emulated the evolutionary process of haloadaptation with natural and designed halophilic polypeptides and applied novel nuclear magnetic resonance (NMR) methodology to study the different mechanisms contributing to protein haloadaptation at a per residue level. Our analysis of an extensive set of NMR observables, determined over several proteins, allowed us to disentangle the synergistic contributions of protein haloadaptation: cation exclusion and electrostatic repulsion between negatively charged residues destabilize the denatured state ensemble while cumulative weak cation-protein interactions stabilize the folded conformations.

Structure of the Bacillus anthracis Sortase A Enzyme Bound to Its Sorting Signal: A FLEXIBLE AMINO-TERMINAL APPENDAGE MODULATES SUBSTRATE ACCESS

The endospore forming bacterium Bacillus anthracis causes lethal anthrax disease in humans and animals. The ability of this pathogen to replicate within macrophages is dependent upon the display of bacterial surface proteins attached to the cell wall by the B. anthracis Sortase A ((Ba)SrtA) enzyme. Previously, we discovered that the class A (Ba)SrtA sortase contains a unique N-terminal appendage that wraps around the body of the protein to contact the active site of the enzyme. To gain insight into its function, we determined the NMR structure of (Ba)SrtA bound to a LPXTG sorting signal analog. The structure, combined with dynamics, kinetics, and whole cell protein display data suggest that the N terminus modulates substrate access to the enzyme. We propose that it may increase the efficiency of protein display by reducing the unproductive hydrolytic cleavage of enzyme-protein covalent intermediates that form during the cell wall anchoring reaction. Notably, a key active site loop (β7/β8 loop) undergoes a disordered to ordered transition upon binding the sorting signal, potentially facilitating recognition of lipid II.

Balancing between affinity and speed in target DNA search by zinc-finger proteins via modulation of dynamic conformational ensemble

Although engineering of transcription factors and DNA-modifying enzymes has drawn substantial attention for artificial gene regulation and genome editing, most efforts focus on affinity and specificity of the DNA-binding proteins, typically overlooking the kinetic properties of these proteins. However, a simplistic pursuit of high affinity can lead to kinetically deficient proteins that spend too much time at nonspecific sites before reaching their targets on DNA. We demonstrate that structural dynamic knowledge of the DNA-scanning process allows for kinetically and thermodynamically balanced engineering of DNA-binding proteins. Our current study of the zinc-finger protein Egr-1 (also known as Zif268) and its nuclease derivatives reveals kinetic and thermodynamic roles of the dynamic conformational equilibrium between two modes during the DNA-scanning process: one mode suitable for search and the other for recognition. By mutagenesis, we were able to shift this equilibrium, as confirmed by NMR spectroscopy. Using fluorescence and biochemical assays as well as computational simulations, we analyzed how the shifts of the conformational equilibrium influence binding affinity, target search kinetics, and efficiency in displacing other proteins from the target sites. A shift toward the recognition mode caused an increase in affinity for DNA and a decrease in search efficiency. In contrast, a shift toward the search mode caused a decrease in affinity and an increase in search efficiency. This accelerated site-specific DNA cleavage by the zinc-finger nuclease, without enhancing off-target cleavage. Our study shows that appropriate modulation of the dynamic conformational ensemble can greatly improve zinc-finger technology, which has used Egr-1 (Zif268) as a major scaffold for engineering.

A multiscale approach to simulating the conformational properties of unbound multi-C₂H₂ zinc finger proteins

The conformational properties of unbound multi-Cys2 His2 (mC2H2) zinc finger proteins, in which zinc finger domains are connected by flexible linkers, are studied by a multiscale approach. Three methods on different length scales are utilized. First, atomic detail molecular dynamics simulations of one zinc finger and its adjacent flexible linker confirmed that the zinc finger is more rigid than the flexible linker. Second, the end-to-end distance distributions of mC2H2 zinc finger proteins are computed using an efficient atomistic pivoting algorithm, which only takes excluded volume interactions into consideration. The end-to-end distance distribution gradually changes its profile, from left-tailed to right-tailed, as the number of zinc fingers increases. This is explained by using a worm-like chain model. For proteins of a few zinc fingers, an effective bending constraint favors an extended conformation. Only for proteins containing more than nine zinc fingers, is a somewhat compacted conformation preferred. Third, a mesoscale model is modified to study both the local and the global conformational properties of multi-C2H2 zinc finger proteins. Simulations of the CCCTC-binding factor (CTCF), an important mC2H2 zinc finger protein for genome spatial organization, are presented.

Conformational plasticity surrounding the active site of NADH oxidase from Thermus thermophilus

Biotechnological applications of enzymes can involve the use of these molecules under nonphysiological conditions. Thus, it is of interest to understand how environmental variables affect protein structure and dynamics and how this ultimately modulates enzyme function. NADH oxidase (NOX) from Thermus thermophilus exemplifies how enzyme activity can be tuned by reaction conditions, such as temperature, cofactor substitution, and the addition of cosolutes. This enzyme catalyzes the oxidation of reduced NAD(P)H to NAD(P)(+) with the concurrent reduction of O2 to H2O2, with relevance to biosensing applications. It is thermophilic, with an optimum temperature of approximately 65°C and sevenfold lower activity at 25°C. Moderate concentrations (≈1M) of urea and other chaotropes increase NOX activity by up to a factor of 2.5 at room temperature. Furthermore, it is a flavoprotein that accepts either FMN or the much larger FAD as cofactor. We have used nuclear magnetic resonance (NMR) titration and (15)N spin relaxation experiments together with isothermal titration calorimetry to study how NOX structure and dynamics are affected by changes in temperature, the addition of urea and the substitution of the FMN cofactor with FAD. The majority of signals from NOX are quite insensitive to changes in temperature, cosolute addition, and cofactor substitution. However, a small cluster of residues surrounding the active site shows significant changes. These residues are implicated in coupling changes in the solution conditions of the enzyme to changes in catalytic activity.

BioJazz: in silico evolution of cellular networks with unbounded complexity using rule-based modeling

Systems biologists aim to decipher the structure and dynamics of signaling and regulatory networks underpinning cellular responses; synthetic biologists can use this insight to alter existing networks or engineer de novo ones. Both tasks will benefit from an understanding of which structural and dynamic features of networks can emerge from evolutionary processes, through which intermediary steps these arise, and whether they embody general design principles. As natural evolution at the level of network dynamics is difficult to study, in silico evolution of network models can provide important insights. However, current tools used for in silico evolution of network dynamics are limited to ad hoc computer simulations and models. Here we introduce BioJazz, an extendable, user-friendly tool for simulating the evolution of dynamic biochemical networks. Unlike previous tools for in silico evolution, BioJazz allows for the evolution of cellular networks with unbounded complexity by combining rule-based modeling with an encoding of networks that is akin to a genome. We show that BioJazz can be used to implement biologically realistic selective pressures and allows exploration of the space of network architectures and dynamics that implement prescribed physiological functions. BioJazz is provided as an open-source tool to facilitate its further development and use. Source code and user manuals are available at: http://oss-lab.github.io/biojazz and http://osslab.lifesci.warwick.ac.uk/BioJazz.aspx.

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.

The PRE-Derived NMR Model of the 38.8-kDa Tri-Domain IsdH Protein from Staphylococcus aureus Suggests That It Adaptively Recognizes Human Hemoglobin

Staphylococcus aureus is a medically important bacterial pathogen that, during infections, acquires iron from human hemoglobin (Hb). It uses two closely related iron-regulated surface determinant (Isd) proteins to capture and extract the oxidized form of heme (hemin) from Hb, IsdH and IsdB. Both receptors rapidly extract hemin using a conserved tri-domain unit consisting of two NEAT (near iron transporter) domains connected by a helical linker domain. To gain insight into the mechanism of extraction, we used NMR to investigate the structure and dynamics of the 38.8-kDa tri-domain IsdH protein (IsdH(N2N3), A326-D660 with a Y642A mutation that prevents hemin binding). The structure was modeled using long-range paramagnetic relaxation enhancement (PRE) distance restraints, dihedral angle, small-angle X-ray scattering, residual dipolar coupling and inter-domain NOE nuclear Overhauser effect data. The receptor adopts an extended conformation wherein the linker and N3 domains pack against each other via a hydrophobic interface. In contrast, the N2 domain contacts the linker domain via a hydrophilic interface and, based on NMR relaxation data, undergoes inter-domain motions enabling it to reorient with respect to the body of the protein. Ensemble calculations were used to estimate the range of N2 domain positions compatible with the PRE data. A comparison of the Hb-free and Hb-bound forms reveals that Hb binding alters the positioning of the N2 domain. We propose that binding occurs through a combination of conformational selection and induced-fit mechanisms that may promote hemin release from Hb by altering the position of its F helix.

NMR studies of the dynamics of high-spin nitrophorins: comparative studies of NP4 and NP2 at close to physiological pH

The β-barrel nitrophorin (NP) heme proteins are found in the saliva of the blood-sucking insect Rhodnius prolixus, which synthesizes and stores nitric oxide (NO) in the salivary glands. NO is bound to iron of the NPs and is released by dilution and an increase in pH when the insect spits its saliva into the tissues of a victim, to aid in obtaining a blood meal. In the adult insect, there are four nitrophorins, NP1–NP4, which have sequence similarities in two pairs, NP1 and NP4 (90% identical) and NP2 and NP3 (80% identical). The available crystal structures of NP4 have been used to propose that pH-dependent changes in the conformation of two loops between adjacent β-strands at the front opening of the protein, the A–B and G–H loops, determine the rate of NO release. At pH 7.3, NP4 releases NO 17 times faster than NP2 does. In this work, the aqua complexes of NP4 and NP2 have been investigated by nuclear magnetic resonance (NMR) relaxation measurements to probe the pico- to nanosecond and micro- to millisecond time scale motions at two pH values, 6.5 and 7.3. It is found that NP4-OH₂ is fairly rigid and only residues in the loop regions show dynamics at pH 6.5; at pH 7.3, much more dynamics of the loops and most of the β-strands are observed while the α-helices remain fairly rigid. In comparison, NP2-OH₂ shows much less dynamics, albeit somewhat more than that of the previously reported NP2-NO complex [Muthu, D., Berry, R. E., Zhang, H., and Walker, F. A. (2013) Biochemistry 52, 7910–7925]. The reasons for this major difference between NP4 and NP2 are discussed.

Local isotropic diffusion approximation for coupled internal and overall molecular motions in NMR spin relaxation

The present work demonstrates that NMR spin relaxation rate constants for molecules interconverting between states with different diffusion tensors can be modeled theoretically by combining orientational correlation functions for exchanging spherical molecules with locally isotropic approximations for the diffusion anisotropic tensors. The resulting expressions are validated by comparison with correlation functions obtained by Monte Carlo simulations and are accurate for moderate degrees of diffusion anisotropy typically encountered in investigations of globular proteins. The results are complementary to an elegant, but more complex, formalism that is accurate for all degrees of diffusion anisotropy [Ryabov, Y.; Clore, G. M.; Schwieters, C. D. J. Chem. Phys. 2012, 136, 034108].

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.

NMR approaches for structural analysis of multidomain proteins and complexes in solution

NMR spectroscopy is a key method for studying the structure and dynamics of (large) multidomain proteins and complexes in solution. It plays a unique role in integrated structural biology approaches as especially information about conformational dynamics can be readily obtained at residue resolution. Here, we review NMR techniques for such studies focusing on state-of-the-art tools and practical aspects. An efficient approach for determining the quaternary structure of multidomain complexes starts from the structures of individual domains or subunits. The arrangement of the domains/subunits within the complex is then defined based on NMR measurements that provide information about the domain interfaces combined with (long-range) distance and orientational restraints. Aspects discussed include sample preparation, specific isotope labeling and spin labeling; determination of binding interfaces and domain/subunit arrangements from chemical shift perturbations (CSP), nuclear Overhauser effects (NOEs), isotope editing/filtering, cross-saturation, and differential line broadening; and based on paramagnetic relaxation enhancements (PRE) using covalent and soluble spin labels. Finally, the utility of complementary methods such as small-angle X-ray or neutron scattering (SAXS, SANS), electron paramagnetic resonance (EPR) or fluorescence spectroscopy techniques is discussed. The applications of NMR techniques are illustrated with studies of challenging (high molecular weight) protein complexes.

The role of dynamics in modulating ligand exchange in intracellular lipid binding proteins

Lipids are essential for many biological processes and crucial in the pathogenesis of several diseases. Intracellular lipid-binding proteins (iLBPs) provide mobile hydrophobic binding sites that allow hydrophobic or amphipathic lipid molecules to penetrate into and across aqueous layers. Thus iLBPs mediate the lipid transport within the cell and participate to a spectrum of tissue-specific pathways involved in lipid homeostasis. Structural studies have shown that iLBPs' binding sites are inaccessible from the bulk, implying that substrate binding should involve a conformational change able to produce a ligand entry portal. Many studies have been reported in the last two decades on iLBPs indicating that their dynamics play a pivotal role in regulating ligand binding and targeted release. The ensemble of reported data has not been reviewed until today. This review is thus intended to summarize and possibly generalize the results up to now described, providing a picture which could help to identify the missing notions necessary to improve our understanding of the role of dynamics in iLBPs' molecular recognition. Such notions would clarify the chemistry of lipid binding to iLBPs and set the basis for the development of new drugs.

Nonlinear excitations match correlated motions unveiled by NMR in proteins: a new perspective on allosteric cross-talk

In this paper we propose a novel theoretical framework for interpreting long-range dynamical correlations unveiled in proteins through NMR measurements. The theoretical rationale relies on the hypothesis that correlated motions in proteins may be reconstructed as large-scale, collective modes sustained by long-lived nonlinear vibrations known as discrete breathers (DB) localized at key, hot-spot sites. DBs are spatially localized modes, whose nonlinear nature hinders resonant coupling with the normal modes, thus conferring them long lifetimes as compared to normal modes. DBs have been predicted to exist in proteins, localized at few hot-spot residues typically within the stiffest portions of the structure. We compute DB modes analytically in the framework of the nonlinear network model, showing that the displacement patterns of many DBs localized at key sites match to a remarkable extent the experimentally uncovered correlation blueprint. The computed dispersion relations prove that it is physically possible for some of these DBs to be excited out of thermal fluctuations at room temperature. Based on our calculations, we speculate that transient energy redistribution among the vibrational modes in a protein might favor the emergence of DB-like bursts of long-lived energy at hot-spot sites with lifetimes in the ns range, able to sustain critical, function-encoding correlated motions. More generally, our calculations provide a novel quantitative tool to predict fold-spanning dynamical pathways of correlated residues that may be central to allosteric cross-talk in proteins.

Deriving quantitative dynamics information for proteins and RNAs using ROTDIF with a graphical user interface

To facilitate rigorous analysis of molecular motions in proteins, DNA, and RNA, we present a new version of ROTDIF, a program for determining the overall rotational diffusion tensor from single- or multiple-field nuclear magnetic resonance relaxation data. We introduce four major features that expand the program's versatility and usability. The first feature is the ability to analyze, separately or together, (13)C and/or (15)N relaxation data collected at a single or multiple fields. A significant improvement in the accuracy compared to direct analysis of R2/R1 ratios, especially critical for analysis of (13)C relaxation data, is achieved by subtracting high-frequency contributions to relaxation rates. The second new feature is an improved method for computing the rotational diffusion tensor in the presence of biased errors, such as large conformational exchange contributions, that significantly enhances the accuracy of the computation. The third new feature is the integration of the domain alignment and docking module for relaxation-based structure determination of multi-domain systems. Finally, to improve accessibility to all the program features, we introduced a graphical user interface that simplifies and speeds up the analysis of the data. Written in Java, the new ROTDIF can run on virtually any computer platform. In addition, the new ROTDIF achieves an order of magnitude speedup over the previous version by implementing a more efficient deterministic minimization algorithm. We not only demonstrate the improvement in accuracy and speed of the new algorithm for synthetic and experimental (13)C and (15)N relaxation data for several proteins and nucleic acids, but also show that careful analysis required especially for characterizing RNA dynamics allowed us to uncover subtle conformational changes in RNA as a function of temperature that were opaque to previous analysis.

NMR studies of the dynamics of nitrophorin 2 bound to nitric oxide

The Rhodnius nitrophorins are β-barrel proteins of the lipocalin fold with a heme protruding from the open end of the barrel. They are found in the saliva of the blood-sucking insect Rhodnius prolixus, which synthesizes and stores nitric oxide (NO) in the salivary glands, where NO is bound to iron. NO is released by dilution and an increase in pH when the insect spits its saliva into the tissues of a victim, to aid in obtaining a blood meal. In the adult insect, there are four nitrophorins, NP1-NP4. At pH 7.3, NP4 releases NO 17 times faster than NP2 does, as measured by stopped-flow kinetics. A number of crystal structures of the least abundant protein, NP4, are available. These structures have been used to propose that two loops between adjacent β-strands at the front opening of the protein, the A-B and G-H loops, determine the rate of NO release. To learn how the protein loops contribute to the release of NO for each of the nitrophorins, the dynamics of these proteins are being studied in our laboratory. In this work, the NP2-NO complex has been investigated by nuclear magnetic resonance relaxation measurements to probe the picosecond-to-nanosecond and microsecond-to-millisecond time scale motions at three pH values, 5.0, 6.5, and 7.3. It is found that at pH 5.0 and 6.5, the NP2-NO complex is rigid and only a few residues in the loop regions show dynamics, while at pH 7.3, somewhat more dynamics, particularly of the A-B loop, are observed. Comparison to other lipocalins shows that all are relatively rigid, and that the dynamics of lipocalins in general are much more subtle than those of mainly α-helical proteins.

Structural and dynamical characterization of the Miz-1 zinc fingers 5-8 by solution-state NMR

Myc-interacting zinc finger protein-1 (Miz-1) is a BTB/POZ transcription factor that activates the transcription of cytostatic genes, such as p15(INK4B) or p21(CIP1). The C-terminus of Miz-1 contains 13 consensus C2H2 zinc finger domains (ZF). ZFs 1-4 have been shown to interact with SMAD3/4, while the remaining ZFs are expected to bind the promoters of target genes. We have noted unusual features in ZF 5 and the linker between ZFs 5 and 6. Indeed, a glutamate is found instead of the conserved basic residue two positions before the second zinc-coordinating histidine on the ZF 5 helix, and the linker sequence is DTDKE in place of the classical TGEKP sequence. In a canonical ββα fold, such unusual primary structure elements should cause severe electrostatic repulsions. In this context, we have characterized the structure and the dynamics of a Miz-1 construct comprising ZFs 5-8 (Miz 5-8) by solution-state NMR. Whilst ZFs 5, 7 and 8 were shown to adopt the classical ββα fold for C2H2 ZFs, the number of long-range NOEs was insufficient to define a classical fold for ZF 6. We show by using (15)N-relaxation dispersion experiments that this lack of NOEs is due to the presence of extensive motions on the μs-ms timescale. Since this negatively charged region would have to be located near the phosphodiester backbone in a DNA complex, we propose that in addition to promoting conformational searches, it could serve as a hinge region to keep ZFs 1-4 away from DNA.

NMR methods for structural studies of large monomeric and multimeric proteins

NMR structural studies of large monomeric and multimeric proteins face distinct challenges. In large monomeric proteins, the common occurrence of frequency degeneracies between residues impedes unambiguous assignment of NMR signals. To overcome this barrier, nonuniform sampling (NUS) is used to measure spectra with optimal resolution within reasonable time, new correlation maps resolve previous impasses in assignment strategies, and novel selective methyl labeling schemes provide additional structural probes without cluttering NMR spectra. These advances push the limits of NMR studies of large monomeric proteins. Large multimeric and multidomain proteins are studied by NMR when individual components can also be studied by NMR and have known structures. The structural properties of large assemblies are obtained by identifying binding surfaces, by orienting domains, and employing limited distance constraints. Segmental labeling and the combination of NMR with other methods have helped popularize NMR studies of such systems.

Measuring dynamic and kinetic information in the previously inaccessible supra-τ(c) window of nanoseconds to microseconds by solution NMR spectroscopy

Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful tool that has enabled experimentalists to characterize molecular dynamics and kinetics spanning a wide range of time-scales from picoseconds to days. This review focuses on addressing the previously inaccessible supra-tc window (defined as τ(c) < supra-τ(c) < 40 μs; in which tc is the overall tumbling time of a molecule) from the perspective of local inter-nuclear vector dynamics extracted from residual dipolar couplings (RDCs) and from the perspective of conformational exchange captured by relaxation dispersion measurements (RD). The goal of the first section is to present a detailed analysis of how to extract protein dynamics encoded in RDCs and how to relate this information to protein functionality within the previously inaccessible supra-τ(c) window. In the second section, the current state of the art for RD is analyzed, as well as the considerable progress toward pushing the sensitivity of RD further into the supra-τ(c) scale by up to a factor of two (motion up to 25 μs). From the data obtained with these techniques and methodology, the importance of the supra-τ(c) scale for protein function and molecular recognition is becoming increasingly clearer as the connection between motion on the supra-τ(c) scale and protein functionality from the experimental side is further strengthened with results from molecular dynamics simulations.

Determining interdomain structure and dynamics of a retroviral capsid protein in the presence of oligomerization: implication for structural transition in capsid assembly

Capsid (CA) proteins from all retroviruses, including HIV-1, are structurally homologous dual-domain helical proteins. They form a capsid lattice composed of unitary symmetric CA hexamers. X-ray crystallography has shown that within each hexamer a monomeric CA adopts a single conformation, where most helices are parallel to the symmetry axis. In solution, large differences in averaged NMR spin relaxation rates for the two domains were observed, suggesting they are dynamically independent. One relevant question for the capsid assembly remains: whether the interdomain conformer within a hexamer unit needs to be induced or pre-exists within the conformational space of a monomeric CA. The latter seems more consistent with the relaxation data. However, possible CA protein oligomerization and the structure of each domain will affect relaxation measurements and data interpretation. This study, using CA proteins from equine infectious anemia virus (EIAV) as an example, demonstrates a linear extrapolation approach to obtain backbone (15)N spin relaxation time ratios T1/T2 for a monomeric EIAV-CA in the presence of oligomerization equilibrium. The interdomain motion turns out to be limited. The large difference in the domain averaged <T1/T2> for a CA monomer is a consequence of the orthogonal distributions of helices in the two domains. The new monomeric interdomain conformation in solution is significantly different from that in CA hexamer. Therefore, if capsid assembly follows a nucleation-propagation process, the interdomain conformational change might be a key step during the nucleation, as the configuration in hexagonal assembly is never formed by diffusion of its two domains in solution.