In silico protein design by combinatorial assembly of protein building blocks

Utilizing concepts of protein building blocks, we propose a de novo computational algorithm that is similar to combinatorial shuffling experiments. Our goal is to engineer new naturally occurring folds with low homology to existing proteins. A selected protein is first partitioned into its building blocks based on their compactness, degree of isolation from the rest of the structure, and hydrophobicity. Next, the protein building blocks are substituted by fragments taken from other proteins with overall low sequence identity, but with a similar hydrophobic/hydrophilic pattern and a high structural similarity. These criteria ensure that the designed protein has a similar fold, low sequence identity, and a good hydrophobic core compared with its native counterpart. Here, we have selected two proteins for engineering, protein G B1 domain and ubiquitin. The two engineered proteins share approximately 20% and approximately 25% amino acid sequence identities with their native counterparts, respectively. The stabilities of the engineered proteins are tested by explicit water molecular dynamics simulations. The algorithm implements a strategy of designing a protein using relatively stable fragments, with a high population time. Here, we have selected the fragments by searching for local minima along the polypeptide chain using the protein building block model. Such an approach provides a new method for engineering new proteins with similar folds and low homology.     

Bayesian analysis of experimental epidemics of foot-and-mouth disease

We investigate the transmission dynamics of a certain type of foot-and-mouth disease (FMD) virus under experimental conditions. Previous analyses of experimental data from FMD outbreaks in non-homogeneously mixing populations of sheep have suggested a decline in viraemic level through serial passage of the virus, but these do not take into account possible variation in the length of the chain of viral transmission for each animal, which is implicit in the non-observed transmission process. We consider a susceptible-exposed-infectious-removed non-Markovian compartmental model for partially observed epidemic processes, and we employ powerful methodology (Markov chain Monte Carlo) for statistical inference, to address epidemiological issues under a Bayesian framework that accounts for all available information and associated uncertainty in a coherent approach. The analysis allows us to investigate the posterior distribution of the hidden transmission history of the epidemic, and thus to determine the effect of the length of the infection chain on the recorded viraemic levels, based on the posterior distribution of a p-value. Parameter estimates of the epidemiological characteristics of the disease are also obtained. The results reveal a possible decline in viraemia in one of the two experimental outbreaks. Our model also suggests that individual infectivity is related to the level of viraemia.     

The Stanford Microarray Database: data access and quality assessment tools

The Stanford Microarray Database (SMD; http://genome-www.stanford.edu/microarray/) serves as a microarray research database for Stanford investigators and their collaborators. In addition, SMD functions as a resource for the entire scientific community, by making freely available all of its source code and providing full public access to data published by SMD users, along with many tools to explore and analyze those data. SMD currently provides public access to data from 3500 microarrays, including data from 85 publications, and this total is increasing rapidly. In this article, we describe some of SMD's newer tools for accessing public data, assessing data quality and for data analysis.     

Targeting alphavbeta3 and alpha5beta1 for gene delivery to proliferating VSMCs: synergistic effect of TGF-beta1

TGF-beta1 levels increase after vascular injury and promote vascular smooth muscle cell (VSMC) proliferation. We define a nonviral gene delivery system that targets alphavbeta3 and alpha5beta1 integrins that are expressed on proliferating VSMCs and strongly induced by TGF-beta1. A 15-amino acid RGDNP-containing peptide from American Pit Viper venom was linked to a Lys(16) peptide as vector (molossin vector) and complexed with Lipofectamine or fusogenic peptide for delivery of luciferase or beta-galactosidase reporter genes to primary cultures of human, rabbit, and rat VSMCs. Preincubation of VSMCs with TGF-beta1 for 24 h, but not with PDGF-BB, interferon-gamma, TNF-alpha, nor PMA, increased alphavbeta3 and alpha5beta1 expressions on VSMCs and enhanced gene delivery of molossin vector. Thus beta-galactosidase activity increased from 35 +/- 5% (controls) to 75 +/- 5% after TGF-beta1 treatment, and luciferase activity increased fourfold over control values. Potential use of this system in vessel bypass surgery was examined in an ex vivo rat aortic organ culture model after endothelial damage. Molossin vector system delivered beta-galactosidase to VSMCs in the vessel wall that remained for up to 12 days posttransfection. The molossin vector system, when combined with TGF-beta1, enhances gene delivery to proliferating VSMCs and might have clinical applications for certain vasculoproliferative diseases.     

DeltaBAFF,an alternate splice isoform that regulates receptor binding and biopresentation of the B cell survival cytokine,BAFF

The tumor necrosis family member BAFF is limiting for the survival of follicular B lymphocytes, but excessive BAFF signaling can lead to autoimmunity, suggesting that its activity must be tightly regulated. We have identified a conserved alternate splice isoform of BAFF, called deltaBAFF, which lacks 57 nt encoding the A-A1 loop and is co-expressed with BAFF in many mouse and human myeloid cells. Mouse deltaBAFF appears on the plasma membrane, but unlike BAFF it is inefficiently released by proteolysis. DeltaBAFF can associate with BAFF in heteromultimers and diminish BAFF bioactivity and release. Thus, alternative splicing of the BAFF gene suppresses BAFF B cell stimulatory function in several ways, and deltaBAFF may promote other functions as well.     

Ghrelin can bind to a species of high density lipoprotein associated with paraoxonase

Ghrelin is a 28-residue peptide hormone that is principally released from the stomach during fasting and prior to eating. Two forms are present in human plasma: the unmodified peptide and a less abundant acylated version, in which octanoic acid is attached to the third residue, a serine, via an ester linkage. The acylated form of ghrelin acts as a ligand for the growth hormone secretagogue receptor and can stimulate the release of growth hormone from the pituitary gland. It also initiates behavioral and metabolic adaptations to fasting. Here we show that an immobilized form of ghrelin specifically binds a species of high density lipoprotein associated with the plasma esterase, paraoxonase, and clusterin. Both free ghrelin and paraoxon, a substrate for paraoxonase, can inhibit this interaction. An endogenous species of ghrelin is found to co-purify with high density lipoprotein during density gradient centrifugation and subsequent gel filtration. This interaction links the orexigenic peptide hormone ghrelin to lipid transport and metabolism. Furthermore, the interaction of the esterified hormone ghrelin with a species of HDL containing an esterase suggests a possible mechanism for the conversion of ghrelin to des-acyl ghrelin.     

Effect of four-alpha-helix bundle cavity size on volatile anesthetic binding energetics

Currently, it is thought that inhalational anesthetics cause anesthesia by binding to ligand-gated ion channels. This is being investigated using four-alpha-helix bundles, small water-soluble analogues of the transmembrane domains of the "natural" receptor proteins. The study presented here specifically investigates how multiple alanine-to-valine substitutions (which each decrease the volume of the internal binding cavity by 38 A(3)) affect structure, stability, and anesthetic binding affinity of the four-alpha-helix bundles. Structure remains essentially unchanged when up to four alanine residues are changed to valine. However, stability increases as the number of these substitutions is increased. Anesthetic binding affinities are also affected. Halothane binds to the four-alpha-helix bundle variants with 0, 1, and 2 substitutions with equivalent affinities but binds to the variants with 3 and 4 more tightly. The same order of binding affinities was observed for chloroform, although for a particular variant, chloroform was bound less tightly. The observed differences in binding affinities may be explained in terms of a modulation of van der Waals and hydrophobic interactions between ligand and receptor. These, in turn, could result from increased four-alpha-helix bundle binding cavity hydrophobicity, a decrease in cavity size, or improved ligand/receptor shape complementarity.     

Production of recombinant proteins using multiple-copy gene integration in high-cell-density fermentations of Ralstonia eutropha

We have previously reported the development of a novel protein expression system based on Ralstonia eutropha. In this study we report on the influence of gene copynumber on recombinant protein expression in R. eutropha. We compare recombinant gene stability and expression levels of chromosomal integration with a plasmid-based expression system. Single, double, and triple copies of a gene encoding organophosphohydrolase (OPH), an enzyme prone to inclusion-body formation in E. coli, were integrated into the R. eutropha chromosome. A linear increase between the concentration of soluble, active OPH and gene copynumber was found. Using a triple-copy integrant, we were able to produce approximately 4.3 g/L of OPH in a high-cell-density fermentation. This represents the highest titer reported to date for this enzyme, and is approximately 30 times greater than expression levels reported in E. coli.     

Hydrolysis of biological peptides by human angiotensin-converting enzyme-related carboxypeptidase

Human angiotensin-converting enzyme-related carboxypeptidase (ACE2) is a zinc metalloprotease whose closest homolog is angiotensin I-converting enzyme. To begin to elucidate the physiological role of ACE2, ACE2 was purified, and its catalytic activity was characterized. ACE2 proteolytic activity has a pH optimum of 6.5 and is enhanced by monovalent anions, which is consistent with the activity of ACE. ACE2 activity is increased approximately 10-fold by Cl(-) and F(-) but is unaffected by Br(-). ACE2 was screened for hydrolytic activity against a panel of 126 biological peptides, using liquid chromatography-mass spectrometry detection. Eleven of the peptides were hydrolyzed by ACE2, and in each case, the proteolytic activity resulted in removal of the C-terminal residue only. ACE2 hydrolyzes three of the peptides with high catalytic efficiency: angiotensin II () (k(cat)/K(m) = 1.9 x 10(6) m(-1) s(-1)), apelin-13 (k(cat)/K(m) = 2.1 x 10(6) m(-1) s(-1)), and dynorphin A 1-13 (k(cat)/K(m) = 3.1 x 10(6) m(-1) s(-1)). The ACE2 catalytic efficiency is 400-fold higher with angiotensin II () as a substrate than with angiotensin I (). ACE2 also efficiently hydrolyzes des-Arg(9)-bradykinin (k(cat)/K(m) = 1.3 x 10(5) m(-1) s(-1)), but it does not hydrolyze bradykinin. An alignment of the ACE2 peptide substrates reveals a consensus sequence of: Pro-X((1-3 residues))-Pro-Hydrophobic, where hydrolysis occurs between proline and the hydrophobic amino acid.