Genomic distribution and heterogeneity of MocR-like transcriptional factors containing a domain belonging to the superfamily of the pyridoxal-5'-phosphate dependent enzymes of fold type I.

Bacterial proteins belonging to the MocR/GabR family are chimeric proteins incorporating a short N-terminal helix-turn-helix containing domain with DNA-binding properties, and a long C-terminal domain belonging to the superfamily of the pyridoxal-5'-phosphate enzymes of fold type I. The first purpose of this report is to give an overview of the distribution of these factors among the different taxonomical bacterial divisions and to determine the degree of conservation of the main structural features of the PLP binding domain. Complete proteomes of bacteria phyla were scanned with a hidden Markov model representative of the MocR family. Results indicate that presence of MocR factors is heterogeneous even within the single bacterial phylum: some species miss completely the factors, while others possess one or even more regulators. Absence of MocR factors is distinctive of some phyla such as Chlamydiae. The genomic distribution of MocR is, as expected, highly correlated to the size of the genome. At variance, phyla missing MocR regulators generally are characterized by compact genomes, of the order of 1.0-2.0 Mb, such as the case of Mollicutes or Chlamydiae. Apparently, the minimum genome size compatible with the presence of MocR genes is around 2.0-2.5 Mb. Conservation of the residues corresponding to those involved in the interaction with the cofactor pyridoxal-5'-phosphate in the homologous 2-aminoadipate aminotransferase, was analyzed in the multiple sequence alignments of MocR within each phyla considered. In the vast majority of cases, residues are conserved or conservatively replaced. This result suggests that, in most cases, MocR factors preserve at least ability to bind the cofactor and very likely some catalytic abilities.

Breast Tumor Malignancy Classification using Smartphone Compression-induced Sensing System and Deformation Index Ratio.

Smartphone-based compression-induced sensing system uses the light diffusion pattern to characterize the early-stage breast tumor noninvasively. The system is built on a smartphone and cloud platform to capture, transfer, and interface with the user. The compressed tissue's deforming pattern creates distinctive tactile images due to the size and hardness of the tumor. From the compression-induced images, we estimate the size of the tumor using projection analysis and the tumor's malignancy using the tissue deformation index ratio. Deformation index ratio is based on the changes of a healthy region over the tumorous region. By using the projection analysis, the human patient tumor size estimation resulted in 52.3% of the average error. For a small number (seven) of the feasibility test, the tumor's malignancy was classified based on the deformation index ratio with 67.0% of sensitivity and 100% specificity.

The primary structure of mitochondrial aspartate aminotransferase from human heart.

The complete amino acid sequence of the mitochondrial aspartate aminotransferase (L-aspartate:2-oxoglutarate aminotransferase, EC 2.6.1.1) from human heart has been determined based mainly on analysis of peptides obtained by digestion with trypsin and by chemical cleavage with cyanogen bromide. Comparison of the sequence with those of the isotopic isoenzymes from pig, rat and chicken showed 27, 29 and 55 differences, respectively, out of a total of 401 amino acid residues. Evidence for structural microheterogeneity at position 317 has also been obtained.

Analysis of insertions/deletions in protein structures.

An analysis of insertions and deletions (indels) occurring in a databank of multiple sequence alignments based on protein tertiary structure is reported. Indels prefer to be short (1 to 5 residues). The average intervening sequence length between them versus the percentage of residue identity in pairwise alignments shows an exponential behaviour, suggesting a stochastic process such that nearly every loop in an ancestral structure is a possible target for indels during evolution. The results also suggest a limit to the average size of indels accommodated by protein structures. The preferred indel conformations are reverse turn and coil as are the preferred conformations at the indel edges (N- and C-terminal sides). Interruptions in helices and strands were observed as very rare events.

Similarity between pyridoxal/pyridoxamine phosphate-dependent enzymes involved in dideoxy and deoxyaminosugar biosynthesis and other pyridoxal phosphate enzymes.

A multiple sequence alignment among aspartate aminotransferase, dialkylglycine decarboxylase, and serine hydroxymethyltransferase (DAS) was used for profile databank search. The DAS profile could detect similarities to other pyridoxal or pyridoxamine phosphate-dependent enzymes, like several gene products involved in dideoxysugar and deoxyaminosugar synthesis. The alignment among DAS and such gene products shows the conservation of aspartate 222 and lysine 258, which, in aspartate aminotransferase, interacts with the N1 of the coenzyme pyridine ring and forms the internal Schiff base, respectively. The lysine is replaced by histidine in the pyridoxamine phosphate-dependent gene products. The alignment indicates also that the region encompassing the coenzyme binding site is the most conserved.

Active site model for gamma-aminobutyrate aminotransferase explains substrate specificity and inhibitor reactivities.

A homology model for the pig isozyme of the pyridoxal phosphate-dependent enzyme gamma-aminobutyrate (GABA) aminotransferase has been built based mainly on the structure of dialkylglycine decarboxylase and on a multiple sequence alignment of 28 evolutionarily related enzymes. The proposed active site structure is presented and analyzed. Hypothetical structures for external aldimine intermediates explain several characteristics of the enzyme. In the GABA external aldimine model, the pro-S proton at C4 of GABA, which abstracted in the 1,3-azaallylic rearrangement interconverting the aldimine and ketimine intermediates, is oriented perpendicular to the plane of the pyridoxal phosphate ring. Lys 329 is in close proximity and is probably the general base catalyst for the proton transfer reaction. The carboxylate group of GABA interacts with Arg 192 and Lys 203, which determine the specificity of the enzyme for monocarboxylic omega-amino acids such as GABA. In the proposed structure for the L-glutamate external aldimine, the alpha-carboxylate interacts with Arg 445. Glu 265 is proposed to interact with this same arginine in the GABA external aldimine, enabling the enzyme to act on omega-amino acids in one half-reaction and on alpha-amino acids in the other. The reactivities of inhibitors are well explained by the proposed active site structure. The R and S isomers of beta-substituted phenyl and p-chlorophenyl GABA would bind in very different modes due to differential steric interactions, with the reactive S isomer leaving the orientation of the GABA moiety relatively unperturbed compared to that of the natural substrate. In our model, only the reactive S isomer of the mechanism-based inhibitor vinyl-GABA, an effective anti-epileptic drug known clinically as Vigabatrin, would orient the scissile C4-H bond perpendicular to the coenzyme ring plane and present the proton to Lys 329, the proposed general base catalyst of the reaction. The R isomer would direct the vinyl group toward Lys 329 and the C4-H bond toward Arg 445. The active site model presented provides a basis for site-directed mutagenesis and drug design experiments.

BSTI, a trypsin inhibitor from skin secretions of Bombina bombina related to protease inhibitors of nematodes.

From skin secretions of the European frog Bombina bombina, a new peptide has been isolated that contains 60 amino acids, including 10 cysteine residues. Its sequence was determined by automated Edman degradation and confirmed by analysis of the cDNA encoding the precursor. A search in the databanks demonstrated that the pattern of cysteine residues in this skin peptide is similar to the ones found in protease inhibitors from Ascaris and in a segment of human von Willebrand factor. The 3D structure of the trypsin inhibitor from Ascaris suum could be used as a template to build a model of the amphibian peptide. In addition, we have demonstrated that this constituent of skin secretion is indeed an inhibitor of trypsin and thrombin, with K(i) values in the range of 0.1 to 1 microM. The new peptide was thus named BSTI for Bombina skin trypsin/thrombin inhibitor.

The structure of serine hydroxymethyltransferase as modeled by homology and validated by site-directed mutagenesis.

We describe a model for the three-dimensional structure of E. coli serine hydroxymethyltransferase based on its sequence homology with other PLP enzymes of the alpha-family and whose tertiary structures are known. The model suggests that certain amino acid residues at the putative active site of the enzyme can adopt specific roles in the catalytic mechanism. These proposals were supported by analysis of the properties of a number of site-directed mutants. New active site features are also proposed for further experimental testing.

Similarity between serine hydroxymethyltransferase and other pyridoxal phosphate-dependent enzymes.

A structural homology of the pyridoxal-5'-phosphate (PLP)-dependent enzyme serine hydroxymethyltransferase (SHMT) with aspartate aminotransferase (AAT) is proposed. Although the two sequences are very dissimilar, a reasonable alignment was obtained using the profile analysis method. Sequences of AAT and dialkylglycine decarboxylase (DGD), for which crystal structure data are available, have been aligned on the basis of their structure superposition. A profile was then calculated and SHMT sequence aligned to it. Three of the four residues conserved in all aminotransferases (including the PLP-binding lysine) are matched. A profile search with DGD-AAT-SHMT profile is more selective and sensitive than individual sequence profiles for PLP-dependent enzyme detection. Potential homologies with the eryC1 gene product involved in erythromycin biosynthesis and with amino acid decarboxylases were observed. Homology with AAT will be used as a guideline for planning site-directed mutagenesis experiments on SHMT.

Selective inhibition of human mast cell tryptase by gabexate mesylate, an antiproteinase drug.

Gabexate mesylate is a non-antigenic synthetic inhibitor of trypsin-like serine proteinases that is therapeutically used in the treatment of pancreatitis and disseminated intravascular coagulation and as a regional anticoagulant for hemodialysis. Considering the structural similarity between gabexate mesylate and arginine-based inhibitors of trypsin-like serine proteinases, the effect of gabexate mesylate on human and bovine mast cell tryptase action was investigated. Values of the inhibition constant (K(i)) for gabexate mesylate binding to human and bovine tryptase were 3.4 x 10(-9) M and 1.8 x 10(-7) M (at pH 7.4 and 37.0 degrees ), respectively. Furthermore, gabexate mesylate inhibited the fibrinogenolytic activity of human tryptase. On the basis of the available x-ray crystal structure of human tryptase, the possible binding mode of gabexate mesylate to human and bovine tryptase was analyzed. Human tryptase inhibition by gabexate mesylate may account for the reported prevention of inflammation, erosion, and ulceration of skin and mucosae.

A simple microcomputer program for predicting the secondary structure of proteins.

An Applesoft BASIC program is presented which applies the Chou-Fasman and Garnier-Osguthorpe-Robson rules for predicting the secondary structure of proteins. This program plots the result of predictions and shows the charge distribution and the consensus prediction. The program is run on an Apple IIe and the results are printed by an Imagewriter printer.

Structural analysis of toxic volkensin, a type 2 ribosome inactivating protein from Adenia volkensii Harm (kilyambiti plant): molecular modeling and surface analysis by computational methods and limited proteolysis.

Volkensin, isolated from Adenia volkensii, is one of the most toxic type 2 ribosome-inactivating protein (RIP), exerting its biological function by inhibiting protein synthesis. Despite the high sequence identity with type 2 RIPs, including ricin, volkensin shows interesting peculiar properties. In this work a computational model building of volkensin was performed. The volkensin electrostatic potential charge distribution, the hydrophobic profile and the surface topology analyses were also carried out to aid the understanding of structure-function relationships of this potent toxin. Volkensin surface topology was probed by applying a limited proteolysis approach with the aim to gain insights into volkensin conformational features.

PROTEUS: a suite of programs for prediction of structural features of proteins using an Apple IIe.

We have implemented several algorithms, developed by various authors for predicting structural features of proteins from their primary structure, on an Apple IIe and collected them in a suite, named PROTEUS. This suite incorporates: (i) methods for predicting secondary structure; (ii) the algorithm for computing the hydropathy profile using one out of five available sets of parameters; (iii) the algorithms for calculating the hydrophobic moment plot; and (iv) for performing the amphipathic analysis using one out of four available sets of parameters. The suite has a utility program for storing on a disk the sequence to be analysed. As an example, we applied some of the methods included in PROTEUS to predict the structure of a mitochondrial leader peptide. The results suggest the occurrence of structural features possibly related to the import of proteins into mitochondria.

A data bank merging related protein structures and sequences.

A data collection which merges protein structural and sequence information is described. Structural superpositions amongst proteins with similar main-chain fold were performed or collected from the literature. Sequences taken from the protein primary structure databases were associated with the multiple structural alignments providing they were at least 50% homologous in residue identity to one of the structural sequences and at least 50% of the structural sequence residues were alignable. Such restrictions allow reasonable confidence that the primary sequences share the conformation of the tertiary structural templates, except in the less conserved loop regions. Multiple structural superpositions were collected for 38 familial groups containing a total of 209 tertiary structures; 45 structures had no superposable mates and were used individually. Other information is also provided as main-chain and side-chain conformational angles, secondary structural assignments and the like. Wedding the primary and tertiary structural data resulted in an 8-fold increase of data bank sequence entries over those associated with the known three-dimensional architectures alone.

Conservation of amphipathic conformations in multiple protein structural alignments.

Protein amphipathic conformations, mainly alpha-helices and beta-strands, are believed to play an important role in protein folding, stability and function. The most popular method for characterizing such structures is the hydrophobic moment. We have analyzed the distribution of hydrophobic moment characteristics (peak magnitude, amphipathic indices and characteristic frequency) in a data bank containing several families of distant sequences multiply aligned by structural superposition. Sequence fragments were classified according to alpha-helix, beta-strand, non-alpha and non-beta conformations. This data bank provided an enhanced sample space compared with those previously reported in the literature. Precautions were taken to reduce over-representation of homologous sequences. Approximately 50% of all individual alpha-helices showed a hydrophobic moment peak in the expected position of the periodicity spectrum while only 38% of individual beta-strands fell in the expected range. False positives account for a surprisingly large 14 and 36% of the non-alpha and non-beta samples respectively. Conservation of hydrophobic moment characteristics and mainly the hydrophobic peak position in the expected periodicity range was examined in the multiple alignments of the distant sequences. Helices tend to conserve more frequently their hydrophobic moment than any other conformation and yet only 13% of all helical segments display such conservation in three-quarters or more of the familial sequences; the similar observation for beta-strands was even lower at 9%. Nonetheless, strongly hydrophobic positions within the structural segments were more conserved than expected.

Structural adaptation of enzymes to low temperatures.

A systematic comparative analysis of 21 psychrophilic enzymes belonging to different structural families from prokaryotic and eukaryotic organisms is reported. The sequences of these enzymes were multiply aligned to 427 homologous proteins from mesophiles and thermophiles. The net flux of amino acid exchanges from meso/thermophilic to psychrophilic enzymes was measured. To assign the observed preferred exchanges to different structural environments, such as secondary structure, solvent accessibility and subunit interfaces, homology modeling was utilized to predict the secondary structure and accessibility of amino acid residues for the psychrophilic enzymes for which no experimental three-dimensional structure is available. Our results show a clear tendency for the charged residues Arg and Glu to be replaced at exposed sites on alpha-helices by Lys and Ala, respectively, in the direction from 'hot' to 'cold' enzymes. Val is replaced by Ala at buried regions in alpha-helices. Compositional analysis of psychrophilic enzymes shows a significant increase in Ala and Asn and a decrease in Arg at exposed sites. Buried sites in beta-strands tend to be depleted of VAL: Possible implications of the observed structural variations for protein stability and engineering are discussed.

Easy method to predict solvent accessibility from multiple protein sequence alignments.

An easy and uncomplicated method to predict the solvent accessibility state of a site in a multiple protein sequence alignment is described. The approach is based on amino acid exchange and compositional preference matrices for each of three accessibility states: buried, exposed, and intermediate. Calculations utilized a modified version of the 3D_ali databank, a collection of multiple sequence alignments anchored through protein tertiary structural superpositions. The technique achieves the same accuracy as much more complex methods and thus provides such advantages as computational affordability, facile updating, and easily understood residue substitution patterns useful to biochemists involved in protein engineering, design, and structural prediction. The program is available from the authors; and, due to its simplicity, the algorithm can be readily implemented on any system. For a given alignment site, a hand calculation can yield a comparative prediction.

Quantification of secondary structure prediction improvement using multiple alignments.

The use of multiple sequence alignments for secondary structure predictions is analysed. Seven different protein families, containing only sequences of known structure, were considered to provide a range of alignment and prediction conditions. Using alignments obtained by spatial superposition of main chain atoms in known tertiary protein structures allowed a mean of 8% in secondary structure prediction accuracy, when compared to those obtained from the individual sequences. Substitution of these alignments by those determined directly from an automated sequence alignment algorithm showed variations in the prediction accuracy which correlated with the quality of the multiple alignments and distance of the primary sequence. Secondary structure predictions can be reliably improved using alignments from an automatic alignment procedure with a mean increase of 6.8%, giving an overall prediction accuracy of 68.5%, if there is a minimum of 25% sequence identity between all sequences in a family.

The primary structure of rabbit liver cytosolic serine hydroxymethyltransferase.

The complete amino acid sequence of cytosolic serine hydroxymethyltransferase from rabbit liver was determined. The sequence was determined from analysis of peptides isolated from tryptic and cyanogen bromide cleavages of the enzyme. Special procedures were used to isolate and sequence the C-terminal and blocked N-terminal peptides. Each of the four identical subunits of the enzyme consists of 483 residues. The sequence could be easily aligned with the sequence of Escherichia coli serine hydroxymethyltransferase. The primary structural homology between the rabbit and E. coli enzymes is about 42%. The importance of the primary and predicted secondary structural homology between the two enzymes is discussed.

Serine hydroxymethyltransferase: origin of substrate specificity.

All forms of serine hydroxymethyltransferase, for which a primary structure is known, have five threonine residues near the active-site lysyl residue (K229) that forms the internal aldimine with pyridoxal phosphate. For Escherichia coli serine hydroxymethyltransferase each of these threonine residues has been changed to an alanine residue. The resulting five mutant enzymes were purified and characterized with respect to kinetic and spectral properties. The mutant enzymes T224A and T227A showed no significant changes in kinetic and spectral properties compared to the wild-type enzyme. The T225A and T230A enzymes exhibited differences in Km and kcat values but exhibited the same spectral properties as the wild-type enzyme. The four threonine residues at positions 224, 225, 227, and 230 do not play a critical role in the mechanism of the enzyme. The T226A enzyme had nearly normal affinity for substrates and coenzymes but had only 3% of the catalytic activity of the wild-type enzyme. The spectrum of the T226A enzyme in the presence of amino acid substrates showed a large absorption maximum at 343 nm with only a small absorption band at 425 nm, unlike the wild-type enzyme whose enzyme-substrate complexes absorb at 425 nm. Rapid reaction studies showed that when amino acid substrates and substrate analogues were added to the T226A enzyme, the internal aldimine absorbing at 422 nm was rapidly converted to a complex absorbing at 343 nm in a second-order process. This was followed by a very slow first-order formation of a complex absorbing at 425 nm.(ABSTRACT TRUNCATED AT 250 WORDS)