Using Amino Acid Correlation and Community Detection Algorithms to Identify Functional Determinants in Protein Families

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Amino acid sequence and crystal structure of BaP1, a metalloproteinase from Bothrops asper snake venom that exerts multiple tissue-damaging activities

BaP1 is a 22.7-kD P-I-type zinc-dependent metalloproteinase isolated from the venom of the snake Bothrops asper, a medically relevant species in Central America. This enzyme exerts multiple tissue-damaging activities, including hemorrhage, myonecrosis, dermonecrosis, blistering, and edema. BaP1 is a single chain of 202 amino acids that shows highest sequence identity with metalloproteinases isolated front the venoms of snakes of the subfamily Crotalinae. It has six Cys residues involved in three disulfide bridges (Cys 117-Cys 197, Cys 159-Cys 181, Cys 157-Cys 164). It has the consensus sequence H(142)E(143)XXH(146)XXGXXH(152), as well as the sequence C164I165M166, which characterize the metzincin superfamily of metalloproteinases. The active-site cleft separates a major subdomain (residues 1-152), comprising four a-helices and a five-stranded beta-sheet, from the minor subdomain, which is formed by a single a-helix and several loops. The catalytic zinc ion is coordinated by the N-epsilon2 nitrogen atoms of His 142, His 146, and His 152, in addition to a solvent water molecule, which in turn is bound to Glu 143. Several conserved residues contribute to the formation of the hydrophobic pocket, and Met 166 serves as a hydrophobic base for the active-site groups. Sequence and structural comparisons of hemorrhagic and nonhemorrhagic P-I metalloproteinases from snake venoms revealed differences in several regions. In particular, the loop comprising residues 153 to 176 has marked structural differences between metalloproteinases with very different hemorrhagic activities. Because this region lies in close proximity to the active-site microenvironment, it may influence the interaction of these enzymes with physiologically relevant substrates in the extracellular matrix.

Structural studies of prephenate dehydratase from Mycobacterium tuberculosis H37Rv by SAXS, ultracentrifugation, and computational analysis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Structural studies of shikimate 5-dehydrogenase from Mycobacterium tuberculosis

Tuberculosis (TB) remains the leading cause of mortality due to a single bacterial pathogen, Mycobacterium tuberculosis. The reemergence of TB as a potential public health threat, the high susceptibility of human immunodeficiency virus-infected persons to the disease, the proliferation of multi-drug-resistant strains (MDR-TB) and, more recently, of extensively drug resistant isolates (XDR-TB) have created a need for the development of new antimycobacterial agents. Amongst the several proteins and/or enzymes to be studied as potential targets to develop novel drugs against M. tuberculosis, the enzymes of the shikimate pathway are attractive targets because they are essential in algae, higher plants, bacteria, and fungi, but absent from mammals. The mycobacterial shikimate pathway leads to the biosynthesis of chorismate, which is a precursor of aromatic amino acids, naphthoquinones, menaquinones, and mycobactins. Here we report the structural studies by homology modeling and circular dichroism spectroscopy of the shikimate dehydrogenase from M. tuberculosis (MtSDH), which catalyses the fourth step of the shikimate pathway. Our structural models show that the MtSDH has similar structure to other shikimate dehydrogenase structures previously reported either in presence or absence of NADP, despite the low amino acid sequence identity. The circular dichroism spectra corroborate the secondary structure content observed in the MtSDH models developed. The enzyme was stable up to 50 degrees C presenting a cooperative unfolding profile with the midpoint of the unfolding temperature value of similar to 63-64 degrees C, as observed in the unfolding experiment followed by circular dichroism. Our MtSDH structural models and circular dichroism data showed small conformational changes induced by NADP binding. We hope that the data presented here will assist the rational design of antitubercular agents.

Crystallization and preliminary X-ray diffraction analysis of an l-amino-acid oxidase from Bothrops jararacussu venom

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

FUNCTIONAL ALPHA-TROPOMYOSIN PRODUCED IN ESCHERICHIA-COLI - A DIPEPTIDE EXTENSION CAN SUBSTITUTE THE AMINO-TERMINAL ACETYL GROUP

Unlike the muscle protein, alpha-tropomyosin expressed in Escherichia coli does not bind actin, does not exhibit head-to-tail polymerization, and does not inhibit actomyosin ATPase activity in the absence of troponin. The only chemical difference between recombinant and muscle tropomyosins is that the first methionine is not acetylated in the recombinant protein (Hitchcock-DeGregori, S. E., and Heald, R. W. (1987) J. Biol. Chem. 262, 9730-9735). We expressed three fusion tropomyosins in E. coli with 2, 3, and 17 amino acids fused to its amino terminus. Ah three fusions restored actin binding, head-to-tail polymerization, and the capacity to inhibit the actomyosin ATPase to these unacetylated tropomyosins. Unlike larger fusions, the small fusions of 2 and 3 amino acids do not interfere with regulatory function. Therefore the presence of a fused dipeptide at the amino terminus of unacetylated tropomyosin is sufficient to replace the function of the N-acetyl group present in muscle tropomyosin. A structural interpretation for the function of the acetyl group, based on our results and the coiled coil structure of tropomyosin, is presented.

Mutational analyses of human eIF5A-1-identification of amino acid residues critical for eIF5A activity and hypusine modification

The eukaryotic translation initiation factor 5A (eIF5A) is the only protein that contains hypusine [N-epsilon-(4-amino-2-hydroxybutyl)lysine], which is required for its activity. Hypusine is formed by post-translational modification of one specific lysine (Lys50 for human eIF5A) by deoxyhypusine synthase and deoxyhypusine hydroxylase. To investigate the features of eIF5A required for its activity, we generated 49 mutations in human eIF5A-1, with a single amino acid substitution at the highly conserved residues or with N-terminal or C-terminal truncations, and tested mutant proteins in complementing the growth of a Saccharomyces cerevisiae eIF5A null strain. Growth-supporting activity was abolished in only a few mutant eIF5As (K47D, G49A, K50A, K50D, K50I, K50R, G52A and K55A), with substitutions at or near the hypusine modification site or with truncation of 21 amino acids from either the N-terminus or C-terminus. The inactivity of the Lys50 substitution proteins is obviously due to lack of deoxyhypusine modification. In contrast, K47D and G49A were effective substrates for deoxyhypusine synthase, yet failed to support growth, suggesting critical roles of Lys47 and Gly49 in eIF5A activity, possibly in its interaction with effector(s). By use of a UBHY-R strain harboring genetically engineered unstable eIF5A, we present evidence for the primary function of eIF5A in protein synthesis. When selected eIF5A mutant proteins were tested for their activity in protein synthesis, a close correlation was observed between their ability to enhance protein synthesis and growth, lending further support for a central role of eIF5A in translation.

The Mycobacterium tuberculosis Rv2540c DNA sequence encodes a bifunctional chorismate synthase

Background. The emergence of multi- and extensively-drug resistant Mycobacterium tuberculosis strains has created an urgent need for new agents to treat tuberculosis (TB). The enzymes of shikimate pathway are attractive targets to the development of antitubercular agents because it is essential for M. tuberculosis and is absent from humans. Chorismate synthase (CS) is the seventh enzyme of this route and catalyzes the NADH- and FMN-dependent synthesis of chorismate, a precursor of aromatic amino acids, naphthoquinones, menaquinones, and mycobactins. Although the M. tuberculosis Rv2540c (aroF) sequence has been annotated to encode a chorismate synthase, there has been no report on its correct assignment and functional characterization of its protein product. Results. In the present work, we describe DNA amplification of aroF-encoded CS from M. tuberculosis (MtCS), molecular cloning, protein expression, and purification to homogeneity. N-terminal amino acid sequencing, mass spectrometry and gel filtration chromatography were employed to determine identity, subunit molecular weight and oligomeric state in solution of homogeneous recombinant MtCS. The bifunctionality of MtCS was determined by measurements of both chorismate synthase and NADH:FMN oxidoreductase activities. The flavin reductase activity was characterized, showing the existence of a complex between FMN ox and MtCS. FMNox and NADH equilibrium binding was measured. Primary deuterium, solvent and multiple kinetic isotope effects are described and suggest distinct steps for hydride and proton transfers, with the former being more rate-limiting. Conclusion. This is the first report showing that a bacterial CS is bifunctional. Primary deuterium kinetic isotope effects show that C4-proS hydrogen is being transferred during the reduction of FMNox by NADH and that hydride transfer contributes significantly to the rate-limiting step of FMN reduction reaction. Solvent kinetic isotope effects and proton inventory results indicate that proton transfer from solvent partially limits the rate of FMN reduction and that a single proton transfer gives rise to the observed solvent isotope effect. Multiple isotope effects suggest a stepwise mechanism for the reduction of FMNox. The results on enzyme kinetics described here provide evidence for the mode of action of MtCS and should thus pave the way for the rational design of antitubercular agents. © 2008 Ely et al; licensee BioMed Central Ltd.

Positive selection results in frequent reversible amino acid replacements in the G protein gene of human respiratory syncytial virus

Human respiratory syncytial virus (HRSV) is the major cause of lower respiratory tract infections in children under 5 years of age and the elderly, causing annual disease outbreaks during the fall and winter. Multiple lineages of the HRSVA and HRSVB serotypes co-circulate within a single outbreak and display a strongly temporal pattern of genetic variation, with a replacement of dominant genotypes occurring during consecutive years. In the present study we utilized phylogenetic methods to detect and map sites subject to adaptive evolution in the G protein of HRSVA and HRSVB. A total of 29 and 23 amino acid sites were found to be putatively positively selected in HRSVA and HRSVB, respectively. Several of these sites defined genotypes and lineages within genotypes in both groups, and correlated well with epitopes previously described in group A. Remarkably, 18 of these positively selected tended to revert in time to a previous codon state, producing a flipflop phylogenetic pattern. Such frequent evolutionary reversals in HRSV are indicative of a combination of frequent positive selection, reflecting the changing immune status of the human population, and a limited repertoire of functionally viable amino acids at specific amino acid sites.

Structural insights into selectivity and cofactor binding in snake venom l-amino acid oxidases

l-Amino acid oxidases (LAAOs) are flavoenzymes that catalytically deaminate l-amino acids to corresponding α-keto acids with the concomitant production of ammonia (NH 3) and hydrogen peroxide (H 2O 2). Particularly, snake venom LAAOs have been attracted much attention due to their diverse clinical and biological effects, interfering on human coagulation factors and being cytotoxic against some pathogenic bacteria and Leishmania ssp. In this work, a new LAAO from Bothrops jararacussu venom (BjsuLAAO) was purified, functionally characterized and its structure determined by X-ray crystallography at 3.1å resolution. BjsuLAAO showed high catalytic specificity for aromatic and aliphatic large side-chain amino acids. Comparative structural analysis with prokaryotic LAAOs, which exhibit low specificity, indicates the importance of the active-site volume in modulating enzyme selectivity. Surprisingly, the flavin adenine dinucleotide (FAD) cofactor was found in a different orientation canonically described for both prokaryotic and eukaryotic LAAOs. In this new conformational state, the adenosyl group is flipped towards the 62-71 loop, being stabilized by several hydrogen-bond interactions, which is equally stable to the classical binding mode. © 2012 Elsevier Inc.

Characterization and fatty acids profile of the oils from Amazon nuts and walnuts: Characterization and fatty acids profile of the oilseeds

Purpose: The purpose of this paper is to characterize oils extracted from Bertholletia excelsa, Lecythis pisonis, Dipteryx lacunifera, Carya illinoensis and Juglans regia, regarding their characterization and fatty acid profile. Design/methodology/approach: The oils were extracted from oilseeds by cold pressing and physico-chemical characterization was performed by using standard methods for oils and fats. The oxidative stability and fatty acid profile also were determined. Findings: According to the results, the physico-chemical properties of oils from nuts and walnuts were comparable to those of good quality conventional oils. The oil seeds are a good source of unsaturated fatty acids, especially oleic and linoleic acids. Research limitations/implications: Implies the identification of fatty acid profile and physico-chemical properties of oils extracted from nuts and walnuts, and to prevent certain types of diseases. Originality/value: The paper identifies a new source of essential fatty acids extracted from oilseeds. © Emerald Group Publishing Limited.

Development of sustainable precision farming systems for swine: Estimating realtime individual amino acid requirements in growing-finishing pigs

The objective of this study was to develop and evaluate a mathematical model used to estimate the daily amino acid requirements of individual growing-finishing pigs. The model includes empirical and mechanistic model components. The empirical component estimates daily feed intake (DFI), BW, and daily gain (DG) based on individual pig information collected in real time. Based on DFI, BW, and DG estimates, the mechanistic component uses classic factorial equations to estimate the optimal concentration of amino acids that must be offered to each pig to meet its requirements. The model was evaluated with data from a study that investigated the effect of feeding pigs with a 3-phase or daily multiphase system. The DFI and BW values measured in this study were compared with those estimated by the empirical component of the model. The coherence of the values estimated by the mechanistic component was evaluated by analyzing if it followed a normal pattern of requirements. Lastly, the proposed model was evaluated by comparing its estimates with those generated by the existing growth model (InraPorc). The precision of the proposed model and InraPorc in estimating DFI and BW was evaluated through the mean absolute error. The empirical component results indicated that the DFI and BW trajectories of individual pigs fed ad libitum could be predicted 1 d (DFI) or 7 d (BW) ahead with the average mean absolute error of 12.45 and 1.85%, respectively. The average mean absolute error obtained with the InraPorc for the average individual of the population was 14.72% for DFI and 5.38% for BW. Major differences were observed when estimates from InraPorc were compared with individual observations. The proposed model, however, was effective in tracking the change in DFI and BW for each individual pig. The mechanistic model component estimated the optimal standardized ileal digestible Lys to NE ratio with reasonable between animal (average CV = 7%) and overtime (average CV = 14%) variation. Thus, the amino acid requirements estimated by model are animal- and time-dependent and follow, in real time, the individual DFI and BW growth patterns. The proposed model can follow the average feed intake and feed weight trajectory of each individual pig in real time with good accuracy. Based on these trajectories and using classical factorial equations, the model makes it possible to estimate dynamically the AA requirements of each animal, taking into account the intake and growth changes of the animal. © 2012 American Society of Animal Science. All rights reserved.

Purification, characterization, and specificity determination of a new serine protease secreted by penicillium waksmanii

The purpose of this work was to purify a protease from Penicillium waksmanii and to determine its biochemical characteristics and specificity. The extracellular protease isolated that was produced by P. waksmanii is a serine protease that is essential for the reproduction and growth of the fungus. The protease isolated showed 32 kDa, and has optimal activity at pH 8.0 and 35 C towards the substrate Abz-KLRSSKQ-EDDnp. The protease is active in the presence of CaCl2, KCl, and BaCl, and partially inhibited by CuCl2, CoCl2 and totally inhibited by AlCl3 and LiCl. In the presence of 1 M urea, the protease remains 50 % active. The activity of the protease increases 60 % when it is exposed to 0.4 % nonionic surfactant-Triton X-100 and loses 10 % activity in the presence of 0.4 % Tween-80. Using fluorescence resonance energy transfer analysis, the protease showed the most specificity for the peptide Abz-KIRSSKQ-EDDnp with k cat/K m of 10,666 mM-1 s-1, followed by the peptide Abz-GLRSSKQ-EDDnp with a k cat/K m of 7,500 mM -1 s-1. Basic and acidic side chain-containing amino acids performed best at subsite S1. Subsites S2, S3, S′ 2, and S′ 1, S ′ 3 showed a preference for binding for amino acids with hydrophobic and basic amino acid side chain, respectively. High values of k cat/K m were observed for the subsites S2, S3, and S′ 2. The sequence of the N-terminus (ANVVQSNVPSWGLARLSSKKTGTTDYTYD) showed high similarity to the fungi Penicillium citrinum and Penicillium chrysogenum, with 89 % of identity at the amino acid level. © 2012 Springer Science+Business Media New York.

EIF5A dimerizes not only in vitro but also in vivo and its molecular envelope is similar to the EF-P monomer

The protein eukaryotic initiation factor 5A (eIF5A) is highly conserved among archaea and eukaryotes, but not in bacteria. Bacteria have the elongation factor P (EF-P), which is structurally and functionally related to eIF5A. eIF5A is essential for cell viability and the only protein known to contain the amino acid residue hypusine, formed by post-translational modification of a specific lysine residue. Although eIF5A was initially identified as a translation initiation factor, recent studies strongly support a function for eIF5A in the elongation step of translation. However, the mode of action of eIF5A is still unknown. Here, we analyzed the oligomeric state of yeast eIF5A. First, by using size-exclusion chromatography, we showed that this protein exists as a dimer in vitro, independent of the hypusine residue or electrostatic interactions. Protein-protein interaction assays demonstrated that eIF5A can form oligomers in vitro and in vivo, in an RNA-dependent manner, but independent of the hypusine residue or the ribosome. Finally, small-angle X-ray scattering (SAXS) experiments confirmed that eIF5A behaves as a stable dimer in solution. Moreover, the molecular envelope determined from the SAXS data shows that the eIF5A dimer is L-shaped and superimposable on the tRNAPhe tertiary structure, analogously to the EF-P monomer. © 2012 Springer-Verlag.