Low-resolution Structural Models Research Articles

The low-resolution structural models of hepatitis C virus RNA subdomain 5BSL3.2 and its distal complex with domain 3'X point to conserved regulatory mechanisms within the Flaviviridae family.

Subdomain 5BSL3.2 of hepatitis C virus RNA lies at the core of a network of distal RNA–RNA contacts that connect the 5′ and 3′ regions of the viral genome and regulate the translation and replication stages of the viral cycle. Using small-angle X-ray scattering and NMR spectroscopy experiments, we have determined at low resolution the structural models of this subdomain and its distal complex with domain 3′X, located at the 3′-terminus of the viral RNA chain. 5BSL3.2 adopts a characteristic ‘L’ shape in solution, whereas the 5BSL3.2–3′X distal complex forms a highly unusual ‘Y’-shaped kissing junction that blocks the dimer linkage sequence of domain 3′X and promotes translation. The structure of this complex may impede an effective association of the viral polymerase with 5BSL3.2 and 3′X to start negative-strand RNA synthesis, contributing to explain the likely mechanism used by these sequences to regulate viral replication and translation. In addition, sequence and shape features of 5BSL3.2 are present in functional RNA motifs of flaviviruses, suggesting conserved regulatory processes within the Flaviviridae family.

Deriving RNA topological structure from SAXS.

Structures of well-folded RNA molecules can be determined with atomic resolution by either X-ray crystallography, cryo-EM, or NMR spectroscopy, but those of conformationally-flexible RNAs often are difficult to study with these methods. However, flexible RNAs have biological relevance and likely represent the majority of the RNA conformational space. Due to the high electron density of the phosphate-sugar backbone, RNA is very sensitive to small-angle X-ray scattering (SAXS), and SAXS data can be recorded with sub-μM concentrations and under near-physiological solution conditions without the need for labeling. For these reasons, SAXS has significant advantages over other techniques for obtaining global structural information of flexible RNAs in the form of molecular envelopes or low-resolution topological structural models. The SAXS-derived information is extremely valuable for bridging secondary structure data, often determined by other techniques, with a three-dimensional structure description. In this chapter, we present a detailed account of the principle, algorithms, and experimental and computational protocols for topological structure determination of RNA molecules in solution. To illustrate the applications of the methodology, we provide several case studies that cover a broad spectrum of the RNA conformational landscape.

Predicting the Shapes of Protein Complexes through Collision Cross Section Measurements and Database Searches

In structural biology, collision cross sections (CCSs) from ion mobility mass spectrometry (IM-MS) measurements are routinely compared to computationally or experimentally derived protein structures. Here, we investigate whether CCS data can inform about the shape of a protein in the absence of specific reference structures. Analysis of the proteins in the CCS database shows that protein complexes with low apparent densities are structurally more diverse than those with a high apparent density. Although assigning protein shapes purely on CCS data is not possible, we find that we can distinguish oblate- and prolate-shaped protein complexes by using the CCS, molecular weight, and oligomeric states to mine the Protein Data Bank (PDB) for potentially similar protein structures. Furthermore, comparing the CCS of a ferritin cage to the solution structures in the PDB reveals significant deviations caused by structural collapse in the gas phase. We then apply the strategy to an integral membrane protein by comparing the shapes of a prokaryotic and a eukaryotic sodium/proton antiporter homologue. We conclude that mining the PDB with IM-MS data is a time-effective way to derive low-resolution structural models.

Structural properties of [2Fe-2S] ISCA2-IBA57: a complex of the mitochondrial iron-sulfur cluster assembly machinery

In mitochondria, a complex protein machinery is devoted to the maturation of iron-sulfur cluster proteins. Structural information on the last steps of the machinery, which involve ISCA1, ISCA2 and IBA57 proteins, needs to be acquired in order to define how these proteins cooperate each other. We report here the use of an integrative approach, utilizing information from small-angle X-ray scattering (SAXS) and bioinformatics-driven docking prediction, to determine a low-resolution structural model of the human mitochondrial [2Fe-2S]2+ ISCA2-IBA57 complex. In the applied experimental conditions, all the data converge to a structural organization of dimer of dimers for the [2Fe-2S]2+ ISCA2-IBA57 complex with ISCA2 providing the homodimerization core interface. The [2Fe-2S] cluster is out of the ISCA2 core while being shared with IBA57 in the dimer. The specific interaction pattern identified from the dimeric [2Fe-2S]2+ ISCA2-IBA57 structural model allowed us to define the molecular grounds of the pathogenic Arg146Trp mutation of IBA57. This finding suggests that the dimeric [2Fe-2S] ISCA2-IBA57 hetero-complex is a physiologically relevant species playing a role in mitochondrial [4Fe-4S] protein biogenesis.

Noninvasive Structural Analysis of Intermediate Species During Fibrillation: An Application of Small-Angle X-Ray Scattering.

Structural investigation of intermediately formed oligomers and pre-fibrillar species is of tremendous importance in order to elucidate the structural principles of fibrillation, and because intermediate species have been suggested as the pathogenic agents in several amyloid diseases. Structural investigations are however greatly complicated by the dynamic changesbetween structural states of very different sizes and life-times. Small angle X-ray scattering (SAXS) is an ideal method to handle this challenge. The method provides information about the fibrillation process (number of species present and their volume fractions) and low-resolution 3-dimensional structural models of individual species, notably also of the intermediately formed, in-process species from undisturbed fibrillation equilibria. Here, we provide a detailed description of the methods used for the measurement and analysis of SAXS data from fibrillating samples, exemplified using data from our own research.

Melting and freezing temperatures of confined Bi nanoparticles over a wide size range

The size dependences of the melting and freezing temperatures, T m and T f, respectively, of spherical Bi nanoparticles embedded in a sodium borate glass were determined by applying a new experimental procedure based on the combined and simultaneous use of small-angle X-ray scattering (SAXS) and wide-angle X-ray scattering (WAXS). This experimental procedure is particularly useful for materials in which a widely polydisperse set of nanoparticles are embedded. The results provide additional and stronger evidence supporting the main previous conclusions: (i) the melting and freezing temperatures both decrease linearly for increasing reciprocal radius (1/R); and (ii) the effect of undercooling is suppressed for Bi nanoparticles with radii smaller than a critical value equal to 1.8 nm. These results confirm a previously proposed low-resolution structural model for Bi nanocrystals below their melting temperature and with radius R > 1.8 nm, which consists of a crystalline core surrounded by a disordered shell. In the present work, a number of samples with different and partially overlapping radius distributions were studied, allowing the determination of T m(R) and T f(R) functions over a wide range of radii (1 < R < 11 nm). Comparison of the experimentally determined T m(R) and T f(R) functions corresponding to different samples indicates good reproducibility of the experimental results. This allowed the verification of the robustness of the experimental procedure based on in situ combined use of SAXS and WAXS for determination of the radius dependence of the melting and freezing temperatures of spherical nanoparticles in dilute solution.

COFACTOR: improved protein function prediction by combining structure, sequence and protein-protein interaction information.

The COFACTOR web server is a unified platform for structure-based multiple-level protein function predictions. By structurally threading low-resolution structural models through the BioLiP library, the COFACTOR server infers three categories of protein functions including gene ontology, enzyme commission and ligand-binding sites from various analogous and homologous function templates. Here, we report recent improvements of the COFACTOR server in the development of new pipelines to infer functional insights from sequence profile alignments and protein–protein interaction networks. Large-scale benchmark tests show that the new hybrid COFACTOR approach significantly improves the function annotation accuracy of the former structure-based pipeline and other state-of-the-art functional annotation methods, particularly for targets that have no close homology templates. The updated COFACTOR server and the template libraries are available at http://zhanglab.ccmb.med.umich.edu/COFACTOR/.

IDPs: Less Disordered and More Ordered than Expected

In this issue of the Biophysical Journal, Solyom et al. (9) provide a telling example of structural preformation in intrinsically disordered proteins (IDPs) and the relevance of distinctly folded and compact conformational substates for molecular recognition processes.

Human mitochondrial Hsp70 (mortalin): shedding light on ATPase activity, interaction with adenosine nucleotides, solution structure and domain organization.

The human mitochondrial Hsp70, also called mortalin, is of considerable importance for mitochondria biogenesis and the correct functioning of the cell machinery. In the mitochondrial matrix, mortalin acts in the importing and folding process of nucleus-encoded proteins. The in vivo deregulation of mortalin expression and/or function has been correlated with age-related diseases and certain cancers due to its interaction with the p53 protein. In spite of its critical biological roles, structural and functional studies on mortalin are limited by its insoluble recombinant production. This study provides the first report of the production of folded and soluble recombinant mortalin when co-expressed with the human Hsp70-escort protein 1, but it is still likely prone to self-association. The monomeric fraction of mortalin presented a slightly elongated shape and basal ATPase activity that is higher than that of its cytoplasmic counterpart Hsp70-1A, suggesting that it was obtained in the functional state. Through small angle X-ray scattering, we assessed the low-resolution structural model of monomeric mortalin that is characterized by an elongated shape. This model adequately accommodated high resolution structures of Hsp70 domains indicating its quality. We also observed that mortalin interacts with adenosine nucleotides with high affinity. Thermally induced unfolding experiments indicated that mortalin is formed by at least two domains and that the transition is sensitive to the presence of adenosine nucleotides and that this process is dependent on the presence of Mg2+ ions. Interestingly, the thermal-induced unfolding assays of mortalin suggested the presence of an aggregation/association event, which was not observed for human Hsp70-1A, and this finding may explain its natural tendency for in vivo aggregation. Our study may contribute to the structural understanding of mortalin as well as to contribute for its recombinant production for antitumor compound screenings.

Protein structure refinement of CASP target proteins using GNEIMO torsional dynamics method.

A longstanding challenge in using computational methods for protein structure prediction is the refinement of low-resolution structural models derived from comparative modeling methods into highly accurate atomistic models useful for detailed structural studies. Previously, we have developed and demonstrated the utility of the internal coordinate molecular dynamics (MD) technique, generalized Newton–Euler inverse mass operator (GNEIMO), for refinement of small proteins. Using GNEIMO, the high-frequency degrees of freedom are frozen and the protein is modeled as a collection of rigid clusters connected by torsional hinges. This physical model allows larger integration time steps and focuses the conformational search in the low frequency torsional degrees of freedom. Here, we have applied GNEIMO with temperature replica exchange to refine low-resolution protein models of 30 proteins taken from the continuous assessment of structure prediction (CASP) competition. We have shown that GNEIMO torsional MD method leads to refinement of up to 1.3 Å in the root-mean-square deviation in coordinates for 30 CASP target proteins without using any experimental data as restraints in performing the GNEIMO simulations. This is in contrast with the unconstrained all-atom Cartesian MD method performed under the same conditions, where refinement requires the use of restraints during the simulations.

Structural analysis of influenza A virus matrix protein M1 and its self-assemblies at low pH.

Influenza A virus matrix protein M1 is one of the most important and abundant proteins in the virus particles broadly involved in essential processes of the viral life cycle. The absence of high-resolution data on the full-length M1 makes the structural investigation of the intact protein particularly important. We employed synchrotron small-angle X-ray scattering (SAXS), analytical ultracentrifugation and atomic force microscopy (AFM) to study the structure of M1 at acidic pH. The low-resolution structural models built from the SAXS data reveal a structurally anisotropic M1 molecule consisting of a compact NM-fragment and an extended and partially flexible C-terminal domain. The M1 monomers co-exist in solution with a small fraction of large clusters that have a layered architecture similar to that observed in the authentic influenza virions. AFM analysis on a lipid-like negatively charged surface reveals that M1 forms ordered stripes correlating well with the clusters observed by SAXS. The free NM-domain is monomeric in acidic solution with the overall structure similar to that observed in previously determined crystal structures. The NM-domain does not spontaneously self assemble supporting the key role of the C-terminus of M1 in the formation of supramolecular structures. Our results suggest that the flexibility of the C-terminus is an essential feature, which may be responsible for the multi-functionality of the entire protein. In particular, this flexibility could allow M1 to structurally organise the viral membrane to maintain the integrity and the shape of the intact influenza virus.

Characterization of the human dynein light chain Rp3 and its use as a non-viral gene delivery vector

Dynein light chains mediate the interaction between the cargo and the dynein motor complex during retrograde microtubule-mediated transport in eukaryotic cells. In this study, we expressed and characterized the recombinant human dynein light chain Rp3 and developed a modified variant harboring an N-terminal DNA-binding domain (Rp3-Db). Our approach aimed to explore the retrograde cell machinery based on dynein to enhance plasmid DNA (pDNA) traffic along the cytosol toward the nucleus. In the context of non-viral gene delivery, Rp3-Db is expected to simultaneously interact with DNA and dynein, thereby enabling a more rapid and efficient transport of the genetic material across the cytoplasm. We successfully purified recombinant Rp3 and obtained a low-resolution structural model using small-angle X-ray scattering. Additionally, we observed that Rp3 is a homodimer under reducing conditions and remains stable over a broad pH range. The ability of Rp3 to interact with the dynein intermediate chain in vitro was also observed, indicating that the recombinant Rp3 is correctly folded and functional. Finally, Rp3-Db was successfully expressed and purified and exhibited the ability to interact with pDNA and mediate the transfection of cultured HeLa cells. Rp3-Db was also capable of interacting in vitro with dynein intermediate chains, indicating that the addition of the N-terminal DNA-binding domain does not compromise its function. The transfection level observed for Rp3-Db is far superior than that reported for protamine and is comparable to that of the cationic lipid Lipofectamine™. This report presents an initial characterization of a non-viral delivery vector based on the dynein light chain Rp3 and demonstrates the potential use of modified human light chains as gene delivery vectors.

Structural Studies of the Apo and Ca2+-Bound States of the Human BK (SLO1) Channel Gating Ring in Solution

The gating ring (GR) regulates the activity of large-conductance voltage- and Ca^(2+)-activated K^+ channels (BK) by interacting with intracellular signaling molecules. To understand the operation of this biological sensor under physiological conditions, we performed Small-Angle X-ray Scattering (SAXS) analysis, at beamline 4-2 at the Stanford Synchrotron Radiation laboratory. SAXS measurements of the purified GR were performed in the absence or in the presence of 35 μM free Ca^(2+), found to be a saturating concentration in previous work. The quality of the circularly-averaged scattering data was evaluated with Guinier analysis, while the ATSAS software suite was used to derive structural information. The radius of gyration (R_g) and maximum interparticle distance (D_(max)) of the apo GR were 48.65±1.372 A and 185 A, respectively. These values are comparable to data obtained from crystal structure of GR (3NAF), where the envelope R_g, calculated with CRYSOL, is 45.55 A, and its diameter 155.6 A. Ca^(2+)-bound GR shows a decrease in R_g to 42.77±1.058 A and D_(max) to 160 A, demonstrating the structural response of GR to Ca^(2+). Low-resolution structural models of the GR were generated from the experimental scattering pattern using DAMMIN. The Ca^(2+)-bound GR revealed notable changes in both flexible and assembly interfaces of the superstructure's constituent RCK1 (Regulator of Conductance for K^+) and RCK2 domains. Since the structural changes are resolved under physiologically-relevant conditions, we speculate that they represent the molecular transitions that initiate the Ca^(2+)-induced activation of human BK channels.

Probing subunit-subunit interactions in the yeast vacuolar ATPase by peptide arrays.

BackgroundVacuolar (H+)-ATPase (V-ATPase; V1Vo-ATPase) is a large multisubunit enzyme complex found in the endomembrane system of all eukaryotic cells where its proton pumping action serves to acidify subcellular organelles. In the plasma membrane of certain specialized tissues, V-ATPase functions to pump protons from the cytoplasm into the extracellular space. The activity of the V-ATPase is regulated by a reversible dissociation mechanism that involves breaking and re-forming of protein-protein interactions in the V1-ATPase - Vo-proton channel interface. The mechanism responsible for regulated V-ATPase dissociation is poorly understood, largely due to a lack of detailed knowledge of the molecular interactions that are responsible for the structural and functional link between the soluble ATPase and membrane bound proton channel domains.Methodology/Principal FindingsTo gain insight into where some of the stator subunits of the V-ATPase associate with each other, we have developed peptide arrays from the primary sequences of V-ATPase subunits. By probing the peptide arrays with individually expressed V-ATPase subunits, we have identified several key interactions involving stator subunits E, G, C, H and the N-terminal domain of the membrane bound a subunit.ConclusionsThe subunit-peptide interactions identified from the peptide arrays complement low resolution structural models of the eukaryotic vacuolar ATPase obtained from transmission electron microscopy. The subunit-subunit interaction data are discussed in context of our current model of reversible enzyme dissociation.

Dynamic and Thermodynamic Signatures of Native and Non-Native Protein States with Application to the Improvement of Protein Structures.

Accurate knowledge of the 3D structural ensemble of proteins is important for understanding of their biological function. We report here the application of microsecond all-atom molecular dynamics (MD) simulations in explicit solvent for the improvement of the quality of low-resolution structures obtained by protein structure prediction (decoys). Seventy MD simulations of ∼1 μs average duration were performed on 13 different protein systems starting from X-ray crystal structures and decoys. Their behavior can be divided into three groups: 22 trajectories converged toward the native state, 27 trajectories displayed a quasi-equilibrium by populating mainly a single non-native free energy basin, and 21 trajectories drifted away from their initial decoy structure transiently visiting multiple free energy minima. To determine whether the native structure can be identified among non-native ensembles, the free energy was determined for each basin by the MM/GBSA method together with the von Mises entropy estimator in dihedral angle space. For the proteins studied here, it is found that the ensembles belonging to free energy basins with the lowest free energies and the longest residence times are most native-like. The results demonstrate that explicit solvent microsecond MD simulations using the latest generation of protein force fields and free energy metrics are sufficiently accurate to permit positive identification of native state ensembles against low-resolution structural models and decoys. The approach can be applied to the direct refinement of predicted or experimental low-resolution protein structures.

Recognizing Protein-Ligand Binding Sites by Global Structural Alignment and Local Geometry Refinement

Proteins perform functions through interacting with other molecules. However, structural details for most of the protein-ligand interactions are unknown. We present a comparative approach (COFACTOR) to recognize functional sites of protein-ligand interactions using low-resolution protein structural models, based on a global-to-local sequence and structural comparison algorithm. COFACTOR was tested on 501 proteins, which harbor 582 natural and drug-like ligand molecules. Starting from I-TASSER structure predictions, the method successfully identifies ligand-binding pocket locations for 65% of apo receptors with an average distance error 2Å. The average precision of binding-residue assignments is 46% and 137% higher than that by FINDSITE and ConCavity. In CASP9, COFACTOR achieveda binding-site prediction precision 72% and Matthews correlation coefficient 0.69 for 31 blind test proteins, which was significantly higher than all other participating methods. These data demonstrate the power of structure-based approaches to protein-ligand interaction predictions applicable for genome-wide structural and functional annotations.

Structural and functional analysis of the symmetrical Type I restriction endonuclease R.EcoR124INT

Type I restriction-modification (RM) systems are comprised of two multi-subunit enzymes, the methyltransferase (∼160 kDa), responsible for methylation of DNA, and the restriction endonuclease (∼400 kDa), responsible for DNA cleavage. Both enzymes share a number of subunits. An engineered RM system, EcoR124INT, based on the N-terminal domain of the specificity subunit of EcoR124I was constructed that recognises the symmetrical sequence GAAN7TTC and is active as a methyltransferase. Here, we investigate the restriction endonuclease activity of R. EcoR124INT in vitro and the subunit assembly of the multi-subunit enzyme. Finally, using small-angle neutron scattering and selective deuteration, we present a low-resolution structural model of the endonuclease and locate the motor subunits within the multi-subunit enzyme. We show that the covalent linkage between the two target recognition domains of the specificity subunit is not required for subunit assembly or enzyme activity, and discuss the implications for the evolution of Type I enzymes.

Validation of Membrane Protein Topology Models by Oxidative Labeling and Mass Spectrometry

Computer-assisted topology predictions are widely used to build low-resolution structural models of integral membrane proteins (IMPs). Experimental validation of these models by traditional methods is labor intensive and requires modifications that might alter the IMP native conformation. This work employs oxidative labeling coupled with mass spectrometry (MS) as a validation tool for computer-generated topology models. ·OH exposure introduces oxidative modifications in solvent-accessible regions, whereas buried segments (e.g., transmembrane helices) are non-oxidizable. The Escherichia coli protein WaaL (O-antigen ligase) is predicted to have 12 transmembrane helices and a large extramembrane domain (Pérez et al., Mol. Microbiol. 2008, 70, 1424). Tryptic digestion and LC-MS/MS were used to map the oxidative labeling behavior of WaaL. Met and Cys exhibit high intrinsic reactivities with ·OH, making them sensitive probes for solvent accessibility assays. Overall, the oxidation pattern of these residues is consistent with the originally proposed WaaL topology. One residue (M151), however, undergoes partial oxidation despite being predicted to reside within a transmembrane helix. Using an improved computer algorithm, a slightly modified topology model was generated that places M151 closer to the membrane interface. On the basis of the labeling data, it is concluded that the refined model more accurately reflects the actual topology of WaaL. We propose that the combination of oxidative labeling and MS represents a useful strategy for assessing the accuracy of IMP topology predictions, supplementing data obtained in traditional biochemical assays. In the future, it might be possible to incorporate oxidative labeling data directly as constraints in topology prediction algorithms.

Statistical mechanical modeling of RNA folding: from free energy landscape to tertiary structural prediction.

In spite of the success of computational methods for predicting RNA secondary structure, the problem of predicting RNA tertiary structure folding remains. Low-resolution structural models show promise as they allow for rigorous statistical mechanical computation for the conformational entropies, free energies, and the coarse-grained structures of tertiary folds. Molecular dynamics refinement of coarse-grained structures leads to all-atom 3D structures. Modeling based on statistical mechanics principles also has the unique advantage of predicting the full free energy landscape, including local minima and the global free energy minimum. The energy landscapes combined with the 3D structures form the basis for quantitative predictions of RNA functions. In this chapter, we present an overview of statistical mechanical models for RNA folding and then focus on a recently developed RNA statistical mechanical model -- the Vfold model. The main emphasis is placed on the physics underpinning the models, the computational strategies, and the connections to RNA biology.

Predicting nucleic acid binding interfaces from structural models of proteins

The function of DNA- and RNA-binding proteins can be inferred from the characterization and accurate prediction of their binding interfaces. However, the main pitfall of various structure-based methods for predicting nucleic acid binding function is that they are all limited to a relatively small number of proteins for which high-resolution three-dimensional structures are available. In this study, we developed a pipeline for extracting functional electrostatic patches from surfaces of protein structural models, obtained using the I-TASSER protein structure predictor. The largest positive patches are extracted from the protein surface using the patchfinder algorithm. We show that functional electrostatic patches extracted from an ensemble of structural models highly overlap the patches extracted from high-resolution structures. Furthermore, by testing our pipeline on a set of 55 known nucleic acid binding proteins for which I-TASSER produces high-quality models, we show that the method accurately identifies the nucleic acids binding interface on structural models of proteins. Employing a combined patch approach we show that patches extracted from an ensemble of models better predicts the real nucleic acid binding interfaces compared with patches extracted from independent models. Overall, these results suggest that combining information from a collection of low-resolution structural models could be a valuable approach for functional annotation. We suggest that our method will be further applicable for predicting other functional surfaces of proteins with unknown structure.