Comparative analysis of the Borrelia garinii genome.

Three members of the genus Borrelia (B.burgdorferi, B.garinii, B.afzelii) cause tick-borne borreliosis. Depending on the Borrelia species involved, the borreliosis differs in its clinical symptoms. Comparative genomics opens up a way to elucidate the underlying differences in Borrelia species. We analysed a low redundancy whole-genome shotgun (WGS) assembly of a B.garinii strain isolated from a patient with neuroborreliosis in comparison to the B.burgdorferi genome. This analysis reveals that most of the chromosome is conserved (92.7% identity on DNA as well as on amino acid level) in the two species, and no chromosomal rearrangement or larger insertions/deletions could be observed. Furthermore, two collinear plasmids (lp54 and cp26) seem to belong to the basic genome inventory of Borrelia species. These three collinear parts of the Borrelia genome encode 861 genes, which are orthologous in the two species examined. The majority of the genetic information of the other plasmids of B.burgdorferii is also present in B.garinii although orthology is not easy to define due to a high redundancy of the plasmid fraction. Yet, we did not find counterparts of the B.burgdorferi plasmids lp36 and lp38 or their respective gene repertoire in the B.garinii genome. Thus, phenotypic differences between the two species could be attributable to the presence or absence of these two plasmids as well as to the potentially positively selected genes.

More hydrogen bonds for the (structural) biologist.

Why does a given protein structure form and why is this structure stable? These fundamental biochemical questions remain fascinating and challenging problems because the physical bases of the forces that govern protein structure, stability and folding are still not well understood. Now, a general concept of hydrogen bonding in proteins is emerging. This concept involves not only N-H and O-H donor groups, but also C-H, and not only N and O as acceptor groups, but also pi-systems. We postulate that the incorporation of the entirety of these interactions leads to a more complete description of the problem, and that this could provide new perspectives and possibly new answers.

Toxicology of simple and complex mixtures.

Humans are exposed to mixtures of chemicals, rather than to individual chemicals. From a public health point of view, it is most relevant to answer the question of whether or not the components in a mixture interact in a way that results in an increase in their overall effect compared with the sum of the effects of the individual components. In this article, options for the hazard identification and risk assessment of simple and complex chemical mixtures will be discussed. In addition, key research needed to continue the development of hazard characterization of chemical mixtures will be described. Clearly, more collaboration among toxicologists, model developers and pharmacologists will be necessary.

Quantum-chemical analysis of C-H...O and C-H...N interactions in RNA base pairs--H-bond versus anti-H-bond pattern.

Geometries and interaction energies of unusual UU and AA base pairs with one standard hydrogen bond (H-bond) and additional C-H...O or C-H...N contacts have been determined by quantum-chemical methods taking into account electron correlation. Whereas the C-H bond length in the UU C-H...O contact increases upon complex formation (H-bond pattern), the C-H bond of the AA C-H....N interaction is shortened (anti-H-bond pattern). The same properties are found for model complexes between U or A and formaldehyde that have intermolecular C-H...acceptor contacts but no standard H-bonds. Both the C-H...acceptor H-bond and anti-H-bond interactions are attractive. A possible influence of the donor CH group charge distribution on the interaction pattern is discussed.

C-H...pi-interactions in proteins.

A non-redundant set of 1154 protein structures from the Protein Data Bank was examined with respect to close interactions between C-H-donor and pi-acceptor groups. A total of 31,087 interactions were found to satisfy our selection criteria. Their geometric parameters suggest that these interactions can be classified as weak hydrogen bonds.A set of 12 interaction classes were defined based on the division of the donors into three groups and the acceptors into four groups. These classes were examined separately, and the respective interactions described in detail in each class. Most prominent were interactions between aliphatic C-H donors and aromatic pi-acceptors and interactions between aromatic C-H donors and aromatic pi-acceptors. About three-quarters of the Trp-rings, half of all Phe and Tyr-rings and a quarter of all His-rings were found to be involved as acceptors in C-H...pi-interactions. On the donor side, a preference for aromatic C-H groups was observed, but also for the aliphatic side-chains of the long, extended amino acid residues Lys, Arg and Met, and the Pro ring. The average distance between the C-donor and the center-of-mass of the pi-acceptor was observed to be significantly longer in the 174 protein structures determined at >2.5 A resolution. Also, the distribution is significantly wider. This resolution dependence suggests that the force fields commonly used for the refinement of protein structures may not be adequate. C-H...pi-interactions involving aromatic groups either as donor or as acceptor groups are found mostly in the interior of the protein. The more hydrophilic the participating groups are, the closer to the surface are the interactions located. About 40 % of all C-H...pi-interactions occur between amino acid residue side-chains that are separated by nine or less residues in sequence. Dependent on the interaction class, different preferences for secondary structure, residue type and side-chain conformation were observed. It is likely that the C-H...pi-interactions contribute significantly to the overall stability of a protein.

Molecular dynamics simulation reveals conformational switching of water-mediated uracil-cytosine base-pairs in an RNA duplex.

A 4 ns molecular dynamics simulation of an RNA duplex (r-GGACUUCGGUCC)(2 )in solution with Na+ and Cl- as counterions was performed. The X-ray structure of this duplex includes two water-mediated uracil-cytosine pairs. In contrast to the other base-pairs in the duplex the water-mediated pairs switch between different conformations. One conformation corresponds to the geometry of the water-mediated UC pairs in the duplex X-ray structure with water acting both as hydrogen-bond donor and acceptor. Another conformation is close to that of a water-mediated UC base-pair found in the X-ray structure of the 23 S rRNA sarcin/ricin domain. In this case the oxygen of the water molecule is linked to two-base donor sites. For a very short time also a direct UC base-pair and a further conformation that is similar to the one found in the RNA duplex structure but exhibits an increased H3(U)...N3(C) distance is observed. Water molecules with unusually long residence times are involved in the water-mediated conformations. These results indicate that the dynamic behaviour of the water-mediated UC base-pairs differs from that of the duplex Watson-Crick and non-canonical guanine-uracil pairs with two or three direct hydrogen bonds. The conformational variability and increased flexibility has to be taken into account when considering these base-pairs as RNA building blocks and as recognition motifs.

Beyond nucleic acid base pairs: from triads to heptads.

Hydrogen-bonded base pairs are an important determinant of nucleic acid structure and function. However, other interactions such as base-base stacking, base-backbone, and backbone-backbone interactions as well as effects exerted by the solvent and by metal or NH(4)(+) ions also have to be taken into account. In addition, hydrogen-bonded base complexes involving more than two bases can occur. With the rapidly increasing number and structural diversity of nucleic acid structures known at atomic detail higher-order hydrogen-bonded base complexes, base polyads, have attracted much interest. This review provides an overview on the occurrence of base polyads in nucleic acid structures and describes computational studies on these nucleic acid building blocks.

A molecular dynamics simulation study of coaxial stacking in RNA.

We report on unrestrained molecular dynamics simulations of an RNA tetramer binding to a tetra-nucleotide overhang at the 5'-end of an RNA hairpin (nicked structure) and of the corresponding continuous hairpin with Na+ as counterions. The simulations lead to stable structures and in this way a structural model for the coaxially stacked RNA hairpin is generated. The stacking interface in the coaxially stacked nicked hairpin structure is characterized by a reduced twist and shift and a slightly increased propeller twist as compared to the continuous system. This leads to an increased overlap between C22 and G23 in the stacking interface of the nicked structure. In the simulations the continuous RNA hairpin has an almost straight helical axis. On the other hand, the corresponding axis for the nicked structure exhibits a marked kink of 39 degrees. The stacking interface exhibits no increased flexibility as compared to the corresponding base pair step in the continuous structure.

The IMB Jena Image Library of biological macromolecules.

The IMB Jena Image Library of Biological Macro-molecules (http://www. imb-jena.de/IMAGE.html ) is aimed at a better dissemination of information on three-dimensional biopolymer structures with an emphasis on visualization and analysis. It provides access to all structure entries deposited at the Protein Data Bank (PDB) and Nucleic Acid Database (NDB). By combining automatic and manual processing it is possible to keep pace with the rapidly growing number of known biopolymer structures and to provide, for selected entries, information not available from automatic procedures. Each entry page contains basic information on the structure, various visualization and analysis tools as well as links to other databases. The visualization techniques adopted include static mono/stereo raster or vector graphics representations, virtual reality modeling (VRML), RasMol/Chime scripts and Java applets. A helix and bending analysis tool provides consistent information on about 750 DNA and RNA duplex structures. Access to metal-containing PDB entries is possible via the Periodic Table of Elements. Finally, general information on amino acids, cis -peptide bonds, structural elements in proteins, base pairs, nucleic acid model conformations and experimental methods for biopolymer structure determination is provided.

A molecular dynamics simulation of the flavin mononucleotide-RNA aptamer complex.

We report on an unrestrained molecular dynamics simulation of the flavin mononucleotide (FMN)-RNA aptamer. The simulated average structure maintains both cross-strand and intermolecular FMN-RNA nuclear Overhauser effects from the nmr experiments and has all qualitative features of the nmr structure including the G10-U12-A25 base triple and the A13-G24, A8-G28, and G9-G27 mismatches. However, the relative orientation of the hairpin loop to the remaining part of the molecule differs from the nmr structure. The simulation predicts that the flexible phosphoglycerol part of FMN moves toward G27 and forms hydrogen bonds. There are structurally long-lived water molecules in the FMN binding pocket forming hydrogen bonds within FMN and between FMN and RNA. In addition, long-lived water is found bridging primarily RNA backbone atoms. A general feature of the environment of long-lived "structural" water is at least two and in most cases three or four potential acceptor atoms. The 2'-OH group of RNA usually acts as an acceptor in interactions with the solvent. There are almost no intrastrand O2'H(n) vertical ...O4'(n + 1) hydrogen bonds within the RNA backbone. In the standard case the preferred orientation of the 2'-OH hydrogen atoms is approximately toward O3' of the same nucleotide. However, a relatively large number of conformations with the backbone torsional angle gamma in the trans orientation is found. A survey of all experimental RNA x-ray structures shows that this backbone conformation occurs but is less frequent than found in the simulation. Experimental nmr RNA aptamer structures have a higher fraction of this conformation as compared to the x-ray structures. The backbone conformation of nucleotide n + 1 with the torsional angle gamma in the trans orientation leads to a relatively short distance between 2'-OH(n) and O5'(n + 1), enabling hydrogen-bond formation. In this case the preferred orientation of the 2'-OH hydrogen atom is approximately toward O5'(n + 1). We find two relatively short and dynamically stable types of backbone-backbone next-neighbor contacts, namely C2'(H)(n) vertical ...O4'(n + 1) and C5'(H)(n + 1) vertical ...O2'(n). These interactions may affect both backbone rigidity and thermodynamic stability of RNA helical structures.

CombiTool--a new computer program for analyzing combination experiments with biologically active agents.

CombiTool is a new computer program for the analysis of combination effects of biologically active agents. It performs model calculations and an analysis of experimental combination effects for two or three agents according to both the Bliss independence and the Loewe additivity criteria. Zero interaction response surfaces are calculated from single-agent dose-response relations and compared to experimental combination data. The calculation of response surfaces for Loewe additivity is based on a new approach which combines the implicit definition equation in terms of doses alone with single-agent dose-response relations. The simultaneous analysis of experimental data according to both Loewe additivity and Bliss independence within one program can hopefully contribute to a better understanding of the meaning and limits of the two criteria. CombiTool has a built-in graphics facility which allows the direct visualization of the response surfaces or the corresponding contour plots and the experimental data.

Molecular modelling and footprinting studies of DNA minor groove binders: bisquaternary ammonium heterocyclic compounds.

We report new quantitative footprinting data which reveal differences in binding constants of bisquaternary ammonium heterocyclic compounds (BQA) with AT-rich DNA sites depending on the ligand structure and on the size and sequence of the DNA binding site. In an attempt to understand the dependence of binding affinity on the ligand structure we have performed quantum-chemical AM1 calculations on the BQA compounds and on subunits to explore the conformational space and to calculate the electronic and structural features of individual ligand conformations. Due to the properties of the rotatable backbone bonds, there is a large number of possible conformations with almost equal energy. We present a new method for the calculation of the radius of curvature of molecular structures. Assuming that strong binders should have a shape complementary to the DNA minor groove, this measure is used to select the optimum conformations for DNA-drug binding. The approach yields the correct ligand conformation for SN6999, for which an X-ray DNA-drug structure is known. The curvature of the optimum conformations of all ligands is compared with the experimental binding constants. A correlation is found between curvature and binding constant provided other structural factors do not vary. Therefore, we conclude that within structurally similar BQA compounds the extent of curvature is the relevant quantity which modulates the binding affinity.

Parallel dose-response curves in combination experiments.

A possible experimental design for combination experiments is to compare the dose-response curve of a single agent with the corresponding curve of the same agent using either a fixed amount of a second one or a fixed dose ratio. No interaction is then often defined by a parallel shift of these curves. We have performed a systematic study for various types of dose-response relations both for the dose-additivity (Loewe additivity) and for the independence (Bliss independence) criteria for defining zero interaction. Parallelism between dose-response curves of a single agent and those of the same agent in the presence of a fixed amount of another one is found for the Loewe-additivity criterion for linear dose-response relations. For nonlinear relations, one has to differentiate between effect parallelism (parallel shift on the effect scale) and dose parallelism (parallel shift on the dose scale). In the case of Loewe additivity, zero-interaction dose parallelism is found for power, Weibull, median-effect and logistic dose-response relations, given that special parameter relationships are fulfilled. The mechanistic model of competitive interaction exhibits dose parallelism but not effect parallelism for Loewe additivity. Bliss independence and Loewe additivity lead to identical results for exponential dose-response curves. This is the only case for which dose parallelism was found for Bliss independence. Parallelism between single-agent dose-response relations and Loewe additivity mixture relations is found for examples with a fixed dose-ratio design. However, this is again not a general property of the design adopted but holds only if special conditions are fulfilled. The comparison of combination dose-response curves with single-agent relations has to be performed taking into account both potency and shape parameters. The results of this analysis lead to the conclusion that parallelism between zero interaction combination and single-agent dose-response relations is found only for special cases and cannot be used as a general criterion for defining zero-interaction in combined-action assessment even if the correct potency shift is taken into account.

Computer analysis of phytochrome sequences and reevaluation of the phytochrome secondary structure by Fourier transform infrared spectroscopy.

A repertoire of various methods of computer sequence analysis was applied to phytochromes in order to gain new insights into their structure and function. A statistical analysis of 23 complete phytochrome sequences revealed regions of non-random amino acid composition, which are supposed to be of particular structural or functional importance. All phytochromes other than phyD and phyE from Arabidopsis have at least one such region at the N-terminus between residues 2 and 35. A sequence similarity search of current databases indicated striking homologies between all phytochromes and a hypothetical 84.2-kDa protein from the cyanobacterium Synechocystis. Furthermore, scanning the phytochrome sequences for the occurrence of patterns defined in the PROSITE database detected the signature of the WD repeats of the beta-transducin family within the functionally important 623-779 region (sequence numbering of phyA from Avena) in a number of phytochromes. A multiple sequence alignment performed with 23 complete phytochrome sequences is made available via the IMB Jena World-Wide Web server (http://www.imb-jena.de/PHYTO.html). It can be used as a working tool for future theoretical and experimental studies. Based on the multiple alignment striking sequence differences between phytochromes A and B were detected directly at the N-terminal end, where all phytochromes B have an additional stretch of 15-42 amino acids. There is also a variety of positions with totally conserved but different amino acids in phytochromes A and B. Most of these changes are found in the sequence segment 150-200. It is, therefore, suggested that this region might be of importance in determining the photosensory specificity of the two phytochromes. The secondary structure prediction based on the multiple alignment resulted in a small but significant beta-sheet content. This finding is confirmed by a reevaluation of the secondary structure using FTIR spectroscopy.

Quantum-chemical ab initio study on the adenine-difluorotoluene complex--a mimic for the adenine-thymine base pair.

Recent experiments have shown that difluorotoluene (F), a nonpolar isostere for thymine (T), codes efficiently and specifically for adenine (A) in DNA replication. F has almost the same shape as thymine but it is unable to form conventional hydrogen bonds with adenine. Therefore, it has been claimed that not hydrogen bonding but shape complementary may be important for the selection of the correct bases by DNA-replicating enzymes. In order to gain deeper insight into structure, charge distribution and energetics of the A-F and A-T base pairs we have performed quantum-chemical ab initio and density functional calculations at the HF, MP2 and B3LYP levels. The interaction energy of the A-F complex amounts to -3.8 kcal/mol (MP2) and is thus substantially smaller than typical ab initio interaction energies for Watson-Crick or non-canonical base pairs. The A-T and A-F complexes are planar and their overall geometries are similar (root-mean-square deviation: 0.4 A). The calculated donor acceptor atom distances in A-T are in good agreement with the experimental mean values obtained from an analysis of 21 high resolution DNA structures. However, A-F shows a base pair opening as compared to A-T. Even though the interaction energy in the A-F base pair is small, the distances for the N6-H...F and N1...H-C3 contacts are still below the sum of the van-der-Waals radii, which means that the interaction is not governed by van-der-Waals forces alone. If the experimental findings can be confirmed, then our results indicate that DNA polymerase is able to retain high fidelity with base pairs of much smaller interaction energies than found for the conventional Watson-Crick and non-canonical base pairs.

Zero interaction response surfaces for combined-action assessment.

Combined-action assessment requires first the calculation of combination effects expected for the case of no interaction. In a second step the experimental combination effects have to be compared to the expected ones in order to classify the combination as synergistic or antagonistic. Adopting the response surface methodology, a new procedure for the calculation of zero interaction response surfaces which can be applied both to various criteria for defining zero interaction and to various types of dose-response relations is proposed. The approach requires only information from single-agent dose-response relations, but provides nevertheless a comprehensive overview over the complete dose range under study. It can replace the tedious isobolographic analysis if the dose-additivity criterion is adopted.

HBexplore--a new tool for identifying and analysing hydrogen bonding patterns in biological macromolecules.

The program HBexplore is a new tool for identifying and analysing H-bonding patterns in biological macromolecules. It selects all potential H-bonds according to geometrical criteria. The H-bond table can then be subjected to further automatic or interactive analysis tools. These tools include the calculation of mean values and distributions of geometrical H-bond parameters for parts of a single structure, for complete single structures and for structure sets, the classification of each H-bond according to the participation of backbone, side chain or base, ligand and water parts of nucleic acids or proteins, identification of Watson-Crick nucleotide pairs and of H-bonded pairs of equal nucleotides, the calculation of the mean number of H-bonds per residue, and of the fraction of potential donor and acceptor atoms involved in H-bonds. HBexplore further generates automatically a H-bond residue interaction table. This table lists for all residues of the structure the other residues, ligands or water molecules directly connected via a H-bond. By means of a binary tree search algorithm, this table is then converted into a H-bond cluster table. Clusters are understood here as an uninterrupted network of H-bonded residues. For nucleic acids, secondary structures and tertiary interactions are automatically derived from the H-bonding pattern. HBexplore is applied to two example RNA structures: a pseudoknot and a hairpin. It provides a comprehensive listing of individual H-bonds and statistical information for larger structure sets. In addition, it can identify interesting new H-bond motifs. One example is a pentanucleotide base-base H-bond interaction motif in the RNA pseudoknot. HBexplore is intended to contribute both to the elucidation of general principles of the architecture of biological macromolecules, and to the prediction and refinement of single structures.

Image library of biological macromolecules.

An Image Library of Biological Macromolecules is described, which contains image and text files related to structures of biological macromolecules. Currently, the Library has approximately 3000 image files of approximately 300 structures of biological macromolecules whose coordinates are available in the Protein Data Bank and in the Nucleic Acid Database. The entries include all RNA structures, approximately 70 DNA structures, 150 proteins and a few carbohydrates. The Library contains further images of amino acids, of standard and modified nucleotides and of nucleic acid model structures. Each entry consists of an annotation file with bibliographic and sequence information and possibly comments, of a color-coded distance plot and of structure images. Almost all of the images are available both in a mono and in a stereo representation. Standard procedures for generating these images were strictly avoided. Therefore, mixed rendering, coloring and labeling techniques were used extensively. Since May 1995 the Library has a growing division of images in the new Virtual Reality Modeling Language (VRML) format. The Image Library of Biological Macromolecules can be accessed via the World-Wide Web (http://www.imb-jena.de/IMAGE.html). There is a large number of structures determined by experimental and/or modeling techniques which are not intended to be included into the Protein Data Bank or Nucleic Acid Database for some reason. The Image Library could be a repository of these structures and of images of these and other structures of biological macromolecules including structures which are not known at atomic detail. Authors who are willing to make available images or coordinates to the scientific community via the Image Library of Biological Macromolecules are requested to contact the author.

Zero interaction response surfaces, interaction functions and difference response surfaces for combinations of biologically active agents.

In the field of combination experiments there is wide-spread confusion over definitions, terminology and methods for the evaluation of interaction between biologically active agents. According to our view the widely used isobole approach is the method of choice. In this contribution it is shown how the combination of the classical isobole approach with response surface modeling and computer graphics leads to powerful new methods for the assessment of interaction of biologically active agents. In particular, zero interaction response surfaces, difference response surfaces and interaction functions are proposed. Zero interaction response surfaces represent surfaces which display zero interaction in the whole dose range. Difference response surfaces display the difference between an actual response surface and the corresponding zero interaction response surface. Interaction functions are a generalization of the index of interaction, which describe the dose dependence of this quantity.