Machine learning methods for prediction of disulphide bonding states of cysteine residues in proteins

The goal of this thesis work is to develop a computational method based on machine learning techniques for predicting disulfide-bonding states of cysteine residues in proteins, which is a sub-problem of a bigger and yet unsolved problem of protein structure prediction. Improvement in the prediction of disulfide bonding states of cysteine residues will help in putting a constraint in the three dimensional (3D) space of the respective protein structure, and thus will eventually help in the prediction of 3D structure of proteins. Results of this work will have direct implications in site-directed mutational studies of proteins, proteins engineering and the problem of protein folding. We have used a combination of Artificial Neural Network (ANN) and Hidden Markov Model (HMM), the so-called Hidden Neural Network (HNN) as a machine learning technique to develop our prediction method. By using different global and local features of proteins (specifically profiles, parity of cysteine residues, average cysteine conservation, correlated mutation, sub-cellular localization, and signal peptide) as inputs and considering Eukaryotes and Prokaryotes separately we have reached to a remarkable accuracy of 94% on cysteine basis for both Eukaryotic and Prokaryotic datasets, and an accuracy of 90% and 93% on protein basis for Eukaryotic dataset and Prokaryotic dataset respectively. These accuracies are best so far ever reached by any existing prediction methods, and thus our prediction method has outperformed all the previously developed approaches and therefore is more reliable. Most interesting part of this thesis work is the differences in the prediction performances of Eukaryotes and Prokaryotes at the basic level of input coding when ‘profile’ information was given as input to our prediction method. And one of the reasons for this we discover is the difference in the amino acid composition of the local environment of bonded and free cysteine residues in Eukaryotes and Prokaryotes. Eukaryotic bonded cysteine examples have a ‘symmetric-cysteine-rich’ environment, where as Prokaryotic bonded examples lack it.

Myelin basic protein kinase activity in tomato leaves is induced systemically by wounding and increases in response to systemin and oligosaccharide elicitors

In response to wounding, a 48-kDa myelin basic protein (MBP) kinase is activated within 2 min, both locally and systemically, in leaves of young tomato plants. The activating signal is able to pass through a steam girdle on the stem, indicating that it moves through the xylem and does not require intact phloem tissue. A 48-kDa MBP kinase is also activated by the 18-amino acid polypeptide systemin, a potent wound signal for the synthesis of systemic wound response proteins (swrps). The kinase activation by systemin is strongly inhibited by a systemin analog having a Thr-17 → Ala-17 substitution, which is a powerful antagonist of systemin activation of swrp genes. A 48-kDa MBP kinase activity also increases in response to polygalacturonic acid and chitosan but not in response to jasmonic acid or phytodienoic acid. In def1, a mutant tomato line having a defective octadecanoid pathway, the 48-kDa MBP kinase is activated by wounding and systemin as in the wild-type plants. This indicates that MBP kinase functions between the perception of primary signals and the DEF1 gene product. In response to wounding, the MBP kinase is phosphorylated on phosphotyrosine residues, indicating a relationship to the mitogen-activated protein kinase family of protein kinases.

Neuronal type information encoded in the basic-helix–loop–helix domain of proneural genes

The proneural genes encode basic-helix–loop–helix (bHLH) proteins and promote the formation of distinct types of sensory organs. In Drosophila, two sets of proneural genes, atonal (ato) and members of the achaete–scute complex (ASC), are required for the formation of chordotonal (ch) organs and external sensory (es) organs, respectively. We assayed the production of sensory organs in transgenic flies expressing chimeric genes of ato and scute (sc), a member of ASC, and found that the information that specifies ch organs resides in the bHLH domain of ato; chimeras containing the b domain of ato and the HLH domain of sc also induced ch organ formation, but to a lesser extent than those containing the bHLH domain of ato. The b domains of ato and sc differ in seven residues. Mutations of these seven residues in the b domain of ato suggest that most or perhaps all of these residues are required for induction of ch organs. None of these seven residues is predicted to contact DNA directly by computer simulation using the structure of the myogenic factor MyoD as a model, implying that interaction of ato with other cofactors is likely to be involved in neuronal type specification.

Structural basis of DNA bending and oriented heterodimer binding by the basic leucine zipper domains of Fos and Jun

Interactions among transcription factors that bind to separate sequence elements require bending of the intervening DNA and juxtaposition of interacting molecular surfaces in an appropriate orientation. Here, we examine the effects of single amino acid substitutions adjacent to the basic regions of Fos and Jun as well as changes in sequences flanking the AP-1 site on DNA bending. Substitution of charged amino acid residues at positions adjacent to the basic DNA-binding domains of Fos and Jun altered DNA bending. The change in DNA bending was directly proportional to the change in net charge for all heterodimeric combinations between these proteins. Fos and Jun induced distinct DNA bends at different binding sites. Exchange of a single base pair outside of the region contacted in the x-ray crystal structure altered DNA bending. Substitution of base pairs flanking the AP-1 site had converse effects on the opposite directions of DNA bending induced by homodimers and heterodimers. These results suggest that Fos and Jun induce DNA bending in part through electrostatic interactions between amino acid residues adjacent to the basic region and base pairs flanking the AP-1 site. DNA bending by Fos and Jun at inverted binding sites indicated that heterodimers bind to the AP-1 site in a preferred orientation. Mutation of a conserved arginine within the basic regions of Fos and transversion of the central C:G base pair in the AP-1 site to G:C had complementary effects on the orientation of heterodimer binding and DNA bending. The conformational variability of the Fos–Jun–AP-1 complex may contribute to its functional versatility at different promoters.

Identification of surface residues mediating tissue factor binding and catalytic function of the serine protease factor VIIa

Factor VIIa (VIIa), the serine protease that initiates the coagulation pathways, is catalytically activated upon binding to its cell surface receptor and cofactor tissue factor (TF). This study provides a comprehensive analysis of the functional surface of VIIa by alanine scanning mutagenesis of 112 residues. Residue side chains were defined which contribute to TF binding and factor X hydrolysis. Energetically important binding contacts at the interface with TF were identified in the first epidermal growth factor domain of VIIa (Gln-64, Ile-69, Phe-71, Arg-79) and in the protease domain (Arg-277, Met-306, Asp-309). The observed energetic defects are in good agreement with the corresponding residues in TF, suggesting that the VIIa light chain plays a prominent role in high affinity binding of cofactor. Mutation of protease domain interface residues indicated that TF allosterically influences the active site of VIIa. Stabilization of a labile zymogen to enzyme transition could explain the activating effect of TF on VIIa catalytic function. Residues important for factor X hydrolysis were found in three regions of the protease domain: (i) specificity determinants in the catalytic cleft and adjacent loops, (ii) an exosite near the TF binding site, and (iii) a large electronegative exosite which is in a position analogous to the basic exosite I of thrombin. TF regions involved in factor X activation are positioned on the same face of the TF·VIIa complex as the two exosites identified on the protease domain surface, providing evidence for an extended interaction of TF·VIIa with macromolecular substrate.

Physical reasons for the unusual α-helix stabilization afforded by charged or neutral polar residues in alanine-rich peptides

We have carried out conformational energy calculations on alanine-based copolymers with the sequence Ac-AAAAAXAAAA-NH2 in water, where X stands for lysine or glutamine, to identify the underlying source of stability of alanine-based polypeptides containing charged or highly soluble polar residues in the absence of charge–charge interactions. The results indicate that ionizable or neutral polar residues introduced into the sequence to make them soluble sequester the water away from the CO and NH groups of the backbone, thereby enabling them to form internal hydrogen bonds. This solvation effect dictates the conformational preference and, hence, modifies the conformational propensity of alanine residues. Even though we carried out simulations for specific amino acid sequences, our results provide an understanding of some of the basic principles that govern the process of folding of these short sequences independently of the kind of residues introduced to make them soluble. In addition, we have investigated through simulations the effect of the bulk dielectric constant on the conformational preferences of these peptides. Extensive conformational Monte Carlo searches on terminally blocked 10-mer and 16-mer homopolymers of alanine in the absence of salt were carried out assuming values for the dielectric constant of the solvent ɛ of 80, 40, and 2. Our simulations show a clear tendency of these oligopeptides to augment the α-helix content as the bulk dielectric constant of the solvent is lowered. This behavior is due mainly to a loss of exposure of the CO and NH groups to the aqueous solvent. Experimental evidence indicates that the helical propensity of the amino acids in water shows a dramatic increase on addition of certain alcohols, such us trifluoroethanol. Our results provide a possible explanation of the mechanism by which alcohol/water mixtures affect the free energy of helical alanine oligopeptides relative to nonhelical ones.

Clusters of charged residues in protein three-dimensional structures.

Statistically significant charge clusters (basic, acidic, or of mixed charge) in tertiary protein structures are identified by new methods from a large representative collection of protein structures. About 10% of protein structures show at least one charge cluster, mostly of mixed type involving about equally anionic and cationic residues. Positive charge clusters are very rare. Negative (or histidine-acidic) charge clusters often coordinate calcium, or magnesium or zinc ions [e.g., thermolysin (PDB code: 3tln), mannose-binding protein (2msb), aminopeptidase (1amp)]. Mixed-charge clusters are prominent at interchain contacts where they stabilize quaternary protein formation [e.g., glutathione S-transferase (2gst), catalase (8act), and fructose-1,6-bisphosphate aldolase (1fba)]. They are also involved in protein-protein interaction and in substrate binding. For example, the mixed-charge cluster of aspartate carbamoyl-transferase (8atc) envelops the aspartate carbonyl substrate in a flexible manner (alternating tense and relaxed states) where charge associations can vary from weak to strong. Other proteins with charge clusters include the P450 cytochrome family (BM-3, Terp, Cam), several flavocytochromes, neuraminidase, hemagglutinin, the photosynthetic reaction center, and annexin. In each case in Table 2 we discuss the possible role of the charge clusters with respect to protein structure and function.

In vivo selection of basic region-leucine zipper proteins with altered DNA-binding specificities.

A transcription interference assay was used to generate mutant basic region-leucine zipper proteins with altered DNA-binding specificities. A library of mutants of a CCAAT/enhancer binding protein was constructed by randomizing five DNA-contacting amino acids in the basic region Asn-18, Ala-15, Val-14, Ser-11, and Arg-10. These mutants were then selected for their ability to bind mutant recognition sequences containing substitutions at the 2 and 3 positions of the wild-type sequence 5'-A5T4T3G2C1G1'C2'A3A4'T5'-3'. Mutants containing the sequence Leu-18Tyr-15Xaa-14Tyr-11Arg-10, in which four of the five contact residues are altered, were identified that recognize the palindromic sequence 5'-ATCYCGY'GAT-3' (Xaa = asparagine when Y = G; Xaa = methionine when Y = A). Moreover, in a selection against the sequence 5'-ATTACGTAAT-3', mutants were obtained containing substitutions not only in the basic region but also in the hinge region between the basic and leucine zipper regions. The mutant proteins showed high specificity in a functional transcription interference assay. A model for the interaction of these mutants with the target DNA sequences is discussed.

Structure of human T-cell receptors specific for an immunodominant myelin basic protein peptide: positioning of T-cell receptors on HLA-DR2/peptide complexes.

T-cell receptors (TCRs) recognize peptide bound within the relatively conserved structural framework of major histocompatibility complex (MHC) class I or class II molecules but can discriminate between closely related MHC molecules. The structural basis for the specificity of ternary complex formation by the TCR and MHC/peptide complexes was examined for myelin basic protein (MBP)-specific T-cell clones restricted by different DR2 subtypes. Conserved features of this system allowed a model for positioning of the TCR on DR2/peptide complexes to be developed: (i) The DR2 subtypes that presented the immunodominant MBP peptide differed only at a few polymorphic positions of the DR beta chain. (ii) TCR recognition of a polymorphic residue on the helical portion of the DR beta chain (position DR beta 67) was important in determining the MHC restriction. (iii) The TCR variable region (V) alpha 3.1 gene segment was used by all of the T-cell clones. TCR V beta usage was more diverse but correlated with the MHC restriction--i.e., with the polymorphic DR beta chains. (iv) Two clones with conserved TCR alpha chains but different TCR beta chains had a different MHC restriction but a similar peptide specificity. The difference in MHC restriction between these T-cell clones appeared due to recognition of a cluster of polymorphic DR beta-chain residues (DR beta 67-71). MBP-(85-99)-specific TCRs therefore appeared to be positioned on the DR2/peptide complex such that the TCR beta chain contacted the polymorphic DR beta-chain helix while the conserved TCR alpha chain contacted the nonpolymorphic DR alpha chain.

Amino Acid Residues Implicated in the Interaction of Melanocortin Ligands and their Receptors: A Study of MC2R Selectivity

Author: Kristopher D. Veo Title: Amino acid residues implicated in the interaction of Melanocortin ligands and their receptors: A study of MC2R selectivity Advisor: Dr. Robert M. Dores Degree Date: August 2009 ABSTRACT Melanocortin receptor ligand selectivity has been a question not easily answered. The inability to functionally express melanocortin 2 receptor (MC2R) has inhibited the study of why MC2R is only stimulated by ACTH, a melanocortin hormone. With the recent discovery of the MC2R accessory protein (MRAP), creating a heterologous system is now feasible. Using a general cell line like CHO-K1 cells, which do not express endogenous MCRs, we were able to create a heterologous expression system and test the selectivity of MC2R using analog variants of ACTH(1-24). Our results indicate an amino acid requirement in the C-terminal portion of ACTH(1-24) for activation, which supports the 2-step method of activation hypothesized for MC2R. This site, the tetra basic cleavage site, when altered does not stimulate cAMP production and does not compete with ACTH(1-24) for binding. We also demonstrate the potential for a non-mammalian MC2R system in cloning full length Silurana tropicalis MC2R and completed localization studies with this system with MRAP using CHO-K1 cells.

Role of proline residues in the expression and function of the human noradrenaline transporter

The aim was to investigate the roles of proline residues in extracellular loop 2 (P172, P183, P188 and P209) and transmembrane domains 2, 5, 11 and 12 (P108, P270, P526, P551, P552 and P570) in determining noradrenaline transporter (NET) expression and function. Mutants of human NET with these residues mutated to alanine were pharmacologically characterized. Mutation of P108, P270 and P526 disrupted cell surface expression, from [H-3]nisoxetine binding and confocal microscopy data. Mutations of P526, P551 and P570 reduced transporter turnover (V-max of [H-3]noradrenaline uptake/B-max of [H-3]nisoxetine binding) by 1.5-1.7-fold compared with wild-type NET, so these residues might be involved in conformational changes associated with substrate translocation. Conversely, mutations of P172, P183, P188 and P209 increased V-max/B-max by 2-3-fold compared with wild-type, indicating that the presence of these proline residues limits turnover of the NET. The mutations had few effects on apparent affinities of substrates or affinities of inhibitors, except decreases in inhibitor affinities after mutations of the P270 and P570 residues, and increases after mutation of the P526 residue. Hence, proline residues in extracellular loop 2 and in transmembrane domains have a range of roles in determining expression and function of the NET.

The Basic Reproduction Number Obtained From Jacobian And Next Generation Matrices - A Case Study Of Dengue Transmission Modelling.

The basic reproduction number is a key parameter in mathematical modelling of transmissible diseases. From the stability analysis of the disease free equilibrium, by applying Routh-Hurwitz criteria, a threshold is obtained, which is called the basic reproduction number. However, the application of spectral radius theory on the next generation matrix provides a different expression for the basic reproduction number, that is, the square root of the previously found formula. If the spectral radius of the next generation matrix is defined as the geometric mean of partial reproduction numbers, however the product of these partial numbers is the basic reproduction number, then both methods provide the same expression. In order to show this statement, dengue transmission modelling incorporating or not the transovarian transmission is considered as a case study. Also tuberculosis transmission and sexually transmitted infection modellings are taken as further examples.

[access To Basic Health Services And Associated Factors: A Population-based Study].

This study sought to identify factors involved in access to the services of a basic health unit. It is a cross-sectional, population-based study involving 101 randomly-selected families residing in the area covered by the health unit. An adult resident of each household was interviewed. The response variable was whether or not the resident frequented the health unit if he/she or anyone in the family required assistance to resolve a health issue. The independent variables investigated were service provision aspects, demographic and socio-economic characteristics, individual habits, morbidities and use of the health unit. In addition to descriptive and univariate analysis, logistic regression was applied in the multivariate analysis. The results show that access to the basic health unit is associated with the treatment received previously (OR = 3,224) with accessibility (OR = 0,146) and micro-area of residence (OR = 10,918). These findings suggest that access is related to the impressions created by the care received at the health unit and is based on experiences with the service, but can also be strongly modulated by individual aspects and factors related to the territory.

Phtx-ii A Basic Myotoxic Phospholipase A2 From Porthidium Hyoprora Snake Venom, Pharmacological Characterization And Amino Acid Sequence By Mass Spectrometry.

A monomeric basic PLA2 (PhTX-II) of 14149.08 Da molecular weight was purified to homogeneity from Porthidium hyoprora venom. Amino acid sequence by in tandem mass spectrometry revealed that PhTX-II belongs to Asp49 PLA2 enzyme class and displays conserved domains as the catalytic network, Ca2+-binding loop and the hydrophobic channel of access to the catalytic site, reflected in the high catalytic activity displayed by the enzyme. Moreover, PhTX-II PLA2 showed an allosteric behavior and its enzymatic activity was dependent on Ca2+. Examination of PhTX-II PLA2 by CD spectroscopy indicated a high content of alpha-helical structures, similar to the known structure of secreted phospholipase IIA group suggesting a similar folding. PhTX-II PLA2 causes neuromuscular blockade in avian neuromuscular preparations with a significant direct action on skeletal muscle function, as well as, induced local edema and myotoxicity, in mice. The treatment of PhTX-II by BPB resulted in complete loss of their catalytic activity that was accompanied by loss of their edematogenic effect. On the other hand, enzymatic activity of PhTX-II contributes to this neuromuscular blockade and local myotoxicity is dependent not only on enzymatic activity. These results show that PhTX-II is a myotoxic Asp49 PLA2 that contributes with toxic actions caused by P. hyoprora venom.

Motivating Medical Students To Learn Basic Science Concepts Using Chronic Myeloid Leukemia As An Integration Theme.

To report on the use of chronic myeloid leukemia as a theme of basic clinical integration for first year medical students to motivate and enable in-depth understanding of the basic sciences of the future physician. During the past thirteen years we have reviewed and updated the curriculum of the medical school of the Universidade Estadual de Campinas. The main objective of the new curriculum is to teach the students how to learn to learn. Since then, a case of chronic myeloid leukemia has been introduced to first year medical students and discussed in horizontal integration with all themes taught during a molecular and cell biology course. Cell structure and components, protein, chromosomes, gene organization, proliferation, cell cycle, apoptosis, signaling and so on are all themes approached during this course. At the end of every topic approached, the students prepare in advance the corresponding topic of clinical cases chosen randomly during the class, which are then presented by them. During the final class, a paper regarding mutations in the abl gene that cause resistance to tyrosine kinase inhibitors is discussed. After each class, three tests are solved in an interactive evaluation. The course has been successful since its beginning, 13 years ago. Great motivation of those who participated in the course was observed. There were less than 20% absences in the classes. At least three (and as many as nine) students every year were interested in starting research training in the field of hematology. At the end of each class, an interactive evaluation was performed and more than 70% of the answers were correct in each evaluation. Moreover, for the final evaluation, the students summarized, in a written report, the molecular and therapeutic basis of chronic myeloid leukemia, with scores ranging from 0 to 10. Considering all 13 years, a median of 78% of the class scored above 5 (min 74%-max 85%), and a median of 67% scored above 7. Chronic myeloid leukemia is an excellent example of a disease that can be used for clinical basic integration as this disorder involves well known protein, cytogenetic and cell function abnormalities, has well-defined diagnostic strategies and a target oriented therapy.