Voronoi binding site models: Calculation of binding modes and influence of drug binding data accuracy

A new and accurate method for calculating the geometrically allowed modes of binding of a ligand molecule to a Voronoi site model is reported. It is shown that the feasibility of the binding of a group of atoms to a Voronoi site reduces to a simple set of linear and quadratic inequalities and quadratic equalities which can be solved by minimization of a simple function. Newton's numerical method of solution coupled to a line search proved to be successful. Moreover, we have developed efficient molecular and site data bases to discard quickly infeasible binding modes without time-consuming numerical calculation. The method is tested with a data set consisting of the binding constants for a series of biphenyls binding to prealbumin. After determination of the conformation space of the molecules and proposal of a Voronoi site geometry, the geometrically feasible modes are calculated and the energy interaction parameters determined to fit the observed binding energies to the site within experimental error ranges. We actually allowed these ranges to vary in order to study the influence of their broadness on the site geometry and found that as they increase, one can first model the receptor as a three-region site then as a single region site, but never as a two-region site.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/38280/1/540100509_ftp.pd

Temporal and Spatial Impact of Human Cadaver Decomposition on Soil Bacterial and Arthropod Community Structure and Function

As vertebrate carrion decomposes, there is a release of nutrient-rich fluids into theunderlying soil, which can impact associated biological community structure andfunction. How these changes alter soil biogeochemical cycles is relatively unknown and may prove useful in the identification of carrion decomposition islands that have long lasting, focal ecological effects. This study investigated the spatial (0, 1, and 5 m) and temporal (3–732 days) dynamics of human cadaver decomposition on soil bacterial and arthropod community structure and microbial function. We observed strong evidence of a predictable response to cadaver decomposition that varies over space for soil bacterial and arthropod community structure, carbon (C) mineralization and microbial substrate utilization patterns. In the presence of a cadaver (i.e., 0 m samples), the relative abundance of Bacteroidetes and Firmicutes was greater, while the relative abundance of Acidobacteria, Chloroflexi, Gemmatimonadetes, and Verrucomicrobia was lower when compared to samples at 1 and 5 m. Micro-arthropods were more abundant (15 to 17-fold) in soils collected at 0 m compared to either 1 or 5 m, but overall, micro-arthropod community composition was unrelated to either bacterial community composition or function. Bacterial community structure and microbial function also exhibited temporal relationships, whereas arthropod community structure did not. Cumulative precipitation was more effective in predicting temporal variations in bacterial abundance and microbial activity than accumulated degree days. In the presence of the cadaver (i.e., 0 m samples), the relative abundance of Actinobacteria increased significantly with cumulative precipitation. Furthermore, soil bacterial communities and C mineralization were sensitive to the introduction of human cadavers as they diverged from baseline levels and did not recover completely in approximately 2 years. These data are valuable for understanding ecosystem function surrounding carrion decomposition islands and can be applicable to environmental bio-monitoring and forensic sciences

The measurement of molecular diversity by receptor site interaction simulation

The assembly of large compound libraries for the purpose of screening against various receptor targets to identify chemical leads for drug discovery programs has created a need for methods to measure the molecular diversity of such libraries. The method described here, for which we propose the acronym RESIS (for Receptor Site Interaction Simulation), relates directly to this use. A database is built of three-dimensional representations of the compounds in the library and a set of three-point three-dimensional theoretical receptor sites is generated based on putative hydrophobic and polar interactions. A series of flexible, three-dimensional searches is then performed over the database, using each of the theoretical sites as the basis for one such search. The resulting pattern of hits across the grid of theoretical receptor sites provides a measure of the molecular diversity of the compound library. This can be conveniently displayed as a density map which provides a readily comprehensible visual impression of the library diversity characteristics. A library of 7500 drug compounds derived from the CIPSLINEPC databases was characterized with respect to molecular diversity using the RESIS method. Some specific uses for the information obtained from application of the method are discussed. A comparison was made of the results from the RESIS method with those from a recently published two-dimensional approach for assessing molecular diversity using sets of compounds from the Maybridge database (MAY).Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/42964/1/10822_2004_Article_165958.pd

Statistical Properties of Contact Maps

A contact map is a simple representation of the structure of proteins and other chain-like macromolecules. This representation is quite amenable to numerical studies of folding. We show that the number of contact maps corresponding to the possible configurations of a polypeptide chain of N amino acids, represented by (N-1)-step self avoiding walks on a lattice, grows exponentially with N for all dimensions D>1. We carry out exact enumerations in D=2 on the square and triangular lattices for walks of up to 20 steps and investigate various statistical properties of contact maps corresponding to such walks. We also study the exact statistics of contact maps generated by walks on a ladder.Comment: Latex file, 15 pages, 12 eps figures. To appear on Phys. Rev.

Molecular dynamics simulation of polymer helix formation using rigid-link methods

Molecular dynamics simulations are used to study structure formation in simple model polymer chains that are subject to excluded volume and torsional interactions. The changing conformations exhibited by chains of different lengths under gradual cooling are followed until each reaches a state from which no further change is possible. The interactions are chosen so that the true ground state is a helix, and a high proportion of simulation runs succeed in reaching this state; the fraction that manage to form defect-free helices is a function of both chain length and cooling rate. In order to demonstrate behavior analogous to the formation of protein tertiary structure, additional attractive interactions are introduced into the model, leading to the appearance of aligned, antiparallel helix pairs. The simulations employ a computational approach that deals directly with the internal coordinates in a recursive manner; this representation is able to maintain constant bond lengths and angles without the necessity of treating them as an algebraic constraint problem supplementary to the equations of motion.Comment: 15 pages, 14 figure

Semiclassical Mechanics of the Wigner 6j-Symbol

The semiclassical mechanics of the Wigner 6j-symbol is examined from the standpoint of WKB theory for multidimensional, integrable systems, to explore the geometrical issues surrounding the Ponzano-Regge formula. The relations among the methods of Roberts and others for deriving the Ponzano-Regge formula are discussed, and a new approach, based on the recoupling of four angular momenta, is presented. A generalization of the Yutsis-type of spin network is developed for this purpose. Special attention is devoted to symplectic reduction, the reduced phase space of the 6j-symbol (the 2-sphere of Kapovich and Millson), and the reduction of Poisson bracket expressions for semiclassical amplitudes. General principles for the semiclassical study of arbitrary spin networks are laid down; some of these were used in our recent derivation of the asymptotic formula for the Wigner 9j-symbol.Comment: 64 pages, 50 figure

Computation of protein geometry and its applications: Packing and function prediction

This chapter discusses geometric models of biomolecules and geometric constructs, including the union of ball model, the weigthed Voronoi diagram, the weighted Delaunay triangulation, and the alpha shapes. These geometric constructs enable fast and analytical computaton of shapes of biomoleculres (including features such as voids and pockets) and metric properties (such as area and volume). The algorithms of Delaunay triangulation, computation of voids and pockets, as well volume/area computation are also described. In addition, applications in packing analysis of protein structures and protein function prediction are also discussed.Comment: 32 pages, 9 figure