Probing membrane protein unfolding with pulse proteolysis

Steric trapping strategy for studying the folding of helical membrane proteins

Elucidating the folding energy landscape of membrane proteins is essential to the understanding of the proteins' stabilizing forces, folding mechanisms, biogenesis, and quality control. This is not a trivial task because the reversible control of folding is inherently difficult in a lipid bilayer environment. Recently, novel methods have been developed, each of which has a unique strength in investigating specific aspects of membrane protein folding. Among such methods, steric trapping is a versatile strategy allowing a reversible control of membrane protein folding with minimal perturbation of native protein-water and protein-lipid interactions. In a nutshell, steric trapping exploits the coupling of spontaneous denaturation of a doubly biotinylated protein to the simultaneous binding of bulky monovalent streptavidin molecules. This strategy has been evolved to investigate key elements of membrane protein folding such as thermodynamic stability, spontaneous denaturation rates, conformational features of the denatured states, and cooperativity of stabilizing interactions. In this review, we describe the critical methodological advancement, limitation, and outlook of the steric trapping strategy.

Fundamentals to function: Quantitative and scalable approaches for measuring protein stability

Folding a linear chain of amino acids into a three-dimensional protein is a complex physical process that ultimately confers an impressive range of diverse functions. Although recent advances have driven significant progress in predicting three-dimensional protein structures from sequence, proteins are not static molecules. Rather, they exist as complex conformational ensembles defined by energy landscapes spanning the space of sequence and conditions. Quantitatively mapping the physical parameters that dictate these landscapes and protein stability is therefore critical to develop models that are capable of predicting how mutations alter function of proteins in disease and informing the design of proteins with desired functions. Here, we review the approaches that are used to quantify protein stability at a variety of scales, from returning multiple thermodynamic and kinetic measurements for a single protein sequence to yielding indirect insights into folding across a vast sequence space. The physical parameters derived from these approaches will provide a foundation for models that extend beyond the structural prediction to capture the complexity of conformational ensembles and, ultimately, their function.

Identify the interactions between phytochemicals and proteins in the complicated food matrix

Phytochemicals usually mix with food proteins in our regular diet. Unexpected interactions may lead to changes in bioaccessibility, bioactivity, and bioavailability of phytochemicals. However, our understanding of these interactions between phytochemical and food proteins is limited because of the experimental restrictions. Here, we used pulse-proteolysis to conduct the unfolding equilibrium and dose-dependent experiments on the food proteins for the first time. The interaction between epigallocatechin gallate (EGCG) and caseins was identified in the complicated food matrix, whole milk. Another food proteome, soymilk, was also optimized for identifying the binding targets of EGCG and caffeine. Among the identified interactions, the mixing of milk with coffee generates the most prominent masking effect of 46.61 ± 3.86% relative to the calculated antioxidant capacity. Our results demonstrated that pulse proteolysis is applicable for identifying the interactions between phytochemicals and proteins in the complicated food matrix.

Solid-state NMR spectroscopy based atomistic view of a membrane protein unfolding pathway

Membrane protein folding, structure, and function strongly depend on a cell membrane environment, yet detailed characterization of folding within a lipid bilayer is challenging. Studies of reversible unfolding yield valuable information on the energetics of folding and on the hierarchy of interactions contributing to protein stability. Here, we devise a methodology that combines hydrogen-deuterium (H/D) exchange and solid-state NMR (SSNMR) to follow membrane protein unfolding in lipid membranes at atomic resolution through detecting changes in the protein water-accessible surface, and concurrently monitoring the reversibility of unfolding. We obtain atomistic description of the reversible part of a thermally induced unfolding pathway of a seven-helical photoreceptor. The pathway is visualized through SSNMR-detected snapshots of H/D exchange patterns as a function of temperature, revealing the unfolding intermediate and its stabilizing factors. Our approach is transferable to other membrane proteins, and opens additional ways to characterize their unfolding and stabilizing interactions with atomic resolution.

Kinetic stability of membrane proteins

Although membrane proteins constitute an important class of biomolecules involved in key cellular processes, study of the thermodynamic and kinetic stability of their structures is far behind that of soluble proteins. It is known that many membrane proteins become unstable when removed by detergent extraction from the lipid environment. In addition, most of them undergo irreversible denaturation, even under mild experimental conditions. This process was found to be associated with partial unfolding of the polypeptide chain exposing hydrophobic regions to water, and it was proposed that the formation of kinetically trapped conformations could be involved. In this review, we will describe some of the efforts toward understanding the irreversible inactivation of membrane proteins. Furthermore, its modulation by phospholipids, ligands, and temperature will be herein discussed.

Time, space, and disorder in the expanding proteome universe

Proteins are highly dynamic entities. Their myriad functions require specific structures, but proteins' dynamic nature ranges all the way from the local mobility of their amino acid constituents to mobility within and well beyond single cells. A truly comprehensive view of the dynamic structural proteome includes: (i) alternative sequences, (ii) alternative conformations, (iii) alternative interactions with a range of biomolecules, (iv) cellular localizations, (v) alternative behaviors in different cell types. While these aspects have traditionally been explored one protein at a time, we highlight recently emerging global approaches that accelerate comprehensive insights into these facets of the dynamic nature of protein structure. Computational tools that integrate and expand on multiple orthogonal data types promise to enable the transition from a disjointed list of static snapshots to a structurally explicit understanding of the dynamics of cellular mechanisms.

Bioanalytical approaches to assess the proteolytic stability of therapeutic fusion proteins

Therapeutic fusion proteins (TFPs) are designed to improve the therapeutic profile of an endogenous protein or protein fragment with a limited dose frequency providing the desired pharmacological activity in vivo. Fusion of a therapeutic protein to a half-life extension or targeting domain can improve the disposition of the molecule or introduce a novel mechanism of action. Prolonged exposure and altered biodistribution of an endogenous protein through fusion technology increases the potential for local protein unfolding during circulation increasing the chance for partial proteolysis of the therapeutic domain. Characterizing the proteolytic liabilities of a TFP can guide engineering efforts to inhibit or hinder partial proteolysis. This review focuses on considerations and techniques for evaluating the stability of a TFP both in vivo and in vitro.

Residue-Dependent Thermodynamic Cost and Barrel Plasticity Balances Activity in the PhoPQ-Activated Enzyme PagP of Salmonella typhimurium

PagP is an eight-stranded transmembrane β-barrel enzyme indispensable for lipid A palmitoylation in Gram-negative bacteria. The severity of infection by pathogens, including Salmonella, Legionella, and Bordetella, and resistance to antimicrobial peptides, relies on lipid A remodeling by PagP, rendering PagP a sought-after drug target. Despite a conserved sequence, more robust palmitoylation of lipid A is observed in Salmonella typhimurium compared to Escherichia coli, a possible consequence of the differential regulation of PagP expression and/or specific activity. Work here identifies molecular signatures that demarcate thermodynamic stability and variances in catalytic efficiency between S. typhimurium (PagP-St) and E. coli (PagP-Ec) transmembrane PagP barrel variants. We demonstrate that Salmonella PagP displays a 2-fold destabilization of the barrel, while achieving 15-20 magnitude higher lipase efficiency, through subtle alterations of lipid-facing residues distal from the active site. We find that catalytic properties of these homologues are retained across different lipid environments such as micelles, vesicles, and natural extracts. By comparing thermodynamic stability with activity of selectively designed mutants, we conclude that activity-stability trade-offs can be influenced by factors secluded from the catalytic region. Our results provide a compelling correlation of the primary protein structure with enzymatic activity, barrel thermodynamic stability, and scaffold plasticity. Our analysis can open avenues for the development of potent pharmaceuticals against salmonellosis.

Bacteriorhodopsin folds through a poorly organized transition state

The folding mechanisms of helical membrane proteins remain largely uncharted. Here we characterize the kinetics of bacteriorhodopsin folding and employ φ-value analysis to explore the folding transition state. First, we developed and confirmed a kinetic model that allowed us to assess the rate of folding from SDS-denatured bacteriorhodopsin (bRU) and provides accurate thermodynamic information even under influence of retinal hydrolysis. Next, we obtained reliable φ-values for 16 mutants of bacteriorhodopsin with good coverage across the protein. Every φ-value was less than 0.4, indicating the transition state is not uniquely structured. We suggest that the transition state is a loosely organized ensemble of conformations.

The safety dance: biophysics of membrane protein folding and misfolding in a cellular context

Most biological processes require the production and degradation of proteins, a task that weighs heavily on the cell. Mutations that compromise the conformational stability of proteins place both specific and general burdens on cellular protein homeostasis (proteostasis) in ways that contribute to numerous diseases. Efforts to elucidate the chain of molecular events responsible for diseases of protein folding address one of the foremost challenges in biomedical science. However, relatively little is known about the processes by which mutations prompt the misfolding of α-helical membrane proteins, which rely on an intricate network of cellular machinery to acquire and maintain their functional structures within cellular membranes. In this review, we summarize the current understanding of the physical principles that guide membrane protein biogenesis and folding in the context of mammalian cells. Additionally, we explore how pathogenic mutations that influence biogenesis may differ from those that disrupt folding and assembly, as well as how this may relate to disease mechanisms and therapeutic intervention. These perspectives indicate an imperative for the use of information from structural, cellular, and biochemical studies of membrane proteins in the design of novel therapeutics and in personalized medicine.

Energetics-based methods for protein folding and stability measurements

Over the past 15 years, a series of energetics-based techniques have been developed for the thermodynamic analysis of protein folding and stability. These techniques include Stability of Unpurified Proteins from Rates of amide H/D Exchange (SUPREX), pulse proteolysis, Stability of Proteins from Rates of Oxidation (SPROX), slow histidine H/D exchange, lysine amidination, and quantitative cysteine reactivity (QCR). The above techniques, which are the subject of this review, all utilize chemical or enzymatic modification reactions to probe the chemical denaturant- or temperature-induced equilibrium unfolding properties of proteins and protein-ligand complexes. They employ various mass spectrometry-, sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE)-, and optical spectroscopy-based readouts that are particularly advantageous for high-throughput and in some cases multiplexed analyses. This has created the opportunity to use protein folding and stability measurements in new applications such as in high-throughput screening projects to identify novel protein ligands and in mode-of-action studies to identify protein targets of a particular ligand.

Measuring membrane protein stability under native conditions

The thermodynamic stability of proteins is typically measured at high denaturant concentrations and then extrapolated back to zero denaturant conditions to obtain unfolding free energies under native conditions. For membrane proteins, the extrapolations are fraught with considerable uncertainty as the denaturants may have complex effects on the membrane or micellar structure. We therefore sought to measure stability under native conditions, using a method that does not perturb the properties of the membrane or membrane mimetics. We use a technique called steric trapping to measure the thermodynamic stability of bacteriorhodopsin in bicelles and micelles. We find that bacteriorhodopsin has a high thermodynamic stability, with an unfolding free energy of ∼11 kcal/mol in dimyristoyl phosphatidylcholine bicelles. Nevertheless, the stability is much lower than predicted by extrapolation of measurements made at high denaturant concentrations. We investigated the discrepancy and found that unfolding free energy is not linear with denaturant concentration. Apparently, long extrapolations of helical membrane protein unfolding free energies must be treated with caution. Steric trapping, however, provides a method for making these measurements.

Folding the proteome

Protein folding is an essential prerequisite for protein function and hence cell function. Kinetic and thermodynamic studies of small proteins that refold reversibly were essential for developing our current understanding of the fundamentals of protein folding mechanisms. However, we still lack sufficient understanding to accurately predict protein structures from sequences, or the effects of disease-causing mutations. To date, model proteins selected for folding studies represent only a small fraction of the complexity of the proteome and are unlikely to exhibit the breadth of folding mechanisms used in vivo. We are in urgent need of new methods - both theoretical and experimental - that can quantify the folding behavior of a truly broad set of proteins under in vivo conditions. Such a shift in focus will provide a more comprehensive framework from which to understand the connections between protein folding, the molecular basis of disease, and cell function and evolution.

Reversible folding of human peripheral myelin protein 22, a tetraspan membrane protein

Misfolding of the α-helical membrane protein peripheral myelin protein 22 (PMP22) has been implicated in the pathogenesis of the common neurodegenerative disease known as Charcot-Marie-Tooth disease (CMTD) and also several other related peripheral neuropathies. Emerging evidence suggests that the propensity of PMP22 to misfold in the cell may be due to an intrinsic lack of conformational stability. Therefore, quantitative studies of the conformational equilibrium of PMP22 are needed to gain insight into the molecular basis of CMTD. In this work, we have investigated the folding and unfolding of wild type (WT) human PMP22 in mixed micelles. Both kinetic and thermodynamic measurements demonstrate that the denaturation of PMP22 by n-lauroyl sarcosine (LS) in dodecylphosphocholine (DPC) micelles is reversible. Assessment of the conformational equilibrium indicates that a significant fraction of unfolded PMP22 persists even in the absence of the denaturing detergent. However, we find the stability of PMP22 is increased by glycerol, which facilitates quantitation of thermodynamic parameters. To our knowledge, this work represents the first report of reversible unfolding of a eukaryotic multispan membrane protein. The results indicate that WT PMP22 possesses minimal conformational stability in micelles, which parallels its poor folding efficiency in the endoplasmic reticulum. Folding equilibrium measurements for PMP22 in micelles may provide an approach to assess the effects of cellular metabolites or potential therapeutic agents on its stability. Furthermore, these results pave the way for future investigation of the effects of pathogenic mutations on the conformational equilibrium of PMP22.

Simplified proteomics approach to discover protein-ligand interactions

Identifying targets of biologically active small molecules is an essential but still challenging task in drug research and chemical genetics. Energetics-based target identification is an approach that utilizes the change in the conformational stabilities of proteins upon ligand binding in order to identify target proteins. Different from traditional affinity-based capture approaches, energetics-based methods do not require any labeling or immobilization of the test molecule. Here, we report a surprisingly simple version of energetics-based target identification, which only requires ion exchange chromatography, SDS PAGE, and minimal use of mass spectrometry. The complexity of a proteome is reduced through fractionation by ion exchange chromatography. Urea-induced unfolding of proteins in each fraction is then monitored by the significant increase in proteolytic susceptibility upon unfolding in the presence and the absence of a ligand. Proteins showing a different degree of unfolding with the ligand are identified by SDS PAGE followed by mass spectrometry. Using this approach, we identified ATP-binding proteins in the Escherichia coli proteome. In addition to known ATP-binding proteins, we also identified a number of proteins that were not previously known to interact with ATP. To validate one such finding, we cloned and purified phosphoglyceromutase, which was not previously known to bind ATP, and confirmed that ATP indeed stabilizes this protein. The combination of fractionation and pulse proteolysis offers an opportunity to investigate protein-drug or protein-metabolite interactions on a proteomic scale with minimal instrumentation and without modification of a molecule of interest.

Cystathionine beta-synthase mutants exhibit changes in protein unfolding: conformational analysis of misfolded variants in crude cell extracts

Protein misfolding has been proposed to be a common pathogenic mechanism in many inborn errors of metabolism including cystathionine β-synthase (CBS) deficiency. In this work, we describe the structural properties of nine CBS mutants that represent a common molecular pathology in the CBS gene. Using thermolysin in two proteolytic techniques, we examined conformation of these mutants directly in crude cell extracts after expression in E. coli. Proteolysis with thermolysin under native conditions appeared to be a useful technique even for very unstable mutant proteins, whereas pulse proteolysis in a urea gradient had limited values for the study of the majority of CBS mutants due to their instability. Mutants in the active core had either slightly increased unfolding (p.A114V, p.E302K and p.G307S) or extensive unfolding with decreased stability (p.H65R, p.T191M, p.I278T and p.R369C). The extent of the unfolding inversely correlated with the previously determined degree of tetrameric assembly and with the catalytic activity. In contrast, mutants bearing aminoacid substitutions in the C-terminal regulatory domain (p.R439Q and p.D444N) had increased global stability with decreased flexibility. This study shows that proteolytic techniques can reveal conformational abnormalities even for CBS mutants that have activity and/or a degree of assembly similar to the wild-type enzyme. We present here a methodological strategy that may be used in cell lysates to evaluate properties of proteins that tend to misfold and aggregate and that may be important for conformational studies of disease-causing mutations in the field of inborn errors of metabolism.

Transmembrane and extramembrane contributions to membrane protein thermal stability: studies with the NaChBac sodium channel

The thermal stabilities of the extramembranous and transmembranous regions of the bacterial voltage-gated sodium channel NaChBac have been characterised using thermal-melt synchrotron radiation circular dichroism (SRCD) spectroscopy. A series of constructs, ranging from the full-length protein containing both the C-terminal cytoplasmic and the transmembranous domains, to proteins with decreasing amounts of the cytoplasmic domain, were examined in order to separately define the roles of these two types of domains in the stability and processes of unfolding of a membrane protein. The sensitivity of the SRCD measurements over a wide range of wavelengths and temperatures has meant that subtle but reproducible conformational changes could be detected with accuracy. The residues in the C-terminal extramembranous domain were highly susceptible to thermal denaturation, but for the most part the transmembrane residues were not thermally-labile and retained their helical character even at very elevated temperatures. The process of thermal unfolding involved an initial irreversible unfolding of the highly labile distal extramembranous C-terminal helical region, which was accompanied by a reversible unfolding of a small number of helical residues in the transmembrane domain. This was then followed by the irreversible unfolding of a limited number of additional transmembrane helical residues at greatly elevated temperatures. Hence this study has been able to determine the different contributions and roles of the transmembrane and extramembrane residues in the processes of thermal denaturation of this multipass integral membrane protein.

Knot formation in newly translated proteins is spontaneous and accelerated by chaperonins

Topological knots are found in a considerable number of protein structures, but it is not clear how they knot and fold within the cellular environment. We investigated the behavior of knotted protein molecules as they are first synthesized by the ribosome using a cell-free translation system. We found that newly translated knotted proteins can spontaneously self-tie and do not require the assistance of molecular chaperones to fold correctly to their trefoil-knotted structures. This process is slow but efficient, and we found no evidence of misfolded species. A kinetic analysis indicates that the knotting process is rate limiting, occurs post-translationally, and is specifically and significantly (P < 0.001) accelerated by the GroEL-GroES chaperonin complex. This demonstrates a new active mechanism for this molecular chaperone and suggests that chaperonin-catalyzed knotting probably dominates in vivo. These results explain how knotted protein structures have withstood evolutionary pressures despite their topological complexity.

Revisiting the folding kinetics of bacteriorhodopsin

The elucidation of the physical principles that govern the folding and stability of membrane proteins is one of the greatest challenges in protein science. Several insights into the folding of α-helical membrane proteins have come from the investigation of the conformational equilibrium of H. halobium bacteriorhodopsin (bR) in mixed micelles using SDS as a denaturant. In an effort to confirm that folded bR and SDS-denatured bR reach the same conformational equilibrium, we found that bR folding is significantly slower than has been previously known. Interrogation of the effect of the experimental variables on folding kinetics reveals that the rate of folding is dependent not only on the mole fraction of SDS but also on the molar concentrations of mixed micelle components, a variable that was not controlled in the previous study of bR folding kinetics. Moreover, when the molar concentrations of mixed micelle components are fixed at the concentrations commonly employed for bR equilibrium studies, conformational relaxation in the transition zone is slower than hydrolysis of the retinal Schiff base. As a result, the conformational equilibrium between folded bR and SDS-denatured bR cannot be achieved under the conventional condition. Our finding suggests that the molar concentrations of mixed micelle components are important experimental variables in the investigation of the kinetics and thermodynamics of bR folding and should be accounted for to ensure the accurate assessment of the conformational equilibrium of bR without the interference of retinal hydrolysis.

Folding and stability of membrane transport proteins in vitro

Transmembrane transporters are responsible for maintaining a correct internal cellular environment. The inherent flexibility of transporters together with their hydrophobic environment means that they are challenging to study in vitro, but recently significant progress been made. This review will focus on in vitro stability and folding studies of transmembrane alpha helical transporters, including reversible folding systems and thermal denaturation. The successful re-assembly of a small number of ATP binding cassette transporters is also described as this is a significant step forward in terms of understanding the folding and assembly of these more complex, multi-subunit proteins. The studies on transporters discussed here represent substantial advances for membrane protein studies as well as for research into protein folding. The work demonstrates that large flexible hydrophobic proteins are within reach of in vitro folding studies, thus holding promise for furthering knowledge on the structure, function and biogenesis of ubiquitous membrane transporter families. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Spatial structure peculiarities of influenza A virus matrix M1 protein in an acidic solution that simulates the internal lysosomal medium

The structure of the C-terminal domain of the influenza virus A matrix M1 protein, for which X-ray diffraction data were still missing, was studied in acidic solution. Matrix M1 protein was bombarded with thermally-activated tritium atoms, and the resulting intramolecular distribution of the tritium label was analyzed to assess the steric accessibility of the amino acid residues in this protein. This technique revealed that interdomain loops and the C-terminal domain of the protein are the most accessible to labeling with tritium atoms. A model of the spatial arrangement of the C-terminal domain of matrix M1 protein was generated using rosetta software adjusted to the data obtained by tritium planigraphy experiments. This model suggests that the C-terminal domain is an almost flat layer with a three-α-helical structure. To explain the high level of tritium label incorporation into the C-terminal domain of the M1 protein in an acidic solution, we also used independent experimental approaches (CD spectroscopy, limited proteolysis and MALDI-TOF MS analysis of the proteolysis products, dynamic light scattering and analytical ultracentrifugation), as well as multiple computational algorithms, to analyse the intrinsic protein disorder. Taken together, the results obtained in the present study indicate that the C-terminal domain is weakly structured. We hypothesize that the specific 3D structural peculiarities of the M1 protein revealed in acidic pH solution allow the protein greater structural flexibility and enable it to interact effectively with the components of the host cell.

Computational design of membrane proteins

This article reviews the recent successes of computational protein design techniques applied to integral membrane proteins. This emerging area is still handicapped by significant difficulties in the experimental characterization of the stability and structure of the designed proteins. Nevertheless, by focusing on oligomeric complexes of single-span transmembrane (TM) peptides with detectable activity, the computational design of membrane proteins has already produced very exciting results. The 'take-home message' is that optimization of van der Waals packing and hydrogen bonding (both 'canonical' and weak Cα-H⋯O bonds) can produce functional structures of remarkable stability and specificity in the membrane.

Kinetic folding mechanism of an integral membrane protein examined by pulsed oxidative labeling and mass spectrometry

We report the application of pulsed oxidative labeling for deciphering the folding mechanism of a membrane protein. SDS-denatured bacteriorhodopsin (BR) was refolded by mixing with bicelles in the presence of free retinal. At various time points (20 ms to 1 day), the protein was exposed to a microsecond ·OH pulse that induces oxidative modifications at solvent-accessible methionine side chains. The extent of labeling was determined by mass spectrometry. These measurements were complemented by stopped-flow spectroscopy. Major time-dependent changes in solvent accessibility were detected for M20 (helix A) and M118 (helix D). Our kinetic data indicate a sequential folding mechanism, consistent with models previously suggested by others on the basis of optical data. Yet, ·OH labeling provides additional structural insights. An initial folding intermediate I(1) gets populated within 20 ms, concomitantly with formation of helix A. Subsequent structural consolidation leads to a transient species I(2). Noncovalent retinal binding to I(2) induces folding of helix D, thereby generating an intermediate I(R). In the absence of retinal, the latter transition does not take place. Hence, formation of helix D depends on retinal binding, whereas this is not the case for helix A. As the cofactor settles deeper into its binding pocket, a final transient species I(R) is generated. This intermediate converts into native BR within minutes by formation of the retinal-K216 Schiff base linkage. The combination of pulsed covalent labeling and optical spectroscopy employed here should also be suitable for exploring the folding mechanisms of other membrane proteins.

Energetics-based discovery of protein-ligand interactions on a proteomic scale

Biochemical functions of proteins in cells frequently involve interactions with various ligands. Proteomic methods for the identification of proteins that interact with specific ligands such as metabolites, signaling molecules, and drugs are valuable in investigating the regulatory mechanisms of cellular metabolism, annotating proteins with unknown functions, and elucidating pharmacological mechanisms. Here we report an energetics-based target identification method in which target proteins in a cell lysate are identified by exploiting the effect of ligand binding on their stabilities. Urea-induced unfolding of proteins in cell lysates is probed by a short pulse of proteolysis, and the effect of a ligand on the amount of folded protein remaining is monitored on a proteomic scale. As proof of principle, we identified proteins that interact with ATP in the Escherichia coli proteome. Literature and database mining confirmed that a majority of the identified proteins are indeed ATP-binding proteins. Four identified proteins that were previously not known to interact with ATP were cloned and expressed to validate the result. Except for one protein, the effects of ATP on urea-induced unfolding were confirmed. Analyses of the protein sequences and structure models were also employed to predict potential ATP binding sites in the identified proteins. Our results demonstrate that this energetics-based target identification approach is a facile method to identify proteins that interact with specific ligands on a proteomic scale.