The hexameric structure of the human mitochondrial replicative helicase Twinkle

The mitochondrial replicative helicase Twinkle is involved in strand separation at the replication fork of mitochondrial DNA (mtDNA). Twinkle malfunction is associated with rare diseases that include late onset mitochondrial myopathies, neuromuscular disorders and fatal infantile mtDNA depletion syndrome. We examined its 3D structure by electron microscopy (EM) and small angle X-ray scattering (SAXS) and built the corresponding atomic models, which gave insight into the first molecular architecture of a full-length SF4 helicase that includes an N-terminal zinc-binding domain (ZBD), an intermediate RNA polymerase domain (RPD) and a RecA-like hexamerization C-terminal domain (CTD). The EM model of Twinkle reveals a hexameric two-layered ring comprising the ZBDs and RPDs in one layer and the CTDs in another. In the hexamer, contacts in trans with adjacent subunits occur between ZBDs and RPDs, and between RPDs and CTDs. The ZBDs show important structural heterogeneity. In solution, the scattering data are compatible with a mixture of extended hexa- and heptameric models in variable conformations. Overall, our structural data show a complex network of dynamic interactions that reconciles with the structural flexibility required for helicase activity.

pyDockSAXS: protein-protein complex structure by SAXS and computational docking

Structural characterization of protein–protein interactions at molecular level is essential to understand biological processes and identify new therapeutic opportunities. However, atomic resolution structural techniques cannot keep pace with current advances in interactomics. Low-resolution structural techniques, such as small-angle X-ray scattering (SAXS), can be applied at larger scale, but they miss atomic details. For efficient application to protein–protein complexes, low-resolution information can be combined with theoretical methods that provide energetic description and atomic details of the interactions. Here we present the pyDockSAXS web server (http://life.bsc.es/pid/pydocksaxs) that provides an automatic pipeline for modeling the structure of a protein–protein complex from SAXS data. The method uses FTDOCK to generate rigid-body docking models that are subsequently evaluated by a combination of pyDock energy-based scoring function and their capacity to describe SAXS data. The only required input files are structural models for the interacting partners and a SAXS curve. The server automatically provides a series of structural models for the complex, sorted by the pyDockSAXS scoring function. The user can also upload a previously computed set of docking poses, which opens the possibility to filter the docking solutions by potential interface residues or symmetry restraints. The server is freely available to all users without restriction.

p15PAF is an intrinsically disordered protein with nonrandom structural preferences at sites of interaction with other proteins

We present to our knowledge the first structural characterization of the proliferating-cell-nuclear-antigen-associated factor p15(PAF), showing that it is monomeric and intrinsically disordered in solution but has nonrandom conformational preferences at sites of protein-protein interactions. p15(PAF) is a 12 kDa nuclear protein that acts as a regulator of DNA repair during DNA replication. The p15(PAF) gene is overexpressed in several types of human cancer. The nearly complete NMR backbone assignment of p15(PAF) allowed us to measure 86 N-H(N) residual dipolar couplings. Our residual dipolar coupling analysis reveals nonrandom conformational preferences in distinct regions, including the proliferating-cell-nuclear-antigen-interacting protein motif (PIP-box) and the KEN-box (recognized by the ubiquitin ligase that targets p15(PAF) for degradation). In accordance with these findings, analysis of the (15)N R2 relaxation rates shows a relatively reduced mobility for the residues in these regions. The agreement between the experimental small angle x-ray scattering curve of p15(PAF) and that computed from a statistical coil ensemble corrected for the presence of local secondary structural elements further validates our structural model for p15(PAF). The coincidence of these transiently structured regions with protein-protein interaction and posttranslational modification sites suggests a possible role for these structures as molecular recognition elements for p15(PAF).

Multiple stable conformations account for reversible concentration-dependent oligomerization and autoinhibition of a metamorphic metallopeptidase

Molecular plasticity controls enzymatic activity: the native fold of a protein in a given environment is normally unique and at a global free-energy minimum. Some proteins, however, spontaneously undergo substantial fold switching to reversibly transit between defined conformers, the "metamorphic" proteins. Here, we present a minimal metamorphic, selective, and specific caseinolytic metallopeptidase, selecase, which reversibly transits between several different states of defined three-dimensional structure, which are associated with loss of enzymatic activity due to autoinhibition. The latter is triggered by sequestering the competent conformation in incompetent but structured dimers, tetramers, and octamers. This system, which is compatible with a discrete multifunnel energy landscape, affords a switch that provides a reversible mechanism of control of catalytic activity unique in nature.

Conformational transitions in human translin enable nucleic acid binding

Translin is a highly conserved RNA- and DNA-binding protein that plays essential roles in eukaryotic cells. Human translin functions as an octamer, but in the octameric crystallographic structure, the residues responsible for nucleic acid binding are not accessible. Moreover, electron microscopy data reveal very different octameric configurations. Consequently, the functional assembly and the mechanism of nucleic acid binding by the protein remain unclear. Here, we present an integrative study combining small-angle X-ray scattering (SAXS), site-directed mutagenesis, biochemical analysis and computational techniques to address these questions. Our data indicate a significant conformational heterogeneity for translin in solution, formed by a lesser-populated compact octameric state resembling the previously solved X-ray structure, and a highly populated open octameric state that had not been previously identified. On the other hand, our SAXS data and computational analyses of translin in complex with the RNA oligonucleotide (GU)₁₂ show that the internal cavity found in the octameric assemblies can accommodate different nucleic acid conformations. According to this model, the nucleic acid binding residues become accessible for binding, which facilitates the entrance of the nucleic acids into the cavity. Our data thus provide a structural basis for the functions that translin performs in RNA metabolism and transport.

Flexible-meccano: a tool for the generation of explicit ensemble descriptions of intrinsically disordered proteins and their associated experimental observables

MOTIVATION

Intrinsically disordered proteins (IDPs) represent a significant fraction of the human proteome. The classical structure function paradigm that has successfully underpinned our understanding of molecular biology breaks down when considering proteins that have no stable tertiary structure in their functional form. One convenient approach is to describe the protein in terms of an equilibrium of rapidly inter-converting conformers. Currently, tools to generate such ensemble descriptions are extremely rare, and poorly adapted to the prediction of experimental data.

RESULTS

We present flexible-meccano-a highly efficient algorithm that generates ensembles of molecules, on the basis of amino acid-specific conformational potentials and volume exclusion. Conformational sampling depends uniquely on the primary sequence, with the possibility of introducing additional local or long-range conformational propensities at an amino acid-specific resolution. The algorithm can also be used to calculate expected values of experimental parameters measured at atomic or molecular resolution, such as nuclear magnetic resonance (NMR) and small angle scattering, respectively. We envisage that flexible-meccano will be useful for researchers who wish to compare experimental data with those expected from a fully disordered protein, researchers who see experimental evidence of deviation from 'random coil' behaviour in their protein, or researchers who are interested in working with a broad ensemble of conformers representing the flexibility of the IDP of interest.

AVAILABILITY

A fully documented multi-platform executable is provided, with examples, at http://www.ibs.fr/science-213/scientific-output/software/flexible-meccano/

CONTACT

martin.blackledge@ibs.fr.

Measurement and analysis of NMR residual dipolar couplings for the study of intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are predicted to represent a significant fraction of all functional proteins. Their inherent plasticity allows them to sample more efficiently their surroundings and thereby increase the probability of interaction with one or several different biological partners. Due to their high flexibility, IDPs cannot be represented by a single, three-dimensional structure; rather, an ensemble description can be invoked, where the protein is assumed to interconvert between different conformations. This chapter focuses on the use of NMR spectroscopy to characterize the dynamic behavior of IDPs, in particular residual dipolar couplings, that provide highly sensitive tools for the study of intrinsic structural propensity and conformational transitions accompanying protein function.

Indirect DNA readout by an H-NS related protein: structure of the DNA complex of the C-terminal domain of Ler

Ler, a member of the H-NS protein family, is the master regulator of the LEE pathogenicity island in virulent Escherichia coli strains. Here, we determined the structure of a complex between the DNA-binding domain of Ler (CT-Ler) and a 15-mer DNA duplex. CT-Ler recognizes a preexisting structural pattern in the DNA minor groove formed by two consecutive regions which are narrower and wider, respectively, compared with standard B-DNA. The compressed region, associated with an AT-tract, is sensed by the side chain of Arg90, whose mutation abolishes the capacity of Ler to bind DNA. The expanded groove allows the approach of the loop in which Arg90 is located. This is the first report of an experimental structure of a DNA complex that includes a protein belonging to the H-NS family. The indirect readout mechanism not only explains the capacity of H-NS and other H-NS family members to modulate the expression of a large number of genes but also the origin of the specificity displayed by Ler. Our results point to a general mechanism by which horizontally acquired genes may be specifically recognized by members of the H-NS family.

Pathogenic Escherichia coli strains and other enterobacteria carry genes acquired from other bacteria by a process known as horizontal gene transfer. Proper regulation of the genes that are expressed in a given moment is crucial for the success of the bacteria. The protein H-NS is a global regulator that binds DNA and maintains a large number of genes silent until they are required, for example, to sustain the bacteria's colonization of a new host. Ler is a member of the H-NS family that competes with H-NS to activate the expression of a group of horizontally acquired genes that encode for a molecular machine used by E. coli to infect human cells. Ler and H-NS share a similar DNA-binding domain and can bind to different DNA sequences. Here, we present the structure of a complex between the DNA-binding domain of Ler and a natural DNA fragment. This structure reveals that Ler recognizes specific DNA shapes, explaining its capacity to regulate genes with different sequences. A single arginine residue is key for the recognition of a DNA narrow minor groove, which is one of, though not the only, hallmarks of the DNA shapes that are recognized by H-NS and Ler.

Structural analysis of intrinsically disordered proteins by small-angle X-ray scattering

Small-angle scattering of X-rays (SAXS) is an established method to study the overall structure and structural transitions of biological macromolecules in solution. For folded proteins, the technique provides three-dimensional low resolution structures ab initio or it can be used to drive rigid-body modeling. SAXS is also a powerful tool for the quantitative analysis of flexible systems, including intrinsically disordered proteins (IDPs), and is highly complementary to the high resolution methods of X-ray crystallography and NMR. Here we present the basic principles of SAXS and review the main approaches to the characterization of IDPs and flexible multidomain proteins using SAXS. Together with the standard approaches based on the analysis of overall parameters, a recently developed Ensemble Optimization Method (EOM) is now available. The latter method allows for the co-existence of multiple protein conformations in solution compatible with the scattering data. Analysis of the selected ensembles provides quantitative information about flexibility and also offers insights into structural features. Examples of the use of SAXS and combined approaches with NMR, X-ray crystallography, and computational methods to characterize completely or partially disordered proteins are presented.

Protein loop compaction and the origin of the effect of arginine and glutamic acid mixtures on solubility, stability and transient oligomerization of proteins

Addition of a 50 mM mixture of L: -arginine and L: -glutamic acid (RE) is extensively used to improve protein solubility and stability, although the origin of the effect is not well understood. We present Small Angle X-ray Scattering (SAXS) and Nuclear Magnetic Resonance (NMR) results showing that RE induces protein compaction by collapsing flexible loops on the protein core. This is suggested to be a general mechanism preventing aggregation and improving resistance to proteases and to originate from the polyelectrolyte nature of RE. Molecular polyelectrolyte mixtures are expected to display long range correlation effects according to dressed interaction site theory. We hypothesize that perturbation of the RE solution by dissolved proteins is proportional to the volume occupied by the protein. As a consequence, loop collapse, minimizing the effective protein volume, is favored in the presence of RE.

Structural analysis of an equilibrium folding intermediate in the apoflavodoxin native ensemble by small-angle X-ray scattering

Intermediate conformations are crucial to our understanding of how proteins fold into their native structures and become functional. Conventional spectroscopic measurements of thermal denaturation transitions allow the detection of equilibrium intermediates but often provide little structural detail; thus, application of more informative techniques is required. Here we used small-angle X-ray scattering (SAXS) to study the thermal denaturation of four variants of Anabaena PCC 7119 flavodoxin, including the wild-type apo and holo forms, and two mutants, E20K/E72K and F98N. Denaturation was monitored from changes in SAXS descriptors. Although the starting and final points of the denaturation were similar for the flavodoxin variants tested, substantial differences in the unfolding pathway were apparent between them. In agreement with calorimetric data, analysis of the SAXS data sets indicated a three-state unfolding equilibrium for wild-type apoflavodoxin, a two-state equilibrium for the F98N mutant, and increased thermostability of the E20K/E72K mutant and holoflavodoxin. Although the apoflavodoxin intermediate consistently appeared mixed with significant amounts of either native or unfolded conformations, its SAXS profile was derived from the deconvolution of the temperature-dependent SAXS data set. The apoflavodoxin thermal intermediate was structurally close to the native state but less compact, thereby indicating incipient unfolding. The residues that foster denaturation were explored by an ensemble of equilibrium ϕ-value restrained molecular dynamics. These simulations pointed to residues located in the cofactor and partner-protein recognition regions as the initial sites of denaturation and suggest a conformational adaptation as the mechanism of action in apoflavodoxin.

Design and structure of an equilibrium protein folding intermediate: a hint into dynamical regions of proteins

Partly unfolded protein conformations close to the native state may play important roles in protein function and in protein misfolding. Structural analyses of such conformations which are essential for their fully physicochemical understanding are complicated by their characteristic low populations at equilibrium. We stabilize here with a single mutation the equilibrium intermediate of apoflavodoxin thermal unfolding and determine its solution structure by NMR. It consists of a large native region identical with that observed in the X-ray structure of the wild-type protein plus an unfolded region. Small-angle X-ray scattering analysis indicates that the calculated ensemble of structures is consistent with the actual degree of expansion of the intermediate. The unfolded region encompasses discontinuous sequence segments that cluster in the 3D structure of the native protein forming the FMN cofactor binding loops and the binding site of a variety of partner proteins. Analysis of the apoflavodoxin inner interfaces reveals that those becoming destabilized in the intermediate are more polar than other inner interfaces of the protein. Natively folded proteins contain hydrophobic cores formed by the packing of hydrophobic surfaces, while natively unfolded proteins are rich in polar residues. The structure of the apoflavodoxin thermal intermediate suggests that the regions of natively folded proteins that are easily responsive to thermal activation may contain cores of intermediate hydrophobicity.

A self-consistent description of the conformational behavior of chemically denatured proteins from NMR and small angle scattering

Characterization of the conformational properties of unfolded proteins is essential for understanding the mechanisms of protein folding and misfolding. This information is also fundamental to determining the relationship between flexibility and function in the highly diverse families of intrinsically disordered proteins. Here we present a self-consistent model of conformational sampling of chemically denatured proteins in agreement with experimental data reporting on long-range distance distributions in unfolded proteins using small-angle x-ray scattering and nuclear magnetic resonance pulse-field gradient-based measurements. We find that standard statistical coil models, selected from folded protein databases with secondary structural elements removed, need to be refined to correct backbone dihedral angle sampling of denatured proteins, although they appear to be appropriate for intrinsically disordered proteins. For denatured proteins, pervasive increases in the sampling of more-extended regions of Ramachandran space {50 degrees Collapse

Key Words Collapse

MESH Headings

• Computer Simulation
• Databases, Protein
• Models, Chemical
• Nuclear Magnetic Resonance, Biomolecular
• Protein Conformation
• Protein Denaturation
• Protein Folding
• Protein Structure, Secondary
• Scattering, Small Angle

Collapse

Grants Collapse

Weak oligomerization of low-molecular-weight protein tyrosine phosphatase is conserved from mammals to bacteria

The well-characterized self-association of a mammalian low-molecular-weight protein tyrosine phosphatase (lmwPTP) produces inactive oligomers that are in equilibrium with active monomers. A role of the inactive oligomers as supramolecular proenzymes has been suggested. The oligomerization equilibrium of YwlE, a lmwPTP from Bacillus subtilis, was studied by NMR. Chemical shift data and NMR relaxation confirm that dimerization takes place through the enzyme's active site, and is fully equivalent to the dimerization previously characterized in a eukaryotic low-molecular-weight phosphatase, with similarly large dissociation constants. The similarity between the oligomerization of prokaryotic and eukaryotic phosphatases extends beyond the dimer and involves higher order oligomers detected by NMR relaxation analysis at high protein concentrations. The conservation across different kingdoms of life suggests a physiological role for lmwPTP oligomerization in spite of the weak association observed in vitro. Structural data suggest that substrate modulation of the oligomerization equilibrium could be a regulatory mechanism leading to the generation of signaling pulses. The presence of a phenylalanine residue in the dimerization site of YwlE, replacing a tyrosine residue conserved in all eukaryotic lmwPTPs, demonstrates that lmwPTP regulation by oligomerization can be independent from tyrosine phosphorylation.

Low-resolution structures of transient protein-protein complexes using small-angle X-ray scattering

The determination of the three-dimensional structure of a weak protein-protein complex in solution using small-angle X-ray scattering requires the deconvolution of its contribution from those of other components coexisting in equilibrium. Using the oligomerization equilibrium of low molecular weight phosphatase (lmwPTP) as a model system, we show computationally and experimentally that the individual low-resolution structures of monomeric and dimeric lmwPTP can be determined from a small number of SAXS curves using the multivariate curve resolution with alternating least squares (MCR-ALS) algorithm. The dimeric complex represents no more than 15% of the macromolecules in the most concentrated sample. The derived structures are in good agreement with the crystallographic ones and the dissociation constant matches the one measured by NMR. These results demonstrate the power of SAXS, in combination with MCR-ALS, to study transient biomolecular complexes. The limits of the method were explored using a three-species model that describes the oligomerization of lmwPTP at higher concentrations.

Structural characterization of unphosphorylated STAT5a oligomerization equilibrium in solution by small-angle X-ray scattering

Signal transducer and activator of transcription (STAT) proteins play a crucial role in the activation of gene transcription in response to extracellular stimuli. The regulation and activity of these proteins require a complex rearrangement of the domains. According to the established models, based on crystallographic data, STATs convert from a basal antiparallel inactive dimer into a parallel active one following phosphorylation. The simultaneous analysis of small-angle X-ray scattering data measured at different concentrations of unphosphorylated human STAT5a core domain unambiguously identifies the simultaneous presence of a monomer and a dimer. The dimer is the minor species but could be structurally characterized by SAXS in the presence of the monomer using appropriate computational tools and shown to correspond to the antiparallel assembly. The equilibrium is governed by a moderate dissociation constant of K(d) approximately 90 microM. Integration of these results with previous knowledge of the N-terminal domain structure and dissociation constants allows the modeling of the full-length protein. A complex network of intermolecular interactions of low or medium affinity is suggested. These contacts can be eventually formed or broken to trigger the dramatic modifications in the dimeric arrangement needed for STAT regulation and activity.

Structural characterization of the natively unfolded N-terminal domain of human c-Src kinase: insights into the role of phosphorylation of the unique domain

The N-terminal regions of the members of Src family of non-receptor protein tyrosine kinases are intrinsically unfolded and contain the maximum sequence divergence among them. In this study, we have addressed the structural characterization by nuclear magnetic resonance of this region of 84 residues that encompasses the SH4 and the unique domains (USrc) of the human c-Src. With this aim, the backbone assignment was performed using (13)C-detected experiments that overcome the spectral resolution problems and the large number of prolines that are typical for intrinsically unfolded proteins. The analysis of the residual dipolar couplings measured for the USrc indicates the presence of a low populated helical structure in the 60-75 region. No long-range contacts between remote fragments of the chain were detected with paramagnetic relaxation enhancement experiments. The structural characterization was extended to two different phosphorylation states of USrc that encompassed three different phosphorylated sites, Ser17, Thr37, and Ser75. The structural and conformational changes upon phosphorylation were monitored through chemical shift perturbations and residual dipolar couplings, indicating that modifications occur at local level and no global rearrangements were apparent. These results suggest a scenario where phosphorylation induces a global electrostatic perturbation that could be involved in the membrane unbinding of c-Src and that could be related with the localization of the enzyme. These observations suggest the unique domain of Src kinases as a source of selectivity and reinforce the relevant role of intrinsically disordered proteins in biological processes.

ProtSA: a web application for calculating sequence specific protein solvent accessibilities in the unfolded ensemble

Background

The stability of proteins is governed by the heat capacity, enthalpy and entropy changes of folding, which are strongly correlated to the change in solvent accessible surface area experienced by the polypeptide. While the surface exposed in the folded state can be easily determined, accessibilities for the unfolded state at the atomic level cannot be obtained experimentally and are typically estimated using simplistic models of the unfolded ensemble. A web application providing realistic accessibilities of the unfolded ensemble of a given protein at the atomic level will prove useful.

Results

ProtSA, a web application that calculates sequence-specific solvent accessibilities of the unfolded state ensembles of proteins has been developed and made freely available to the scientific community. The input is the amino acid sequence of the protein of interest. ProtSA follows a previously published calculation protocol which uses the Flexible-Meccano algorithm to generate unfolded conformations representative of the unfolded ensemble of the protein, and uses the exact analytical software ALPHASURF to calculate atom solvent accessibilities, which are averaged on the ensemble.

Conclusion

ProtSA is a novel tool for the researcher investigating protein folding energetics. The sequence specific atom accessibilities provided by ProtSA will allow obtaining better estimates of the contribution of the hydrophobic effect to the free energy of folding, will help to refine existing parameterizations of protein folding energetics, and will be useful to understand the influence of point mutations on protein stability.