Structural model of the CopA copper ATPase of Enterococcus hirae based on chemical cross-linking.

The CopA copper ATPase of Enterococcus hirae belongs to the family of heavy metal pumping CPx-type ATPases and shares 43% sequence similarity with the human Menkes and Wilson copper ATPases. Due to a lack of suitable protein crystals, only partial three-dimensional structures have so far been obtained for this family of ion pumps. We present a structural model of CopA derived by combining topological information obtained by intramolecular cross-linking with molecular modeling. Purified CopA was cross-linked with different bivalent reagents, followed by tryptic digestion and identification of cross-linked peptides by mass spectrometry. The structural proximity of tryptic fragments provided information about the structural arrangement of the hydrophilic protein domains, which was integrated into a three-dimensional model of CopA. Comparative modeling of CopA was guided by the sequence similarity to the calcium ATPase of the sarcoplasmic reticulum, Serca1, for which detailed structures are available. In addition, known partial structures of CPx-ATPase homologous to CopA were used as modeling templates. A docking approach was used to predict the orientation of the heavy metal binding domain of CopA relative to the core structure, which was verified by distance constraints derived from cross-links. The overall structural model of CopA resembles the Serca1 structure, but reveals distinctive features of CPx-type ATPases. A prominent feature is the positioning of the heavy metal binding domain. It features an orientation of the Cu binding ligands which is appropriate for the interaction with Cu-loaded metallochaperones in solution. Moreover, a novel model of the architecture of the intramembranous Cu binding sites could be derived.

The collaboratory for MS3D: a new cyberinfrastructure for the structural elucidation of biological macromolecules and their assemblies using mass spectrometry-based approaches.

Modern biomedical research is evolving with the rapid growth of diverse data types, biophysical characterization methods, computational tools and extensive collaboration among researchers spanning various communities and having complementary backgrounds and expertise. Collaborating researchers are increasingly dependent on shared data and tools made available by other investigators with common interests, thus forming communities that transcend the traditional boundaries of the single research laboratory or institution. Barriers, however, remain to the formation of these virtual communities, usually due to the steep learning curve associated with becoming familiar with new tools, or with the difficulties associated with transferring data between tools. Recognizing the need for shared reference data and analysis tools, we are developing an integrated knowledge environment that supports productive interactions among researchers. Here we report on our current collaborative environment, which focuses on bringing together structural biologists working in the area of mass spectrometric based methods for the analysis of tertiary and quaternary macromolecular structures (MS3D) called the Collaboratory for MS3D (C-MS3D). C-MS3D is a Web-portal designed to provide collaborators with a shared work environment that integrates data storage and management with data analysis tools. Files are stored and archived along with pertinent meta data in such a way as to allow file handling to be tracked (data provenance) and data files to be searched using keywords and modification dates. While at this time the portal is designed around a specific application, the shared work environment is a general approach to building collaborative work groups. The goal of this is to not only provide a common data sharing and archiving system, but also to assist in the building of new collaborations and to spur the development of new tools and technologies.

Partial acetylation of lysine residues improves intraprotein cross-linking.

Intramolecular cross-linking coupled with mass spectrometric identification of cross-linked amino acids is a rapid method for elucidating low-resolution protein tertiary structures or fold families. However, previous cross-linking studies on model proteins, such as cytochrome c and ribonuclease A, identified a limited number of peptide cross-links that are biased toward only a few of the potentially reactive lysine residues. Here, we report an approach to improve the diversity of intramolecular protein cross-linking starting with a systematic quantitation of the reactivity of lysine residues of a model protein, bovine cytochrome c. Relative lysine reactivities among the 18 lysine residues of cytochrome c were determined by the ratio of d0 and acetyl-d3 groups at each lysine after partial acetylation with sulfosuccinimidyl acetate followed by denaturation and quantitative acetylation of remaining unmodified lysines with acetic-d6 anhydride. These lysine reactivities were then compared with theoretically derived pKa and relative solvent accessibility surface values. To ascertain if partial N-acetylation of the most reactive lysine residues prior to cross-linking can redirect and increase the observable Lys-Lys cross-links, partially acetylated bovine cytochrome c was cross-linked with the amine-specific, bis-functional reagent, bis(sulfosuccinimidyl)suberate. After proteolysis and mass spectrometry analysis, partial acetylation was shown to significantly increase the number of observable peptides containing Lys-Lys cross-links, shifting the pattern from the most reactive lysine residues to less reactive ones. More importantly, these additional cross-linked peptides contained novel Lys-Lys cross-link information not seen in the non-acetylated protein and provided additional distance constraints that were consistent with the crystal structure and facilitated the identification of the proper protein fold.

Structure and dynamics of dark-state bovine rhodopsin revealed by chemical cross-linking and high-resolution mass spectrometry.

Recent work using chemical cross-linking to define interresidue distance constraints in proteins has shown that these constraints are useful for testing tertiary structural models. We applied this approach to the G-protein-coupled receptor bovine rhodopsin in its native membrane using lysine- and cysteine-targeted bifunctional cross-linking reagents. Cross-linked proteolytic peptides of rhodopsin were identified by combined liquid chromatography and FT-ICR mass spectrometry with automated data-reduction and assignment software. Tandem mass spectrometry was used to verify cross-link assignments and locate the exact sites of cross-link attachment. Cross-links were observed to form between 10 pairs of residues in dark-state rhodopsin. For each pair, cross-linkers with a range of linker lengths were tested to determine an experimental distance-of-closest-approach (DCA) between reactive side-chain atoms. In all, 28 cross-links were identified using seven different cross-linking reagents. Molecular mechanics procedures were applied to published crystal structure data to calculate energetically achievable theoretical DCAs between reactive atoms without altering the position of the protein backbone. Experimentally measured DCAs are generally in good agreement with the theoretical DCAs. However, a cross-link between C316 and K325 in the C-terminal region cannot be rationalized by DCA simulations and suggests that backbone reorientation relative to the crystal coordinates occurs on the timescale of cross-linking reactions. Biochemical and spectroscopic data from other studies have found that the C-terminal region is highly mobile in solution and not fully represented by X-ray crystallography data. Our results show that chemical cross-linking can provide reliable three-dimensional structural information and insight into local conformational dynamics in a membrane protein.

"Zero-length" cross-linking in solid state as an approach for analysis of protein-protein interactions.

We have developed a new approach for the analysis of interacting interfaces in protein complexes and protein quaternary structure based on cross-linking in the solid state. Protein complexes are freeze-dried under vacuum, and cross-links are introduced in the solid phase by dehydrating the protein in a nonaqueous solvent creating peptide bonds between amino and carboxyl groups of the interacting peptides. Cross-linked proteins are digested into peptides with trypsin in both H2(16)O and H(2)18O and then readily distinguished in mass spectra by characteristic 8 atomic mass unit (amu) shifts reflecting incorporation of two 18O atoms into each C terminus of proteolytic peptides. Computer analysis of mass spectrometry (MS) and MS/MS data is used to identify the cross-linked peptides. We demonstrated specificity and reproducibility of our method by cross-linking homo-oligomeric protein complexes of glutathione-S-transferase (GST) from Schistosoma japonicum alone or in a mixture of many other proteins. Identified cross-links were predominantly of amide origin, but six esters and thioesters were also found. The cross-linked peptides were validated against the GST monomer and dimer X-ray structures and by experimental (MS/MS) analyses. Some of the identified cross-links matched interacting peptides in the native 3D structure of GST, indicating that the structure of GST and its oligomeric complex remained primarily intact after freeze-drying. The pattern of oligomeric GST obtained in solid state was the same as that obtained in solution by Ru (II) Bpy(3)2+ catalyzed, oxidative "zero-length" cross-linking, confirming that it is feasible to use our strategy for analyzing the molecular interfaces of interacting proteins or peptides.

Influence of crosslinker identity and position on gas-phase dissociation of Lys-Lys crosslinked peptides.

A systematic study of the dissociation patterns of crosslinked peptides analyzed by tandem mass spectrometry is reported. A series of 11-mer peptides was designed around either a polyalanine or polyglycine scaffold with arginine at the C terminus. One or two lysine residues were included at various locations within the peptides to effect inter- or intra-molecular crosslinking, respectively. Crosslinked species were generated with four commonly used amine-specific chemical crosslinking reagents: disuccinimidyl suberate (DSS), disuccinimidyl tartarate (DST), dithiobis(succinimidylpropionate) (DSP), and disuccinimidyl glutarate (DSG). The influence of precursor charge state, location of crosslink, and specific crosslinking reagent on the MS/MS dissociation pattern was examined. Observed trends in the dissociation patterns obtained for these species will allow for improvements to software used in the automated interpretation of crosslinked peptide MS/MS data.

PostDOCK: a structural, empirical approach to scoring protein ligand complexes.

In this work we introduce a postprocessing filter (PostDOCK) that distinguishes true binding ligand-protein complexes from docking artifacts (that are created by DOCK 4.0.1). PostDOCK is a pattern recognition system that relies on (1) a database of complexes, (2) biochemical descriptors of those complexes, and (3) machine learning tools. We use the protein databank (PDB) as the structural database of complexes and create diverse training and validation sets from it based on the "families of structurally similar proteins" (FSSP) hierarchy. For the biochemical descriptors, we consider terms from the DOCK score, empirical scoring, and buried solvent accessible surface area. For the machine-learners, we use a random forest classifier and logistic regression. Our results were obtained on a test set of 44 structurally diverse protein targets. Our highest performing descriptor combinations obtained approximately 19-fold enrichment (39 of 44 binding complexes were correctly identified, while only allowing 2 of 44 decoy complexes), and our best overall accuracy was 92%.

Unambiguous assignment of intramolecular chemical cross-links in modified mammalian membrane proteins by Fourier transform-tandem mass spectrometry.

Fourier transform tandem mass spectrometry (FT-MS/MS) can be used to unambiguously assign intramolecular chemical cross-links to specific amino acid residues even when two or more possible cross-linking sites are adjacent in the cross-linked protein. Bovine rhodopsin (Rho) in its dark-adapted state was intramolecularly cross-linked with lysine-cysteine (K-C) or lysine-lysine (K-K) cross-linkers to obtain interatomic distance information. Large, multiply charged, cross-linked peptide ions containing adjacent lysines, corresponding to Rho(50-86) (K(66) or K(67)) cross-linked to Rho(310-317) (C(316)) or Rho(318-348) (K(325) or K(339)), were fragmented by collision-induced dissociation (CID), infrared multiphoton dissociation (IRMPD), and electron capture dissociation (ECD). Complementary sequence-specific information was obtained by combining cross-link assignments; however, only ECD revealed full palmitoylation of adjacent cysteines (C(322) and C(323)) and cross-linking of K(67) (and not K(66)) to C(316), K(325), and K(339). ECD spectra contained crucial c- and z-ions resulting from cleavage of the bond between K(66) and K(67). To our knowledge, this work also presents the first demonstration that ECD can be used to characterize S-linked fatty acid acylation on cysteines. The comprehensive fragmentation of large peptides by CID, IRMPD, and particularly ECD, in conjunction with the high resolution and mass accuracy of FT-MS/MS, is shown to be a valuable means of characterizing mammalian membrane proteins with both chemical and posttranslational modifications.

Probing alpha-crystallin structure using chemical cross-linkers and mass spectrometry.

PURPOSE: Alternatives to X-ray crystallographic techniques are needed to probe the structure of the hetero-oligomeric lens protein alpha-crystallin. We utilized mass spectrometry for 3 dimensional analysis (MS3D) to study the quaternary structural characteristics of this important lens protein and molecular chaperone. METHODS: We have employed two types of chemical cross-linkers to probe key protein-protein and protein-solvent interactions of alpha-crystallin using MS3D. Native alpha-crystallin was exposed to 3,3'-dithiobis[sulfosuccinimidyl propionate] (DTSSP) and the common fixative, formaldehyde. The reaction products were denatured and enriched in cross-linked and modified species using size exclusion chromatography. Tryptic digests of these fractions were purified using reverse phase HPLC and analyzed by both electrospray and matrix assisted laser desorption mass spectrometry. Comprehensive spectra for each C18 fraction were screened for ions with mass unique to each chemical treatment and candidate sequences matching the experimental data were assigned using MS3D "Links" and "ASAP" software. Selected ions were sequenced by collision induced dissociation. RESULTS: Peptides including residues 164-175 of alphaB-crystallin and residues 1-99 of alphaA-crystallin were modified by formaldehyde and partially hydrolyzed DTSSP. Peptides containing modified lysines 11, 78, and 99 of alphaA-crystallin were sequenced and the modified amino acids identified. In addition, ions corresponding to intramolecular and/or intermolecular cross-links were assigned a sequence based on two criteria. First, the mass values observed were unique to a single cross-linking experiment and were not present in a control where no cross-linker was utilized. Second, two unique ions detected from different cross-linking experiments were correlated in that the structures assigned to the masses were equivalent apart from the structure of the cross-linker. One such correlation was found involving lysine121, within the "highly conserved alpha-crystallin domain" of alphaB-crystallin, cross-linked to either lysine11 or lysine99 of alphaA-crystallin. Another two independent correlations involving lysine72 of alphaB-crystallin were found that indicate cross-linking of two subunits of alphaB-crystallin through this same residue. CONCLUSIONS: Sequences of peptides modified by partially hydrolyzed DTSSP and formaldehyde provide experimental evidence for models of alpha-crystallin quaternary structure that suggest a similar tertiary fold for both alphaA-crystallin and alphaB-crystallin. Analogous to multiple phosphorylations along the N-terminus of alphaB-crystallin, our data indicate that the same region of alphaA-crystallin, up to and including lysine99 is also relatively accessible to modification despite its hydrophobicity. Mass correlation between experiments using different reagents suggests that cross-linking occurred between N-termini of adjacent subunits of alphaB-crystallin in the native complex in support of the amphiphilic, toroidal, or "open micelle" models. In addition, multiple cross-links involving lysine121 of the so called "dimer interface" region within the "highly conserved alpha-crystallin domain" indicate that this region is a site of inter-subunit contacts in the native context.

Optimal bundling of transmembrane helices using sparse distance constraints.

We present a two-step approach to modeling the transmembrane spanning helical bundles of integral membrane proteins using only sparse distance constraints, such as those derived from chemical cross-linking, dipolar EPR and FRET experiments. In Step 1, using an algorithm, we developed, the conformational space of membrane protein folds matching a set of distance constraints is explored to provide initial structures for local conformational searches. In Step 2, these structures refined against a custom penalty function that incorporates both measures derived from statistical analysis of solved membrane protein structures and distance constraints obtained from experiments. We begin by describing the statistical analysis of the solved membrane protein structures from which the theoretical portion of the penalty function was derived. We then describe the penalty function, and, using a set of six test cases, demonstrate that it is capable of distinguishing helical bundles that are close to the native bundle from those that are far from the native bundle. Finally, using a set of only 27 distance constraints extracted from the literature, we show that our method successfully recovers the structure of dark-adapted rhodopsin to within 3.2 A of the crystal structure.

Top-down characterization of nucleic acids modified by structural probes using high-resolution tandem mass spectrometry and automated data interpretation.

A top-down approach based on sustained off-resonance irradiation collision-induced dissociation (SORI-CID) has been implemented on an electrospray ionization (ESI) Fourier transform mass spectrometer (FTMS) to characterize nucleic acid substrates modified by structural probes. Solvent accessibility reagents, such as dimethyl sulfate (DMS), 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate (CMCT), and beta-ethoxy-alpha-ketobutyraldehyde (kethoxal, KT) are widely employed to reveal the position of single- vs double-stranded regions and obtain the footprint of bound proteins onto nucleic acids structures. Established methods require end-labeling of the nucleic acid constructs, probe-specific chemistry to produce strand cleavage at the modified nucleotides, and analysis by polyacrylamide gel electrophoresis to determine the position of the susceptible sites. However, these labor-intensive procedures can be avoided when mass spectrometry is used to identify the probe-induced modifications from their characteristic mass signatures. In particular, ESI-FTMS can be directly employed to monitor the conditions of probe application to avoid excessive alkylation, which could induce unwanted distortion or defolding of the substrate of interest. The sequence position of the covalent modifications can be subsequently obtained from classic tandem techniques, which allow for the analysis of individual target adducts present in complex reaction mixtures with no need for separation techniques. Selection and activation by SORI-CID has been employed to reveal the position of adducts in nucleic acid substrates in excess of 6 kDa. The stability of the different covalent modifications under SORI-CID conditions was investigated. Multiple stages of isolation and activation were employed in MS(n)() experiments to obtain the desired sequence information whenever the adduct stability was not particularly favorable, and SORI-CID induced the facile loss of the modified base. A new program called MS2Links was developed for the automated reduction and interpretation of fragmentation data obtained from modified nucleic acids. Based on an algorithm that searches for plausible isotopic patterns, the data reduction module is capable of discriminating legitimate signals from noise spikes of comparable intensity. The fragment identification module calculates the monoisotopic mass of ion products expected from a certain sequence and user-defined covalent modifications, which are finally matched with the signals selected by the data reduction program. Considering that MS2Links can generate similar fragment libraries for peptides and their covalent conjugates with other peptides or nucleic acids, this program provides an integrated platform for the structural investigation of protein-nucleic acid complexes based on cross-linking strategies and top-down ESI-FTMS.

A top-down method for the determination of residue-specific solvent accessibility in proteins.

We present a method employing top-down Fourier transform mass spectrometry (FTMS) for the rapid profiling of amino acid side-chain reactivity. The reactivity of side-chain groups can be used to infer residue-specific solvent accessibility and can also be used in the same way as H/D exchange reactions to probe protein structure and interactions. We probed the reactivity of the N-terminal and epsilon-lysine amino groups of ubiquitin by reaction with N-hydroxysuccinimidyl acetate (NHSAc), which specifically acetylates primary amines. Using a hybrid Q-FTMS instrument, we observed several series of multiply acetylated ubiquitin ions that varied with the NHSAc:protein stoichiometry. We isolated and fragmented each member of the series of acetylated ubiquitin ions in the front end of the instrument and measured the fragment ion masses in the FTMS analyzer cell to determine which residue positions were modified. As we increased the NHSAc:protein stoichiometric ratio, identification of the fragments from native protein and protein with successively increasing modification allowed the assignment of the complete order of reactivity of the primary amino groups in ubiquitin (Met 1 approximately Lys 6 approximately Lys 48 approximately Lys 63>Lys 33>Lys 11>Lys 27, Lys 29). These results are in excellent agreement with the reactivity expected from other studies and predicted from the known crystal structure of ubiquitin. The top-down approach eliminates the need for proteolytic digestion, high-performance liquid chromatographic separations and all other chemical steps except the labeling reaction, making it rapid and amenable to automation using small quantities of protein.

MS2Assign, automated assignment and nomenclature of tandem mass spectra of chemically crosslinked peptides.

In a previous report (Young et al., Proc. Natl. Acad. Sci. U.S.A. 2000, 97, 5802-5806), we provided a proof-of-principle for fold recognition of proteins using a homobifunctional amine-specific chemical crosslinking reagent in combination with mass spectrometry analysis and homology modeling. In this current work, we propose a systematic nomenclature to describe the types of peptides that are generated after proteolysis of crosslinked proteins, their fragmentation by tandem mass spectrometry, and an automated algorithm for MS/MS spectral assignment called "MS2Assign." Several examples are provided from crosslinked peptides and proteins including HIV-integrase, cytochrome c, ribonuclease A, myoglobin, cytidine 5-monophosphate N-acetylneuraminic acid synthetase, and the peptide thymopentin. Tandem mass spectra were obtained from various crosslinked peptides using post source decay MALDI-TOF and collision induced dissociation on a quadrupole-TOF instrument, along with their automated interpretation using MS2Assign. A variety of possible outcomes are described and categorized according to the number of modified lysines and/or peptide chains involved, as well as the presence of singly modified (dead-end) lysine residues. In addition, the proteolysis and chromatographic conditions necessary for optimized crosslinked peptide recovery are presented.

A top down approach to protein structural studies using chemical cross-linking and Fourier transform mass spectrometry.

Mass spectrometric analysis of wild-type proteins that have been covalently modified by bifunctional cross-linking reagents and then digested proteolytically can be used to obtain low-resolution distance constraints, which can be useful for protein structure determination. Limitations of this approach include time-consuming separation steps, such as the separation of internally cross-linked protein monomers from covalent dimers, and a susceptibility to artifacts due to low levels of natural and man-made peptide modifications that can be mistaken for cross-linked species. The results presented here show that when a crude cross-linked protein mixture is injected into an electrospray ionization Fourier transform mass spectrometry (ESI-FTMS) instrument, the cross-link positions can be localized by fragmentation and mass spectrometry on the 'gas-phase purified' singly internally cross-linked monomer. Our results show that reaction of ubiquitin with the homobifunctional lysine-lysine cross-linking reagent dissuccinimidyl suberate (DSS) resulted in two cross-links consistent with the known ubiquitin tertiary structure (K6-K11 and K48-K63). Because no protein or peptide chemistry steps are needed, other than the initial cross-linking, this new top down approach appears well suited for high-throughput experiments with multiple cross-linkers and reaction conditions. Published in 2002 by John Wiley & Sons, Ltd.

A top-down approach to protein structure studies using chemical cross-linking and Fourier transform mass spectrometry.

In a preliminary communication we described a top-down approach to the determination of chemical cross-link location in proteins using Fourier transform mass spectrometry (FT-MS). We have since extended the approach to use a series of homobifunctional cross-linkers with the same reactive functional groups, but different cross-linker arm lengths. Correlating cross-linking data across a series of related linkers allows the distance constraint derived from a cross-link between two reactive side chains to be determined more accurately and increases the confidence in the assignment of the cross-links. In ubiquitin, there are seven lysines with primary amino groups and the amino terminus. Disuccinimidyl suberate (DSS, cross-linker arm length = 11.4 A), disuccinimidyl glutarate (DSG, cross-linker arm length = 7.5 A) and disuccinimidyl tartrate (DST, cross- linker arm length = 5.8 A) are homobifunctional cross-linking reagents that react specifically with primary amines. Using tandem mass spectrometry (MS/MS) on the singly, internally cross-linked precursor ion of ubiquitin, we found cross-links with DSS and DSG between the amino terminus and Lys 6, between Lys 6 and Lys 11, and between Lys 63 and Lys 48. Using disuccinimidyl tartrate (DST), the shortest cross-linker in the series, only the cross-links between the amino terminus and Lys 6, and between Lys 6 and Lys 11 were observed. The observed cross-links are consistent with the crystal structure of ubiquitin, if the lysine side chains and the amino terminus are assumed to have considerable flexibility. In a separate study, we probed the reactivity of the primary amino groups in ubiquitin using the amino acetylating reagent, N-hydroxy succinimidyl acetate (NHSAc), and a top-down approach to localize the acetylated lysine residues. The reactivity order obtained in that study (M1 approximate, equals K6 approximate, equals K48 approximate, equals K63) > K33 > K11 > (K27, K29), shows that the cross-link first formed in ubiquitin by reaction with DSS and DSG occurs between the most reactive residues.

Computer-assisted mutagenesis of ecotin to engineer its secondary binding site for urokinase inhibition.

Inhibitors of urokinase-type plasminogen activator (uPA) were selected in vitro from two ecotin phage-display libraries to study the effect on binding of amino acid substitutions at critical positions 108, 110, 112, and 113 within the 100s loop (RNKL, respectively, in wild type ecotin). The first, a focused library, was the result of a computation-assisted approach using the three-dimensional structure of the ecotin-trypsin complex to guide the modeling of amino acid substitutions predicted to increase affinity for uPA. The second, a complete library, allowed for all substitutions at the above identified positions. The consensus sequences selected from the focused, and complete libraries were RRWS and R(R/N)QL, respectively. Inhibition constant determinations showed ecotin variants containing these sequences to be similarly potent (K(i) = 1-2 nm). These substitutions were combined with previously identified substitutions in another critical region of ecotin. One of these combinations (D70R/M84R/RRQL) is the tightest (K(i) = 50 pm) ecotin variant inhibitor of uPA. The blending of combinatorial methods and computer algorithms designed to predict stronger binders has allowed us to obtain protein derivatives that exhibit greatly increased affinity for a predetermined target. This technology can be applied to select for enhanced binding interactions at protein-protein interfaces and accelerate the process of protease inhibitor development.