Probing the structure of the Caulobacter crescentus ribosome with chemical labeling and mass spectrometry

ETD-Cleavable Linker for Confident Cross-linked Peptide Identifications

Peptide cross-links formed using the homobifunctional-linker diethyl suberthioimidate (DEST) are shown to be ETD-cleavable. DEST has a spacer arm consisting of a 6-carbon alkyl chain and it cleaves at the amidino groups created upon reaction with primary amines. In ETD MS² spectra, DEST cross-links can be recognized based on mass pairs consisting of peptide-NH₂^• and peptide+linker+NH₃ ions, and backbone cleavages are more equally distributed over the two constituent peptides compared with collisional activation. Dead ends that are often challenging to distinguish from cross-links are diagnosed by intense reporter ions. ETD mass pairs can be used in MS³ experiments to confirm cross-link identifications. These features provide a simple but reliable approach to identify cross-links that should facilitate studies of protein complexes.

Covalent Surface Modification of Prions: A Mass Spectrometry-Based Means of Detecting Distinctive Structural Features of Prion Strains

Prions (PrP(Sc)) are molecular pathogens that are able to convert the isosequential normal cellular prion protein (PrP(C)) into a prion. The only demonstrated difference between PrP(C) and PrP(Sc) is conformational: they are isoforms. A given host can be infected by more than one kind or strain of prion. Five strains of hamster-adapted scrapie [Sc237 (=263K), drowsy, 139H, 22AH, and 22CH] and recombinant PrP were reacted with five different concentrations (0, 1, 5, 10, and 20 mM) of reagent (N-hydroxysuccinimide ester of acetic acid) that acetylates lysines. The extent of lysine acetylation was quantitated by mass spectrometry. The lysines in rPrP react similarly. The lysines in the strains react differently from one another in a given strain and react differently when strains are compared. Lysines in the C-terminal region of prions have different strain-dependent reactivity. The results are consistent with a recently proposed model for the structure of a prion. This model proposes that prions are composed of a four-rung β-solenoid structure comprised of four β-sheets that are joined by loops and turns of amino acids. Variation in the amino acid composition of the loops and β-sheet structures is thought to result in different strains of prions.

Dynamics of ribosomal protein S1 on a bacterial ribosome with cross-linking and mass spectrometry

Ribosomal protein S1 has been shown to be a significant effector of prokaryotic translation. The protein is in fact capable of efficiently initiating translation, regardless of the presence of a Shine-Dalgarno sequence in mRNA. Structural insights into this process have remained elusive, as S1 is recalcitrant to traditional techniques of structural analysis, such as x-ray crystallography. Through the application of protein cross-linking and high resolution mass spectrometry, we have detailed the ribosomal binding site of S1 and have observed evidence of its dynamics. Our results support a previous hypothesis that S1 acts as the mRNA catching arm of the prokaryotic ribosome. We also demonstrate that in solution the major domains of the 30S subunit are remarkably flexible, capable of moving 30-50Å with respect to one another.

Discovering free energy basins for macromolecular systems via guided multiscale simulation

An approach for the automated discovery of low free energy states of macromolecular systems is presented. The method does not involve delineating the entire free energy landscape but proceeds in a sequential free energy minimizing state discovery; i.e., it first discovers one low free energy state and then automatically seeks a distinct neighboring one. These states and the associated ensembles of atomistic configurations are characterized by coarse-grained variables capturing the large-scale structure of the system. A key facet of our approach is the identification of such coarse-grained variables. Evolution of these variables is governed by Langevin dynamics driven by thermal-average forces and mediated by diffusivities, both of which are constructed by an ensemble of short molecular dynamics runs. In the present approach, the thermal-average forces are modified to account for the entropy changes following from our knowledge of the free energy basins already discovered. Such forces guide the system away from the known free energy minima, over free energy barriers, and to a new one. The theory is demonstrated for lactoferrin, known to have multiple energy-minimizing structures. The approach is validated using experimental structures and traditional molecular dynamics. The method can be generalized to enable the interpretation of nanocharacterization data (e.g., ion mobility-mass spectrometry, atomic force microscopy, chemical labeling, and nanopore measurements).

Chemical reactivity of brome mosaic virus capsid protein

Viral particles are biological machines that have evolved to package, protect, and deliver the viral genome into the host via regulated conformational changes of virions. We have developed a procedure to modify lysine residues with S-methylthioacetimidate across the pH range from 5.5 to 8.5. Lysine residues that are not completely modified are involved in tertiary or quaternary structural interactions, and their extent of modification can be quantified as a function of pH. This procedure was applied to the pH-dependent structural transitions of brome mosaic virus (BMV). As the reaction pH increases from 5.5 to 8.5, the average number of modified lysine residues in the BMV capsid protein increases from 6 to 12, correlating well with the known pH-dependent swelling behavior of BMV virions. The extent of reaction of each of the capsid protein's lysine residues has been quantified at eight pH values using coupled liquid chromatography-tandem mass spectrometry. Each lysine can be assigned to one of three structural classes identified by inspection of the BMV virion crystal structure. Several lysine residues display reactivity that indicates their involvement in dynamic interactions that are not obvious in the crystal structure. The influence of several capsid protein mutants on the pH-dependent structural transition of BMV has also been investigated. Mutant H75Q exhibits an altered swelling transition accompanying solution pH increases. The H75Q capsids show increased reactivity at lysine residues 64 and 130, residues distal from the dimer interface occupied by H75, across the entire pH range.

Ratiometric pulse-chase amidination mass spectrometry as a probe of biomolecular complex formation

Selective chemical modification of protein side chains coupled with mass spectrometry is often most informative when used to compare residue-specific reactivities in a number of functional states or macromolecular complexes. Herein, we develop ratiometric pulse-chase amidination mass spectrometry (rPAm-MS) as a site-specific probe of lysine reactivities at equilibrium using the Cu(I)-sensing repressor CsoR from Bacillus subtilis as a model system. CsoR in various allosteric states was reacted with S-methyl thioacetimidate (SMTA) for pulse time, t, and chased with excess of S-methyl thiopropionimidate (SMTP) (Δ = 14 amu), quenched and digested with chymotrypsin or Glu-C protease, and peptides were quantified by high-resolution matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry and/or liquid chromatography electrospray ionization tandem mass spectrometry (LC-ESI-MS/MS). We show that the reactivities of individual lysines from peptides containing up to three Lys residues are readily quantified using this method. New insights into operator DNA binding and the Cu(I)-mediated structural transition in the tetrameric copper sensor CsoR are also obtained.

Multiscale simulation of microbe structure and dynamics

A multiscale mathematical and computational approach is developed that captures the hierarchical organization of a microbe. It is found that a natural perspective for understanding a microbe is in terms of a hierarchy of variables at various levels of resolution. This hierarchy starts with the N -atom description and terminates with order parameters characterizing a whole microbe. This conceptual framework is used to guide the analysis of the Liouville equation for the probability density of the positions and momenta of the N atoms constituting the microbe and its environment. Using multiscale mathematical techniques, we derive equations for the co-evolution of the order parameters and the probability density of the N-atom state. This approach yields a rigorous way to transfer information between variables on different space-time scales. It elucidates the interplay between equilibrium and far-from-equilibrium processes underlying microbial behavior. It also provides framework for using coarse-grained nanocharacterization data to guide microbial simulation. It enables a methodical search for free-energy minimizing structures, many of which are typically supported by the set of macromolecules and membranes constituting a given microbe. This suite of capabilities provides a natural framework for arriving at a fundamental understanding of microbial behavior, the analysis of nanocharacterization data, and the computer-aided design of nanostructures for biotechnical and medical purposes. Selected features of the methodology are demonstrated using our multiscale bionanosystem simulator DeductiveMultiscaleSimulator. Systems used to demonstrate the approach are structural transitions in the cowpea chlorotic mosaic virus, RNA of satellite tobacco mosaic virus, virus-like particles related to human papillomavirus, and iron-binding protein lactoferrin.

Structural analysis of a prokaryotic ribosome using a novel amidinating cross-linker and mass spectrometry

The structure of the Escherichia coli ribosome, a 2.5 MDa ribonucleoprotein complex containing more than 50 proteins, was probed using the novel amidinating cross-linker diethyl suberthioimidate (DEST) and mass spectrometry. Peptide cross-links derived from this complex structure were identified at high confidence (FDR 0.8%) from precursor mass measurements and collision-induced dissociation (CID) fragmentation spectra. The acquired cross-linking data were found to be in excellent agreement with the crystal structure of the E. coli ribosome. DEST cross-links are particularly amenable to strong cation exchange (SCX) chromatography, facilitating a large-scale analysis. SCX enrichment and fractionation were shown to increase the number of cross-link spectra matches in our analysis 10-fold. Evidence is presented that these techniques can be used to study complex interactomes.

Novel amidinating cross-linker for facilitating analyses of protein structures and interactions

A novel bifunctional thioimidate cross-linking reagent (diethyl suberthioimidate) that modifies amines without sacrificing their native basicity is developed. Intermolecular cross-linking of neurotensin and intramolecular cross-linking of cytochrome c under physiological conditions is investigated with this reagent. Because it does not perturb the electrostatic properties of a protein, it is unlikely to lead to artifactual conclusions about native protein structure. The interpeptide cross-links formed with this reagent are easily separated from other tryptic fragments using strong cation exchange chromatography, and they have a readily identified mass spectrometric signature. The use of this novel amidinating protein cross-linking reagent holds great promise for efficient, large-scale structural analysis of complex systems.

All-atom multiscale simulation of cowpea chlorotic mottle virus capsid swelling

An all-atom multiscale computational modeling approach, molecular dynamics/order parameter extrapolation (MD/OPX), has recently been developed for simulating large bionanosystems. It accelerates MD simulations and addresses rapid atomistic fluctuations and slowly varying nanoscale dynamics of bionanosystems simultaneously. With modules added to account for water molecules and ions, MD/OPX is applied to simulate the swelling of cowpea chlorotic mottle virus (CCMV) capsid solvated in a host medium in this study. Simulation results show that the N-terminal arms of capsid proteins undergo large deviations from the initial configurations with their length extended quickly during the early stage of capsid swelling. The capsid swelling is a symmetry-breaking process involving local initiation and front propagation. The capsid swelling rate is approximately 0.25 nm/ns (npn) during the early stage of the simulation, and propagation of the structural transition across the capsid is roughly 0.6 npn. The system conditions that affect swelling of the capsid are analyzed. Prospects for creating a phase diagram for CCMV capsid swelling and using predictions to guide experiments are discussed.

De novo sequencing of tryptic peptides derived from Deinococcus radiodurans ribosomal proteins using 157 nm photodissociation MALDI TOF/TOF mass spectrometry

Vacuum ultraviolet photodissociation of peptide ions in a matrix assisted laser desorption ionization (MALDI) tandem time-of-flight (TOF) mass spectrometer is used to characterize peptide mixtures derived from Deinococcus radiodurans ribosomal proteins. Tryptic peptides from 52 proteins were separated by reverse-phase liquid chromatography and spotted onto a MALDI plate. From 192 sample spots, 492 peptide ions were isolated, fragmented by both photodissociation and postsource decay (PSD), and then de novo sequenced. Three-hundred seventy-two peptides yielded sequences with 5 or more amino acids. Homology searches of these sequences against the whole bacterial proteome identified 49 ribosomal proteins, 45 of which matched with two or more peptides. Peptide de novo sequencing identified slightly more proteins than conventional database searches using Mascot and was particularly advantageous in identifying unexpected peptide modifications. In the present analysis, 52 peptide modifications were identified by de novo sequencing, most of which were not recognized by database searches.

'Fixed charge' chemical derivatization and data dependant multistage tandem mass spectrometry for mapping protein surface residue accessibility

Protein surface accessible residues play an important role in protein folding, protein-protein interactions and protein-ligand binding. However, a common problem associated with the use of selective chemical labeling methods for mapping protein solvent accessible residues is that when a complicated peptide mixture resulting from a large protein or protein complex is analyzed, the modified peptides may be difficult to identify and characterize amongst the largely unmodified peptide population (i.e., the 'needle in a haystack' problem). To address this challenge, we describe here the development of a strategy involving the synthesis and application of a novel 'fixed charge' sulfonium ion containing lysine-specific protein modification reagent, S,S'-dimethylthiobutanoylhydroxysuccinimide ester (DMBNHS), coupled with capillary HPLC-ESI-MS, automated CID-MS/MS, and data-dependant neutral loss mode MS(3) in an ion trap mass spectrometer, to map the surface accessible lysine residues in a small model protein, cellular retinoic acid binding protein II (CRABP II). After reaction with different reagent:protein ratios and digestion with Glu-C, modified peptides are selectively identified and the number of modifications within each peptide are determined by CID-MS/MS, via the exclusive neutral loss(es) of dimethylsulfide, independently of the amino acid composition and precursor ion charge state (i.e., proton mobility) of the peptide. The observation of these characteristic neutral losses are then used to automatically 'trigger' the acquisition of an MS(3) spectrum to allow the peptide sequence and the site(s) of modification to be characterized. Using this approach, the experimentally determined relative solvent accessibilities of the lysine residues were found to show good agreement with the known solution structure of CRABP II.

B. subtilis ribosomal proteins: structural homology and post-translational modifications

Ribosomal proteins of the model gram-positive bacterium B. subtilis 168 were extensively characterized in a proteomic study. Mass spectra of the 52 proteins expected to be constitutive components of the 70S ribosome were recorded. Peptide MS/MS analysis with an average sequence coverage of 85% supported the identification of these proteins and facilitated the unambiguous assignment of post-translational modifications, including the methylation of S7, L11, and L16 and the N-terminal acetylation of S9. In addition, the high degree of structural homology between B. subtilis and other eubacterial ribosomal proteins was demonstrated through chemical labeling with S-methylthioacetimidate. One striking difference from previous characterizations of bacterial ribosomal proteins is that dozens of protein masses were found to be in error and not easily accounted for by post-translational modifications. This, in turn, led us to discover an inordinate number of sequencing errors in the reference genome of B. subtilis 168. We have found that these errors have been corrected in a recently revised version of the genome.

Structural characterization of the conformational change in calbindin-D28k upon calcium binding using differential surface modification analyzed by mass spectrometry

Calbindin-D28k is a calcium binding protein with six EF hand domains. Calbindin-D28k is unique in that it functions as both a calcium buffer and a sensor protein. It is found in many tissues, including brain, pancreas, kidney, and intestine, playing important roles in each. Calbindin-D28k is known to bind four calcium ions and upon calcium binding undergoes a conformational change. The structure of apo calbindin-D28k is in an ordered state, transitioning into a disordered state as calcium is bound. Once fully loaded with four calcium ions, it again takes on an ordered state. The solution structure of disulfide-reduced holo-calbindin-D28k has been determined by NMR, while the structure of apo calbindin-D28k has yet to be determined. Differential surface modification of lysine and histidine residues analyzed by mass spectrometry has been used in this study to identify, for the first time, the specific regions of calbindin-D28k undergoing conformational changes between the holo and apo states. Using differential surface modification in combination with mass spectrometry, EF hands 1 and 4 as well as the linkers before EF hand 1 and the linkers between EF hands 4 and 5 and EF hands 5 and 6 were identified as regions of conformational change between apo and holo calbindin-D28k. Under the experimental conditions employed, EF hands 2 and 6, which are known not to bind calcium, were unaffected in either form. EF hand 2 is highly accessible; however, EF hand 6 was determined not to be surface accessible in either form. Previous research has identified a disulfide bond between cysteines 94 and 100 in the holo state. Until now, it was unknown whether this bond also exists in the apo form. Our data confirm the presence of the disulfide bond between cysteines 94 and 100 in the holo form and indicate that there is predominantly no disulfide bond between these residues in the apoprotein.

Probing protein structure by amino acid-specific covalent labeling and mass spectrometry

For many years, amino acid-specific covalent labeling has been a valuable tool to study protein structure and protein interactions, especially for systems that are difficult to study by other means. These covalent labeling methods typically map protein structure and interactions by measuring the differential reactivity of amino acid side chains. The reactivity of amino acids in proteins generally depends on the accessibility of the side chain to the reagent, the inherent reactivity of the label and the reactivity of the amino acid side chain. Peptide mass mapping with ESI- or MALDI-MS and peptide sequencing with tandem MS are typically employed to identify modification sites to provide site-specific structural information. In this review, we describe the reagents that are most commonly used in these residue-specific modification reactions, details about the proper use of these covalent labeling reagents, and information about the specific biochemical problems that have been addressed with covalent labeling strategies.

Limited proteolysis analysis of the ribosome is affected by subunit association

Our understanding of the structural organization of ribosome assembly intermediates, in particular those intermediates that result from misfolding leading to their eventual degradation within the cell, is limited because of the lack of methods available to characterize assembly intermediate structures. Because conventional structural approaches, such as NMR, X-ray crystallography, and cryo-EM, are not ideally suited to characterize the structural organization of these flexible and sometimes heterogeneous assembly intermediates, we have set out to develop an approach combining limited proteolysis with matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) that might be applicable to ribonucleoprotein complexes as large as the ribosome. This study focuses on the limited proteolysis behavior of appropriately assembled ribosome subunits. Isolated subunits were analyzed using limited proteolysis and MALDI-MS and the results were compared with previous data obtained from 70S ribosomes. Generally, ribosomal proteins were found to be more stable in 70S ribosomes than in their isolated subunits, consistent with a reduction in conformational flexibility on subunit assembly. This approach demonstrates that limited proteolysis combined with MALDI-MS can reveal structural changes to ribosomes on subunit assembly or disassembly, and provides the appropriate benchmark data from 30S, 50S, and 70S proteins to enable studies of ribosome assembly intermediates. (c) 2009 Wiley Periodicals, Inc. Biopolymers 91: 410-422, 2009.

Frameshifting RNA pseudoknots: structure and mechanism

Programmed ribosomal frameshifting (PRF) is one of the multiple translational recoding processes that fundamentally alters triplet decoding of the messenger RNA by the elongating ribosome. The ability of the ribosome to change translational reading frames in the −1 direction (−1 PRF) is employed by many positive strand RNA viruses, including economically important plant viruses and many human pathogens, such as retroviruses, e.g., HIV-1, and coronaviruses, e.g., the causative agent of severe acute respiratory syndrome (SARS), in order to properly express their genomes. −1 PRF is programmed by a bipartite signal embedded in the mRNA and includes a heptanucleotide “slip site” over which the paused ribosome “backs up” by one nucleotide, and a downstream stimulatory element, either an RNA pseudoknot or a very stable RNA stem–loop. These two elements are separated by six to eight nucleotides, a distance that places the 5′ edge of the downstream stimulatory element in direct contact with the mRNA entry channel of the 30S ribosomal subunit. The precise mechanism by which the downstream RNA stimulates −1 PRF by the translocating ribosome remains unclear. This review summarizes the recent structural and biophysical studies of RNA pseudoknots and places this work in the context of our evolving mechanistic understanding of translation elongation. Support for the hypothesis that the downstream stimulatory element provides a kinetic barrier to the ribosome-mediated unfolding is discussed.

Stochastic dynamics of bionanosystems: Multiscale analysis and specialized ensembles

An approach for simulating bionanosystems such as viruses and ribosomes is presented. This calibration-free approach is based on an all-atom description for bionanosystems, a universal interatomic force field, and a multiscale perspective. The supramillion-atom nature of these bionanosystems prohibits the use of a direct molecular dynamics approach for phenomena such as viral structural transitions or self-assembly that develop over milliseconds or longer. A key element of these multiscale systems is the cross-talk between, and consequent strong coupling of processes over many scales in space and time. Thus, overall nanoscale features of these systems control the relative probability of atomistic fluctuations, while the latter mediate the average forces and diffusion coefficients that induce the dynamics of these nanoscale features. This feedback loop is overlooked in typical coarse-grained methods. We elucidate the role of interscale cross-talk and overcome bionanosystem simulation difficulties with (1) automated construction of order parameters (OPs) describing suprananometer scale structural features, (2) construction of OP-dependent ensembles describing the statistical properties of atomistic variables that ultimately contribute to the entropies driving the dynamics of the OPs, and (3) the derivation of a rigorous equation for the stochastic dynamics of the OPs. As the OPs capture hydrodynamic modes in the host medium, "long-time tails" in the correlation functions yielding the generalized diffusion coefficients do not emerge. Since the atomic-scale features of the system are treated statistically, several ensembles are constructed that reflect various experimental conditions. Attention is paid to the proper use of the Gibbs hypothesized equivalence of long-time and ensemble averages to accommodate the varying experimental conditions. The theory provides a basis for a practical, quantitative bionanosystem modeling approach that preserves the cross-talk between the atomic and nanoscale features. A method for integrating information from nanotechnical experimental data in the derivation of equations of stochastic OP dynamics is also introduced.

Biophysical studies of bacterial ribosome assembly

The assembly of the bacterial ribosome involves the association of over 50 proteins to 3 large RNA molecules, and it represents a major metabolic activity for rapidly growing bacteria. The availability of atomic structures of the ribosome and the application of biochemical and biophysical methods have led to rapid progress in understanding the mechanistic details of ribosome assembly. The basic steps required to assemble a ribosome are outlined, and the contributions of mass spectrometry, computational methods, and RNA-folding studies in understanding these steps are detailed. This complex process takes place with both sequential and parallel processing that is coordinated to ensure efficient and complete assembly of ribosomes to meet the demands of cell growth.

A top-down/bottom-up study of the ribosomal proteins of Caulobacter crescentus

Ribosomes from the Gram-negative alpha-proteobacterium Caulobacter crescentus were isolated using standard methods. Proteins were separated using a two-dimensional liquid chromatographic system that allowed the analysis of whole proteins by direct coupling to an ESI-QTOF mass spectrometer and of proteolytic digests by a number of mass spectrometric methods. The masses of 53 of 54 ribosomal proteins were directly measured. Protein identifications and proposed post-translational modifications were supported by proteolysis with trypsin, endoprotease Glu-C, and exoproteases carboxypeptidases Y and P. Tryptic peptide mass maps show an average sequence coverage of 62%, and carboxypeptidase C-terminal sequence tagging provided unambiguous identification of the small, highly basic proteins of the large subunit. C. crescentus presents some post-translational modifications that are similar to those of Escherichia coli (e.g., N-terminal acetylation of S9 and S18) along with some unique variations, such as a near absence of L7 and extensive modification of L11. The comprehensive description of this organism's ribosomal proteome provides a foundation for the study of ribosome structure, dependence of post-translational modifications on growth conditions, and the evolution of subcellular organelles.