Transimulation - protein biosynthesis web service.

PubMed

Siwiak, Marlena; Zielenkiewicz, Piotr

2013-01-01

Although translation is the key step during gene expression, it remains poorly characterized at the level of individual genes. For this reason, we developed Transimulation - a web service measuring translational activity of genes in three model organisms: Escherichia coli, Saccharomyces cerevisiae and Homo sapiens. The calculations are based on our previous computational model of translation and experimental data sets. Transimulation quantifies mean translation initiation and elongation time (expressed in SI units), and the number of proteins produced per transcript. It also approximates the number of ribosomes that typically occupy a transcript during translation, and simulates their propagation. The simulation of ribosomes' movement is interactive and allows modifying the coding sequence on the fly. It also enables uploading any coding sequence and simulating its translation in one of three model organisms. In such a case, ribosomes propagate according to mean codon elongation times of the host organism, which may prove useful for heterologous expression. Transimulation was used to examine evolutionary conservation of translational parameters of orthologous genes. Transimulation may be accessed at http://nexus.ibb.waw.pl/Transimulation (requires Java version 1.7 or higher). Its manual and source code, distributed under the GPL-2.0 license, is freely available at the website.

PDB to AMPL Conversion

DOE Office of Scientific and Technical Information (OSTI.GOV)

Anna Johnston, SNL 9215

2002-09-01

PDB to AMPL Conversion was written to convert protein data base files to AMPL files. The protein data bases on the internet contain a wealth of information about the structue and makeup of proteins. Each file contains information derived by one or more experiments and contains information on how the experiment waw performed, the amino acid building blocks of each chain, and often the three-dimensional structure of the protein extracted from the experiments. The way a protein folds determines much about its function. Thus, studying the three-dimensional structure of the protein is of great interest. Analysing the contact maps ismore » one way to examine the structure. A contact map is a graph which has a linear back bone of amino acids for nodes (i.e., adjacent amino acids are always connected) and vertices between non-adjacent nodes if they are close enough to be considered in contact. If the graphs are similar then the folds of the protein and their function should also be similar. This software extracts the contact maps from a protein data base file and puts in into AMPL data format. This format is designed for use in AMPL, a programming language for simplifying linear programming formulations.« less

DNAtraffic--a new database for systems biology of DNA dynamics during the cell life.

PubMed

Kuchta, Krzysztof; Barszcz, Daniela; Grzesiuk, Elzbieta; Pomorski, Pawel; Krwawicz, Joanna

2012-01-01

DNAtraffic (http://dnatraffic.ibb.waw.pl/) is dedicated to be a unique comprehensive and richly annotated database of genome dynamics during the cell life. It contains extensive data on the nomenclature, ontology, structure and function of proteins related to the DNA integrity mechanisms such as chromatin remodeling, histone modifications, DNA repair and damage response from eight organisms: Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Escherichia coli and Arabidopsis thaliana. DNAtraffic contains comprehensive information on the diseases related to the assembled human proteins. DNAtraffic is richly annotated in the systemic information on the nomenclature, chemistry and structure of DNA damage and their sources, including environmental agents or commonly used drugs targeting nucleic acids and/or proteins involved in the maintenance of genome stability. One of the DNAtraffic database aim is to create the first platform of the combinatorial complexity of DNA network analysis. Database includes illustrations of pathways, damage, proteins and drugs. Since DNAtraffic is designed to cover a broad spectrum of scientific disciplines, it has to be extensively linked to numerous external data sources. Our database represents the result of the manual annotation work aimed at making the DNAtraffic much more useful for a wide range of systems biology applications.

DNAtraffic—a new database for systems biology of DNA dynamics during the cell life

PubMed Central

Kuchta, Krzysztof; Barszcz, Daniela; Grzesiuk, Elzbieta; Pomorski, Pawel; Krwawicz, Joanna

2012-01-01

DNAtraffic (http://dnatraffic.ibb.waw.pl/) is dedicated to be a unique comprehensive and richly annotated database of genome dynamics during the cell life. It contains extensive data on the nomenclature, ontology, structure and function of proteins related to the DNA integrity mechanisms such as chromatin remodeling, histone modifications, DNA repair and damage response from eight organisms: Homo sapiens, Mus musculus, Drosophila melanogaster, Caenorhabditis elegans, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Escherichia coli and Arabidopsis thaliana. DNAtraffic contains comprehensive information on the diseases related to the assembled human proteins. DNAtraffic is richly annotated in the systemic information on the nomenclature, chemistry and structure of DNA damage and their sources, including environmental agents or commonly used drugs targeting nucleic acids and/or proteins involved in the maintenance of genome stability. One of the DNAtraffic database aim is to create the first platform of the combinatorial complexity of DNA network analysis. Database includes illustrations of pathways, damage, proteins and drugs. Since DNAtraffic is designed to cover a broad spectrum of scientific disciplines, it has to be extensively linked to numerous external data sources. Our database represents the result of the manual annotation work aimed at making the DNAtraffic much more useful for a wide range of systems biology applications. PMID:22110027

Overlapping clusters for distributed computation.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mirrokni, Vahab; Andersen, Reid; Gleich, David F.

2010-11-01

Scalable, distributed algorithms must address communication problems. We investigate overlapping clusters, or vertex partitions that intersect, for graph computations. This setup stores more of the graph than required but then affords the ease of implementation of vertex partitioned algorithms. Our hope is that this technique allows us to reduce communication in a computation on a distributed graph. The motivation above draws on recent work in communication avoiding algorithms. Mohiyuddin et al. (SC09) design a matrix-powers kernel that gives rise to an overlapping partition. Fritzsche et al. (CSC2009) develop an overlapping clustering for a Schwarz method. Both techniques extend an initialmore » partitioning with overlap. Our procedure generates overlap directly. Indeed, Schwarz methods are commonly used to capitalize on overlap. Elsewhere, overlapping communities (Ahn et al, Nature 2009; Mishra et al. WAW2007) are now a popular model of structure in social networks. These have long been studied in statistics (Cole and Wishart, CompJ 1970). We present two types of results: (i) an estimated swapping probability {rho}{infinity}; and (ii) the communication volume of a parallel PageRank solution (link-following {alpha} = 0.85) using an additive Schwarz method. The volume ratio is the amount of extra storage for the overlap (2 means we store the graph twice). Below, as the ratio increases, the swapping probability and PageRank communication volume decreases.« less

Determination of paraldehyde by gas chromatography in whole blood from children.

PubMed

Githiga, Isaiah M; Muchohi, Simon N; Ogutu, Bernhards R; Newton, Charles R J C; Otieno, Godfrey O; Gitau, Evelyn N; Kokwaro, Gilbert O

2004-06-15

A rapid, sensitive and selective gas chromatographic method with flame ionization detection was developed for the determination of paraldehyde in small blood samples taken from children. Whole blood samples (300 microl) collected in a 3 ml Wheaton glass sample vial were spiked with acetone (internal standard: 15 ng) followed by addition of concentrated hydrochloric acid. The mixture was heated in the sealed airtight sample vial in a water bath (96 Celsius; 5 min) to depolymerize paraldehyde to acetaldehyde. A 2 ml aliquot of the headspace was analyzed by gas chromatography with flame ionization detector using a stainless steel column (3 m x 4 mm i.d.) packed with 10% Carbowax 20 M/ 2% KOH on 80/100 Chromosorb WAW. Calibration curves were linear from 1.0-20 microg (r2>0.99). The limit of detection was 1.5 microg/ml, while relative mean recoveries at 2 and 18 microg were 105.6 +/- 8.4 and 101.2 +/- 5.9%, respectively (n = 10 for each level). Intra- and inter-assay relative standard deviations at 2, 10 and 18 microg were <15%. There was no interference from other drugs concurrently used in children with severe malaria, such as anticonvulsants (diazepam, phenytoin, phenobarbitone), antipyretics/analgesics (paracetamol and salicylate), antibiotics (gentamicin, chloramphenicol, benzyl penicillin) and antimalarials (chloroquine, quinine, proguanil, cycloguanil, pyrimethamine and sulfadoxine). The method was successfully applied for pharmacokinetic studies of paraldehyde in children with convulsions associated with severe malaria.

The presence of bromuconazole fungicide pollutant in organic waste anaerobic fermentation

NASA Astrophysics Data System (ADS)

Hariyadi, H. R.

2017-03-01

The presence of bromuconazole fungicide pollutant in organic waste anaerobic fermentation was carried out as well as the influence phenol and benzoate, and biodegradation of bromuconazole. Bromuconazole is a fungicide effective against Ascomycetes, Basidiomycetes and fungi imperfecti in cereals, grapes, top fruits and vegetables. It is also effective against Alternaria and Fusarium sp. The remaining fungicide in leaves might contaminates landfill. One month of organic waste added with bromuconazole was anaerobically incubated in 500 mL bottles at 30°C without shaking in dark room. High-Performance Liquid Chromatography (HPLC) with UV detector and a 100 RP 185μm Lichrosphere column was used to determine bromuconazole concentration. Methane content was determined by Gas Chromatography (GC) method equipped with a flame ionization detector and a metal column packed with 5% neopentyl glycol sebacate and 1% H3PO4 on Chromosorb W-AW (mesh 80-100). After incubation for 225 days, bromuconazole of 200 mg/L inhibited the production of methane (99.5 mM) significantly, but did not inhibit the production of volatile fatty acids. The addition of 100 mg/L phenol or 146 mg/L benzoate increased the production of methane, 143 mM and 135.2 mM, respectively compared with control (121.8 mM). In anaerobic conditions, the presence of toxic pollutants such as fungicide bromuconazole in landfills sites may cause further problems with the accumulation of volatile fatty acids in leachate. Further study to determine the threshold, the presence of bromconazole in low concentration (less than 200 mg/L) on the methane production is recommended.

Ionosphere Waves Service (IWS) - a problem-oriented tool in ionosphere and Space Weather research produced by POPDAT project

NASA Astrophysics Data System (ADS)

Ferencz, Csaba; Lizunov, Georgii; Crespon, François; Price, Ivan; Bankov, Ludmil; Przepiórka, Dorota; Brieß, Klaus; Dudkin, Denis; Girenko, Andrey; Korepanov, Valery; Kuzmych, Andrii; Skorokhod, Tetiana; Marinov, Pencho; Piankova, Olena; Rothkaehl, Hanna; Shtus, Tetyana; Steinbach, Péter; Lichtenberger, János; Sterenharz, Arnold; Vassileva, Any

2014-05-01

In the frame of the FP7 POPDAT project the Ionosphere Waves Service (IWS) has been developed and opened for public access by ionosphere experts. IWS is forming a database, derived from archived ionospheric wave records to assist the ionosphere and Space Weather research, and to answer the following questions: How can the data of earlier ionospheric missions be reprocessed with current algorithms to gain more profitable results? How could the scientific community be provided with a new insight on wave processes that take place in the ionosphere? The answer is a specific and unique data mining service accessing a collection of topical catalogs that characterize a huge number of recorded occurrences of Whistler-like Electromagnetic Wave Phenomena, Atmosphere Gravity Waves, and Traveling Ionosphere Disturbances. IWS online service (http://popdat.cbk.waw.pl) offers end users to query optional set of predefined wave phenomena, their detailed characteristics. These were collected by target specific event detection algorithms in selected satellite records during database buildup phase. Result of performed wave processing thus represents useful information on statistical or comparative investigations of wave types, listed in a detailed catalog of ionospheric wave phenomena. The IWS provides wave event characteristics, extracted by specific software systems from data records of the selected satellite missions. The end-user can access targets by making specific searches and use statistical modules within the service in their field of interest. Therefore the IWS opens a new way in ionosphere and Space Weather research. The scientific applications covered by IWS concern beyond Space Weather also other fields like earthquake precursors, ionosphere climatology, geomagnetic storms, troposphere-ionosphere energy transfer, and trans-ionosphere link perturbations.

Mechanical interaction between volcanic systems in Libya

NASA Astrophysics Data System (ADS)

Elshaafi, Abdelsalam; Gudmundsson, Agust

2018-01-01

The spatial distributions of monogenetic volcanoes, primarily volcanic craters, within the four principal volcanic provinces of Libya are examined and presented on a volcano-density map. Six main volcanic clusters have been identified, referred to as volcanic systems. Remarkably, the Al Haruj (AHVP) and Nuqay (NVP) volcanic provinces have double-peak volcano-density distributions, while the Gharyan (GVP) and As Sawda (SVP) volcanic provinces have single-peak volcano-density distributions. We interpret each volcano-density peak as corresponding to a separate volcanic system, so that there is a total of six systems in these four provinces. There was an overlap in volcanic activity in these provinces with at least three simultaneously active. We propose that each of the 6 volcanic systems was/is supplied with magma from a large sill-like reservoir - similar in lateral dimensions to the systems/clusters themselves. Numerical results show zones of high tensile and shear stresses between the reservoirs that coincide roughly with the main swarms of extension (dykes and volcanic fissures) and shear (faults) fractures in the areas. The most recent volcanic eruptions in Libya fall within the modelled high-stress concentration zones, primarily eruptions in the volcano Waw an Namus and the Holocene Al Mashaqaq lava flow. There are no known eruptions in Libya in historical time, but some or all the volcanic systems may have had one or more arrested historical dyke injections. In particular, part of the recurrent seismic events in the Hun Graben in the northwest Libya may be related to dyke propagation and arrest. If some of the inferred magma reservoirs are still fluid, as is likely, they pose earthquake and volcanic hazards to parts of Libya, particularly to the city of Gharyan and Zallah town, as well as to many oil-field operations.

Warsaw Catalogue of cometary orbits: 119 near-parabolic comets

NASA Astrophysics Data System (ADS)

Królikowska, Małgorzata

2014-07-01

Context. The dynamical evolution of near-parabolic comets strongly depends on the starting values of the orbital elements derived from the positional observations. In addition, when drawing conclusions about the origin of these objects, it is crucial to control the uncertainties of orbital elements at each stage of the dynamical evolution. Aims: I apply a completely homogeneous approach to determine the cometary orbits and their uncertainties. The resulting catalogue is suitable for the investigation of the origin and future of near-parabolic comets. Methods: First, osculating orbits were determined on the basis of positional data. Second, the dynamical calculations were performed backwards and forwards up to 250 au from the Sun to derive original and future barycentric orbits for each comet. In the present investigation of dynamical evolution, the numerical calculations for a given object start from the swarm of virtual comets constructed using the previously determined osculating (nominal) orbit. In this way, the uncertainties of orbital elements were derived at the end of numerical calculations. Results: Homogeneous sets of orbital elements for osculating, original and future orbits are given. The catalogue of 119 cometary orbits constitutes about 70 per cent of all the first class so-called Oort spike comets discovered during the period 1801-2010 and about 90 per cent of those discovered in 1951-2010, for which observations were completed at the end of 2013. Non-gravitational (NG) orbits are derived for 45 comets, including asymmetric NG solution for six of them. Additionally, the new method for cometary orbit-quality assessment is applied for all these objects. The catalogue is available at http://ssdp.cbk.waw.pl/LPCs and also at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr (ftp://130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/567/A126

Ionospheric research for space weather service support

NASA Astrophysics Data System (ADS)

Stanislawska, Iwona; Gulyaeva, Tamara; Dziak-Jankowska, Beata

2016-07-01

Knowledge of the behavior of the ionosphere is very important for space weather services. A wide variety of ground based and satellite existing and future systems (communications, radar, surveillance, intelligence gathering, satellite operation, etc) is affected by the ionosphere. There are the needs for reliable and efficient support for such systems against natural hazard and minimalization of the risk failure. The joint research Project on the 'Ionospheric Weather' of IZMIRAN and SRC PAS is aimed to provide on-line the ionospheric parameters characterizing the space weather in the ionosphere. It is devoted to science, techniques and to more application oriented areas of ionospheric investigation in order to support space weather services. The studies based on data mining philosophy increasing the knowledge of ionospheric physical properties, modelling capabilities and gain applications of various procedures in ionospheric monitoring and forecasting were concerned. In the framework of the joint Project the novel techniques for data analysis, the original system of the ionospheric disturbance indices and their implementation for the ionosphere and the ionospheric radio wave propagation are developed since 1997. Data of ionosonde measurements and results of their forecasting for the ionospheric observatories network, the regional maps and global ionospheric maps of total electron content from the navigational satellite system (GNSS) observations, the global maps of the F2 layer peak parameters (foF2, hmF2) and W-index of the ionospheric variability are provided at the web pages of SRC PAS and IZMIRAN. The data processing systems include analysis and forecast of geomagnetic indices ap and kp and new eta index applied for the ionosphere forecasting. For the first time in the world the new products of the W-index maps analysis are provided in Catalogues of the ionospheric storms and sub-storms and their association with the global geomagnetic Dst storms is investigated. The products of the Project web sites at http://www.cbk.waw.pl/rwc and http://www.izmiran.ru/services/iweather are widely used in scientific investigations and numerous applications by the telecommunication and navigation operators and users whose number at the web sites is growing substantially from month to month.

Final Report - "CO2 Sequestration in Cell Biomass of Chlorobium Thiosulfatophilum"

DOE Office of Scientific and Technical Information (OSTI.GOV)

James L. Gaddy, PhD; Ching-Whan Ko, PhD

2009-05-04

World carbon dioxide emissions from the combustion of fossil fuels have increased at a rate of about 3 percent per year during the last 40 years to over 24 billion tons today. While a number of methods have been proposed and are under study for dealing with the carbon dioxide problem, all have advantages as well as disadvantages which limit their application. The anaerobic bacterium Chlorobium thiosulfatophilum uses hydrogen sulfide and carbon dioxide to produce elemental sulfur and cell biomass. The overall objective of this project is to develop a commercial process for the biological sequestration of carbon dioxide andmore » simultaneous conversion of hydrogen sulfide to elemental sulfur. The Phase I study successfully demonstrated the technical feasibility of utilizing this bacterium for carbon dioxide sequestration and hydrogen sulfide conversion to elemental sulfur by utilizing the bacterium in continuous reactor studies. Phase II studies involved an advanced research and development to develop the engineering and scale-up parameters for commercialization of the technology. Tasks include culture isolation and optimization studies, further continuous reactor studies, light delivery systems, high pressure studies, process scale-up, a market analysis and economic projections. A number of anaerobic and aerobic microorgansims, both non-photosynthetic and photosynthetic, were examined to find those with the fastest rates for detailed study to continuous culture experiments. C. thiosulfatophilum was selected for study to anaerobically produce sulfur and Thiomicrospira crunogena waws selected for study to produce sulfate non-photosynthetically. Optimal conditions for growth, H2S and CO2 comparison, supplying light and separating sulfur were defined. The design and economic projections show that light supply for photosynthetic reactions is far too expensive, even when solar systems are considered. However, the aerobic non-photosynthetic reaction to produce sulfate with T. crunogena produces a reasonable return when treating a sour gas stream of 120 million SCFD containing 2.5 percent H2S. In this case, the primary source of revenue is from desulfurization of the gas stream. While the technology has significant application in sequestering carbon dioxide in cell biomass or single cell proten (SCP), perhaps the most immediate application is in desulfurizing LGNG or other gas streams. This biological approach is a viable economical alternative to existing hydrogen sulfide removal technology, and is not sensitive to the presence of hydrocarbons which act as catalyst poisons.« less

Restitic or not? Insights from trace element content and crystal - Structure of spinels in African mantle xenoliths

NASA Astrophysics Data System (ADS)

Lenaz, Davide; Musco, Maria Elena; Petrelli, Maurizio; Caldeira, Rita; De Min, Angelo; Marzoli, Andrea; Mata, Joao; Perugini, Diego; Princivalle, Francesco; Boumehdi, Moulay Ahmed; Bensaid, Idris Ali Ahmadi; Youbi, Nasrrddine

2017-05-01

The lithospheric architecture of Africa consists of several Archean cratons and smaller cratonic fragments, stitched together and flanked by polycyclic fold belts. Here we investigate the structure and chemistry of spinels from lithospheric mantle xenoliths from distinct tectonic settings, i.e. from the Saharan metacraton in Libya (Waw-En-Namus) which could show archaic chemical features, Cameroon (Barombi Koto and Nyos Lakes) where the Sub Continental Lithospheric Mantle was modified during the Pan-African event and fluxed by asthenospheric melts of the Tertiary Cameroon Volcanic Line and Morocco (Tafraoute, Bou-Ibalrhatene maars) in the Middle Atlas where different metasomatic events have been recorded. From a structural point of view it is to notice that the Libyan spinels can be divided into two groups having different oxygen positional parameter (u > 0.2632 and u < 0.2627, respectively), while those from Cameroon are in between those values as the Moroccan ones already studied by other authors. The intracrystalline closure temperature (Tc) of the here studied spinels is different among the different samples with one Libyan group (LB I) showing Tc in the range 490-640 °C and the other 680-950 °C (LB II). Cameroon and Morocco spinels show a Tc in the range 630-760 °C. About 150 different spinels have been studied for their trace element content and it can be seen that many of them are related to Cr content, while Zn and Co are not and clearly distinguish the occurrences. Differences in the trace element chemistry, in the structural parameters and in the intracrystalline closure temperatures suggest that a different history should be considered for Cameroon, Morocco and LB I and LB II spinels. Even if it was not considered for this purpose, we tentatively used the Fe2 +/Fe3 + vs. TiO2 diagram that discriminate between peridotitic and the so-called "magmatic" spinels, i.e. spinel crystallized from melts. LB I and LB II spinels plot in the peridotitic field while Cameroon and Morocco spinels fall in the magmatic one. Consequently, the xenoliths sampled from a probably juvenile SCLM at the edge of the most important lithospheric roots (i.e. Cameroon and Morocco) apparently have spinels possibly fractionated in situ from percolating melts and do not represent a real spinel-peridotite facies. On the contrary mantle xenoliths from Libya exhibit spinels with peridotitic features compatible with a slow ascent of a mantle diapir (plume).

New catalogue of single-apparition comets discovered in the years 1901-1950. Part I

NASA Astrophysics Data System (ADS)

Królikowska, M.; Sitarski, G.; Pittich, E.; Szutowicz, S.; Ziołkowski, K.; Rickman, H.; Gabryszewski, R.; Rickman, B.

2014-07-01

A new catalogue of cometary orbits derived using a completely homogeneous method of data treatment, accurate methods of numerical integration, and modern model of the Solar System is presented. We constructed a sample of near-parabolic comets from the first half of the twentieth century with original reciprocals of semimajor axes less than 0.000130 au^{-1} in the Marsden and Williams Catalogue of Cometary Orbits (2008, hereafter MW08), i.e., comets of original semimajor axes larger than 7700 au. We found 38 such comets in MW08, where 32 have first-quality orbits (class 1A or 1B) and the remaining 6 have second-quality orbits (2A or 2B). We presented satisfactory non-gravitational (hereafter NG) models for thirteen of the investigated comets. The four main features, distinguishing this catalogue of orbits of single- apparition comets discovered in the early twentieth century from other catalogues of orbits of similarly old objects, are the following. 1. Old cometary positional observations require a very careful analysis. For the purpose of this new catalogue, great emphasis has been placed in collecting sets of observations as complete as possible for the investigated comets. Moreover, for many observations, comet-minus-star-type measurements were also available. This type of data was particularly valuable as the most original measurements of comet positions and has allowed us to recalculate new positions of comets using the PPM star catalogue. 2. Old cometary observations were prepared by observers usually as apparent positions in Right Ascension and Declination or as reduced positions for the epoch of the beginning of the year of a given observation. This was a huge advantage of these data, because this allows us to uniformly take into account all necessary corrections associated with the data reduction to the standard epoch. 3. The osculating orbits of single-apparition comets discovered more than sixty years ago have been formerly determined with very different numerical methods and assumptions on the model of the Solar System, including the number of planets taken into account. This new catalogue changes this situation. We offer a new catalogue of cometary orbits derived using completely homogeneous methods of data treatment, accurate methods of numerical integration, and a modern model of the Solar System. 4. The osculating, original, and future sets of orbits are presented for each catalogue comet. In the case of a comet with detectable NG effects, we give both types of orbit: purely gravitational and non- gravitational. We concluded, however, that all thirteen NG orbital solutions given in the catalogue better represent the actual motions of the investigated comets. Surprisingly, the NG effects were detectable in data for five comets of second-quality-class orbits. Among these five are three comets with hyperbolic original, barycentric GR orbits. This publication will be accompanied by an online catalogue available at ssdp.cbk.waw.pl/LPCs, providing entries to orbital elements of considered comets as well as to full swarms of original and future virtual comets that formed the basis for the further analysis of dynamical evolution.