Analytical shape computation of macromolecules: I. molecular area and volume through alpha shape

The size and shape of macromolecules such as proteins and nucleic acids play an important role in their functions. Prior efforts to quantify these properties have been based on various discretization or tessellation procedures involving analytical or numerical computations. In this article, we present an analytically exact method for computing the metric properties of macromolecules based on the alpha shape theory. This method uses the duality between alpha complex and the weighted Voronoi decomposition of a molecule. We describe the intuitive ideas and concepts behind the alpha shape theory and the algorithm for computing areas and volumes of macromolecules. We apply our method to compute areas and volumes of a number of protein systems. We also discuss several difficulties commonly encountered in molecular shape computations and outline methods to overcome these problems. Proteins 33:1–17, 1998. © 1998 Wiley-Liss, Inc.     

Spatial Point Pattern Analysis by using Voronoi Diagrams and Delaunay Tessellations – A Comparative Study

Given a spatial point pattern, we use various characteristics of its Voronoi diagram and Delaunay tessellation to extract information of the dependence between points. In particular, we use the characteristics to construct statistics for testing complete spatial randomness. It is shown that the minimum angle of a typical Delaunay triangle is sensitive to both regularity and clustering alternatives, whilst the triangle's area or perimeter is more sensitive to clustering than regularity. These statistics are also sensitive to the Baddeley‐Silverman cell process.     

Spatially balanced sampling through the pivotal method

A simple method to select a spatially balanced sample using equal or unequal inclusion probabilities is presented. For populations with spatial trends in the variables of interest, the estimation can be much improved by selecting samples that are well spread over the population. The method can be used for any number of dimensions and can hence also select spatially balanced samples in a space spanned by several auxiliary variables. Analysis and examples indicate that the suggested method achieves a high degree of spatial balance and is therefore efficient for populations with trends.     

Comprehensive mutagenesis of HIV-1 protease: a computational geometry approach

A computational geometry technique based on Delaunay tessellation of protein structure, represented by C(alpha) atoms, is used to study effects of single residue mutations on sequence-structure compatibility in HIV-1 protease. Profiles of residue scores derived from the four-body statistical potential are constructed for all 1881 mutants of the HIV-1 protease monomer and compared with the profile of the wild-type protein. The profiles for an isolated monomer of HIV-1 protease and the identical monomer in a dimeric state with an inhibitor are analyzed to elucidate changes to structural stability. Protease residues shown to undergo the greatest impact are those forming the dimer interface and flap region, as well as those known to be involved in inhibitor binding.     

Water in protein cavities: A procedure to identify internal water and exchange pathways and application to fatty acid-binding protein

A computational approach based on Delaunay triangulation is presented to identify internal water molecules in proteins and to capture pathways of exchange with the bulk. The implemented procedure is computationally efficient and can easily be applied to long molecular dynamics trajectories of protein simulations. In an application to fatty acid-binding protein in apo-form and with bound palmitate, several protein orifices known from crystal structures have been confirmed to be major portals of solvent exchange. Differences between the two forms of the protein are observed and discussed.     

A modelling approach towards epidermal homoeostasis control

In order to grasp the features arising from cellular discreteness and individuality, in large parts of cell tissue modelling agent-based models are favoured. The subclass of off-lattice models allows for a physical motivation of the intercellular interaction rules. We apply an improved version of a previously introduced off-lattice agent-based model to the steady-state flow equilibrium of skin. The dynamics of cells is determined by conservative and drag forces, supplemented with delta-correlated random forces. Cellular adjacency is detected by a weighted Delaunay triangulation. The cell cycle time of keratinocytes is controlled by a diffusible substance provided by the dermis. Its concentration is calculated from a diffusion equation with time-dependent boundary conditions and varying diffusion coefficients. The dynamics of a nutrient is also taken into account by a reaction-diffusion equation. It turns out that the analysed control mechanism suffices to explain several characteristics of epidermal homoeostasis formation. In addition, we examine the question of how in silico melanoma with decreased basal adhesion manage to persist within the steady-state flow equilibrium of the skin. Interestingly, even for melanocyte cell cycle times being substantially shorter than for keratinocytes, tiny stochastic effects can lead to completely different outcomes. The results demonstrate that the understanding of initial states of tumour growth can profit significantly from the application of off-lattice agent-based models in computer simulations.     

Spatial Graphs: Principles and Applications for Habitat Connectivity

ABSTRACT Well-founded methods to assess habitat connectivity are essential to inform land management decisions that include conservation and restoration goals. Indeed, to be able to develop a conservation plan that maintains animal movement through a fragmented landscape, spatial locations of habitat and paths among them need to be represented. Graph-based approaches have been proposed to determine paths among habitats at various scales and dispersal movement distances, and balance data requirements with information content. Conventional graphs, however, do not explicitly maintain geographic reference, reducing communication capacity and utility of other geo-spatial information. We present spatial graphs as a unifying theory for applying graph-based methods in a geographic context. Spatial graphs integrate a geometric reference system that ties patches and paths to specific spatial locations and spatial dimensions. Arguably, the complete graph, with paths between every pair of patches, may be one of the most relevant graphs from an ecosystem perspective, but it poses challenges to compute, process and visualize. We developed Minimum Planar Graphs as a spatial generalization of Delaunay triangulations to provide a reasonable approximation of complete graphs that facilitates visualization and comprehension of the network of connections across landscapes. If, as some authors have suggested, the minimum spanning tree identifies the connectivity “backbone” of a landscape, then the Minimum Planar Graph identifies the connectivity “network”. We applied spatial graphs, and in particular the Minimum Planar Graph, to analyze woodland caribou habitat in Manitoba, Canada to support the establishment of a national park.     

Genetic conservation of Ficus bonijesulapensis R.M. Castro in a dry forest on limestone outcrops

Ficus bonijesulapensis is endemic to a seasonally dry forest on limestone outcrops and it is arranged in disjunct areas in the Cerrado and Caatinga Domains. Species of this genus are considered to be key plants in tropical forests, since they provide resources during periods of scarcity of other resources and, additionally, they contribute in a plant community's restoration. Therefore, the conservation of these species in their natural habitat contributes to the maintenance of the long-term population viability and their genetic diversity. We used nine ISSR primers to analyze the genetic diversity and the genetic spatial patterns of 15 populations of F. bonijesulapensis. We obtained 75 polymorphic bands, the expected heterozygosity (He) was 0.30 and AMOVA showed that most of the genetic diversity was found within populations (77%). The historical gene flow was 1.1 migrants and the Bayesian genetic structure assigned individuals genotypes to eight groups. However, there was not a spatial pattern of genetic variability according to a multivariate correlogram and as confirmed by Mantel's test (r = 0.06, p = 0.68). Eight management units (MU) were proposed and were aiming at the MU's ability to maintain minimally viable populations and a higher genetic variability.     

Mechanical property of a TIM-barrel protein

The mechanical response of a TIM-barrel protein to an applied pressure has been studied. We generated structures under an applied pressure by assuming the volume change to be a linear function of normal mode variables. By Delaunay tessellation, the space occupied by protein atoms is divided uniquely into tetrahedra, whose four vertices correspond to atomic positions. Based on the atoms that define them, the resulting Delaunay tetrahedra are classified as belonging to various secondary structures in the protein. The compressibility of various regions identified with respect to secondary structural elements in this protein is obtained from volume changes of respective regions in two structures with and without an applied pressure. We found that the β barrel region located at the core of the protein is quite soft. The interior of the β barrel, occupied by side chains of β strands, is the softest. The helix, strand, and loop segments themselves are extremely rigid, while the regions existing between these secondary structural elements are soft. These results suggest that the regions between secondary structural elements play an important role in protein dynamics. Another aspect of tetrahedra, referred to as bond distance, is introduced to account for rigidities of the tetrahedra. Bond distance is a measure of separation of the atoms of a tetrahedron in terms of number of bonds along the polypeptide chain or side chains. Tetrahedra with longer bond distances are found to be softer on average. From this behavior, we derive a simple empirical equation, which well describes the compressibilities of various regions. © 1997 Wiley-Liss Inc.