Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “affinity proteomics”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,114 records · Page 62Linked to original sources

[Creation of functional muteins using phage libraries for pharmacoproteomic-based drug discovery and development of DDS].

Tumor necrosis factor-alpha (TNF-alpha) has been expected to be a promising new antitumor agent, but toxic side effects by the systemic administration of TNF-alpha limit its clinical application. In this study, we attempted to improve the therapeutic potency of TNF-alpha by using our protein-drug innovation systems. Among phage libraries displaying various mutant TNF-alphas, we isolated some lysine-deficient super mutant TNF-alphas, typified by mTNF-alpha-K90R, with higher TNF-receptor affinities and stronger bioactivity in vitro, in spite of the importance of lysine residues for trimer formation and receptor binding. mTNF-alpha-K90R showed more than 10 times stronger in vivo antitumor effects and 1.3 times less toxicity than wild-type TNF-alpha (wTNF-alpha). Site-specifically mono-PEGylated mTNF-alpha-K90R (sp-PEG-mTNF-alpha-K90R) at N-terminus showed higher in vitro bioactivity than unmodified wTNF-alpha, whereas randomly mono-PEGylated wTNF-alpha at a lysine residue (ran-PEG-wTNF-alpha) had less than 6% of the bioactivity of wTNF-alpha. The antitumor therapeutic window of sp-PEG-mTNF-alpha-K90R was extended by about 5 times, 60 times and 18 times compared with those of mTNF-alpha-K90R, wTNF-alpha and ran-PEG-wTNF-alpha, respectively. sp-PEG-mTNF-alpha-K90R may, thus, be a potential systemic anti-tumor therapeutic agent. These data suggested that our fusion protein-drug innovation system composed of a creation system of functional mutant proteins based on phage display technique and a site-specific PEGylation system may open up a new avenue to the optimal protein therapy.

Animals↗

Quantitative phosphoproteomics of vasopressin-sensitive renal cells: regulation of aquaporin-2 phosphorylation at two sites.

Protein phosphorylation plays a key role in vasopressin signaling in the renal-collecting duct. Large-scale identification and quantification of phosphorylation events triggered by vasopressin is desirable to gain a comprehensive systems-level understanding of this process. We carried out phosphoproteomic analysis of rat inner medullary collecting duct cells by using a combination of phosphopeptide enrichment by immobilized metal affinity chromatography and phosphorylation site identification by liquid chromatography-mass spectrometry(n) neutral loss scanning. A total of 714 phosphorylation sites on 223 unique phosphoproteins were identified from inner medullary collecting duct samples treated short-term with either calyculin A or vasopressin. A number of proteins involved in cytoskeletal reorganization, vesicle trafficking, and transcriptional regulation were identified. Previously unidentified phosphorylation sites were found for membrane proteins essential to collecting duct physiology, including eight sites among aquaporin-2 (AQP2), aquaporin-4, and urea transporter isoforms A1 and A3. Through label-free quantification of phosphopeptides, we identified a number of proteins that significantly changed phosphorylation state in response to short-term vasopressin treatment: AQP2, Bclaf1, LRRC47, Rgl3, and SAFB2. In the presence of vasopressin, AQP2 monophosphorylated at S256 and diphosphorylated AQP2 (pS256/261) increased in abundance, whereas AQP2 monophosphorylated at S261 decreased, raising the possibility that both sites are involved in vasopressin-dependent AQP2 trafficking. This study reveals the practicality of liquid chromatography-mass spectrometry(n) neutral loss scanning for large-scale identification and quantification of protein phosphorylation in the analysis of cell signaling in a native mammalian system.

Amino Acid Sequence↗

Prediction of binding modes for ligands in the cytochromes P450 and other heme-containing proteins.

The cytochromes P450 (P450s) are a family of heme-containing monooxygenase enzymes involved in a variety of functions, including the metabolism of endogenous and exogenous substances in the human body. During lead optimization, and in drug development, many potential drug candidates are rejected because of the affinity they display for drug-metabolising P450s. Recently, crystal structures of human enzymes involved in drug metabolism have been determined, significantly augmenting the prospect of using structure-based design to modulate the binding and metabolizing properties of compounds against P450 proteins. An important step in the application of structure-based metabolic optimization is the accurate prediction of docking modes in heme binding proteins. In this paper we assess the performance of the docking program GOLD at predicting the binding mode of 45 heme-containing complexes. We achieved success rates of 64% and 57% for Chemscore and Goldscore respectively; these success rates are significantly lower than the value of 79% observed with both scoring functions for the full GOLD validation set. Re-parameterization of metal-acceptor interactions and lipophilicity of planar nitrogen atoms in the scoring functions resulted in a significant increase in the percentage of successful dockings against the heme binding proteins (Chemscore 73%, Goldscore 65%). The modified scoring functions will be useful in docking applications on P450 enzymes and other heme binding proteins.

Algorithms↗

Local modeling of global interactome networks.

MOTIVATION: Systems biology requires accurate models of protein complexes, including physical interactions that assemble and regulate these molecular machines. Yeast two-hybrid (Y2H) and affinity-purification/mass-spectrometry (AP-MS) technologies measure different protein-protein relationships, and issues of completeness, sensitivity and specificity fuel debate over which is best for high-throughput 'interactome' data collection. Static graphs currently used to model Y2H and AP-MS data neglect dynamic and spatial aspects of macromolecular complexes and pleiotropic protein function. RESULTS: We apply the local modeling methodology proposed by Scholtens and Gentleman (2004) to two publicly available datasets and demonstrate its uses, interpretation and limitations. Specifically, we use this technology to address four major issues pertaining to protein-protein networks. (1) We motivate the need to move from static global interactome graphs to local protein complex models. (2) We formally show that accurate local interactome models require both Y2H and AP-MS data, even in idealized situations. (3) We briefly discuss experimental design issues and how bait selection affects interpretability of results. (4) We point to the implications of local modeling for systems biology including functional annotation, new complex prediction, pathway interactivity and coordination with gene-expression data. AVAILABILITY: The local modeling algorithm and all protein complex estimates reported here can be found in the R package apComplex, available at http://www.bioconductor.org CONTACT: dscholtens@northwestern.edu SUPPLEMENTARY INFORMATION: http://daisy.prevmed.northwestern.edu/~denise/pubs/LocalModeling

Algorithms↗

Protein-protein interaction analysis using an affinity peptide tag and hydrophilic polystyrene plate.

A sandwich ELISA method using peptide tags showing a specific affinity to a hydrophilic polystyrene surface (PS-tags), PS 19 composed of RAFIASRRIKRP and KPS19R10 of KRAFIASRRIRRP and a hydrophilic polystyrene (phi-PS) plate was used to analyze protein-protein interactions. An Escherichia coli cysteine synthase complex, in which serine acetyltransferase (SAT) interacts with O-acetylserine sulfhydrylase-A (OASS) was used as a model system. When the interaction was detected by the conventional sandwich ELISA method using a hydrophobic polystyrene (pho-PS) plate, for the exclusive use of ELISA, the signal intensity was barely detectable due to conformational change of the ligand protein, OASS in the adsorbed state. On the contrary, when OASS, genetically fused with PS19 (OASS-PS19) or chemically conjugated with KPS19R10 (OASS-KPS19R10), was immobilized on the phi-PS plate, a high signal intensity was detected. Furthermore, by applying the two-step sandwich ELISA, in which OASS-PS19 or OASS-KPS19R10 formed a complex with SAT in the blocking solution before immobilization on the phi-PS plate, the signal intensity was further increased with a much shorter operational time, because SAT in the blocking solution formed a complex with OASS-PS19 or OASS-KPS19R10 without any steric hindrance.

Enzyme-Linked Immunosorbent Assay↗

Identification of protein-protein interactions using in vivo cross-linking and mass spectrometry.

The purification of protein complexes can be accomplished by different types of affinity chromatography. In a typical immunoaffinity experiment, protein complexes are captured from a cell lysate by an immobilized antibody that recognizes an epitope on one of the known components of the complex. After extensive washing to remove unspecifically bound proteins, the complexes are eluted and analyzed by mass spectrometry (MS). Transient complexes, which are characterized by high dissociation constants, are typically lost by this approach. In the present study, we describe a novel method for identifying transient protein-protein interactions using in vivo cross-linking and MS-based protein identification. Live cells are treated with formaldehyde, which rapidly permeates the cell membrane and generates protein-protein cross-links. Proteins cross-linked to a Myc-tagged protein of interest are copurified by immunoaffinity chromatography and subjected to a procedure which dissociates the cross-linked complexes. After separation by SDS-PAGE, proteins are identified by tandem mass spectrometry. Application of this method enabled the identification of numerous proteins that copurified with a constitutively active form of M-Ras (M-Ras(Q71L)). Among these, we identified the RasGAP-related protein IQGAP1 to be a novel interaction partner of M-Ras(Q71L). This method is applicable to many proteins and will aid in the study of protein-protein interactions.

Animals↗

Environmentally modulated phosphoproteome of photosynthetic membranes in the green alga Chlamydomonas reinhardtii.

Mapping of in vivo protein phosphorylation sites in photosynthetic membranes of the green alga Chlamydomonas reinhardtii revealed that the major environmentally dependent changes in phosphorylation are clustered at the interface between the photosystem II (PSII) core and its light-harvesting antennae (LHCII). The photosynthetic membranes that were isolated form the algal cells exposed to four distinct environmental conditions affecting photosynthesis: (i) dark aerobic, corresponding to photosynthetic State 1; (ii) dark under nitrogen atmosphere, corresponding to photosynthetic State 2; (iii) moderate light; and (iv) high light. The surface-exposed phosphorylated peptides were cleaved from the membrane by trypsin, methyl-esterified, enriched by immobilized metal affinity chromatography, and sequenced by nanospray-quadrupole time-of-flight mass spectrometry. A total of 19 in vivo phosphorylation sites were mapped in the proteins corresponding to 15 genes in C. reinhardtii. Amino-terminal acetylation of seven proteins was concomitantly determined. Sequenced amino termini of six mature LHCII proteins differed from the predicted ones. The State 1-to-State 2 transition induced phosphorylation of the PSII core components D2 and PsbR and quadruple phosphorylation of a minor LHCII antennae subunit, CP29, as well as phosphorylation of constituents of a major LHCII complex, Lhcbm1 and Lhcbm10. Exposure of the algal cells to either moderate or high light caused additional phosphorylation of the D1 and CP43 proteins of the PSII core. The high light treatment led to specific hyperphosphorylation of CP29 at seven distinct residues, phosphorylation of another minor LHCII constituent, CP26, at a single threonine, and double phosphorylation of additional subunits of a major LHCII complex including Lhcbm4, Lhcbm6, Lhcbm9, and Lhcbm11. Environmentally induced protein phosphorylation at the interface of PSII core and the associated antenna proteins, particularly multiple differential phosphorylations of CP29 linker protein, suggests the mechanisms for control of photosynthetic state transitions and for LHCII uncoupling from PSII under high light stress to allow thermal energy dissipation.

Algal Proteins↗

Pro-apoptotic proteins released from the mitochondria regulate the protein composition and caspase-processing activity of the native Apaf-1/caspase-9 apoptosome complex.

The apoptosome is a large caspase-activating ( approximately 700-1400 kDa) complex, which is assembled from Apaf-1 and caspase-9 when cytochrome c is released during mitochondrial-dependent apoptotic cell death. Apaf-1 the core scaffold protein is approximately 135 kDa and contains CARD (caspase recruitment domain), CED-4, and multiple (13) WD40 repeat domains, which can potentially interact with a variety of unknown regulatory proteins. To identify such proteins we activated THP.1 lysates with dATP/cytochrome c and used sucrose density centrifugation and affinity-based methods to purify the apoptosome for analysis by MALDI-TOF mass spectrometry. First, we used a glutathione S-transferase (GST) fusion protein (GST-casp9(1-130)) containing the CARD domain of caspase-9-(1-130), which binds to the CARD domain of Apaf-1 when it is in the apoptosome and blocks recruitment/activation of caspase-9. This affinity-purified apoptosome complex contained only Apaf-1XL and GST-casp9(1-130), demonstrating that the WD40 and CED-4 domains of Apaf-1 do not stably bind other cytosolic proteins. Next we used a monoclonal antibody to caspase-9 to immunopurify the native active apoptosome complex from cell lysates, containing negligible levels of cytochrome c, second mitochondria-derived activator of caspase (Smac), or Omi/HtrA2. This apoptosome complex exhibited low caspase-processing activity and contained four stably associated proteins, namely Apaf-1, pro-p35/34 forms of caspase-9, pro-p20 forms of caspase-3, X-linked inhibitor of apoptosis (XIAP), and cytochrome c, which was only bound transiently to the complex. However, in lysates containing Smac and Omi/HtrA2, the caspase-processing activity of the purified apoptosome complex increased 6-8-fold and contained only Apaf-1 and the p35/p34-processed subunits of caspase-9. During apoptosis, Smac, Omi/HtrA2, and cytochrome c are released simultaneously from mitochondria, and thus it is likely that the functional apoptosome complex in apoptotic cells consists primarily of Apaf-1 and processed caspase-9.

Amino Acid Sequence↗

Comparative serum glycoproteomics using lectin selected sialic acid glycoproteins with mass spectrometric analysis: application to pancreatic cancer serum.

A strategy is developed in this study for identifying sialylated glycoprotein markers in human cancer serum. This method consists of three steps: lectin affinity selection, a liquid separation and characterization of the glycoprotein markers using mass spectrometry. In this work, we use three different lectins (Wheat Germ Agglutinin, (WGA) Elderberry lectin,(SNA), Maackia amurensis lectin, (MAL)) to extract sialylated glycoproteins from normal and cancer serum. Twelve highly abundant proteins are depleted from the serum using an IgY-12 antibody column. The use of the different lectin columns allows one to monitor the distribution of alpha(2,3) and alpha(2,6) linkage type sialylation in cancer serum vs that in normal samples. Extracted glycoproteins are fractionated using NPS-RP-HPLC followed by SDS-PAGE. Target glycoproteins are characterized further using mass spectrometry to elucidate the carbohydrate structure and glycosylation site. We applied this approach to the analysis of sialylated glycoproteins in pancreatic cancer serum. Approximately 130 sialylated glycoproteins are identified using microLC-MS/MS. Sialylated plasma protease C1 inhibitor is identified to be down-regulated in cancer serum. Changes in glycosylation sites in cancer serum are also observed by glycopeptide mapping using microLC-ESI-TOF-MS where the N83 glycosylation of alpha1-antitrypsin is down regulated. In addition, the glycan structures of the altered proteins are assigned using MALDI-QIT-MS. This strategy offers the ability to quantitatively analyze changes in glycoprotein abundance and detect the extent of glycosylation alteration as well as the carbohydrate structure that correlate with cancer.

Case-Control Studies↗

Prolonged selection in aerobic, glucose-limited chemostat cultures of Saccharomyces cerevisiae causes a partial loss of glycolytic capacity.

Prolonged cultivation of Saccharomyces cerevisiae in aerobic, glucose-limited chemostat cultures (dilution rate, 0.10 h(-1)) resulted in a progressive decrease of the residual glucose concentration (from 20 to 8 mg l(-1) after 200 generations). This increase in the affinity for glucose was accompanied by a fivefold decrease of fermentative capacity, and changes in cellular morphology. These phenotypic changes were retained when single-cell isolates from prolonged cultures were used to inoculate fresh chemostat cultures, indicating that genetic changes were involved. Kinetic analysis of glucose transport in an 'evolved' strain revealed a decreased Km, while Vmax was slightly increased relative to the parental strain. Apparently, fermentative capacity in the evolved strain was not controlled by glucose uptake. Instead, enzyme assays in cell extracts of the evolved strain revealed strongly decreased capacities of enzymes in the lower part of glycolysis. This decrease was corroborated by genome-wide transcriptome analysis using DNA microarrays. In aerobic batch cultures on 20 g glucose l(-1), the specific growth rate of the evolved strain was lower than that of the parental strain (0.28 and 0.37 h(-1), respectively). Instead of the characteristic instantaneous production of ethanol that is observed when aerobic, glucose-limited cultures of wild-type S. cerevisiae are exposed to excess glucose, the evolved strain exhibited a delay of approximately 90 min before aerobic ethanol formation set in. This study demonstrates that the effects of selection in glucose-limited chemostat cultures extend beyond glucose-transport kinetics. Although extensive physiological analysis offered insight into the underlying cellular processes, the evolutionary 'driving force' for several of the observed changes remains to be elucidated.

Aerobiosis↗

LC-MS/MS analysis of apical and basolateral plasma membranes of rat renal collecting duct cells.

We used biotinylation and streptavidin affinity chromatography to label and enrich proteins from apical and basolateral membranes of rat kidney inner medullary collecting ducts (IMCDs) prior to LC-MS/MS protein identification. To enrich apical membrane proteins and bound peripheral membrane proteins, IMCDs were perfusion-labeled with primary amine-reactive biotinylation reagents at 2 degrees C using a double barreled pipette. The perfusion-biotinylated proteins and proteins bound to them were isolated with CaptAvidin-agarose beads, separated with SDS-PAGE, and sliced into continuous gel pieces for LC-MS/MS protein identification (LTQ, Thermo Electron Corp.). 17 integral and glycosylphosphatidylinositol (GPI)-linked membrane proteins and 44 non-integral membrane proteins were identified. Immunofluorescence confocal microscopy confirmed ACVRL1, H(+)/K(+)-ATPase alpha1, NHE2, and TauT expression in the IMCDs. Basement membrane and basolateral membrane proteins were biotinylated via incubation of IMCD suspensions with biotinylation reagents on ice. 23 integral and GPI-linked membrane proteins and 134 non-integral membrane proteins were identified. Analyses of non-integral membrane proteins preferentially identified in the perfusion-biotinylated and not in the incubation-biotinylated IMCDs revealed protein kinases, scaffold proteins, SNARE proteins, motor proteins, small GTP-binding proteins, and related proteins that may be involved in vasopressin-stimulated AQP2, UT-A1, and ENaC regulation. A World Wide Web-accessible database was constructed of 222 membrane proteins (integral and GPI-linked) from this study and prior studies.

Animals↗

Study of serum haptoglobin and its glycoforms in the diagnosis of hepatocellular carcinoma: a glycoproteomic approach.

Increased serum haptoglobin concentration and changes in its glycosylation have been reported in certain cancer types. Information for hepatocellular carcinoma (HCC) has not yet been available. In this study, we aimed to carry out a systematic analysis of serum concentrations of haptoglobin (Hp) and its glycoforms in the patients with HCC and noncancer patients only with chronic liver diseases (CLD) and to examine their clinical values. This study was divided into two major parts, (1) measurement of serum Hp concentration, and investigation of its value in the diagnosis of HCC, and (2) quantitative analysis of Hp glycoforms with alpha-2,6-sialylation and/or alpha-1,6-fucosylation by using lectin affinity purification and 2D gel electrophoresis and investigation of their relationships with tumor stage. The concentrations of serum Hp in HCC patients were significantly higher than those in noncancer patients with CLD. With the use of serum concentrations of Hp and alpha-fetoprotein, a logistic regression (LR) model was developed from the training data set and used to classify the validation cases. At a specificity of 95%, the sensitivity for HCC detection was 79%. Comparing serum concentrations of alpha-2,6-sialylated Hp (S-Hp) and alpha-1,6-fucosylated Hp (F-Hp) between HCC and CLD patients suggests that purification of S-Hp and F-Hp could enrich the glycosylation variants associated with HCC. 2D gel analysis of S-Hp and F-Hp identified a total of 18 glycoforms. A unique pattern of Hp glycoforms comprising both hypersialylated fucosylated and hyposialylated fucosylated species was found in the HCC patients. Serum concentrations of these glycoproteins were significantly higher in the patients with advanced tumors, suggesting their tumor-specific nature. We have shown that serum Hp is a potential biomarker in the diagnosis of HCC. The combined use of Hp and AFP could greatly improve the diagnostic accuracy. A unique pattern of Hp glycoforms with altered sialylation and fucosylation is specific to HCC and associated tumor progression.

Adult↗

Identification of beta-tubulin isoforms as tumor antigens in neuroblastoma.

There is currently substantial interest in the identification of human tumor antigens for diagnosis and immunotherapy of cancer. We have implemented a proteomic approach for the identification of tumor proteins that elicit a humoral response in cancer patients, which we have applied to neuroblastoma. Proteins from neuroblastoma tumors and cell lines were separated by two-dimensional PAGE and transferred to poly(vinylidene difluoride) membranes. Sera from 23 newly diagnosed patients with neuroblastoma, from 12 newly diagnosed children with other solid tumors, and from 13 normal individuals were screened for IgG and IgM autoantibodies against neuroblastoma proteins by means of Western blot analysis. Sera from 11 patients with neuroblastoma and from 1 patient with a primitive neuroectodermal tumor, but none of the other controls exhibited IgG-based reactivity against a protein constellation with an estimated Mr 50,000. NH2-terminal sequence and mass spectrometric analysis identified the major constituents of this constellation as beta-tubulin isoforms I and III. The IgG antibodies were additionally characterized to be of the subclass IgG1. Neuroblastoma patient sera that contained anti-beta-tubulin IgG antibodies also contained IgM antibodies specific against the full-length beta-tubulin molecule and against COOH-terminal beta-tubulin cleavage products. Neuroblastoma patient sera that reacted with beta-tubulin I and III isoforms in neuroblastoma tissues did not react with beta-tubulin I and III isoforms found in normal brain tissue. Our findings indicate the occurrence of beta-tubulin peptides in neuroblastoma, which are immunogenic. The occurrence of immunogenic peptides in neuroblastoma may have utility in diagnosis and in immunotherapy of this aggressive childhood tumor.

Adolescent↗

Global landscape of protein complexes in the yeast Saccharomyces cerevisiae.

Identification of protein-protein interactions often provides insight into protein function, and many cellular processes are performed by stable protein complexes. We used tandem affinity purification to process 4,562 different tagged proteins of the yeast Saccharomyces cerevisiae. Each preparation was analysed by both matrix-assisted laser desorption/ionization-time of flight mass spectrometry and liquid chromatography tandem mass spectrometry to increase coverage and accuracy. Machine learning was used to integrate the mass spectrometry scores and assign probabilities to the protein-protein interactions. Among 4,087 different proteins identified with high confidence by mass spectrometry from 2,357 successful purifications, our core data set (median precision of 0.69) comprises 7,123 protein-protein interactions involving 2,708 proteins. A Markov clustering algorithm organized these interactions into 547 protein complexes averaging 4.9 subunits per complex, about half of them absent from the MIPS database, as well as 429 additional interactions between pairs of complexes. The data (all of which are available online) will help future studies on individual proteins as well as functional genomics and systems biology.

Biological Evolution↗

Quantitative peptidomics of pituitary glands from mice deficient in copper transport.

We previously described a method of quantitating levels of peptides in Cpe(fat)/Cpe(fat) mice using affinity chromatography to isolate peptide-processing intermediates and differential isotopic labeling/mass spectrometry. In the present study, we compared two different isotopic labels, acetic anhydride and succinic anhydride for detection and quantitation of peptides in wild type mice. As previously found for acetic anhydride, succinic anhydride efficiently labels all primary amines in various peptides. Of these two reagents, succinic anhydride provides better resolution between the heavy and light peaks of the labelled peptides due to a greater mass difference between the deuterated (heavy) and non-deuterated (light) form of this label (4 Da for succinate, 3 Da for acetate). Using succinic anhydride labeling, the accuracy of measuring 1:1 and 1:2 ratios of peptides in pituitary extracts was within 5% of the theoretical value for most peptides. The accuracy with succinic anhydride is comparable to the accuracy of acetic anhydride and more peptides could be detected and quantitated with succinic anhydride. The two labels were then used to examine pituitary peptides in mice with a defect in copper transport (Atp7a mice) vs wild type mice. Using succinic anhydride, 13 peptides could be detected, 12 of which matched the theoretical mass of known pituitary peptides. Five of the six peptides which contain C-terminal amide groups were significantly decreased in the Atp7a mice relative to wild type mice, whereas only one non-amidated peptide was significantly decreased in Atp7a mice. With acetic anhydride, only five peptides could be quantitated. The three peptides which contain C-terminal amide groups were decreased approximately 30% in the Atp7a mice. The selective decrease in amidated peptides in Atp7a mice is consistent with the copper-requirement of the enzyme that forms C-terminal amides.

Acetic Anhydrides↗

Receptor rigidity and ligand mobility in trypsin-ligand complexes.

The trypsin-like serine proteases comprise a structurally similar family of proteins with a wide diversity of biological functions. Members of this family play roles in digestion, hemostasis, immune responses, and cancer metastasis. Bovine trypsin is an archetypical member of this family that has been extensively characterized both functionally and structurally, and that preferentially hydrolyzes Arg/Lys-Xaa peptide bonds. We have used molecular dynamics (MD) simulations to study bovine trypsin complexed with the two noncovalent small-molecule ligands, benzamidine and tranylcypromine, that have the same hydrogen-bond donating moieties as Arg and Lys side-chains, respectively. Multiple (10) simulations ranging from 1 ns to 2.2 ns, with explicit water molecules and periodic boundary conditions, were performed. The simulations reveal that the trypsin binding pocket residues are relatively rigid regardless of whether there is no ligand, a high-affinity ligand (benzamidine), or a low-affinity ligand (tranylcypromine). The thermal average of the conformations sampled by benzamidine bound to trypsin is planar and consistent with the planar internal geometry of the benzamidine crystallographic model coordinates. However, the most probable bound benzamidine conformations are +/-25 degrees out of plane, implying that the observed X-ray electron density represents an average of densities from two mirror symmetric, nonplanar conformations. Solvated benzamidine has free energy minima at +/-45 degrees , and the induction of a more planar geometry upon binding is associated with approximately 1 kcal/mol of intramolecular strain. Tranylcypromine's hydrogen-bonding pattern in the MD differs substantially from that inferred from the X-ray electron density. Early in simulations of this system, tranylcypromine adopts an alternative binding conformation, changing from the crystallographic conformation, with a direct hydrogen bond between its amino moiety and the backbone oxygen of Gly219, to one having a bridging water molecule. This result is consistently seen with the CHARMM22, Amber, or OPLS-AA force fields. The trypsin-tranylcypromine hydrogen-bonding pattern observed in the simulations also occurs as the crystallographic binding mode of the Lys15 side-chain of bovine pancreatic trypsin inhibitor bound to trypsin. In this latter cocrystal, a bridging crystallographic water does reside between the side-chain's amino group and the trypsin Gly219 backbone oxygen. Furthermore, the trypsin-tranylcypromine simulations sample two different stable noncrystallographic binding poses. These data suggest that some of the electron density ascribed to tranylcypromine in the X-ray model is rather due to a bound water molecule, and that multiple tranylcypromine binding conformations (crystallographic disorder) may be the cause of ambiguous electron density. The combined trypsin-benzamidine and trypsin- tranylcypromine results highlight the ability of simulations to augment protein-ligand complex structural data by deconvoluting the effects of thermal and structural averaging, and by finding energetically optimal ligand and bound water positions for weakly bound ligands.

Animals↗