Ligand-Info,searching for similar small compounds using index profiles

MOTIVATION: The Ligand-Info system is based on the assumption that small molecules with similar structure have similar functional (binding) properties. The developed system enables a fast and sensitive index based search for similar compounds in large databases. Index profiles, constructed by averaging indexes of related molecules are used to increase the specificity of the search. The utilization of index profiles helps to focus on frequent, common features of a family of compounds. RESULTS: A Java-based tool for clustering and scanning of small molecules has been created. The tool can interactively cluster sets of molecules and create index profiles on the user side and automatically download similar molecules from a databases of 250 000 compounds. The results of the application of index profiles demonstrate that the profile based search strategy can increase the quality of the selection process. AVAILABILITY: The system is available at http://Ligand.Info. The application requires the Java Runtime Environment 1.4, which can be automatically installed during the first use on desktop systems, which support it. A standalone version of the program is available from the authors upon request.     

A digital system for generating dynamic sinusoidal gas concentration signals

A computer-controlled gas-mixing system is presented. It is capable of mixing four gases, the concentration of three of which will follow a path to be determined by the user. For our purposes the output O2 fraction is maintained constant and the levels of Ar and N2O vary sinusoidally and independently, with periods between 0.25 and 30 min. A fourth gas, N2 is necessary to make the sum of the individual fractions 100%. The system uses banks of between one and four solenoid valves each linked via a sonic choke to a common mixing chamber. A regime of pulse frequency modulation is employed. All calculations and timing of valve switching are performed by a dedicated microcomputer built for the purpose. The device has been used to provide respiratory gas forcing functions for a program of research in respiratory monitoring.     

PEANUT: Computer graphics program to represent atomic displacement parameters

PEANUT is an easy to use computer graphics program for the visualization and real-time manipulation of the atomic displacement parameters of small molecules. A flexible, dynamic data structure allows the user to compute complicated, nonspherical atomic surfaces and to handle the point group symmetry of the molecules automatically. Pictures with hidden-line removal may be plotted in publication quality on appropriate output devices.     

Software for analysis and manipulation of genetic linkage data.

We present eight computer programs written in the C programming language that are designed to analyze genotypic data and to support existing software used to construct genetic linkage maps. Although each program has a unique purpose, they all share the common goals of affording a greater understanding of genetic linkage data and of automating tasks to make computers more effective tools for map building. The PIC/HET and FAMINFO programs automate calculation of relevant quantities such as heterozygosity, PIC, allele frequencies, and informativeness of markers and pedigrees. PREINPUT simplifies data submissions to the Centre d'Etude du Polymorphisme Humain (CEPH) data base by creating a file with genotype assignments that CEPH's INPUT program would otherwise require to be input manually. INHERIT is a program written specifically for mapping the X chromosome: by assigning a dummy allele to males, in the nonpseudoautosomal region, it eliminates falsely perceived noninheritances in the data set. The remaining four programs complement the previously published genetic linkage mapping software CRI-MAP and LINKAGE. TWOTABLE produces a more readable format for the output of CRI-MAP two-point calculations; UNMERGE is the converse to CRI-MAP's merge option; and GENLINK and LINKGEN automatically convert between the genotypic data file formats required by these packages. All eight applications read input from the same types of data files that are used by CRI-MAP and LINKAGE. Their use has simplified the management of data, has increased knowledge of the content of information in pedigrees, and has reduced the amount of time needed to construct genetic linkage maps of chromosomes.     

Automatic localization of curved wires used in brachytherapy. The COUREP program.

The localization of sources used for brachytherapy has become more difficult since radium sources are replaced by flexible wires such as Iridium 192 wires. We describe the COUREP program which automatically reconstructs the real 3-D arrangement of curved wires from two orthogonal radiographs of the implant. The data used by the program consist mainly of the projections of the wires on the radiographs expressed as the co-ordinates of points selected on these projections individually on each radiograph. The output consists of the 3-D co-ordinates of points located on the wires with optionally a graphic output including the plotting of a perspective view of the implant.     

Computer assisted microdensitometric analysis of footprinting autoradiographic DATA.

A system comprised of a linear scanning microdensitometer interfaced to a personal computer has been developed to facilitate analysis of ligand-DNA footprinting autoradiograms. The system, which can be used to record density and sequence information from autoradiographic films, enables the user to relate the area under an autoradiographic band to the concentration of radiolabeled molecules present in the electrophoresis gel. This report describes the computer program which performs the calculations and discusses the ability of the system to accurately determine oligonucleotide concentration, as a function of band separation, photographic response, and the computational algorithm used to calculate band areas.     

MixtureTree Annotator: A Program for Automatic Colorization and Visual Annotation of MixtureTree

The MixtureTree Annotator, written in JAVA, allows the user to automatically color any phylogenetic tree in Newick format generated from any phylogeny reconstruction program and output the Nexus file. By providing the ability to automatically color the tree by sequence name, the MixtureTree Annotator provides a unique advantage over any other programs which perform a similar function. In addition, the MixtureTree Annotator is the only package that can efficiently annotate the output produced by MixtureTree with mutation information and coalescent time information. In order to visualize the resulting output file, a modified version of FigTree is used. Certain popular methods, which lack good built-in visualization tools, for example, MEGA, Mesquite, PHY-FI, TreeView, treeGraph and Geneious, may give results with human errors due to either manually adding colors to each node or with other limitations, for example only using color based on a number, such as branch length, or by taxonomy. In addition to allowing the user to automatically color any given Newick tree by sequence name, the MixtureTree Annotator is the only method that allows the user to automatically annotate the resulting tree created by the MixtureTree program. The MixtureTree Annotator is fast and easy-to-use, while still allowing the user full control over the coloring and annotating process.     

A vector-based method for drawing RNA secondary structure.

MOTIVATION: To produce a polygonal display of RNA secondary structure with minimal overlap and distortion of structural elements, with minimal search for positioning them, and with minimal user intervention. RESULTS: A new algorithm for automatically drawing RNA secondary structure has been developed. The algorithm represents the direction and space for a structural element using vector and vector space. Two heuristics are used. The first heuristic is concerned with ordering structural elements to be positioned and the second with positioning them in space. The algorithm and a graphical user interface have been implemented in a working program called VizQFolder on IBM PC compatibles. Experimental results demonstrate that VizQFolder is capable of automatically generating nearly overlap-free polygonal displays for long RNA molecules. The only distortion performed to avoid overlap is the rotation of helices, leading to efficient generation of a polygonal display without sacrificing its readability. VizQFolder is not coupled to a specific prediction program of RNA secondary structure, and thus can be used for visualizing secondary structure models obtained by any means. AVAILABILITY: The executable code of VizQFolder is available at http://automation.inha.ac.kr/khan. It can also be obtained from the authors upon request.     

CO: A chemical ontology for identification of functional groups and semantic comparison of small molecules

A novel chemical ontology based on chemical functional groups automatically, objectively assigned by a computer program, was developed to categorize small molecules. It has been applied to PubChem and the small molecule interaction database to demonstrate its utility as a basic pharmacophore search system. Molecules can be compared using a semantic similarity score based on functional group assignments rather than 3D shape, which succeeds in identifying small molecules known to bind a common binding site. This ontology will serve as a powerful tool for searching chemical databases and identifying key functional groups responsible for biological activities.     

Calculation and display of electrostatic potentials

A general methodology is developed for incorporating accurate electrostatic information from ab initio molecular orbital calculations into molecular mechanics calculations. Examples are given of the method applied to simple aromatic organic molecules. A program has been developed for displaying the results of the ab initio calculations on a Silicon Graphics workstation. The technique developed here provides an alternative method for including electrostatic interactions in molecular mechanics calculations and is compared with other methods for determining atomic charges.     

Pseudocontact shifts as constraints for energy minimization and molecular dynamics calculations on solution structures of paramagnetic metalloproteins

The pseudocontact shifts of NMR signals, which arise from the magnetic susceptibility anisotropy of paramagnetic molecules, have been used as structural constraints under the form of a pseudopotential in the SANDER module of the AMBER 4.1 molecular dynamics software package. With this procedure, restrained energy minimization (REM) and restrained molecular dynamics (RMD) calculations can be performed on structural models by using pseudocontact shifts. The structure of the cyanide adduct of the Met80Ala mutant of the yeast iso-1-cytochrome c has been used for successfully testing the calculations. For this protein, a family of structures is available, which was obtained by using NOE and pseudocontact shifts as constraints in a distance geometry program. The structures obtained by REM and RMD calculations with the inclusion of pseudocontact shifts are analyzed. Proteins 29:68–76, 1997. © 1997 Wiley-Liss, Inc.     

Use of the Computer Automated Structure Evaluation program in determining quantitative structure-activity relationships within hallucinogenic phenylalkylamines

The Computer Automated Structure Evaluation (CASE) program has been successfully used to generate automatically and identify molecular fragments relevant to the hallucinogenic activity expressed by some phenylalkylamines. Utilizing these major fragments, Quantitative Structure-Activity Relationship (QSAR) calculations were carried out to obtain an equation which was used for predictions of potencies. Correlations of these major activating/inactivating fragments with the biological activity of the compounds, as well as predictive capabilities of the CASE program, are discussed.     

Bacteriorhodospin: a trans-membrane pump containing alpha-helix.

The probable arrangement of the bacteriorhodopsin molecules in the purple membrane of Halobacterium halobium is in clusters of three, with a 3-fold axis at the centre of each cluster; the axis is at right angles to the plane of the membrane. The proposed arrangement and the results of model calculations together indicate that each protein molecule spans the entire thickness of the membrane. An earlier proposal for the structure had the protein molecules in two layers, and it was symmetric in projection onto the profile-axis. This model is now rejected since it would be difficult to account for the recently discovered function of pumping protons. There remains a discrepancy in that the calculated number of protein molecules in the unit-cell is 3.4 compared to the three expected.The X-ray diffraction patterns from dispersions of the lipids extracted from the red and purple membranes of H. halobium are described.Model calculations are reported, which are based on the bilayer profile calculated for the extracted lipids and on two simple profiles for the protein. The calculations favour a structure for the purple membrane having the lipid molecules in two layers, as in a bilayer, although there may be more of the lipid on one side of the membrane than on the other. Assuming bilayer structure, the diffraction nearest the centre of the oriented pattern suggests that the lipid molecules may be located mainly in a few discrete regions, roughly 20 Å across, between the protein molecules. An uninterrupted monolayer of the lipid on one surface of a sheet of the protein molecules gives poor agreement with the observed profile-diffraction.The X-ray diffraction pattern from the oriented membranes suggested α-helix in the bacteriorhodopsin, and this has been confirmed by recording a 1.5 Å-reflection oriented on the profile-axis. There appear to be at least two segments of α-helix, which are somewhat inclined to one another, and the two may be packed together. Prominent diffraction on the in-plane axis near 10 Å is consistent with the segments lying more or less perpendicular to the plane of the membrane.     

Explicit solvent models in protein pKa calculations.

Continuum methods for calculation of protein electrostatics treat buried and ordered water molecules by one of two approximations; either the dielectric constant of regions containing ordered water molecules is equal to the bulk solvent dielectric constant, or it is equal to the protein dielectric constant though no fixed atoms are used to represent water molecules. A method for calculating the titration behavior of individual residues in proteins has been tested on models of hen egg white lysozyme containing various numbers of explicit water molecules. Water molecules were included based on hydrogen bonding, solvent accessibility, and/or proximity to titrating groups in the protein. Inclusion of water molecules significantly alters the calculated titration behavior of individual titrating sites, shifting calculated pKa values by up to 0.5 pH unit. Our results suggest that approximately one water molecule within hydrogen-bonding distance of each charged group should be included in protein electrostatics calculations.     

Automated quality assurance and patient dosimetry in diagnostic radiology

It has been generally considered that improved methods of quality assurance would reduce the population dose from diagnostic radiology. This paper describes the development of a computerized method of automatically monitoring tube and generator parameters to perform on-line quality assurance, whilst undertaking various patient dosimetry measurements and calculations for each exposure. The method involves interfacing a microcomputer to a microprocessor controlled X-ray generator. Details of the various interfacing methods and modifications to the X-ray unit are given. The instrument enables quality assurance to be performed for every exposure by comparing tube and generator parameters against nominal settings. The software automatically warns the operator of any deviations from accepted limiting values. When a patient is examined, details of the examination and projection are entered into a database. The exposure area product and field size are monitored for each exposure. This data, together with information on tube potential and examination/projection is used to deduce patient entrance skin dose and energy imparted. Doses to individual organs are estimated using normalized organ dose data and a knowledge of tube potential and field size.     

AMIDE: a free software tool for multimodality medical image analysis

Amide's a Medical Image Data Examiner (AMIDE) has been developed as a user-friendly, open-source software tool for displaying and analyzing multimodality volumetric medical images. Central to the package's abilities to simultaneously display multiple data sets (e.g., PET, CT, MRI) and regions of interest is the on-demand data reslicing implemented within the program. Data sets can be freely shifted, rotated, viewed, and analyzed with the program automatically handling interpolation as needed from the original data. Validation has been performed by comparing the output of AMIDE with that of several existing software packages. AMIDE runs on UNIX, Macintosh OS X, and Microsoft Windows platforms, and it is freely available with source code under the terms of the GNU General Public License.     

A visual data flow environment for macromolecular crystallographic computing

This article describes the integration of programs from the widely used CCP4 macromolecular crystallography package into a modern data flow visualization environment (application visualization system [AVS]), which provides a simple graphical user interface, a visual programming paradigm, and a variety of 1-, 2-, and 3-D data visualization tools for the display of graphical information and the results of crystallographic calculations, such as electron density and Patterson maps. The CCP4 suite comprises a number of separate Fortran 77 programs, which communicate via common file formats. Each program is encapsulated into an AVS macro module, and may be linked to others in a data flow network, reflecting the nature of many crystallo-graphic calculations. Named pipes are used to pass input parameters from a graphical user interface to the program module, and also to intercept line printer output, which can be filtered to extract graphical information and significant numerical parameters. These may be passed to downstream modules, permitting calculations to be automated if no user interaction is required, or giving the user the opportunity to make selections in an interactive manner.     

A three-dimensional method for calculating currents induced in bodies by extremely low-frequency electric fields

A user-friendly, numerical program has been developed to permit the calculation of induced currents in modeled bodies of human and infrahuman subjects. The program is based on a charge-simulation method (CSM), and it takes into account the three-dimensional (3-D) character of the extremely-low-frequency (ELF) electric field and of the models to be exposed. The principle of the method is to simulate a 3-D object, for example, an animal model, by a combination of several parts (blocks) having simple geometric forms such as a sphere, a cylinder, or a cone. This approach permits easy preparation of input data on the dimensions of the blocks and their positions in a 3-D arrangement. Other input data, such as the coordinates of the contour points and the imaginary values of charges inside objects, which are necessary in the calculations by CSM, are produced automatically by selecting an appropriate "level" for each block, according to its importance. To simulate parts having irregular shapes, special blocks may be added. In one series of experiments, induced currents were calculated for a baboon model in various postures: standing upright, positioned on four legs, and sitting on the floor. Calculated currents, the total induced current in particular, agreed very well with experimental values. Local currents in parts of the baboon models were more variable, ranging from 5% to 17% of measured values in the case of induced currents in the head. Some problems with this method, such as the effect of the dimensions of blocks or the choice of block levels, are discussed.     

Search for β2 Adrenergic Receptor Ligands by Virtual Screening via Grid Computing and Investigation of Binding Modes by Docking and Molecular Dynamics Simulations

We designed a program called MolGridCal that can be used to screen small molecule database in grid computing on basis of JPPF grid environment. Based on MolGridCal program, we proposed an integrated strategy for virtual screening and binding mode investigation by combining molecular docking, molecular dynamics (MD) simulations and free energy calculations. To test the effectiveness of MolGridCal, we screened potential ligands for β₂ adrenergic receptor (β₂AR) from a database containing 50,000 small molecules. MolGridCal can not only send tasks to the grid server automatically, but also can distribute tasks using the screensaver function. As for the results of virtual screening, the known agonist BI-167107 of β₂AR is ranked among the top 2% of the screened candidates, indicating MolGridCal program can give reasonable results. To further study the binding mode and refine the results of MolGridCal, more accurate docking and scoring methods are used to estimate the binding affinity for the top three molecules (agonist BI-167107, neutral antagonist alprenolol and inverse agonist ICI 118,551). The results indicate agonist BI-167107 has the best binding affinity. MD simulation and free energy calculation are employed to investigate the dynamic interaction mechanism between the ligands and β₂AR. The results show that the agonist BI-167107 also has the lowest binding free energy. This study can provide a new way to perform virtual screening effectively through integrating molecular docking based on grid computing, MD simulations and free energy calculations. The source codes of MolGridCal are freely available at http://molgridcal.codeplex.com.