Thermodynamic analysis of ligand-induced changes in protein thermal unfolding applied to high-throughput determination of ligand affinities with extrinsic fluorescent dyes.

The quantification of protein-ligand interactions is essential for systems biology, drug discovery, and bioengineering. Ligand-induced changes in protein thermal stability provide a general, quantifiable signature of binding and may be monitored with dyes such as Sypro Orange (SO), which increase their fluorescence emission intensities upon interaction with the unfolded protein. This method is an experimentally straightforward, economical, and high-throughput approach for observing thermal melts using commonly available real-time polymerase chain reaction instrumentation. However, quantitative analysis requires careful consideration of the dye-mediated reporting mechanism and the underlying thermodynamic model. We determine affinity constants by analysis of ligand-mediated shifts in melting-temperature midpoint values. Ligand affinity is determined in a ligand titration series from shifts in free energies of stability at a common reference temperature. Thermodynamic parameters are obtained by fitting the inverse first derivative of the experimental signal reporting on thermal denaturation with equations that incorporate linear or nonlinear baseline models. We apply these methods to fit protein melts monitored with SO that exhibit prominent nonlinear post-transition baselines. SO can perturb the equilibria on which it is reporting. We analyze cases in which the ligand binds to both the native and denatured state or to the native state only and cases in which protein:ligand stoichiometry needs to treated explicitly.

Computational Inference of Genome-Wide Protein-DNA Interactions Using High-Throughput Genomic Data

Transcriptional regulation has been studied intensively in recent decades. One important aspect of this regulation is the interaction between regulatory proteins, such as transcription factors (TF) and nucleosomes, and the genome. Different high-throughput techniques have been invented to map these interactions genome-wide, including ChIP-based methods (ChIP-chip, ChIP-seq, etc.), nuclease digestion methods (DNase-seq, MNase-seq, etc.), and others. However, a single experimental technique often only provides partial and noisy information about the whole picture of protein-DNA interactions. Therefore, the overarching goal of this dissertation is to provide computational developments for jointly modeling different experimental datasets to achieve a holistic inference on the protein-DNA interaction landscape.

We first present a computational framework that can incorporate the protein binding information in MNase-seq data into a thermodynamic model of protein-DNA interaction. We use a correlation-based objective function to model the MNase-seq data and a Markov chain Monte Carlo method to maximize the function. Our results show that the inferred protein-DNA interaction landscape is concordant with the MNase-seq data and provides a mechanistic explanation for the experimentally collected MNase-seq fragments. Our framework is flexible and can easily incorporate other data sources. To demonstrate this flexibility, we use prior distributions to integrate experimentally measured protein concentrations.

We also study the ability of DNase-seq data to position nucleosomes. Traditionally, DNase-seq has only been widely used to identify DNase hypersensitive sites, which tend to be open chromatin regulatory regions devoid of nucleosomes. We reveal for the first time that DNase-seq datasets also contain substantial information about nucleosome translational positioning, and that existing DNase-seq data can be used to infer nucleosome positions with high accuracy. We develop a Bayes-factor-based nucleosome scoring method to position nucleosomes using DNase-seq data. Our approach utilizes several effective strategies to extract nucleosome positioning signals from the noisy DNase-seq data, including jointly modeling data points across the nucleosome body and explicitly modeling the quadratic and oscillatory DNase I digestion pattern on nucleosomes. We show that our DNase-seq-based nucleosome map is highly consistent with previous high-resolution maps. We also show that the oscillatory DNase I digestion pattern is useful in revealing the nucleosome rotational context around TF binding sites.

Finally, we present a state-space model (SSM) for jointly modeling different kinds of genomic data to provide an accurate view of the protein-DNA interaction landscape. We also provide an efficient expectation-maximization algorithm to learn model parameters from data. We first show in simulation studies that the SSM can effectively recover underlying true protein binding configurations. We then apply the SSM to model real genomic data (both DNase-seq and MNase-seq data). Through incrementally increasing the types of genomic data in the SSM, we show that different data types can contribute complementary information for the inference of protein binding landscape and that the most accurate inference comes from modeling all available datasets.

This dissertation provides a foundation for future research by taking a step toward the genome-wide inference of protein-DNA interaction landscape through data integration.

De novo design and molecular assembly of a transmembrane diporphyrin-binding protein complex.

The de novo design of membrane proteins remains difficult despite recent advances in understanding the factors that drive membrane protein folding and association. We have designed a membrane protein PRIME (PoRphyrins In MEmbrane) that positions two non-natural iron diphenylporphyrins (Fe(III)DPP's) sufficiently close to provide a multicentered pathway for transmembrane electron transfer. Computational methods previously used for the design of multiporphyrin water-soluble helical proteins were extended to this membrane target. Four helices were arranged in a D(2)-symmetrical bundle to bind two Fe(II/III) diphenylporphyrins in a bis-His geometry further stabilized by second-shell hydrogen bonds. UV-vis absorbance, CD spectroscopy, analytical ultracentrifugation, redox potentiometry, and EPR demonstrate that PRIME binds the cofactor with high affinity and specificity in the expected geometry.

Magnetic resonance water proton relaxation in protein solutions and tissue: T(1rho) dispersion characterization.

BACKGROUND: Image contrast in clinical MRI is often determined by differences in tissue water proton relaxation behavior. However, many aspects of water proton relaxation in complex biological media, such as protein solutions and tissue are not well understood, perhaps due to the limited empirical data. PRINCIPAL FINDINGS: Water proton T(1), T(2), and T(1rho) of protein solutions and tissue were measured systematically under multiple conditions. Crosslinking or aggregation of protein decreased T(2) and T(1rho), but did not change high-field T(1). T(1rho) dispersion profiles were similar for crosslinked protein solutions, myocardial tissue, and cartilage, and exhibited power law behavior with T(1rho)(0) values that closely approximated T(2). The T(1rho) dispersion of mobile protein solutions was flat above 5 kHz, but showed a steep curve below 5 kHz that was sensitive to changes in pH. The T(1rho) dispersion of crosslinked BSA and cartilage in DMSO solvent closely resembled that of water solvent above 5 kHz but showed decreased dispersion below 5 kHz. CONCLUSIONS: Proton exchange is a minor pathway for tissue T(1) and T(1rho) relaxation above 5 kHz. Potential models for relaxation are discussed, however the same molecular mechanism appears to be responsible across 5 decades of frequencies from T(1rho) to T(1).

The Haemophilus influenzae HMW1C protein is a glycosyltransferase that transfers hexose residues to asparagine sites in the HMW1 adhesin.

The Haemophilus influenzae HMW1 adhesin is a high-molecular weight protein that is secreted by the bacterial two-partner secretion pathway and mediates adherence to respiratory epithelium, an essential early step in the pathogenesis of H. influenzae disease. In recent work, we discovered that HMW1 is a glycoprotein and undergoes N-linked glycosylation at multiple asparagine residues with simple hexose units rather than N-acetylated hexose units, revealing an unusual N-glycosidic linkage and suggesting a new glycosyltransferase activity. Glycosylation protects HMW1 against premature degradation during the process of secretion and facilitates HMW1 tethering to the bacterial surface, a prerequisite for HMW1-mediated adherence. In the current study, we establish that the enzyme responsible for glycosylation of HMW1 is a protein called HMW1C, which is encoded by the hmw1 gene cluster and shares homology with a group of bacterial proteins that are generally associated with two-partner secretion systems. In addition, we demonstrate that HMW1C is capable of transferring glucose and galactose to HMW1 and is also able to generate hexose-hexose bonds. Our results define a new family of bacterial glycosyltransferases.

Enhanced myocardial relaxation in vivo in transgenic mice overexpressing the beta2-adrenergic receptor is associated with reduced phospholamban protein.

To assess the effect of targeted myocardial beta-adrenergic receptor (AR) stimulation on relaxation and phospholamban regulation, we studied the physiological and biochemical alterations associated with overexpression of the human beta2-AR gene in transgenic mice. These mice have an approximately 200-fold increase in beta-AR density and a 2-fold increase in basal adenylyl cyclase activity relative to negative littermate controls. Mice were catheterized with a high fidelity micromanometer and hemodynamic recordings were obtained in vivo. Overexpression of the beta2-AR altered parameters of relaxation. At baseline, LV dP/dt(min) and the time constant of LV pressure isovolumic decay (Tau) in the transgenic mice were significantly shorter compared with controls, indicating markedly enhanced myocardial relaxation. Isoproterenol stimulation resulted in shortening of relaxation velocity in control mice but not in the transgenic mice, indicating maximal relaxation in these animals. Immunoblotting analysis revealed a selective decrease in the amount of phospholamban protein, without a significant change in the content for either sarcoplasmic reticulum Ca2+ ATPase or calsequestrin, in the transgenic hearts compared with controls. This study indicates that myocardial relaxation is both markedly enhanced and maximal in these mice and that conditions associated with chronic beta-AR stimulation can result in a selective reduction of phospholamban protein.

G-protein-coupled receptor genes as protooncogenes: constitutively activating mutation of the alpha 1B-adrenergic receptor enhances mitogenesis and tumorigenicity.

The alpha 1B-adrenergic receptor (alpha 1B-ADR) is a member of the G-protein-coupled family of transmembrane receptors. When transfected into Rat-1 and NIH 3T3 fibroblasts, this receptor induces focus formation in an agonist-dependent manner. Focus-derived, transformed fibroblasts exhibit high levels of functional alpha 1B-ADR expression, demonstrate a catecholamine-induced enhancement in the rate of cellular proliferation, and are tumorigenic when injected into nude mice. Induction of neoplastic transformation by the alpha 1B-ADR, therefore, identifies this normal cellular gene as a protooncogene. Mutational alteration of this receptor can lead to activation of this protooncogene, resulting in an enhanced ability of agonist to induce focus formation with a decreased latency and quantitative increase in transformed foci. In contrast to cells expressing the wild-type alpha 1B-ADR, focus formation in "oncomutant"-expressing cell lines appears constitutively activated with the generation of foci in unstimulated cells. Further, these cell lines exhibit near-maximal rates of proliferation even in the absence of catecholamine supplementation. They also demonstrate an enhanced ability for tumor generation in nude mice with a decreased period of latency compared with cells expressing the wild-type receptor. Thus, the alpha 1B-ADR gene can, when overexpressed and activated, function as an oncogene inducing neoplastic transformation. Mutational alteration of this receptor gene can result in the activation of this protooncogene, enhancing its oncogenic potential. These findings suggest that analogous spontaneously occurring mutations in this class of receptor proteins could play a key role in the induction or progression of neoplastic transformation and atherosclerosis.

Beta-adrenergic receptor kinase: identification of a novel protein kinase that phosphorylates the agonist-occupied form of the receptor.

Agonist-promoted desensitization of adenylate cyclase is intimately associated with phosphorylation of the beta-adrenergic receptor in mammalian, avian, and amphibian cells. However, the nature of the protein kinase(s) involved in receptor phosphorylation remains largely unknown. We report here the identification and partial purification of a protein kinase capable of phosphorylating the agonist-occupied form of the purified beta-adrenergic receptor. The enzyme is prepared from a supernatant fraction from high-speed centrifugation of lysed kin- cells, a mutant of S49 lymphoma cells that lacks a functional cAMP-dependent protein kinase. The beta-agonist isoproterenol induces a 5- to 10-fold increase in receptor phosphorylation by this kinase, which is blocked by the antagonist alprenolol. Fractionation of the kin- supernatant on molecular-sieve HPLC and DEAE-Sephacel results in a 50- to 100-fold purified beta-adrenergic receptor kinase preparation that is largely devoid of other protein kinase activities. The kinase activity is insensitive to cAMP, cGMP, cAMP-dependent kinase inhibitor, Ca2+-calmodulin, Ca2+-phospholipid, and phorbol esters and does not phosphorylate general kinase substrates such as casein and histones. Phosphate appears to be incorporated solely into serine residues. The existence of this novel cAMP-independent kinase, which preferentially phosphorylates the agonist-occupied form of the beta-adrenergic receptor, suggests a mechanism that may explain the homologous or agonist-specific form of adenylate cyclase desensitization. It also suggests a general mechanism for regulation of receptor function in which only the agonist-occupied or "active" form of the receptor is a substrate for enzymes inducing covalent modification.

Impaired formation of beta-adrenergic receptor-nucleotide regulatory protein complexes in pseudohypoparathyroidism.

Decreased activity of the guanine nucleotide regulatory protein (N) of the adenylate cyclase system is present in cell membranes of some patients with pseudohypoparathyrodism (PHP-Ia) whereas others have normal activity of N (PHP-Ib). Low N activity in PHP-Ia results in a decrease in hormone (H)-stimulatable adenylate cyclase in various tissues, which might be due to decreased ability to form an agonist-specific high affinity complex composed of H, receptor (R), and N. To test this hypothesis, we compared beta-adrenergic agonist-specific binding properties in erythrocyte membranes from five patients with PHP-Ia (N = 45% of control), five patients with PHP-Ib (N = 97%), and five control subjects. Competition curves that were generated by increasing concentrations of the beta-agonist isoproterenol competing with [125I]pindolol were shallow (slope factors less than 1) and were computer fit to a two-state model with corresponding high and low affinity for the agonist. The agonist competition curves from the PHP-Ia patients were shifted significantly (P less than 0.02) to the right as a result of a significant (P less than 0.01) decrease in the percent of beta-adrenergic receptors in the high affinity state from 64 +/- 22% in PHP-Ib and 56 +/- 5% in controls to 10 +/- 8% in PHP-Ia. The agonist competition curves were computer fit to a "ternary complex" model for the two-step reaction: H + R + N in equilibrium HR + N in equilibrium HRN. The modeling was consistent with a 60% decrease in the functional concentration of N, and was in good agreement with the biochemically determined decrease in erythrocyte N protein activity. These in vitro findings in erythrocytes taken together with the recent observations that in vivo isoproterenol-stimulated adenylate cyclase activity is decreased in patients with PHP (Carlson, H. E., and A. S. Brickman, 1983, J. Clin. Endocrinol. Metab. 56:1323-1326) are consistent with the notion that N is a bifunctional protein interacting with both R and the adenylate cyclase. It may be that in patients with PHP-Ia a single molecular and genetic defect accounts for both decreased HRN formation and decreased adenylate cyclase activity, whereas in PHP-Ib the biochemical lesion(s) appear not to affect HRN complex formation.

Suppression of conformational heterogeneity at a protein-protein interface.

Staphylococcal protein A (SpA) is an important virulence factor from Staphylococcus aureus responsible for the bacterium's evasion of the host immune system. SpA includes five small three-helix-bundle domains that can each bind with high affinity to many host proteins such as antibodies. The interaction between a SpA domain and the Fc fragment of IgG was partially elucidated previously in the crystal structure 1FC2. Although informative, the previous structure was not properly folded and left many substantial questions unanswered, such as a detailed description of the tertiary structure of SpA domains in complex with Fc and the structural changes that take place upon binding. Here we report the 2.3-Å structure of a fully folded SpA domain in complex with Fc. Our structure indicates that there are extensive structural rearrangements necessary for binding Fc, including a general reduction in SpA conformational heterogeneity, freezing out of polyrotameric interfacial residues, and displacement of a SpA side chain by an Fc side chain in a molecular-recognition pocket. Such a loss of conformational heterogeneity upon formation of the protein-protein interface may occur when SpA binds its multiple binding partners. Suppression of conformational heterogeneity may be an important structural paradigm in functionally plastic proteins.

Component-wise gradient boosting and false discovery control in survival analysis with high-dimensional covariates.

MOTIVATION: Technological advances that allow routine identification of high-dimensional risk factors have led to high demand for statistical techniques that enable full utilization of these rich sources of information for genetics studies. Variable selection for censored outcome data as well as control of false discoveries (i.e. inclusion of irrelevant variables) in the presence of high-dimensional predictors present serious challenges. This article develops a computationally feasible method based on boosting and stability selection. Specifically, we modified the component-wise gradient boosting to improve the computational feasibility and introduced random permutation in stability selection for controlling false discoveries. RESULTS: We have proposed a high-dimensional variable selection method by incorporating stability selection to control false discovery. Comparisons between the proposed method and the commonly used univariate and Lasso approaches for variable selection reveal that the proposed method yields fewer false discoveries. The proposed method is applied to study the associations of 2339 common single-nucleotide polymorphisms (SNPs) with overall survival among cutaneous melanoma (CM) patients. The results have confirmed that BRCA2 pathway SNPs are likely to be associated with overall survival, as reported by previous literature. Moreover, we have identified several new Fanconi anemia (FA) pathway SNPs that are likely to modulate survival of CM patients. AVAILABILITY AND IMPLEMENTATION: The related source code and documents are freely available at https://sites.google.com/site/bestumich/issues. CONTACT: yili@umich.edu.

Potential To Streamline Heterologous DNA Prime and NYVAC/Protein Boost HIV Vaccine Regimens in Rhesus Macaques by Employing Improved Antigens.

UNLABELLED: In a follow-up to the modest efficacy observed in the RV144 trial, researchers in the HIV vaccine field seek to substantiate and extend the results by evaluating other poxvirus vectors and combinations with DNA and protein vaccines. Earlier clinical trials (EuroVacc trials 01 to 03) evaluated the immunogenicity of HIV-1 clade C GagPolNef and gp120 antigens delivered via the poxviral vector NYVAC. These showed that a vaccination regimen including DNA-C priming prior to a NYVAC-C boost considerably enhanced vaccine-elicited immune responses compared to those with NYVAC-C alone. Moreover, responses were improved by using three as opposed to two DNA-C primes. In the present study, we assessed in nonhuman primates whether such vaccination regimens can be streamlined further by using fewer and accelerated immunizations and employing a novel generation of improved DNA-C and NYVAC-C vaccine candidates designed for higher expression levels and more balanced immune responses. Three different DNA-C prime/NYVAC-C+ protein boost vaccination regimens were tested in rhesus macaques. All regimens elicited vigorous and well-balanced CD8(+)and CD4(+)T cell responses that were broad and polyfunctional. Very high IgG binding titers, substantial antibody-dependent cellular cytotoxicity (ADCC), and modest antibody-dependent cell-mediated virus inhibition (ADCVI), but very low neutralization activity, were measured after the final immunizations. Overall, immune responses elicited in all three groups were very similar and of greater magnitude, breadth, and quality than those of earlier EuroVacc vaccines. In conclusion, these findings indicate that vaccination schemes can be simplified by using improved antigens and regimens. This may offer a more practical and affordable means to elicit potentially protective immune responses upon vaccination, especially in resource-constrained settings. IMPORTANCE: Within the EuroVacc clinical trials, we previously assessed the immunogenicity of HIV clade C antigens delivered in a DNA prime/NYVAC boost regimen. The trials showed that the DNA prime crucially improved the responses, and three DNA primes with a NYVAC boost appeared to be optimal. Nevertheless, T cell responses were primarily directed toward Env, and humoral responses were modest. The aim of this study was to assess improved antigens for the capacity to elicit more potent and balanced responses in rhesus macaques, even with various simpler immunization regimens. Our results showed that the novel antigens in fact elicited larger numbers of T cells with a polyfunctional profile and a good Env-GagPolNef balance, as well as high-titer and Fc-functional antibody responses. Finally, comparison of the different schedules indicates that a simpler regimen of only two DNA primes and one NYVAC boost in combination with protein may be very efficient, thus showing that the novel antigens allow for easier immunization protocols.