Alexander Grishaev

Structure/Function Implications in a Dynamic Complex of the Intrinsically Disordered Sic1 with the Cdc4 Subunit of an SCF Ubiquitin Ligase

Intrinsically disordered proteins can form highly dynamic complexes with partner proteins. One such dynamic complex involves the intrinsically disordered Sic1 with its partner Cdc4 in regulation of yeast cell cycle progression. Phosphorylation of six N-terminal Sic1 sites leads to equilibrium engagement of each phosphorylation site with the primary binding pocket in Cdc4, the substrate recognition subunit of a ubiquitin ligase. ENSEMBLE calculations using experimental nuclear magnetic resonance and small-angle X-ray scattering data reveal significant transient structure in both phosphorylation states of the isolated ensembles (Sic1 and pSic1) that modulates their electrostatic potential, suggesting a structural basis for the proposed strong contribution of electrostatics to binding. A structural model of the dynamic pSic1-Cdc4 complex demonstrates the spatial arrangements in the ubiquitin ligase complex. These results provide a physical picture of a protein that is predominantly disordered in both its free and bound states, enabling aspects of its structure/function relationship to be elucidated.

Improved fitting of solution X-ray scattering data to macromolecular structures and structural ensembles by explicit water modeling.

A new procedure, AXES, is introduced for fitting small-angle X-ray scattering (SAXS) data to macromolecular structures and ensembles of structures. By using explicit water models to account for the effect of solvent, and by restricting the adjustable fitting parameters to those that dominate experimental uncertainties, including sample/buffer rescaling, detector dark current, and, within a narrow range, hydration layer density, superior fits between experimental high resolution structures and SAXS data are obtained. AXES results are found to be more discriminating than standard Crysol fitting of SAXS data when evaluating poorly or incorrectly modeled protein structures. AXES results for ensembles of structures previously generated for ubiquitin show improved fits over fitting of the individual members of these ensembles, indicating these ensembles capture the dynamic behavior of proteins in solution.

Solution Structure of the 128 kDa Enzyme I Dimer from Escherichia coli and Its 146 kDa Complex with HPr Using Residual Dipolar Couplings and Small- and Wide-Angle X-ray Scattering

The solution structures of free Enzyme I (EI, ∼128 kDa, 575 × 2 residues), the first enzyme in the bacterial phosphotransferase system, and its complex with HPr (∼146 kDa) have been solved using novel methodology that makes use of prior structural knowledge (namely, the structures of the dimeric EIC domain and the isolated EIN domain both free and complexed to HPr), combined with residual dipolar coupling (RDC), small- (SAXS) and wide- (WAXS) angle X-ray scattering and small-angle neutron scattering (SANS) data. The calculational strategy employs conjoined rigid body/torsion/Cartesian simulated annealing, and incorporates improvements in calculating and refining against SAXS/WAXS data that take into account complex molecular shapes in the description of the solvent layer resulting in a better representation of the SAXS/WAXS data. The RDC data orient the symmetrically related EIN domains relative to the C(2) symmetry axis of the EIC dimer, while translational, shape, and size information is provided by SAXS/WAXS. The resulting structures are independently validated by SANS. Comparison of the structures of the free EI and the EI-HPr complex with that of the crystal structure of a trapped phosphorylated EI intermediate reveals large (∼70-90°) hinge body rotations of the two subdomains comprising the EIN domain, as well as of the EIN domain relative to the dimeric EIC domain. These large-scale interdomain motions shed light on the structural transitions that accompany the catalytic cycle of EI.

Structure and Dynamics of Full-Length HIV-1 Capsid Protein in Solution.

The HIV-1 capsid protein plays a crucial role in viral infectivity, assembling into a cone that encloses the viral RNA. In the mature virion, the N-terminal domain of the capsid protein forms hexameric and pentameric rings, while C-terminal domain homodimers connect adjacent N-terminal domain rings to one another. Structures of disulfide-linked hexamer and pentamer assemblies, as well as structures of the isolated domains, have been solved previously. The dimer configuration in C-terminal domain constructs differs in solution (residues 144-231) and crystal (residues 146-231) structures by ∼30°, and it has been postulated that the former connects the hexamers while the latter links pentamers to hexamers. Here we study the structure and dynamics of full-length capsid protein in solution, comprising a mixture of monomeric and dimeric forms in dynamic equilibrium, using ensemble simulated annealing driven by experimental NMR residual dipolar couplings and X-ray scattering data. The complexity of the system necessitated the development of a novel computational framework that should be generally applicable to many other challenging systems that currently escape structural characterization by standard application of mainstream techniques of structural biology. We show that the orientation of the C-terminal domains in dimeric full-length capsid and isolated C-terminal domain constructs is the same in solution, and we obtain a quantitative description of the conformational space sampled by the N-terminal domain relative to the C-terminal domain on the nano- to millisecond time scale. The positional distribution of the N-terminal domain relative to the C-terminal domain is large and modulated by the oligomerization state of the C-terminal domain. We also show that a model of the hexamer/pentamer assembly can be readily generated with a single configuration of the C-terminal domain dimer, and that capsid assembly likely proceeds via conformational selection of sparsely populated configurations of the N-terminal domain within the capsid protein dimer.

CLOUDS, a protocol for deriving a molecular proton density via NMR

We demonstrate the feasibility of computing realistic spatial proton distributions for proteins in solution from experimental NMR nuclear Overhauser effect data only and with minimal assignments. The method, CLOUDS, relies on precise and abundant interproton distance restraints calculated via a relaxation matrix analysis of sets of experimental nuclear Overhauser effect spectroscopy crosspeaks. The MIDGE protocol was adapted for this purpose. A gas of unassigned, unconnected H atoms is condensed into a structured proton distribution (cloud) via a molecular dynamics simulated-annealing scheme in which the internuclear distances and van der Waals repulsive terms are the only active restraints. Proton densities are generated by combining a large number of such clouds, each computed from a different trajectory. After filtering by reference to the cloud closest to the mean, a minimal dispersion proton density (foc) is identified. The latter affords a quasi-continuous hydrogen-only probability distribution that conveys immediate information on the protein surface topology (grooves, protrusions, potential binding site cavities, etc.), directly related to the molecular structure. Feasibility of the method was tested on NMR data measured on two globular protein domains of low regular secondary structure content, the col 2 domain of human matrix metalloproteinase-2 and the kringle 2 domain of human plasminogen, of 60 and 83 amino acid residues, respectively.

The Impact of Hydrogen Bonding on Amide 1H Chemical Shift Anisotropy Studied by Cross-Correlated Relaxation and Liquid Crystal NMR Spectroscopy

Site-specific 1H chemical shift anisotropy (CSA) tensors have been derived for the well-ordered backbone amide moieties in the B3 domain of protein G (GB3). Experimental input data include residual chemical shift anisotropy (RCSA), measured in six mutants that align differently relative to the static magnetic field when dissolved in a liquid crystalline Pf1 suspension, and cross-correlated relaxation rates between the 1HN CSA tensor and either the 1H−15N, the 1H−13C′, or the 1H−13Cα dipolar interactions. Analyses with the assumption that the 1HN CSA tensor is symmetric with respect to the peptide plane (three-parameter fit) or without this premise (five-parameter fit) yield very similar results, confirming the robustness of the experimental input data, and that, to a good approximation, one of the principal components orients orthogonal to the peptide plane. 1HN CSA tensors are found to deviate strongly from axial symmetry, with the most shielded tensor component roughly parallel to the N−H vector, and the least shielded component orthogonal to the peptide plane. DFT calculations on pairs of N-methyl acetamide and acetamide in H-bonded geometries taken from the GB3 X-ray structure correlate with experimental data and indicate that H-bonding effects dominate variations in the 1HN CSA. Using experimentally derived 1HN CSA tensors, the optimal relaxation interference effect needed for narrowest 1HN TROSY line widths is found at ∼1200 MHz.

The second type II module from human matrix metalloproteinase 2: Structure, function and dynamics

Abstract Background: Matrix metalloproteinase 2 (MMP-2, gelatinase A, 72 kDa type IV collagenase) has an important role in extracellular matrix degradation during cell migration and tissue remodeling. It is involved in development, inflammation, wound healing, tumor invasion, metastasis and other physiological and pathological processes. The enzyme cleaves several types of collagen, elastin, fibronectin and laminin. Binding to collagen is mediated by three repeats homologous to fibronectin type II modules, which are inserted in the catalytic domain in proximity to the active site. Results: We have determined the NMR solution structure of the second type II module from human MMP-2 (col-2). The module exhibits a typical type II fold with two short double-stranded antiparallel β sheets and three large loops packed around a cluster of conserved aromatic residues. Backbone amide dynamics, derived from 15 N relaxation experiments, correlate well with solvent accessibility and intramolecular hydrogen bonding. A synthetic peptide with the collagen consensus sequence, (Pro-Pro-Gly) 6 , is shown to interact with the module. Conclusions: Spectral perturbations induced by (Pro-Pro-Gly) 6 binding reveal the region involved in the interaction of col-2 with collagen. The binding surface comprises exposed aromatic residues Phe21, Tyr38, Trp40, Tyr47, Tyr53 and Phe55, and the neighboring Gly33–Gly37 segment.

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.

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

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

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 McClintocks 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.