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)

Antimicrobial activity of an L-amino acid oxidase isolated from Bothrops leucurus snake venom

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Isolation of a new L-amino acid oxidase from Crotalus durissus cascavella venom

A novel L-amino acid oxidase (LAO) (Casca LAO) from Crotalus durissus cascavella venom was purified to a high degree of molecular homogeneity using a combination of molecular exclusion and ion-exchange chromatography system. The purified monomer of LAO presented a molecular mass of 68 kDa and pI estimated in 5.43, which were determined by two-dimensional electrophoresis. The 71st N-terminal amino acid sequence of the LAO from Crotalus durissus cascavella presented a high amino acid sequence similarities with other LAOs from Colloselasma rhosostoma, Crotalus adamanteus, Agkistrodon h. blomhoffi, Agkistrodon h. halys and Trimeresurus stejnegeri. LAO displayed a Michaelis-Menten behavior with a kilometer of 46.7 mu M and an optimum pH for enzymatic activity of 6.5. Casca LAO induced a dose-dependent platelet aggregation, which was abolished by catalase and inhibited by indomethacin and aspirin. These results suggest that the production of H2O2 is involved in subsequent activation of inflammatory enzymes, such as thromboxane. Casca LAO also inhibited the bacterial Growth of Gram-negative (Xanthomonas axonopodis pv passiflorae) and Gram-positive (S. mutans) strains. Electron microscopy assessments of both bacterial strains suggest that the hydrogen peroxide produced by LAO induce bacterial membrane rupture and consequently loss of cytoplasmatic content. This LAO exhibited a high antileishmanic activity against the promastigote of Leishmania amazonensis in vitro, its activity was dependent on the production of hydrogen peroxide, and the 50% inhibitory concentration was estimated in 2.39 mu g/ml. (C) 2005 Elsevier Ltd. All rights reserved.

Eucalyptus ESTs associated with resistance to herbicide inhibitors of aromatic and branched-chain amino acid synthesis

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

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.

Functional and structural characterization of the first prokaryotic member of the L-amino acid transporter (LAT) family: a model for APC transporters

We have identified YkbA from Bacillus subtilis as a novel member of the L-amino acid transporter (LAT) family of amino acid transporters. The protein is approximately 30% identical in amino acid sequence to the light subunits of human heteromeric amino acid transporters. Purified His-tagged YkbA from Escherichia coli membranes reconstituted in proteoliposomes exhibited sodium-independent, obligatory exchange activity for L-serine and L-threonine and also for aromatic amino acids, albeit with less activity. Thus, we propose that YkbA be renamed SteT (Ser/Thr exchanger transporter). Kinetic analysis supports a sequential mechanism of exchange for SteT. Freeze-fracture analysis of purified, functionally active SteT in proteoliposomes, together with blue native polyacrylamide gel electrophoresis and transmission electron microscopy of detergent-solubilized purified SteT, suggest that the transporter exists in a monomeric form. Freeze-fracture analysis showed spherical particles with a diameter of 7.4 nm. Transmission electron microscopy revealed elliptical particles (diameters 6 x 7 nm) with a distinct central depression. To our knowledge, this is the first functional characterization of a prokaryotic member of the LAT family and the first structural data on an APC (amino acids, polyamines, and choline for organocations) transporter. SteT represents an excellent model to study the molecular architecture of the light subunits of heteromeric amino acid transporters and other APC transporters.

Identification of amino acid substitutions with compensational effects in the attachment protein of canine distemper virus.

The hemagglutinin (H) gene of canine distemper virus (CDV) encodes the receptor-binding protein. This protein, together with the fusion (F) protein, is pivotal for infectivity since it contributes to the fusion of the viral envelope with the host cell membrane. Of the two receptors currently known for CDV (nectin-4 and the signaling lymphocyte activation molecule [SLAM]), SLAM is considered the most relevant for host susceptibility. To investigate how evolution might have impacted the host-CDV interaction, we examined the functional properties of a series of missense single nucleotide polymorphisms (SNPs) naturally accumulating within the H-gene sequences during the transition between two distinct but related strains. The two strains, a wild-type strain and a consensus strain, were part of a single continental outbreak in European wildlife and occurred in distinct geographical areas 2 years apart. The deduced amino acid sequence of the two H genes differed at 5 residues. A panel of mutants carrying all the combinations of the SNPs was obtained by site-directed mutagenesis. The selected mutant, wild type, and consensus H proteins were functionally evaluated according to their surface expression, SLAM binding, fusion protein interaction, and cell fusion efficiencies. The results highlight that the most detrimental functional effects are associated with specific sets of SNPs. Strikingly, an efficient compensational system driven by additional SNPs appears to come into play, virtually neutralizing the negative functional effects. This system seems to contribute to the maintenance of the tightly regulated function of the H-gene-encoded attachment protein. Importance: To investigate how evolution might have impacted the host-canine distemper virus (CDV) interaction, we examined the functional properties of naturally occurring single nucleotide polymorphisms (SNPs) in the hemagglutinin gene of two related but distinct strains of CDV. The hemagglutinin gene encodes the attachment protein, which is pivotal for infection. Our results show that few SNPs have a relevant detrimental impact and they generally appear in specific combinations (molecular signatures). These drastic negative changes are neutralized by compensatory mutations, which contribute to maintenance of an overall constant bioactivity of the attachment protein. This compensational mechanism might reflect the reaction of the CDV machinery to the changes occurring in the virus following antigenic variations critical for virulence.

STUDIES ON THE PATTERNS OF NUCLEOTIDE AND AMINO ACID SUBSTITUTION (ELECTROPHORESIS, MUTATION RATE, SELECTION, BASE, COMPOSITION)

Theoretical and empirical studies were conducted on the pattern of nucleotide and amino acid substitution in evolution, taking into account the effects of mutation at the nucleotide level and purifying selection at the amino acid level. A theoretical model for predicting the evolutionary change in electrophoretic mobility of a protein was also developed by using information on the pattern of amino acid substitution. The specific problems studied and the main results obtained are as follows: (1) Estimation of the pattern of nucleotide substitution in DNA nuclear genomes. The pattern of point mutations and nucleotide substitutions among the four different nucleotides are inferred from the evolutionary changes of pseudogenes and functional genes, respectively. Both patterns are non-random, the rate of change varying considerably with nucleotide pair, and that in both cases transitions occur somewhat more frequently than transversions. In protein evolution, substitution occurs more often between amino acids with similar physico-chemical properties than between dissimilar amino acids. (2) Estimation of the pattern of nucleotide substitution in RNA genomes. The majority of mutations in retroviruses accumulate at the reverse transcription stage. Selection at the amino acid level is very weak, and almost non-existent between synonymous codons. The pattern of mutation is very different from that in DNA genomes. Nevertheless, the pattern of purifying selection at the amino acid level is similar to that in DNA genomes, although selection intensity is much weaker. (3) Evaluation of the determinants of molecular evolutionary rates in protein-coding genes. Based on rates of nucleotide substitution for mammalian genes, the rate of amino acid substitution of a protein is determined by its amino acid composition. The content of glycine is shown to correlate strongly and negatively with the rate of substitution. Empirical formulae, called indices of mutability, are developed in order to predict the rate of molecular evolution of a protein from data on its amino acid sequence. (4) Studies on the evolutionary patterns of electrophoretic mobility of proteins. A theoretical model was constructed that predicts the electric charge of a protein at any given pH and its isoelectric point from data on its primary and quaternary structures. Using this model, the evolutionary change in electrophoretic mobilities of different proteins and the expected amount of electrophoretically hidden genetic variation were studied. In the absence of selection for the pI value, proteins will on the average evolve toward a mildly basic pI. (Abstract shortened with permission of author.) ^

A σ32 mutant with a single amino acid change in the highly conserved region 2.2 exhibits reduced core RNA polymerase affinity

σ32, the product of the rpoH gene in Escherichia coli, provides promoter specificity by interacting with core RNAP. Amino acid sequence alignment of σ32 with other sigma factors in the σ70 family has revealed regions of sequence homology. We have investigated the function of the most highly conserved region, 2.2, using purified products of various rpoH alleles. Core RNAP binding analysis by glycerol gradient sedimentation has revealed reduced core RNAP affinity for one of the mutant σ32 proteins, Q80R. This reduced core interaction is exacerbated in the presence of σ70, which competes with σ32 for binding of core RNAP. When a different but more conserved amino acid was introduced at this position by site-directed mutagenesis (Q80N), this mutant sigma factor still displayed a significant reduction in its core RNAP affinity. Based on these results, we conclude that at least one specific amino acid in region 2.2 is involved in core RNAP interaction.

Protective DNA vaccination against organ-specific autoimmunity is highly specific and discriminates between single amino acid substitutions in the peptide autoantigen

DNA vaccines that encode encephalitogenic sequences in tandem can protect from subsequent experimental autoimmune encephalomyelitis induced with the corresponding peptide. The mechanism for this protection and, in particular, if it is specific for the amino acid sequence encoding the vaccine are not known. We show here that a single amino acid exchange in position 79 from serine (nonself) to threonine (self) in myelin basic protein peptide MBP68–85, which is a major encephalitogenic determinant for Lewis rats, dramatically alters the protection. Moreover, vaccines encoding the encephalitogenic sequence MBP68–85 do not protect against the second encephalitogenic sequence MBP89–101 in Lewis rats and vice versa. Thus, protective immunity conferred by DNA vaccination exquisitely discriminates between peptide target autoantigens. No bystander suppression was observed. The exact underlying mechanisms remain elusive because no simple correlation between impact on ex vivo responses and protection against disease were noted.

Mapping rat megalin: the second cluster of ligand binding repeats contains a 46-amino acid pathogenic epitope involved in the formation of immune deposits in Heymann nephritis.

Megalin (gp330), an epithelial endocytic receptor, is a major target antigen of Heymann nephritis (HN), an autoimmune disease in rats. To elucidate the mechanisms of HN, we have mapped a pathogenic epitope in megalin that binds anti-megalin antibodies. We focused our attention on four clusters of cysteine-rich, low density lipoprotein receptor (LDLR) ligand binding repeats in the extracellular domain of megalin because they represent putative ligand binding regions and therefore would be expected to be exposed in vivo and to be able to bind circulating antibodies. Rat megalin cDNA fragments I through IV encoding the first through fourth clusters of ligand-binding repeats, respectively, were expressed in a baculovirus system. All four expression products were detected by immunoblotting with two antisera capable of inducing passive HN (pHN). When antibodies eluted from glomeruli of rats with pHN were used for immunoblotting, only the expression product encoded by fragment II was detected. This indicates that the second cluster of LDLR ligand binding repeats is directly involved in binding anti-megalin antibodies and in the induction of pHN. To narrow the major epitope in this domain, fragment II was used to prepare proteins sequentially truncated from the C- and N-terminal ends by in vitro translation. Analysis of the truncated translation products by immunoprecipitation with anti-megalin IgG revealed that the fifth ligand-binding repeat (amino acids 1160-1205) contains the major epitope recognized. This suggests that a 46-amino acid sequence in the second cluster of LDLR ligand binding repeats contains a major pathogenic epitope that plays a key role in pHN. Identification of this epitope will facilitate studies on the pathogenesis of HN.

D-amino acid residue in a defensin-like peptide from platypus venom: effect on structure and chromatographic properties

The recent discovery that the natriuretic peptide OvCNPb (Ornithorhynchus venom C-type natriuretic peptide B) from platypus (Ornithorynchus anatinus) venom contains a D-amino acid residue suggested that other D-amino-acid-containing peptides might be present in the venom. In the present study, we show that DLP-2 (defensin-like peptide-2), a 42-amino-acid residue polypeptide in the platypus venom, also contains a D-amino acid residue, D-methionine, at position 2, while DLP-4, which has an identical amino acid sequence, has all amino acids in the L-form. These findings were supported further by the detection of isomerase activity in the platypus gland venom extract that converts DLP-4 into DLP-2. In the light of this new information, the tertiary structure of DLP-2 was recalculated using a new structural template with D-Met(2). The structure of DLP-4 was also determined in order to evaluate the effect of a D-amino acid at position 2 on the structure and possibly to explain the large retention time difference observed for the two molecules in reverse-phase HPLC. The solution structures of the DLP-2 and DLP-4 are very similar to each other and to the earlier reported structure of DLP-2, which assumed that all amino acids were in the L-form. Our results suggest that the incorporation of the D-amino acid at position 2 has minimal effect on the overall fold in solution.

Deduced amino acid sequences surrounding the fusion glycoprotein cleavage site and of the carboxyl-terminus of haemagglutinin-neuraminidase protein of the avirulent thermostable vaccine strain I-2 of Newcastle disease virus

A single-tube RT-PCR technique generated a 387 bp or 300 bp cDNA amplicon covering the F-0 cleavage site or the carboxyl (C)-terminus of the HN gene, respectively, of Newcastle disease virus (NDV) strain 1-2. Sequence analysis was used to deduce the amino acid sequences of the cleavage site of F protein and the C-terminus of HN protein, which were then compared with sequences for other NDV strains. The cleavage site of NDV strain 1-2 had a sequence Motif of (112)RKQGRLIG(119), consistent with an avirulent phenotype. Nucleotide sequencing and deduction of amino acids at the C-terminus of HN revealed that strain 1-2 had a 7-amino-acid extension (VEILKDGVREARSSR). This differs from the virulent viruses that caused outbreaks of Newcastle disease in Australia in the 1930s and 1990s, which have HN extensions of 0 and 9 amino acids, respectively. Amino acid sequence analyses of the F and HN genes of strain 1-2 confirmed its avirulent nature and its Australian origin.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Amino acid selective unlabeling for sequence specific resonance assignments in proteins

Sequence specific resonance assignment constitutes an important step towards high-resolution structure determination of proteins by NMR and is aided by selective identification and assignment of amino acid types. The traditional approach to selective labeling yields only the chemical shifts of the particular amino acid being selected and does not help in establishing a link between adjacent residues along the polypeptide chain, which is important for sequential assignments. An alternative approach is the method of amino acid selective `unlabeling' or reverse labeling, which involves selective unlabeling of specific amino acid types against a uniformly C-13/N-15 labeled background. Based on this method, we present a novel approach for sequential assignments in proteins. The method involves a new NMR experiment named, {(CO)-C-12 (i) -N-15 (i+1)}-filtered HSQC, which aids in linking the H-1(N)/N-15 resonances of the selectively unlabeled residue, i, and its C-terminal neighbor, i + 1, in HN-detected double and triple resonance spectra. This leads to the assignment of a tri-peptide segment from the knowledge of the amino acid types of residues: i - 1, i and i + 1, thereby speeding up the sequential assignment process. The method has the advantage of being relatively inexpensive, applicable to H-2 labeled protein and can be coupled with cell-free synthesis and/or automated assignment approaches. A detailed survey involving unlabeling of different amino acid types individually or in pairs reveals that the proposed approach is also robust to misincorporation of N-14 at undesired sites. Taken together, this study represents the first application of selective unlabeling for sequence specific resonance assignments and opens up new avenues to using this methodology in protein structural studies.