Using Helix-coil Models to Study Protein Unfolded States

An abstract of a thesis devoted to using helix-coil models to study unfolded states.\\

Research on polypeptide unfolded states has received much more attention in the last decade or so than it has in the past. Unfolded states are thought to be implicated in various

misfolding diseases and likely play crucial roles in protein folding equilibria and folding rates. Structural characterization of unfolded states has proven to be

much more difficult than the now well established practice of determining the structures of folded proteins. This is largely because many core assumptions underlying

folded structure determination methods are invalid for unfolded states. This has led to a dearth of knowledge concerning the nature of unfolded state conformational

distributions. While many aspects of unfolded state structure are not well known, there does exist a significant body of work stretching back half a century that

has been focused on structural characterization of marginally stable polypeptide systems. This body of work represents an extensive collection of experimental

data and biophysical models associated with describing helix-coil equilibria in polypeptide systems. Much of the work on unfolded states in the last decade has not been devoted

specifically to the improvement of our understanding of helix-coil equilibria, which arguably is the most well characterized of the various conformational equilibria

that likely contribute to unfolded state conformational distributions. This thesis seeks to provide a deeper investigation of helix-coil equilibria using modern

statistical data analysis and biophysical modeling techniques. The studies contained within seek to provide deeper insights and new perspectives on what we presumably

know very well about protein unfolded states. \\

Chapter 1 gives an overview of recent and historical work on studying protein unfolded states. The study of helix-coil equilibria is placed in the context

of the general field of unfolded state research and the basics of helix-coil models are introduced.\\

Chapter 2 introduces the newest incarnation of a sophisticated helix-coil model. State of the art modern statistical techniques are employed to estimate the energies

of various physical interactions that serve to influence helix-coil equilibria. A new Bayesian model selection approach is utilized to test many long-standing

hypotheses concerning the physical nature of the helix-coil transition. Some assumptions made in previous models are shown to be invalid and the new model

exhibits greatly improved predictive performance relative to its predecessor. \\

Chapter 3 introduces a new statistical model that can be used to interpret amide exchange measurements. As amide exchange can serve as a probe for residue-specific

properties of helix-coil ensembles, the new model provides a novel and robust method to use these types of measurements to characterize helix-coil ensembles experimentally

and test the position-specific predictions of helix-coil models. The statistical model is shown to perform exceedingly better than the most commonly used

method for interpreting amide exchange data. The estimates of the model obtained from amide exchange measurements on an example helical peptide

also show a remarkable consistency with the predictions of the helix-coil model. \\

Chapter 4 involves a study of helix-coil ensembles through the enumeration of helix-coil configurations. Aside from providing new insights into helix-coil ensembles,

this chapter also introduces a new method by which helix-coil models can be extended to calculate new types of observables. Future work on this approach could potentially

allow helix-coil models to move into use domains that were previously inaccessible and reserved for other types of unfolded state models that were introduced in chapter 1.

On Provable Algorithms for Determination of Continuous Protein Interdomain Motions from Residual Dipolar Couplings

Dynamics of biomolecules over various spatial and time scales are essential for biological functions such as molecular recognition, catalysis and signaling. However, reconstruction of biomolecular dynamics from experimental observables requires the determination of a conformational probability distribution. Unfortunately, these distributions cannot be fully constrained by the limited information from experiments, making the problem an ill-posed one in the terminology of Hadamard. The ill-posed nature of the problem comes from the fact that it has no unique solution. Multiple or even an infinite number of solutions may exist. To avoid the ill-posed nature, the problem needs to be regularized by making assumptions, which inevitably introduce biases into the result.

Here, I present two continuous probability density function approaches to solve an important inverse problem called the RDC trigonometric moment problem. By focusing on interdomain orientations we reduced the problem to determination of a distribution on the 3D rotational space from residual dipolar couplings (RDCs). We derived an analytical equation that relates alignment tensors of adjacent domains, which serves as the foundation of the two methods. In the first approach, the ill-posed nature of the problem was avoided by introducing a continuous distribution model, which enjoys a smoothness assumption. To find the optimal solution for the distribution, we also designed an efficient branch-and-bound algorithm that exploits the mathematical structure of the analytical solutions. The algorithm is guaranteed to find the distribution that best satisfies the analytical relationship. We observed good performance of the method when tested under various levels of experimental noise and when applied to two protein systems. The second approach avoids the use of any model by employing maximum entropy principles. This 'model-free' approach delivers the least biased result which presents our state of knowledge. In this approach, the solution is an exponential function of Lagrange multipliers. To determine the multipliers, a convex objective function is constructed. Consequently, the maximum entropy solution can be found easily by gradient descent methods. Both algorithms can be applied to biomolecular RDC data in general, including data from RNA and DNA molecules.

Estimate of the (R)USLE rainfall erosivity factor from monthly precipitation data in mainland Spain

The need for continuous recording rain gauges makes it difficult to determine the rainfall erosivity factor (R-factor) of the (R)USLE model in areas without good temporal data coverage. In mainland Spain, the Nature Conservation Institute (ICONA) determined the R-factor at few selected pluviographs, so simple estimates of the R-factor are definitely of great interest. The objectives of this study were: (1) to identify a readily available estimate of the R-factor for mainland Spain; (2) to discuss the applicability of a single (global) estimate based on analysis of regional results; (3) to evaluate the effect of record length on estimate precision and accuracy; and (4) to validate an available regression model developed by ICONA. Four estimators based on monthly precipitation were computed at 74 rainfall stations throughout mainland Spain. The regression analysis conducted at a global level clearly showed that modified Fournier index (MFI) ranked first among all assessed indexes. Applicability of this preliminary global model across mainland Spain was evaluated by analyzing regression results obtained at a regional level. It was found that three contiguous regions of eastern Spain (Catalonia, Valencian Community and Murcia) could have a different rainfall erosivity pattern, so a new regression analysis was conducted by dividing mainland Spain into two areas: Eastern Spain and plateau-lowland area. A comparative analysis concluded that the bi-areal regression model based on MFI for a 10-year record length provided a simple, precise and accurate estimate of the R-factor in mainland Spain. Finally, validation of the regression model proposed by ICONA showed that R-ICONA index overpredicted the R-factor by approximately 19%.

Mapping Protein-DNA Interactions Using ChIP-exo and Illumina-Based Sequencing

Chromatin immunoprecipitation (ChIP) provides a means of enriching DNA associated with transcription factors, histone modifications, and indeed any other proteins for which suitably characterized antibodies are available. Over the years, sequence detection has progressed from quantitative real-time PCR and Southern blotting to microarrays (ChIP-chip) and now high-throughput sequencing (ChIP-seq). This progression has vastly increased the sequence coverage and data volumes generated. This in turn has enabled informaticians to predict the identity of multi-protein complexes on DNA based on the overrepresentation of sequence motifs in DNA enriched by ChIP with a single antibody against a single protein. In the course of the development of high-throughput sequencing, little has changed in the ChIP methodology until recently. In the last three years, a number of modifications have been made to the ChIP protocol with the goal of enhancing the sensitivity of the method and further reducing the levels of nonspecific background sequences in ChIPped samples. In this chapter, we provide a brief commentary on these methodological changes and describe a detailed ChIP-exo method able to generate narrower peaks and greater peak coverage from ChIPped material.

An Ex Vivo Model for Studying Hepatic Schistosomiasis and the Effect of Released Protein from Dying Eggs

BACKGROUND: We report the use of an ex vivo precision cut liver slice (PCLS) mouse model for studying hepatic schistosomiasis. In this system, liver tissue is unfixed, unfrozen, and alive for maintenance in culture and subsequent molecular analysis.

METHODS AND FINDINGS: Using thick naive mouse liver tissue and sterile culture conditions, the addition of soluble egg antigen (SEA) derived from Schistosoma japonicum eggs, followed 4, 24 and 48 hrs time points. Tissue was collected for transcriptional analysis and supernatants collected to quantitate liver enzymes, cytokines and chemokines. No significant hepatotoxicity was demonstrated by supernatant liver enzymes due to the presence of SEA. A proinflammatory response was observed both at the transcriptional level and at the protein level by cytokine and chemokine bead assay. Key genes observed elevated transcription in response to the addition of SEA included: IL1-α and IL1-β, IL6, all associated with inflammation. The recruitment of antigen presenting cells was reflected in increases in transcription of CD40, CCL4 and CSF1. Indications of tissue remodeling were seen in elevated gene expression of various Matrix MetalloProteinases (MMP3, 9, 10, 13) and delayed increases in TIMP1. Collagen deposition was significantly reduced in the presence of SEA as shown in COL1A1 expression by qPCR after 24 hrs culture. Cytokine and chemokine analysis of the culture supernatants confirmed the elevation of proteins including IL6, CCL3, CCL4 and CXCL5.

CONCLUSIONS: This ex vivo model system for the synchronised delivery of parasite antigen to liver tissue provides an insight into the early phase of hepatic schistosomiasis, corresponding with the release of soluble proteins from dying schistosome eggs.

Exchange protein directly activated by cyclic AMP (EPAC) activation reverses neutrophil dysfunction induced by β2-agonists, corticosteroids, and critical illness

Background Neutrophils play a role in the pathogenesis of asthma, chronic obstructive pulmonary disease, and pulmonary infection. Impaired neutrophil phagocytosis predicts hospital-acquired infection. Despite this, remarkably few neutrophil-specific treatments exist. 

Objectives We sought to identify novel pathways for the restoration of effective neutrophil phagocytosis and to activate such pathways effectively in neutrophils from patients with impaired neutrophil phagocytosis. 

Methods Blood neutrophils were isolated from healthy volunteers and patients with impaired neutrophil function. In healthy neutrophils phagocytic impairment was induced experimentally by using β₂-agonists. Inhibitors and activators of cyclic AMP (cAMP)-dependent pathways were used to assess the influence on neutrophil phagocytosis in vitro. 

Results β₂-Agonists and corticosteroids inhibited neutrophil phagocytosis. Impairment of neutrophil phagocytosis by β₂-agonists was associated with significantly reduced RhoA activity. Inhibition of protein kinase A (PKA) restored phagocytosis and RhoA activity, suggesting that cAMP signals through PKA to drive phagocytic impairment. However, cAMP can signal through effectors other than PKA, such as exchange protein directly activated by cyclic AMP (EPAC). An EPAC-activating analog of cAMP (8CPT-2Me-cAMP) reversed neutrophil dysfunction induced by β₂-agonists or corticosteroids but did not increase RhoA activity. 8CPT-2Me-cAMP reversed phagocytic impairment induced by Rho kinase inhibition but was ineffective in the presence of Rap-1 GTPase inhibitors. 8CPT-2Me-cAMP restored function to neutrophils from patients with known acquired impairment of neutrophil phagocytosis. 

Conclusions EPAC activation consistently reverses clinical and experimental impairment of neutrophil phagocytosis. EPAC signals through Rap-1 and bypasses RhoA. EPAC activation represents a novel potential means by which to reverse impaired neutrophil phagocytosis.

Serum Fatty Acid Binding Protein 4 (FABP4) Predicts Pre-eclampsia in Women with Type 1 Diabetes

Objective: To examine the association between fatty acid binding protein 4 (FABP4) and pre-eclampsia risk in women with type 1 diabetes.
Reesearch Design and Methods: Serum FABP4 was measured in 710 women from the Diabetes and Pre-eclampsia Intervention Trial (DAPIT) in early pregnancy and in the second trimester (median 14 and 26 weeks gestation, respectively).
Results: FABP4 was significantly elevated in early pregnancy (geometric mean 15.8 ng/mL [interquartile range 11.6–21.4] vs. 12.7 ng/mL [interquartile range 9.6–17]; P < 0.001) and the second trimester (18.8 ng/mL [interquartile range 13.6–25.8] vs. 14.6 ng/mL [interquartile range 10.8–19.7]; P < 0.001) in women in whom pre-eclampsia later developed. Elevated second-trimester FABP4 level was independently associated with pre-eclampsia (odds ratio 2.87 [95% CI 1.24, 6.68], P = 0.03). The addition of FABP4 to established risk factors significantly improved net reclassification improvement at both time points and integrated discrimination improvement in the second trimester.
Conclusions: Increased second-trimester FABP4 independently predicted pre-eclampsia and significantly improved reclassification and discrimination. FABP4 shows potential as a novel biomarker for pre-eclampsia prediction in women with type 1 diabetes.

Comparison of module detection algorithms in protein networks and investigation of the biological meaning of predicted modules

Background
It is generally acknowledged that a functional understanding of a biological system can only be obtained by an understanding of the collective of molecular interactions in form of biological networks. Protein networks are one particular network type of special importance, because proteins form the functional base units of every biological cell. On a mesoscopic level of protein networks, modules are of significant importance because these building blocks may be the next elementary functional level above individual proteins allowing to gain insight into fundamental organizational principles of biological cells.
Results
In this paper, we provide a comparative analysis of five popular and four novel module detection algorithms. We study these module prediction methods for simulated benchmark networks as well as 10 biological protein interaction networks (PINs). A particular focus of our analysis is placed on the biological meaning of the predicted modules by utilizing the Gene Ontology (GO) database as gold standard for the definition of biological processes. Furthermore, we investigate the robustness of the results by perturbing the PINs simulating in this way our incomplete knowledge of protein networks.
Conclusions
Overall, our study reveals that there is a large heterogeneity among the different module prediction algorithms if one zooms-in the biological level of biological processes in the form of GO terms and all methods are severely affected by a slight perturbation of the networks. However, we also find pathways that are enriched in multiple modules, which could provide important information about the hierarchical organization of the system

The Precipitationshed : Concepts, Methods, and Applications

Human societies are reliant on the functioning of the hydrologic cycle. The atmospheric branch of this cycle, often referred to as moisture recycling in the context of land-to-land exchange, refers to water evaporating, traveling through the atmosphere, and falling out as precipitation. Similar to the surface water cycle that uses the watershed as the unit of analysis, it is also possible to consider a ‘watershed of the sky’ for the atmospheric water cycle. Thus, I explore the precipitationshed - defined as the upwind surface of the Earth that provides evaporation that later falls as precipitation in a specific place. The primary contributions of this dissertation are to (a) introduce the precipitationshed concept, (b) provide a quantitative basis for the study of the precipitationshed, and (c) demonstrate its use in the fields of hydrometeorology, land-use change, social-ecological systems, ecosystem services, and environmental governance. In Paper I, the concept of the precipitationshed is introduced and explored for the first time. The quantification of precipitationshed variability is described in Paper II, and the key finding is that the precipitationsheds for multiple regions are persistent in time and space. Moisture recycling is further described as an ecosystem service in Paper III, to integrate the concept into the existing language of environmental sustainability and management. That is, I identify regions where vegetation more strongly regulates the provision of atmospheric water, as well as the regions that more strongly benefit from this regulation. In Paper IV, the precipitationshed is further explored through the lens of urban reliance on moisture recycling. Using a novel method, I quantify the vulnerability of urban areas to social-ecological changes within their precipitationsheds. In Paper V, I argue that successful moisture recycling governance will require flexible, transboundary institutions that are capable of operating within complex social-ecological systems. I conclude that, in the future, the precipitationshed can be a key tool in addressing the complexity of social-ecological systems. 

Insulin-Like Growth Factor-I (IGF-1), IGF-Binding Protein-3 (IGFBP-3) and mammographic features

Introduction. The IGF system has recently been shown to play an important role in the regulation of breast tumor cell proliferation. However, also breast density is currently considered as the strongest breast cancer risk factor. It is not yet clear whether these factors are interrelated and if and how they are influenced by menopausal status. The purpose of this study was to examine the possible effects of IGF-1 and IGFBP-3 and IGF-1/IGFBP-3 molar ratio on mammographic density stratified by menopausal status. Patients and methods. A group of 341 Italian women were interviewed to collect the following data: family history of breast cancer, reproductive and menstrual factors, breast biopsies, previous administration of hormonal contraceptive therapy, hormone replacement therapy (HRT) in menopause and lifestyle information. A blood sample was drawn for determination of IGF-1, IGFBP-3 levels. IGF-1/ IGFBP-3 molar ratio was then calculated. On the basis of recent mammograms the women were divided into two groups: dense breast (DB) and non-dense breast (NDB). Student’s t-test was employed to assess the association between breast density and plasma level of IGF-1, IGFBP-3 and molar ratio. To assess if this relationship was similar in subgroups of pre- and postmenopausal women, the study population was stratified by menopausal status and Student’s t-test was performed. Finally, multivariate analysis was employed to evaluate if there were confounding factors that might influence the relationship between growth factors and breast density. Results. The analysis of the relationship between mammographic density and plasma level of IGF-1, IGFBP-3 and IGF-1/ IGFBP-3 molar ratio showed that IGF-1 levels and molar ratio varied in the two groups resulting in higher mean values in the DB group (IGF-1: 109.6 versus 96.6 ng/ml; p= 0.001 and molar ratio 29.4 versus 25.5 ng/ml; p= 0.001) whereas IGFBP-3 showed similar values in both groups (DB and NDB). Analysis of plasma level of IGF-1, IGFBP-3 and IGF-1/IGFBP-3 molar ratio compared to breast density after stratification of the study population by menopausal status (premenopausal and postmenopausal) showed that there was no association between the plasma of growth factors and breast density, neither in premenopausal nor in postmenopausal patients. Multivariate analysis showed that only nulliparity, premenopausal status and body mass index (BMI) are determinants of breast density. Conclusions. Our study provides a strong evidence of a crude association between breast density and plasma levels of IGF-1 and molar ratio. On the basis of our results, it is reasonable to assume that the role of IGF-1 and molar ratio in the pathogenesis of breast cancer might be mediated through mammographic density. IGF-1 and molar ratio might thus increase the risk of cancer by increasing mammographic density.

Therapeutic Potential of RNAi Through Endocytotic Methods

The ability to manipulate gene expression promises to be an important tool for the management of infectious diseases and genetic disorders. However, a major limitation to effective delivery of therapeutic RNA to living cells is the cellular toxicity of conventional techniques. Team PANACEA’s research objective was to create new reagents based on a novel small-molecule delivery system that uses a modular recombinant protein vehicle consisting of a specific ligand coupled to a Hepatitis B Virus-derived RNA binding domain (HBV-RBD). Two such recombinant delivery proteins were developed: one composed of Interleukin-8, the other consisting of the Machupo Virus GP1 protein. The ability of these proteins to deliver RNA to cells were then tested. The non-toxic nature of this technology has the potential to overcome limitations of current methods and could provide a platform for the expansion of personalized medicine.

Analysis of Molecular Dynamics Simulations of Protein Folding

Microsecond long Molecular Dynamics (MD) trajectories of biomolecular processes are now possible due to advances in computer technology. Soon, trajectories long enough to probe dynamics over many milliseconds will become available. Since these timescales match the physiological timescales over which many small proteins fold, all atom MD simulations of protein folding are now becoming popular. To distill features of such large folding trajectories, we must develop methods that can both compress trajectory data to enable visualization, and that can yield themselves to further analysis, such as the finding of collective coordinates and reduction of the dynamics. Conventionally, clustering has been the most popular MD trajectory analysis technique, followed by principal component analysis (PCA). Simple clustering used in MD trajectory analysis suffers from various serious drawbacks, namely, (i) it is not data driven, (ii) it is unstable to noise and change in cutoff parameters, and (iii) since it does not take into account interrelationships amongst data points, the separation of data into clusters can often be artificial. Usually, partitions generated by clustering techniques are validated visually, but such validation is not possible for MD trajectories of protein folding, as the underlying structural transitions are not well understood. Rigorous cluster validation techniques may be adapted, but it is more crucial to reduce the dimensions in which MD trajectories reside, while still preserving their salient features. PCA has often been used for dimension reduction and while it is computationally inexpensive, being a linear method, it does not achieve good data compression. In this thesis, I propose a different method, a nonmetric multidimensional scaling (nMDS) technique, which achieves superior data compression by virtue of being nonlinear, and also provides a clear insight into the structural processes underlying MD trajectories. I illustrate the capabilities of nMDS by analyzing three complete villin headpiece folding and six norleucine mutant (NLE) folding trajectories simulated by Freddolino and Schulten [1]. Using these trajectories, I make comparisons between nMDS, PCA and clustering to demonstrate the superiority of nMDS. The three villin headpiece trajectories showed great structural heterogeneity. Apart from a few trivial features like early formation of secondary structure, no commonalities between trajectories were found. There were no units of residues or atoms found moving in concert across the trajectories. A flipping transition, corresponding to the flipping of helix 1 relative to the plane formed by helices 2 and 3 was observed towards the end of the folding process in all trajectories, when nearly all native contacts had been formed. However, the transition occurred through a different series of steps in all trajectories, indicating that it may not be a common transition in villin folding. The trajectories showed competition between local structure formation/hydrophobic collapse and global structure formation in all trajectories. Our analysis on the NLE trajectories confirms the notion that a tight hydrophobic core inhibits correct 3-D rearrangement. Only one of the six NLE trajectories folded, and it showed no flipping transition. All the other trajectories get trapped in hydrophobically collapsed states. The NLE residues were found to be buried deeply into the core, compared to the corresponding lysines in the villin headpiece, thereby making the core tighter and harder to undo for 3-D rearrangement. Our results suggest that the NLE may not be a fast folder as experiments suggest. The tightness of the hydrophobic core may be a very important factor in the folding of larger proteins. It is likely that chaperones like GroEL act to undo the tight hydrophobic core of proteins, after most secondary structure elements have been formed, so that global rearrangement is easier. I conclude by presenting facts about chaperone-protein complexes and propose further directions for the study of protein folding.

Prediction of Phosphorylation Sites from Protein Sequences

Phosphorylation is amongst the most crucial and well-studied post-translational modifications. It is involved in multiple cellular processes which makes phosphorylation prediction vital for understanding protein functions. However, wet-lab techniques are labour and time intensive. Thus, computational tools are required for efficiency. This project aims to provide a novel way to predict phosphorylation sites from protein sequences by adding flexibility and Sezerman Grouping amino acid similarity measure to previous methods, as discovering new protein sequences happens at a greater rate than determining protein structures. The predictor – NOPAY - relies on Support Vector Machines (SVMs) for classification. The features include amino acid encoding, amino acid grouping, predicted secondary structure, predicted protein disorder, predicted protein flexibility, solvent accessibility, hydrophobicity and volume. As a result, we have managed to improve phosphorylation prediction accuracy for Homo sapiens by 3% and 6.1% for Mus musculus. Sensitivity at 99% specificity was also increased by 6% for Homo sapiens and for Mus musculus by 5% on independent test sets. In this study, we have managed to increase phosphorylation prediction accuracy for Homo sapiens and Mus musculus. When there is enough data, future versions of the software may also be able to predict other organisms.

Produção de isolado proteico proveniente de subprodutos da indústria de aves em diferentes escalas

A indústria de aves brasileira destaca-se economicamente, onde um total de 12,3 milhões de toneladas foi produzido no país em 2013. Esta produção em larga escala gera considerável volume de subprodutos, chegando até 35% da ave viva. Tais resíduos são convertidos, por processos tradicionais, em produtos de baixo valor comercial, como por exemplo, farinhas. O processo de variação de pH constitui um importante processo alternativo de obtenção de proteínas com melhores características funcionais e nutricionais. Estudar as variáveis do processo, efetuando aumento dimensional, é fundamental para aplicação das tecnologias desenvolvidas no laboratório e posterior definição final de processos industriais. A produção de isolados proteicos seria uma tecnologia atraente no aproveitamento de subprodutos da indústria de frango, convertendo-os em uma ótima fonte proteica, agregando valor ao produto obtido. Este trabalho teve por objetivo produzir isolados proteicos em diferentes escalas, utilizando subprodutos não comestíveis da indústria de frango. Foi estudada a solubilização das proteínas da matéria-prima (MP) para definir pHs de solubilização e de precipitação isoelétrica. A curva apontou um pH alcalino de 11,0 para etapa de solubilização e de 5,25 para etapa de precipitação proteica. As proteínas obtidas foram caracterizadas quanto sua composição proximal, índice de acidez (IA), índice de peróxidos (IP) e substâncias reativas ao ácido tiobarbitúrico (TBARS) além de propriedades funcionais de solubilidade, capacidade de retenção de água (CRA) e capacidade de retenção de óleo (CRO); e nutricionais de digestibilidade proteica. Comparativamente foram analisadas farinhas de vísceras comerciais nos mesmos parâmetros. Um aumento de escala do processo foi realizado e avaliado pelas mesmas respostas do produto da escala laboratorial. Foi obtido um teor proteico de 82 e 85% em escala laboratorial e aumento de escala, respectivamente, e também uma redução lipídica de 75%, e de cinzas de 85%, em relação à MP. A composição proximal das farinhas analisadas ficou entre 67-72% para proteína bruta, 17-22% para lipídios e 9-15% para cinzas. O IA, apresentou valores de 2,2 e 3,1 meq/g de isolado e de 1,6 a 2,0 meq/g de farinha. Já para IP, obteve-se valores de 0,003 a 0,005 meq/g de isolado e de 0,002 a 0,049 meq/g de farinha. Os índices de TBARS apontaram valores de 0,081 e 0,214 mg MA/g de isolado e 0,041 a 0,128 mg MA/g de farinha. A solubilidade das proteínas do isolado apontou 84 e 81% em pH 3 e 11 respectivamente e de 5% em pH 5, já para farinhas variaram de 22 a 31% em pH de 3 a 11. A CRA obtida no isolado foi 3,1 a 16,5 g água/g de proteína e de 3,8 a 10,9 g água/g de proteína nas farinhas. A CRO ficou em 4,2 mL de óleo/g de proteína do isolados e 2,6 mL de óleo/g de proteína da farinhas. Os isolados proteicos apresentaram 92 e 95% de digestibilidade das proteínas, em comparação aos 84% das farinhas comerciais. Os índices acumulados e apresentados neste trabalho concluíram que foi possível aumentar a escala do processo de variação de pH, sem perder qualidade nos índices físico-químicos e de digestibilidade proteica.

Effects of early nutrition on the wnt and bone morphogenic protein pathways in the developing intestine

Infant formula is consumed by the majority of infants in the United States for at least part of the first year of life. Infant formula lacks many of the bioactive compounds that are naturally occurring in breast milk. Because of this, there has been an increased interest by the companies that manufacture infant formula to include additives that would potentially allow formula to more closely mimic breast milk activity. One such ingredient currently being added to infant formula is prebiotics. Prebiotics are non-digestible food ingredients that beneficially affect the host by selectively stimulating the growth of specific healthful bacteria in the colon. It is speculated that prebiotics replicate the activity of breast milk oligosaccharides, which through the production of butyrate by intestinal microbiota, may interact with the Wnt/BMP pathways. The Wnt/BMP pathways regulate intestinal stem cells, which determine the growth, development and maintenance of the intestine. Therefore, the objective of this study was to explore the effects that the addition of prebiotics to formula have on the regulation of the Wnt/BMP pathways when fed to neonatal piglets, a model commonly used in the study of infant nutrition. Piglets (n=5) were randomized into sow-reared (SR), fed control formula (F), or fed formula with added prebiotics (F+P). Fructooligosaccharides (FOS) (2 g/L) and polydextrose (PDX) (2 g/L) were chosen as the prebiotics for this study, because this combination had been less studied than other combinations. Ileum and ascending colon were collected at 7 and 14 days-of-age. Dry matter content, pH, and short chain fatty acid (SCFA) content was measured. The mRNA expression of β-catenin, sFRP3, sFRP4, frizzled 6, DKK1 (Wnt pathway), gremlin (BMP pathway), TNF-a, HNF-4α and osteopontin (OPN) was measured by RT-qPCR. Piglets fed the F+P diet had greater acetate concentration and lower pH in the ileum at day 14 and in the colon at day 7 and day 14 than F piglets. Butyrate concentrations were highest in SR with F+P not differing from F in ileum at day 14 and colon at day 7 and day 14. Effects of age were seen in all genes, with the exception of OPN, sFRP-3 and sFRP-4. On day 7, no effect of diet was observed in the ileum, however, mRNA expression of DKK1 and frizzled 6 were greater in F+P than SR (p≤0.05). On day 14, gremlin expression was lower and OPN was greater in the ileum of SR piglets compared to F and F+P. Also on day 14, HNF-4α mRNA expression was greater in both ileum and colon of F+P piglets and sFRP3 mRNA expression was greater in the colon than F or SR . In summary, differences were observed between gene expression of F+P and SR piglet intestines, but the supplementation of 2 g/L scFOS and 2 g/L PDX to formula did not shift expression of genes in the Wnt/BMP pathways to be more similar to SR than F. As the Wnt/BMP pathway is known to exist in a gradient along the crypt-villus axis, with Wnt expression dominating in the crypt region and BMP expression dominating in the villi, it was possible that pooling whole tissue reduced our ability to detect treatment effects that would be concentrated in either region. A method was therefore developed to remove intestinal epithelial cells along the villus-to-crypt axis. Twenty-five-day-old F and SR piglets were euthanized and ileal tissue was collected and placed in a dissociation buffer in a shaking water bath. Exfoliated cells were removed at increasing time points from 5 to 100 minutes in order to remove cells along the villus-to-crypt axis. After the final incubation, remaining mucosal tissue was removed using a sterile glass microscope slide and pooled with the final exfoliated cell isolation. After each cell collection, a section of tissue was fixed in formalin for histomorphological examination. Expression of genes in the Wnt/BMP pathways, along with crypt marker genes (CDK5 and v-myb), were measured in both whole ileal tissue, pooled epithelial cells, and separate epithelial cell isolations from the same piglet. The expression of β-catenin, HNF-4α, TNF-α, TGF-β and the crypt marker v-myb matched the expected villus-to-crypt pattern in cells collected after 10 (incubation 1), 30 (incubation 2) and 60 (incubation 3) minutes. However, expression of expression in cells collected after 100 minutes (incubation 4) was variable, which may be due to the fact that crypt cells were not efficiently removed and the presence of unwanted non-epithelial tissue. Gremlin, OPN, DKK1, sFRP3 and sFRP4 expression was not statistically different along the villus-to-crypt axis. Frizzled 6 and CDK5 did not express as we had predicted, with expression highest towards the villi. In summary, the epithelial cell collection method used was not entirely successful. While much of the gene data suggests that cells were removed along the villus-to-crypt axis through the first three incubations, the last incubation, which involved scraping the tissue, removed non-epithelial components of the mucosa, while leaving the crypts intact. In conclusion, the addition of 2 g/L PDX and 2 g/L scFOS did not cause gene expression of the Wnt/BMP pathways to mirror either F or SR expression. New isolation methods to extract cells along the crypt-villus axis should be considered, including the use of a laser capture microdissection. While this combination of prebiotics did not yield the intended effects, future research should be done on other combinations, such as the inclusion of galactooligosaccharides (GOS), which is commonly added to food products including infant formula.