The Ligand-Free State of the TPP Riboswitch: A Partially Folded RNA Structure

Riboswitches are elements of mRNA that regulate gene expression by undergoing structural changes upon binding of small ligands. Although the structures of several riboswitches have been solved with their ligands bound, the ligand-free states of only a few riboswitches have been characterized. The ligand-free state is as important for the functionality of the riboswitch as the ligand-bound form, but the ligand-free state is often a partially folded structure of the RNA, with conformational heterogeneity that makes it particularly challenging to study. Here, we present models of the ligand-free state of a thiamine pyrophosphate riboswitch that are derived from a combination of complementary experimental and computational modeling approaches. We obtain a global picture of the molecule using small-angle X-ray scattering data and use an RNA structure modeling software, MC-Sym, to fit local structural details to these data on an atomic scale. We have used two different approaches to obtaining these models. Our first approach develops a model of the RNA from the structures of its constituent junction fragments in isolation. The second approach treats the RNA as a single entity, without bias from the structure of its individual constituents. We find that both approaches give similar models for the ligand-free form, but the ligand-bound models differ for the two approaches, and only the models from the second approach agree with the ligand-bound structure known previously from X-ray crystallography. Our models provide a picture of the conformational changes that may occur in the riboswitch upon binding of its ligand. Our results also demonstrate the power of combining experimental small-angle X-ray scattering data with theoretical structure prediction tools in the determination of RNA structures beyond riboswitches.     

Conformational interdomain flexibility in a bacterial α-isopropylmalate synthase is necessary for leucine biosynthesis

α-Isopropylmalate synthase (IPMS) catalyzes the first step in leucine (Leu) biosynthesis and is allosterically regulated by the pathway end product, Leu. IPMS is a dimeric enzyme with each chain consisting of catalytic, accessory, and regulatory domains, with the accessory and regulatory domains of each chain sitting adjacent to the catalytic domain of the other chain. The IPMS crystal structure shows significant asymmetry because of different relative domain conformations in each chain. Owing to the challenges posed by the dynamic and asymmetric structures of IPMS enzymes, the molecular details of their catalytic and allosteric mechanisms are not fully understood. In this study, we have investigated the allosteric feedback mechanism of the IPMS enzyme from the bacterium that causes meningitis, Neisseria meningitidis (NmeIPMS). By combining molecular dynamics simulations with small-angle X-ray scattering, mutagenesis, and heterodimer generation, we demonstrate that Leu-bound NmeIPMS is in a rigid conformational state stabilized by asymmetric interdomain polar interactions. Furthermore, we found removing these polar interactions by mutagenesis impaired the allosteric response without compromising Leu binding. Our results suggest that the allosteric inhibition of NmeIPMS is achieved by restricting the flexibility of the accessory and regulatory domains, demonstrating that significant conformational flexibility is required for catalysis.     

Conformational differences among solution structures of the type Ialpha, IIalpha and IIbeta protein kinase A regulatory subunit homodimers: role of the linker regions

The regulatory (R) subunits of the cAMP-dependent protein kinase (protein kinase A or PKA) are multi-domain proteins responsible for conferring cAMP-dependence and localizing PKA to specific subcellular locations. There are four isoforms of the R subunit in mammals that are similar in molecular mass and domain organization, but clearly serve different biological functions. Although high-resolution structures are available for the cAMP-binding domains and dimerization/docking domains of two isoforms, there are no high-resolution structures of any of the intact R subunit homodimer isoforms. The results of small-angle X-ray scattering studies presented here indicate that the RIalpha, RIIalpha, and RIIbeta homodimers differ markedly in overall shape, despite extensive sequence homology and similar molecular masses. The RIIalpha and RIIbeta homodimers have very extended, rod-like shapes, whereas the RIalpha homodimer likely has a compact Y-shape. Based on a comparison of the R subunit sequences, we predict that the linker regions are the likely cause of these large differences in shape among the isoforms. In addition, we show that cAMP binding does not cause large conformational changes in type Ialpha or IIalpha R subunit homodimers, suggesting that the activation of PKA by cAMP involves only local conformational changes in the R subunits.     

Structural Characterization of Protein-Protein Complexes by Integrating Computational Docking with Small-angle Scattering Data

X-ray crystallography and NMR can provide detailed structural information of protein-protein complexes, but technical problems make their application challenging in the high-throughput regime. Other methods such as small-angle X-ray scattering (SAXS) are more promising for large-scale application, but at the cost of lower resolution, which is a problem that can be solved by complementing SAXS data with theoretical simulations. Here, we propose a novel strategy that combines SAXS data and accurate protein-protein docking simulations. The approach has been benchmarked on a large pool of known structures with synthetic SAXS data, and on three experimental examples. The combined approach (pyDockSAXS) provided a significantly better success rate (43% for the top 10 predictions) than either of the two methods alone. Further analysis of the influence of different docking parameters made it possible to increase the success rates for specific cases, and to define guidelines for improving the data-driven protein-protein docking protocols.     

Extended conformational states dominate the Hsp90 chaperone dynamics

The heat shock protein 90 (Hsp90) is a molecular chaperone central to client protein folding and maturation in eukaryotic cells. During its chaperone cycle, Hsp90 undergoes ATPase-coupled large-scale conformational changes between open and closed states, where the N-terminal and middle domains of the protein form a compact dimerized conformation. However, the molecular principles of the switching motion between the open and closed states remain poorly understood. Here we show by integrating atomistic and coarse-grained molecular simulations with small-angle X-ray scattering experiments and NMR spectroscopy data that Hsp90 exhibits rich conformational dynamics modulated by the charged linker, which connects the N-terminal with the middle domain of the protein. We show that the dissociation of these domains is crucial for the conformational flexibility of the open state, with the separation distance controlled by a β-sheet motif next to the linker region. Taken together, our results suggest that the conformational ensemble of Hsp90 comprises highly extended states, which could be functionally crucial for client processing.     

The evidence of large-scale DNA-induced compaction in the mycobacterial chromosomal ParB

The bacterial chromosome trafficking apparatus or the segrosome participates in the mitotic-like segregation of the chromosomes prior to cell division in several bacteria. ParB, which is the parS DNA-binding component of the segrosome, polymerizes on the parS-adjacent chromosome to form a nucleoprotein filament of unknown nature for the segregation function. We combined static light scattering, circular dichroism and small-angle X-ray scattering to present evidence that the apo form of the mycobacterial ParB forms an elongated dimer with intrinsically disordered regions as well as folded domains in solution. A comparison of the solution scattering of the apo and the parS-bound ParBs indicates a rather drastic compaction of the protein upon DNA binding. We propose that this binding-induced conformational transition is priming the ParB for polymerization on the DNA template.     

Molecular Morphology of Eukaryotic Class-1 Translation Termination Factor eRF1 in Solution

The integral structural parameters and the shape of the molecule of human translation termination factor eRF1 were determined from the small-angle X-ray scattering in solution. The molecular shapes were found by bead modeling with nonlinear minimization of the root-mean-square deviation of the calculated from the experimental scattering curve. Comparisons of the small-angle scattering curves computed for atomic-resolution structures of eRF1 with the experimental data on scattering from solution testified that the crystal and the solution conformations are close. In the ribosome, the distance between the eRF1 motifs GGQ and NIKS must be shorter than in crystal or solution (75 versus 100–107 Å). Therefore, like its bacterial counterpart RF2, the eukaryotic eRF1 must change its conformation as it binds to the ribosome. The conformational mobility of eukaryotic and prokaryotic class-1 release factors is another feature making them functionally akin to tRNA.     

Improving small-angle X-ray scattering data for structural analyses of the RNA world

Defining the shape, conformation, or assembly state of an RNA in solution often requires multiple investigative tools ranging from nucleotide analog interference mapping to X-ray crystallography. A key addition to this toolbox is small-angle X-ray scattering (SAXS). SAXS provides direct structural information regarding the size, shape, and flexibility of the particle in solution and has proven powerful for analyses of RNA structures with minimal requirements for sample concentration and volumes. In principle, SAXS can provide reliable data on small and large RNA molecules. In practice, SAXS investigations of RNA samples can show inconsistencies that suggest limitations in the SAXS experimental analyses or problems with the samples. Here, we show through investigations on the SAM-I riboswitch, the Group I intron P4-P6 domain, 30S ribosomal subunit from Sulfolobus solfataricus (30S), brome mosaic virus tRNA-like structure (BMV TLS), Thermotoga maritima asd lysine riboswitch, the recombinant tRNA^val, and yeast tRNA^phe that many problems with SAXS experiments on RNA samples derive from heterogeneity of the folded RNA. Furthermore, we propose and test a general approach to reducing these sample limitations for accurate SAXS analyses of RNA. Together our method and results show that SAXS with synchrotron radiation has great potential to provide accurate RNA shapes, conformations, and assembly states in solution that inform RNA biological functions in fundamental ways.     

Dissection of the ATP-Dependent Conformational Change Cycle of a Group II Chaperonin

Group II chaperonin captures an unfolded protein while in its open conformation and then mediates the folding of the protein during ATP-driven conformational change cycle. In this study, we performed kinetic analyses of the group II chaperonin from a hyperthermophilic archaeon, Thermococcus sp. KS-1 (TKS1-Cpn), by stopped-flow fluorometry and stopped-flow small-angle X-ray scattering to reveal the reaction cycle. Two TKS1-Cpn variants containing a Trp residue at position 265 or position 56 exhibit nearly the same fluorescence kinetics induced by rapid mixing with ATP. Fluorescence started to increase immediately after the start of mixing and reached a maximum at 1–2 s after mixing. Only in the presence of K⁺ that a gradual decrease in fluorescence was observed after the initial peak. Similar results were obtained by stopped-flow small-angle X-ray scattering. A rapid fluorescence increase, which reflects nucleotide binding, was observed for the mutant containing a Trp residue near the ATP binding site (K485W), irrespective of the presence or absence of K⁺. Without K⁺, a small, rapid fluorescence decrease followed the initial increase, and then a gradual decrease was observed. In contrast, with K⁺, a large, rapid fluorescence decrease occurred just after the initial increase, and then the fluorescence gradually increased. Finally, we observed ATP binding signal and also subtle conformational change in an ATPase-deficient mutant with K485W mutation. Based on these results, we propose a reaction cycle model for group II chaperonins.     

Early collapse is not an obligate step in protein folding

The dimensions and secondary structure content of two proteins which fold in a two-state manner are measured within milliseconds of denaturant dilution using synchrotron-based, stopped-flow small-angle X-ray scattering and far-UV circular dichroism spectroscopy. Even upon a jump to strongly native conditions, neither ubiquitin nor common-type acylphosphatase contract prior to the major folding event. Circular dichroism and fluorescence indicate that negligible amounts of secondary and tertiary structures form in the burst phase. Thus, for these two denatured states, collapse and secondary structure formation are not energetically downhill processes even under aqueous, low-denaturant conditions. In addition, water appears to be as good a solvent as that with high concentrations of denaturant, when considering the over-all dimensions of the denatured state. However, the removal of denaturant does subtly alter the distribution of backbone dihedral phi,psi angles, most likely resulting in a shift from the polyproline II region to the helical region of the Ramachandran map. We consider the thermodynamic origins of these behaviors along with implications for folding mechanisms and computer simulations thereof.     

Protein kinase A targeting and activation as seen by small-angle solution scattering

We have studied the solution structures of the multi-functional protein kinase A using small-angle X-ray and neutron scattering and have found a remarkable structural diversity in the different isoforms of this multi-subunit enzyme, in spite of its having high sequence homology and a common domain organization within its sequences. The available high-resolution crystal and NMR structural data for the protein kinase A components have aided in the interpretation of the solution scattering data and enabled us to develop models that bring insights into protein kinase A activation and targeting mechanisms, such as the opening and closing of the catalytic cleft to facilitate substrate binding or inhibition, respectively, and the role of sequence segments that join functional domains in the R subunit in providing a structurally flexible scaffold for interactions with the C subunit and A kinase-anchoring proteins (AKAPs).     

Streptokinase is a flexible multi-domain protein

The structure of streptokinase in solution has been studied by dynamic light scattering, small-angle X-ray scattering and circular dichroism spectroscopy. The Stokes' radius and radius of gyration of the protein monomer are 3.58 nm and 4.03 nm, respectively. The maximum intraparticle distance of the molecule is 14 nm. More than half of the amino acids of the molecule are organized in regular secondary structures. The X-ray scattering curve, the results from dynamic light scattering, and the finding that at least 50% of the amino acid residues are organized in regularly folded secondary structures are consistent with the following structural model. Streptokinase consists of four compact, separately folded, domains linked by mobile segments of the protein chain. The molecule exhibits the conformation of a flexible string-of-beads in solution.     

Modulation of the Pyrococcus abyssi NucS endonuclease activity by replication clamp at functional and structural levels

Pyrococcus abyssi NucS is the founding member of a new family of structure-specific DNA endonucleases that interact with the replication clamp proliferating cell nuclear antigen (PCNA). Using a combination of small angle x-ray scattering and surface plasmon resonance analyses, we demonstrate the formation of a stable complex in solution, in which one molecule of the PabNucS homodimer binds to the outside surface of the PabPCNA homotrimer. Using fluorescent labels, PCNA is shown to increase the binding affinity of NucS toward single-strand/double-strand junctions on 5' and 3' flaps, as well as to modulate the cleavage specificity on the branched DNA structures. Our results indicate that the presence of a single major contact between the PabNucS and PabPCNA proteins, together with the complex-induced DNA bending, facilitate conformational flexibility required for specific cleavage at the single-strand/double-strand DNA junction.     

The active conformation of glutamate synthase and its binding to ferredoxin

Glutamate synthases (GltS) are crucial enzymes in ammonia assimilation in plants and bacteria, where they catalyze the formation of two molecules of L-glutamate from L-glutamine and 2-oxoglutarate. The plant-type ferredoxin-dependent GltS and the functionally homologous alpha subunit of the bacterial NADPH-dependent GltS are complex four-domain monomeric enzymes of 140-165 kDa belonging to the NH(2)-terminal nucleophile family of amidotransferases. The enzymes function through the channeling of ammonia from the N-terminal amidotransferase domain to the FMN-binding domain. Here, we report the X-ray structure of the Synechocystis ferredoxin-dependent GltS with the substrate 2-oxoglutarate and the covalent inhibitor 5-oxo-L-norleucine bound in their physically distinct active sites solved using a new crystal form. The covalent Cys1-5-oxo-L-norleucine adduct mimics the glutamyl-thioester intermediate formed during L-glutamine hydrolysis. Moreover, we determined a high resolution structure of the GltS:2-oxoglutarate complex. These structures represent the enzyme in the active conformation. By comparing these structures with that of GltS alpha subunit and of related enzymes we propose a mechanism for enzyme self-regulation and ammonia channeling between the active sites. X-ray small-angle scattering experiments were performed on solutions containing GltS and its physiological electron donor ferredoxin (Fd). Using the structure of GltS and the newly determined crystal structure of Synechocystis Fd, the scattering experiments clearly showed that GltS forms an equimolar (1:1) complex with Fd. A fundamental consequence of this result is that two Fd molecules bind consecutively to Fd-GltS to yield the reduced FMN cofactor during catalysis.     

Structural characterization of intramolecular Hg(2+) transfer between flexibly linked domains of mercuric ion reductase

The enzyme mercuric ion reductase MerA is the central component of bacterial mercury resistance encoded by the mer operon. Many MerA proteins possess metallochaperone-like N-terminal domains (NmerA) that can transfer Hg²⁺ to the catalytic core domain (Core) for reduction to Hg⁰. These domains are tethered to the homodimeric Core by ∼ 30-residue linkers that are susceptible to proteolysis, the latter of which has prevented characterization of the interactions of NmerA and the Core in the full-length protein. Here, we report purification of homogeneous full-length MerA from the Tn21 mer operon using a fusion protein construct and combine small-angle X-ray scattering and small-angle neutron scattering with molecular dynamics simulation to characterize the structures of full-length wild-type and mutant MerA proteins that mimic the system before and during handoff of Hg²⁺ from NmerA to the Core. The radii of gyration, distance distribution functions, and Kratky plots derived from the small-angle X-ray scattering data are consistent with full-length MerA adopting elongated conformations as a result of flexibility in the linkers to the NmerA domains. The scattering profiles are best reproduced using an ensemble of linker conformations. This flexible attachment of NmerA may facilitate fast and efficient removal of Hg²⁺ from diverse protein substrates. Using a specific mutant of MerA allowed the formation of a metal-mediated interaction between NmerA and the Core and the determination of the position and relative orientation of NmerA to the Core during Hg²⁺ handoff.     

Determination of the Individual Roles of the Linker Residues in the Interdomain Motions of Calmodulin Using NMR Chemical Shifts

Many protein molecules are formed by two or more domains whose structures and dynamics are closely related to their biological functions. It is thus important to develop methods to determine the structural properties of these multidomain proteins. Here, we characterize the interdomain motions in the calcium-bound state of calmodulin (Ca^2 +-CaM) using NMR chemical shifts as replica-averaged structural restraints in molecular dynamics simulations. We find that the conformational fluctuations of the interdomain linker, which are largely responsible for the overall interdomain motions of CaM, can be well described by exploiting the information provided by chemical shifts. We thus identify 10 residues in the interdomain linker region that change their conformations upon substrate binding. Five of these residues (Met76, Lys77, Thr79, Asp80 and Ser81) are highly flexible and cover the range of conformations observed in the substrate-bound state, while the remaining five (Arg74, Lys75, Asp78, Glu82 and Glu83) are much more rigid and do not populate conformations typical of the substrate-bound form. The ensemble of conformations representing the Ca^2 +-CaM state obtained in this study is in good agreement with residual dipolar coupling, paramagnetic resonance enhancement, small-angle X-ray scattering and fluorescence resonance energy transfer measurements, which were not used as restraints in the calculations. These results provide initial evidence that chemical shifts can be used to characterize the conformational fluctuations of multidomain proteins.     

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather than Electrostatic Repulsion

RNA folding occurs via a series of transitions between metastable intermediate states for Mg²⁺ concentrations below those needed to fold the native structure. In general, these folding intermediates are considerably less compact than their respective native states. Our previous work demonstrates that the major equilibrium intermediate of the 154-residue specificity domain (S-domain) of the Bacillus subtilis RNase P RNA is more extended than its native structure. We now investigate two models with falsifiable predictions regarding the origins of the extended intermediate structures in the S-domains of the B. subtilis and the Escherichia coli RNase P RNA that belong to different classes of P RNA and have distinct native structures. The first model explores the contribution of electrostatic repulsion, while the second model probes specific interactions in the core of the folding intermediate. Using small-angle X-ray scattering and Langevin dynamics simulations, we show that electrostatics plays only a minor role, whereas specific interactions largely account for the extended nature of the intermediate. Structural contacts in the core, including a nonnative base pair, help to stabilize the intermediate conformation. We conclude that RNA folding intermediates adopt extended conformations due to short-range, nonnative interactions rather than generic electrostatic repulsion of helical domains. These principles apply to other ribozymes and riboswitches that undergo functionally relevant conformational changes.     

Comparison of the three-dimensional molecular models of bovine submicellar caseins with small-angle X-ray scattering. Influence of protein hydration

To test the applicability of two energy-minimized, three-dimensional structures of the bovine casein submicelle, theoretical small-angle X-ray scattering curves in the presence and absence of water were compared to experimental data. The published method simulates molecular dynamics of proteins in solution by employing adjustable Debye-Waller temperature factors (B factors) for the protein, for the solvent, and for protein-bound water. The programs were first tested upon bovine pancreatic trypsin inhibitor beginning with its known X-ray crystal structure. To approximate the degree of protein hydration previously determined by NMR relaxation experiments (0.01 g water/g protein), 120 water molecules were docked into the large void of the-casein portion of the structure for both the symmetric and asymmetric casein submicelle models. To approximate hydrodynamic hydration (0.244 g water/g protein), 2703 water molecules were added to each of the above structures using the droplet algorithm in the Sybyl molecular modeling package. All structures were then energy-minimized and their solvation energies calculated. Theoretical small-angle X-ray scattering curves were calculated for all unhydrated and hydrated structures and compared with experimentally determined scattering profiles for submicellar casein. Best results were achieved with the 120-bound-water structure for both the symmetric and asymmetric submicelle models. Comparison of results for the protein submicelle models with those for the theoretical and literature values of bovine pancreatic trypsin inhibitor demonstrates the applicability of the methodology.Reference to a brand or firm name does not constitute endorsement by the U.S. Department of Agriculture over others of a similar nature not mentioned.     

Small-angle X-ray scattering studies of 7S and 11S globulins from pea (Pisum sativum)

Small-angle X-ray scattering studies have been conducted on solutions of 11S and 7S globulins isolated from peas (Pisum sativum cv. Filby), and the radii of gyration and molecular weights determined. The general features of the scattering curves were similar to those reported for other seed storage proteins.     

The compact and expanded denatured conformations of apomyoglobin in the methanol-water solvent

We have performed a detailed study of methanol-induced conformational transitions of horse heart apomyoglobin (apoMb) to investigate the existence of the compact and expanded denatured states. A combination of far- and near-ultraviolet circular dichroism, NMR spectroscopy, and small-angle X-ray scattering (SAXS) was used, allowing a phase diagram to be constructed as a function of pH and the methanol concentration. The phase diagram contains four conformational states, the native (N), acid-denatured (U(A)), compact denatured (I(M)), and expanded helical denatured (H) states, and indicates that the compact denatured state (I(M)) is stable under relatively mild denaturing conditions, whereas the expanded denatured states (U(A) and H) are realized under extreme conditions of pH (strong electric repulsion) or alcohol concentration (weak hydrophobic interaction). The results of this study, together with many previous studies in the literature, indicate the general existence of the compact denatured states not only in the salt-pH plane but also in the alcohol-pH plane. Furthermore, to determine the general feature of the H conformation we used several proteins including ubiquitin, ribonuclease A, alpha-lactalbumin, beta-lactoglobulin, and Streptomyces subtilisin inhibitor (SSI) in addition to apoMb. SAXS studies of these proteins in 60% methanol showed that the H states of these all proteins have expanded and nonglobular conformations. The qualitative agreement of the experimental data with computer-simulated Kratky profiles also supports this structural feature of the H state.