Cortisol and behavioral responses of working dogs to environmental challenges.

This paper's primary objective is to analyse the physiological (cortisol) and behavioral responses of military working dogs (MWD). Dogs (N=27) were submitted twice to environmental challenges (challenge 1 and 2, 20 days in-between) composed of social (training), visual (mobile toy car) and auditory (air blast) stimuli. Cortisol levels decreased back to the baseline after the second challenge. The behavioral observations showed that these MWD were more active, and presented less stereotypic behaviors (pacing, manipulation of the environment) during both visual challenges, whereas half low posture was observed during the first but not during the second visual challenge. The present study shows that this group of MWD still has an adaptation capacity to an environmental challenge (return to baseline of the cortisol levels, a higher posture during the second than at the first challenge). These results are encouraging and indicate that the dogs might have a diminished welfare (i.e. stereotypic behaviors), but are not chronically stressed.

Crystallographic analysis of family 11 endo-beta-1,4-xylanase Xyl1 from Streptomyces sp. S38.

Family 11 endo-beta-1,4-xylanases degrade xylan, the main constituent of plant hemicelluloses, and have many potential uses in biotechnology. The structure of Xyl1, a family 11 endo-xylanase from Streptomyces sp. S38, has been solved. The protein crystallized from ammonium sulfate in the trigonal space group P321, with unit-cell parameters a = b = 71.49, c = 130.30 A, gamma = 120.0 degrees. The structure was solved at 2.0 A by X-ray crystallography using the molecular-replacement method and refined to a final R factor of 18.5% (R(free) = 26.9%). Xyl1 has the overall fold characteristic of family 11 xylanases, with two highly twisted beta-sheets defining a long cleft containing the two catalytic residues Glu87 and Glu177.

Structural approach of human MAO-A using fold recognition (threading) techniques.

The major goal of the present work is to further approach the structure of human monoamine oxidase A (MAO-A). A first partial three-dimensional model of human MAO-A has already been established using secondary structure predictions and fold recognition methods [Wouters and Baudoux, 1998]. In this modeled structure, a segment of the sequence (residues 369-393) located near the covalent linkage to the essential flavin cofactor, and potentially involved in the structure of the active site of the protein, could not be modeled. We here propose a possible fold for that segment, based on threading techniques. The identification of regions of the protein potentially involved in its dimerization was also undertaken by studying hydrophobic areas present at the surface of the structure.

Topology prediction of Brucella abortus Omp2b and Omp2a porins after critical assessment of transmembrane beta strands prediction by several secondary structure prediction methods.

In order to propose a reliable model for Brucella porin topology, several structure prediction methods were evaluated in their ability to predict porin topology. Four porins of known structure were selected as test-cases and their secondary structure delineated. The specificity and sensitivity of 11 methods were separately evaluated. Our critical assessment shows that some secondary structure prediction methods (PHD, Dsc, Sopma) originally designed to predict globular protein structure are useful on porin topology prediction. The overall best prediction is obtained by combining these three "generalist" methods with a transmembrane beta strand prediction technique. This "consensus" method was applied to Brucella porins Omp2b and Omp2a, sharing no sequence homology with any other porin. The predicted topology is a 16-stranded antiparallel beta barrel with Omp2a showing a higher number of negatively charged residue in the exposed loops than Omp2b. Experiments are in progress to validate the proposed topology and the functional hypotheses. The ability of the proposed consensus method to predict topology of complex outer membrane protein is briefly discussed.

Purification, cloning, and three-dimensional structure prediction of Micrococcus luteus FAD-containing tyramine oxidase.

The FAD-containing tyramine oxidase enzyme and gene from the Gram (+) bacterium Micrococcus luteus were isolated, and computer prediction was used to propose a preliminary 3D model of the protein. A 2.8-kb Sau3AI fragment containing the structural gene of tyramine oxidase was cloned from a M. luteus genomic DNA library. The 1332 bp gene encodes a protein of 443 amino acids, with a calculated molecular mass of 49.1 kDa. The enzyme was found to be a homodimer with a molecular weight of 49,000. It oxidizes tyramine, adrenaline, 3-hydroxytyramine, dopamine, and noradrenaline, and was reversibly inhibited by FAD-containing monoamine oxidase A and B specific inhibitors. Sequence comparison show that tyramine oxidase is smaller than other FAD-amine oxidases but that it contains well-conserved amino acid residues reported in all other FAD-amine oxidases. A hypothetical three-dimensional structure of tyramine oxidase has also been proposed based on secondary structure predictions, threading, and comparative modeling.

Lys13 plays a crucial role in the functional adaptation of the thermophilic triose-phosphate isomerase from Bacillus stearothermophilus to high temperatures.

The thermophilic triose-phosphate isomerases (TIMs) of Bacillus stearothermophilus (bTIM) and Thermotoga maritima (tTIM) have been found to possess a His12-Lys13 pair instead of the Asn12-Gly13 pair normally present in mesophilic TIMs. His12 in bTIM was proposed to prevent deamidation at high temperature, while the precise role of Lys13 is unknown. To investigate the role of the His12 and Lys13 pair in the enzyme's thermoadaptation, we reintroduced the "mesophilic residues" Asn and Gly into both thermophilic TIMs. Neither double mutant displayed diminished structural stability, but the bTIM double mutant showed drastically reduced catalytic activity. No similar behavior was observed with the tTIM double mutant, suggesting that the presence of the His12 and Lys13 cannot be systematically correlated to thermoadaptation in TIMs. We determined the crystal structure of the bTIM double mutant complexed with 2-phosphoglycolate to 2.4-A resolution. A molecular dynamics simulation showed that upon substitution of Lys13 to Gly an increase of the flexibility of loop 1 is observed, causing an incorrect orientation of the catalytic Lys10. This suggests that Lys13 in bTIM plays a crucial role in the functional adaptation of this enzyme to high temperature. Analysis of bTIM single mutants supports this assumption.

Structure and function prediction of the Brucella abortus P39 protein by comparative modeling with marginal sequence similarities.

A methodology is proposed to solve a difficult modeling problem related to the recently sequenced P39 protein. This sequence shares no similarity with any known 3D structure, but a fold is proposed by several threading tools. The difficulty in aligning the target sequence on one of the proposed template structures is overcome by combining the results of several available prediction methods and by refining a rational consensus between them. In silico validation of the obtained model and a preliminary cross-check with experimental features allow us to state that this borderline prediction is at least reasonable. This model raises relevant hypotheses on the main structural features of the protein and allows the design of site-directed mutations. Knowing the genetic context of the P39 reading frame, we are now able to suggest a function for the P39 protein: it would act as a periplasmic substrate-binding protein.

Mutagenesis of the fructose-6-phosphate-binding site in the 2-kinase domain of 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase.

Multiple alignment of several isozyme sequences of the bifunctional enzyme 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase revealed conserved residues in the 2-kinase domain. Among these residues, three asparagine residues (Asn76, Asn97 and Asn133; numbering refers to the liver isozyme sequence) and three threonine residues (Thr132, Thr134 and Thr135) are located near the fructose 6-phosphate-binding site in the crystal structure of the bifunctional enzyme. The role of these residues in substrate binding and catalysis in the 6-phosphofructo-2-kinase domain has been studied by mutagenesis to alanine. Since the crystal structure of 6-phosphofructo-2-kinase does not contain fructose 6-phosphate, this substrate was docked into the putative binding site by computer modelling, and its interactions with the protein were predicted. Analysis of the mutagenesis-induced changes in kinetic properties and of the substrate-docking model revealed that all these residues are directly or indirectly involved in fructose-6-phosphate binding. All the mutants displayed an increased Km for fructose 6-phosphate (10-200-fold). We propose that Asn133 stabilises Arg138, which itself makes a direct electrostatic bond with the 6-phosphate group of fructose 6-phosphate, that Asn76 interacts with the C3 hydroxyl group of fructose 6-phosphate, that Thr132 makes a hydrogen bond with the C6 oxygen of this substrate, and that Thr134 interacts with two residues involved in fructose-6-phosphate binding, Thr132 and Tyr199. On the other hand, Asn97 and Thr135 play structural roles, by maintaining the structure of the fructose-6-phosphate-binding pocket.

Comparative analysis of seven multiple protein sequence alignment servers: clues to enhance reliability of predictions.

MOTIVATION: The prediction reliability of seven multiple alignment servers currently available on the Internet (ClustalW, MAP, PIMA, Block Maker, MSA, MEME and Match-Box) has been evaluated in terms of power (sensitivity) and confidence (selectivity). Therefore, the alignments obtained have been respectively compared to refined structural alignments for 20 families of related proteins with low levels of identity. RESULTS: Results clearly show that any powerful method remains reliable when the rate of identity falls. For some methods, power and confidence decrease linearly with the rate of identity, while other methods emphasize reliability at the cost of a lower power. Increasing the number of related sequences included in the alignment may either improve or decrease the quality of the predictions substantially. For some methods, the gain in power or in confidence is quite systematic; for others, the effect of the addition of homologous sequences is highly unpredictable. Extracting the consensus between two different methods may increase the overall confidence of the predictions tremendously. Our conclusions induce users of sequence alignment methods on the Internet to select the most suitable technique according to their requirements in terms of selectivity and sensitivity. AVAILABILITY: The aligned sequences of the 20 alignments of structure can be obtained automatically by sending the message 'send: cabios_tests.txt' by e-mail to 'matchbox@biq.fundp.ac.be'. CONTACT: eric.depiereux@fundp.ac.be

Match-Box_server: a multiple sequence alignment tool placing emphasis on reliability.

MOTIVATION: The Match-Box software comprises protein sequence alignment tools based on strict statistical thresholds of similarity between protein segments. The method circumvents the gap penalty requirement: gaps being the result of the alignment and not a governing parameter of the procedure. The reliable conserved regions outlined by Match-Box are particularly relevant for homology modelling of protein structures, prediction of essential residues for site-directed mutagenesis and oligonucleotide design for cloning homologous genes by polymerase chain reaction (PCR). RESULTS: The method produces reliable results, as assessed by tests performed on protein families of known structures and of low sequence similarity. A reliability score is computed in relation to a threshold of similarity progressively raised to extend the aligned regions to their maximal length, up to the significance limit of matching segments. The score obtained at each position is printed below the sequences and allows a discriminant reading of each aligned region. AVAILABILITY: Sequences may be submitted to a Web server at http://www.fundp.ac.be/sciences/biologie/bms/+ ++matchbox_submit.html or sent by e-mail to matchbox/biq.fundp.ac.be (help available by just mailing help).

Modelling the 2-kinase domain of 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase on adenylate kinase.

Simultaneous multiple alignment of available sequences of the bifunctional enzyme 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase revealed several segments of conserved residues in the 2-kinase domain. The sequence of the kinase domain was also compared with proteins of known three-dimensional structure. No similarity was found between the kinase domain of 6-phosphofructo-2-kinase and 6-phosphofructo-1-kinase. This questions the modelling of the 2-kinase domain on bacterial 6-phosphofructo-1-kinase that has previously been proposed [Bazan, Fletterick and Pilkis (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 9642-9646]. However, sequence similarities were found between the 2-kinase domain and several nucleotide-binding proteins, the most similar being adenylate kinase. A structural model of the 2-kinase domain based on adenylate kinase is proposed. It accommodates all the results of site-directed mutagenesis studies carried out to date on residues in the 2-kinase domain. It also allows residues potentially involved in catalysis and/or substrate binding to be predicted.

ANOLEA: a www server to assess protein structures.

ANOLEA (Atomic Non-Local Environment Assessment) is a www server that performs energy calculations at the atomic level in protein structures. The calculations involve the non-local interactions between all the heavy atoms of the twenty standard amino acids in the molecule. The input of the server is a PDB file containing one or more protein chains. The output is an energy profile, which gives an energy value for each amino acid of the protein. High energy zones (HEZs) in the profile correlate with errors or with potential interacting zones of proteins. The output of the server also displays the structure in three dimensions, pointing out the high energy amino acids in the protein. This option requires the CHIME plug-in, which is freely available on Internet and makes possible, in real time, to rotate, translate and change the point of view and presentation of the molecule in three dimensions. Thus, a fast analysis of a protein structure can be done using a personal computer connected to Internet. The server is available at: http:@www.fundp.ac.be/pub/ANOLEA.html.

Identification of residues potentially involved in the interactions between subunits in yeast alcohol dehydrogenases.

The lack of crystal structure for tetrameric yeast alcohol dehydrogenases (ADHs) has precluded, until now, the identification of the residues involved in subunit contacts. In order to address this question, we have characterized the thermal stability and dissociation propensity of native ADH I and ADH II isozymes as well as of several chimeric (ADH I-ADH II) enzymes. Three groups of substitutions affecting the thermostability have been identified among the 24 substitutions observed between isozymes I and II. The first group contains a Cys277-->Ser substitution, located at the interface between subunits in a three-dimensional model of ADH I, based on the crystallographic structure of the dimeric horse liver ADH. In the second group, the Asp236-->Asn substitution is located in the same interaction zone on the model. The stabilizing effect of this substitution can result from the removal of a charge repulsion between subunits. It is shown that the effect of these two groups of substitutions correlates with changes in dissociation propensities. The third group contains the Met168-->Arg substitution that increases the thermal stability, probably by the formation of an additional salt bridge between subunits through the putative interface. These data suggest that at least part of the subunit contacts observed in horse liver ADH are located at homologous positions in yeast ADHs.

Knowledge-based modeling of the D-lactate dehydrogenase three-dimensional structure.

A three-dimensional structure of the NAD-dependent D-lactate dehydrogenase of Lactobacillus bulgaricus is modeled using the structure of the formate dehydrogenase of Pseudomonas sp. as template. Both sequences share only 22% of identical residues. Regions for knowledge-based modeling are defined from the structurally conserved regions predicted by multiple alignment of a set of related protein sequences with low homology. The model of the D-LDH subunit shows, as for the formate dehydrogenase, an alpha/beta structure, with a catalytic domain and a coenzyme binding domain. It points out the catalytic histidine (His-296) and supports the hypothetical catalytic mechanism. It also suggests that the other residues involved in the active site are Arg-235, possibly involved in the binding of the carboxyl group of the pyruvate, and Phe-299, a candidate for stabilizing the methyl group of the substrate.

Prediction of structurally conserved regions of D-specific hydroxy acid dehydrogenases by multiple alignment with formate dehydrogenase.

We propose a multiple alignment of the sequence of formate dehydrogenase with the D-specific 2-hydroxy acid dehydrogenases family. Structurally conserved regions are predicted for those sequences corresponding to important regions of the catalytic and the coenzyme binding domains defined from the known three-dimensional structure of the formate dehydrogenase, namely the nicotinamide binding site (beta D to beta F) and the beta A-loop-alpha B region containing the typical glycine pattern of the adenosine binding site, the catalytic histidine/aspartic acid pair and an arginine probably involved in the interaction with the carboxyl group of the substrate.

Similarities between alanine dehydrogenase and the N-terminal part of pyridine nucleotide transhydrogenase and their possible implication in the virulence mechanism of Mycobacterium tuberculosis.

Recent developments in simultaneous multiple alignment methods of protein sequences allow prediction of structural similarity in related proteins. Alanine dehydrogenase and the N-terminal sequence of pyridine nucleotide transhydrogenase were compared for their sequences. High similarities of sequences were observed especially in their NAD(H)-binding sites. These similarities suggest that antibodies which recognized the alanine dehydrogenase of Mycobacterium tuberculosis can also be directed against the membrane bound pyridine nucleotide transhydrogenase. If this is the case, the virulent property of this pathogen could be linked to its higher synthesis of NADPH necessary for its anabolism.

MATCH-BOX: a fundamentally new algorithm for the simultaneous alignment of several protein sequences.

Original algorithms for simultaneous alignment of protein sequences are presented, including sequence clustering and within- or between-groups multiple alignment. The way of matching similar regions is fundamentally new. Complete matches are formed by segments more similar than expected by random, according to a given probability limit. Any classic or user-defined score matrix can be used to express the similarity between the residues. The algorithm seeks for complete matches common to all the sequences without performing pairwise alignment and regardless of gap weighting. An automatic screening delineates all the similar regions (boxes) that may be defined for a given maximal shift between the sequences. The shift can be large enough to allow the matching of any region of a sequence with any region of another one. It can also be short and used to refine the alignment around anchor points. The algorithm provides the most likely optimal alignment and a comprehensive list of the alignment dilemma. Duality between automatism and interactivity is provided. Depending on the problem complexity, a final alignment is obtained fully automatically or requires some interactive handling to discriminate alternative pathways.

Simultaneous and multivariate alignment of protein sequences: correspondence between physicochemical profiles and structurally conserved regions (SCR).

A general protein sequence alignment methodology for detecting a priori unknown common structural and functional regions is described. The method proposed in this paper is based on two basic requirements for a meaningful alignment. First, each sequence or segment of a sequence is characterized by a multivariate physicochemical profile. Second, the alignment is performed by considering all the sequences simultaneously, and the algorithm detects those regions that form a set of similar profiles. In order to test the structural meaning of the alignment obtained from the sequences, quantitative comparisons are performed with structurally conserved regions (SCR) determined from the X-ray structures of three serine proteases. Results suggest that the limits of the SCR may be predicted from the similarities between the physicochemical profiles of the sequences. The procedures are not completely automated. The final step requires a visual screening of alternative pathways in order to determine an optimal alignment.

Compartmental analysis of steady-state taurocholate transport kinetics by isolated rat hepatocytes.

We used compartmental modeling to describe taurocholate transport by isolated rat liver cells in suspension. Cells are preincubated in the presence of unlabeled taurocholate. When a steady-state for taurocholate is reached, radiolabeled taurocholate is added to the medium and its exchange kinetics between the medium and the cells are followed over time. Because the studies are performed under steady-state conditions, the kinetics can be described by linear compartmental models. We found a closed two-compartment model sufficient to describe the steady-state transport data. Simulations reveal that if the pools of free and bound intracellular taurocholate exchange rapidly, the cells will behave as a single, kinetically homogeneous compartment and intracellular events will not influence the exchange kinetics of taurocholate between the medium and the cells. The two-compartment model was used to study the concentration dependence of taurocholate transport by isolated cells. Steady-state transport rates and taurocholate concentrations in the medium and the cells were calculated using the model equations. Taurocholate influx, accumulation and efflux processes were studied simultaneously by examining the relationship between appropriate combinations of these variables. Application of this approach to study the inhibition of taurocholate transport by taurochenodeoxycholate is illustrated. In conclusion, this method provides a complementary approach to initial rate studies, which are generally used to investigate bile acid transport by isolated cells.

Saturation of hepatic transport of taurocholate in rats in vivo.

A single intravenous injection of [14C]taurocholate was followed up in blood and bile of rats submitted to steady intravenous infusions of taurocholate (TC) at rates of 0.0, 0.5, 1.0, and 1.5 mumol.min-1.100 g body wt-1 for at least 30 min. The transport rate constants and the amounts of TC in different compartments were estimated by weighted least-squares adjustment of a six-compartment model to the experimental data (3 compartments for TC distribution in blood, 2 compartments for liver, and 1 compartment for sinusoidal blood space). The saturation of the TC excretion rate was reached at 0.8 mumol.min-1.100 g body wt-1. It was characterized by a decrease of both the uptake and excretion rate constants, by an increase of the ratio of the amounts of TC in the two intrahepatic compartments (H'/H), and by an intrahepatic TC concentration of approximately 2 mM. When tauroursodeoxycholate (TUDC) was infused at a rate of 0.5 mumol.min-1.100 g body wt-1 together with TC at a rate of 1.5 mumol.min-1.100 g body wt-1, the TC excretion rate increased to 1.2 mumol.min-1.100 g body wt-1, and the excretion rate constant and H'/H decreased toward control values. These results support the hypothesis that the saturation of the transport of TC is due to TC hepatotoxicity and can be reduced by TUDC. Michaelis-Menten parameters, derived from saturation curves for both uptake and excretion steps, closely matched earlier results, thus confirming the good descriptive capacity of the model.