Accurate Prediction of Amide Exchange in the Fast Limit Reveals Thrombin Allostery

Amide hydrogen/deuterium exchange mass spectrometry (HDXMS) of proteins has become extremely popular for identifying ligand-binding sites, protein-protein interactions, intrinsic disorder, and allosteric changes upon protein modification. Such phenomena are revealed when amide exchange is measured in the fast limit, that is, within a few minutes of exchange in deuterated buffer. The HDXMS data have a resolution of the length of peptides and are difficult to interpret because many different phenomena lead to changes in hydrogen/deuterium exchange. We present a quantitative analysis of accelerated molecular dynamics simulations that provides impressive agreement with peptide-length HDXMS data. Comparative analysis of thrombin and a single-point mutant reveals that the simulation analysis can distinguish the subtle differences in exchange due to mutation. In addition, the results provide a deeper understanding of the underlying changes in dynamics revealed by the HDXMS that extend far from the site of mutation.     

Production of large quantities of isotopically labeled protein in Pichia pastoris by fermentation

Heterologous expression in Pichia pastoris has many of the advantages of eukaryotic expression, proper folding and disulfide bond formation, glycosylation, and secretion. Contrary to other eukaryotic systems, protein production from P.pastoris occurs in simple minimal defined media making this system attractive for production of labeled proteins for NMR analysis. P.pastoris is therefore the expression system of choice for NMR of proteins that cannot be refolded from inclusion bodies or that require post-translational modifications for proper folding or function. The yield of expressed proteins from P.pastoris depends critically on growth conditions, and attainment of high cell densities by fermentation has been shown to improve protein yields by 10–100-fold. Unfortunately, the cost of the isotopically enriched fermentation media components, particularly 15NH4OH, is prohibitively high. We report fermentation methods that allow for both 15N- labeling from (15NH4)2SO4 and 13C-labeling from 13C-glucose or 13C-glycerol of proteins produced in Pichia pastoris. Expression of an 83 amino acid fragment of thrombomodulin with two N-linked glycosylation sites shows that fermentation is more cost effective than shake flask growth for isotopic enrichment.     

STING Millennium Suite: integrated software for extensive analyses of 3d structures of proteins and their complexes

Background  

The integration of many aspects of protein/DNA structure analysis is an important requirement for software products in general area of structural bioinformatics. In fact, there are too few software packages on the internet which can be described as successful in this respect. We might say that what is still missing is publicly available, web based software for interactive analysis of the sequence/structure/function of proteins and their complexes with DNA and ligands. Some of existing software packages do have certain level of integration and do offer analysis of several structure related parameters, however not to the extent generally demanded by a user.     

Two different proteins that compete for binding to thrombin have opposite kinetic and thermodynamic profiles

Thrombin binds thrombomodulin (TM) at anion binding exosite 1, an allosteric site far from the thrombin active site. A monoclonal antibody (mAb) has been isolated that competes with TM for binding to thrombin. Complete binding kinetic and thermodynamic profiles for these two protein-protein interactions have been generated. Binding kinetics were measured by Biacore. Although both interactions have similar K(D)s, TM binding is rapid and reversible while binding of the mAb is slow and nearly irreversible. The enthalpic contribution to the DeltaG(bind) was measured by isothermal titration calorimetry and van't Hoff analysis. The contribution to the DeltaG(bind) from electrostatic steering was assessed from the dependence of the k(a) on ionic strength. Release of solvent H(2)O molecules from the interface was assessed by monitoring the decrease in amide solvent accessibility at the interface upon protein-protein binding. The mAb binding is enthalpy driven and has a slow k(d). TM binding appears to be entropy driven and has a fast k(a). The favorable entropy of the thrombin-TM interaction seems to be derived from electrostatic steering and a contribution from solvent release. The two interactions have remarkably different thermodynamic driving forces for competing reactions. The possibility that optimization of binding kinetics for a particular function may be reflected in different thermodynamic driving forces is discussed.     

Solution structure of the smallest cofactor-active fragment of thrombomodulin

A glycosylated fragment of thrombomodulin containing two epidermal growth factor-like domains (TMEGF45) was analyzed by NMR. The 4th-domains structure of this two-domain fragment is similar to that of the individual domain previously determined. The 5th-domain, which has uncrossed disulfide bonds, is not as well determined in the two-domain fragment than the individual domain previously solved. The flexibility of the 5th-domain is consistent with low heteronuclear NOEs. In the individual 5th-domain, Met 388 was disordered, and key thrombin binding residues formed a hydrophobic core. By contrast, in TMEGF45, Met 388 is in the 5th-domain core, positioned by Phe 376 from the 4th-domain. As a result, key thrombin binding residues that were in the core of the individual domain are expelled. Upon thrombin binding, chemical shifts of two residues in the 4th-domain, the three interdomain linker residues, and nearly all of the 5th-domain are perturbed. Thus, TMEGF45 binds thrombin by an induced fit mechanism involving a flexible 5th-domain.     

Granule-bound starch synthase: structure, function, and phylogenetic utility

Interest in the use of low-copy nuclear genes for phylogenetic analyses of plants has grown rapidly, because highly repetitive genes such as those commonly used are limited in number. Furthermore, because low- copy genes are subject to different evolutionary processes than are plastid genes or highly repetitive nuclear markers, they provide a valuable source of independent phylogenetic evidence. The gene for granule-bound starch synthase (GBSSI or waxy) exists in a single copy in nearly all plants examined so far. Our study of GBSSI had three parts: (1) Amino acid sequences were compared across a broad taxonomic range, including grasses, four dicotyledons, and the microbial homologs of GBSSI. Inferred structural information was used to aid in the alignment of these very divergent sequences. The informed alignments highlight amino acids that are conserved across all sequences, and demonstrate that structural motifs can be highly conserved in spite of marked divergence in amino acid sequence. (2) Maximum-likelihood (ML) analyses were used to examine exon sequence evolution throughout grasses. Differences in probabilities among substitution types and marked among-site rate variation contributed to the observed pattern of variation. Of the parameters examined in our set of likelihood models, the inclusion of among-site rate variation following a gamma distribution caused the greatest improvement in likelihood score. (3) We performed cladistic parsimony analyses of GBSSI sequences throughout grasses, within tribes, and within genera to examine the phylogenetic utility of the gene. Introns provide useful information among very closely related species, but quickly become difficult to align among more divergent taxa. Exons are variable enough to provide extensive resolution within the family, but with low bootstrap support. The combined results of amino acid sequence comparisons, maximum-likelihood analyses, and phylogenetic studies underscore factors that might affect phylogenetic reconstruction. In this case, accommodation of the variable rate of evolution among sites might be the first step in maximizing the phylogenetic utility of GBSSI.      

The function of communities in protein interaction networks at multiple scales

Background  

If biology is modular then clusters, or communities, of proteins derived using only protein interaction network structure should define protein modules with similar biological roles. We investigate the link between biological modules and network communities in yeast and its relationship to the scale at which we probe the network.     

Interactions of the NPXY microdomains of the low density lipoprotein receptor‐related protein 1

The low density lipoprotein receptor‐related protein 1 (LRP1) mediates internalization of a large number of proteins and protein–lipid complexes and is widely implicated in Alzheimer's disease. The cytoplasmic domain of LRP1 (LRP1‐CT) can be phosphorylated by activated protein‐tyrosine kinases at two NPXY motifs in LRP1‐CT; Tyr 4507 is readily phosphorylated and must be phosphorylated before phosphorylation of Tyr 4473 occurs. Pull‐down experiments from brain lysate revealed numerous proteins binding to LRP1‐CT, but the results were highly variable. To separate which proteins bind to each NPXY motif and their phosphorylation dependence, each NPXY motif microdomain was prepared in both phosphorylated and non‐phosphorylated forms and used to probe rodent brain extracts for binding proteins. Proteins that bound specifically to the microdomains were identified by LC‐MS/MS, and confirmed by Western blot. Recombinant proteins were then tested for binding to each NPXY motif. The NPXY₄₅₀₇ (membrane distal) was found to interact with a large number of proteins, many of which only bound the tyrosine‐phosphorylated form. This microdomain also bound a significant number of other proteins in the unphosphorylated state. Many of the interactions were later confirmed to be direct with recombinant proteins. The NPXY₄₄₇₃ (membrane proximal) bound many fewer proteins and only to the phosphorylated form.