A systematic method for studying the spatial distribution of water molecules around nucleic acid bases.

A new method to analyze the distribution of water molecules around the bases in DNA is presented. This method relies on the notion of a "hydrated building block," which represents the joint observed hydration around all bases of a particular type, in structures of a particular conformation type. The hydrated building blocks were constructed using atomic coordinates from 40 structures contained in the Nucleic Acid Database. Pseudoelectron densities were calculated for water molecules in each hydrated building block using standard crystallographic procedures. The electron densities were fitted to obtain "average building blocks," which represent bases with waters only at average or probable positions. Both types of building blocks were used to construct models of hydrated DNA oligomers. The essential features of the solvent structure around d(CGCGAATTCGCG)2 in the B form and d(CGCGCG)2 in the Z form were reproduced.     

Time-accurate,parallel, multi-zone,multi-block solver to study the human cardio-vascular system

A parallel, time-accurate flow solver is devised to study the human cardio-vascular system. The solver is capable of dealing with moving boundaries and moving grids. It is designed to handle complex, three-dimensional vascular systems. The computational domain is divided into multiple block subdomains. At each cross section the plane is divided into twelve sub-zones to allow flexibility for handling complex geometries and, if needed, appropriate parallel data partitioning. The unsteady, three-dimensional, incompressible Navier-Stokes equations are solved numerically. A second-order in time and third-order upwind finite volume method for solving time-accurate incompressible flows based on pseudo-compressibility and dual time-stepping technique is used. For parallel execution, the flow domain is partitioned. Communication between the subdomains of the flow on Riken's VPP/700E supercomputer is implemented using MPI message-passing library. A series of numerical simulations of biologically relevant flows is used to validate this code.     

Hierarchical protein folding pathways: a computational study of protein fragments

We have previously presented a building block folding model. The model postulates that protein folding is a hierarchical top-down process. The basic unit from which a fold is constructed, referred to as a hydrophobic folding unit, is the outcome of combinatorial assembly of a set of "building blocks." Results obtained by the computational cutting procedure yield fragments that are in agreement with those obtained experimentally by limited proteolysis. Here we show that as expected, proteins from the same family give very similar building blocks. However, different proteins can also give building blocks that are similar in structure. In such cases the building blocks differ in sequence, stability, contacts with other building blocks, and in their 3D locations in the protein structure. This result, which we have repeatedly observed in many cases, leads us to conclude that while a building block is influenced by its environment, nevertheless, it can be viewed as a stand-alone unit. For small-sized building blocks existing in multiple conformations, interactions with sister building blocks in the protein will increase the population time of the native conformer. With this conclusion in hand, it is possible to develop an algorithm that predicts the building block assignment of a protein sequence whose structure is unknown. Toward this goal, we have created sequentially nonredundant databases of building block sequences. A protein sequence can be aligned against these, in order to be matched to a set of potential building blocks.     

Effects of strychnine on the sodium conductance of the frog node of Ranvier

Strychnine blocks sodium conductance in the frog node of Ranvier. This block was studied by reducing and slowing sodium inactivation with scorpion venom. The block is voltage and time dependent. The more positive the axoplasm the greater the block and the faster the approach to equilibrium. Some evidence is presented suggesting that only open channels can be blocked. The block is reduced by raising external sodium or lithium but not impermeant cations. A quaternary derivative of strychnine was synthesized and found to have the same action only when applied intracellularly. We conclude that strychnine blocks sodium channels by a mechanism analogous to that by which it blocks potassium channels. The potassium channel block had previously been found to be identical to that by tetraethylammonium ion derivatives. In addition, strychnine resembles procaine and its derivatives in both its structure and the mechanism of sodium channel block.     

Comparative analysis of human and bovine papillomaviruses

A method is presented for the analysis and comparison of nucleic acid and protein sequences utilizing all identity blocks (the term "identity block" refers to a set of consecutive matches between two sequences) above a prescribed length. Moreover, such identity blocks are determined for various groupings of amino acids according to chemical, functional, charge, and hydrophobic classifications. Alignment maps based on these classifications and containing all statistically significant identity blocks between two or more sequences are constructed. New theoretical results for determining the expected length of the longest identity block between sequences are also presented and are used, along with permutation procedures, to ascertain the significance of sequence identity blocks. As an example of the type of information that can be obtained, comparison has been made of the complete DNA sequences and the E1, E2, L1, and L2 genes of human and bovine papillomaviruses based on the classification schemes described above.     

Model for the feedback control system of bacterial growth. II. Growth in continuous culture

A mathematical model is developed that describes substrate limited bacterial growth in a continuous culture and that is based upon the conceptual framework elaborated in a previous paper for describing the feedback control system of cell growth [S. Bleecken, (1988). J. theor. Biol. 133, 37.] Central to the theory are the ideas that the limiting substrate is converted into low molecular weight building blocks of macromolecular synthesis which again are converted into biomass (RNA and protein) and that the rates of RNA and protein synthesis are controlled by the intracellular concentration of building blocks. It is shown that a continuous culture can be simulated by two interconnected feedback control systems the actuating signals of which are limiting substrate concentration and the intracellular concentration of building blocks, respectively. Three types of steady-states are found to appear in a continuous culture, besides the well-known stable steady-state of the whole culture there exist two batchlike steady-states of the biotic part of the culture which are metastable. The model is used to analyse the steady-states and their stability properties as well as the dynamic responses of biomass, RNA, protein, building block and substrate concentrations to changes in environmental conditions. Especially the inoculation of a continuous culture and the effects of step changes in dilution rate, inlet substrate concentration and growth temperature are studied in detail. Relations between the growth behaviour of a single cell and that of a continuous culture are derived. The RNA to protein ratio is introduced as a rough measure of the physiological state of cells and it is shown that a cell reacts to environmental changes with a simple pattern of basic responses in growth rate and physiological state. There are reasons to assume that the model presented is the minimal version of a structured model of bacterial growth and represents an optimum compromise between biological relevance and mathematical practicability.     

A balance control model of quiet upright stance based on an optimal control strategy

Models of balance control can aid in understanding the mechanisms by which humans maintain balance. A balance control model of quiet upright stance based on an optimal control strategy is presented here. In this model, the human body was represented by a simple single-segment inverted pendulum during upright stance, and the neural controller was assumed to be an optimal controller that generates ankle control torques according to a certain performance criterion. This performance criterion was defined by several physical quantities relevant to sway. In order to accurately simulate existing experimental data, an optimization procedure was used to specify the set of model parameters to minimize the scalar error between experimental and simulated sway measures. Thirty-two independent simulations were performed for both younger and older adults. The model's capabilities, in terms of reflecting sway behaviors and identifying aging effects, were then analyzed based on the simulation results. The model was able to accurately predict center-of-pressure-based sway measures, and identify potential changes in balance control mechanisms caused by aging. Correlations between sway measures and model parameters are also discussed.     

A statistical framework for haplotype block inference

The existence of haplotype blocks transmitted from parents to offspring has been suggested recently. This has created an interest in the inference of the block structure and length. The motivation is that haplotype blocks that are characterized well will make it relatively easier to quickly map all the genes carrying human diseases. To study the inference of haplotype block systematically, we propose a statistical framework. In this framework, the optimal haplotype block partitioning is formulated as the problem of statistical model selection; missing data can be handled in a standard statistical way; population strata can be implemented; block structure inference/hypothesis testing can be performed; prior knowledge, if present, can be incorporated to perform a Bayesian inference. The algorithm is linear in the number of loci, instead of NP-hard for many such algorithms. We illustrate the applications of our method to both simulated and real data sets.     

A general modeling framework for describing spatially structured population dynamics

Variation in movement across time and space fundamentally shapes the abundance and distribution of populations. Although a variety of approaches model structured population dynamics, they are limited to specific types of spatially structured populations and lack a unifying framework. Here, we propose a unified network‐based framework sufficiently novel in its flexibility to capture a wide variety of spatiotemporal processes including metapopulations and a range of migratory patterns. It can accommodate different kinds of age structures, forms of population growth, dispersal, nomadism and migration, and alternative life‐history strategies. Our objective was to link three general elements common to all spatially structured populations (space, time and movement) under a single mathematical framework. To do this, we adopt a network modeling approach. The spatial structure of a population is represented by a weighted and directed network. Each node and each edge has a set of attributes which vary through time. The dynamics of our network‐based population is modeled with discrete time steps. Using both theoretical and real‐world examples, we show how common elements recur across species with disparate movement strategies and how they can be combined under a unified mathematical framework. We illustrate how metapopulations, various migratory patterns, and nomadism can be represented with this modeling approach. We also apply our network‐based framework to four organisms spanning a wide range of life histories, movement patterns, and carrying capacities. General computer code to implement our framework is provided, which can be applied to almost any spatially structured population. This framework contributes to our theoretical understanding of population dynamics and has practical management applications, including understanding the impact of perturbations on population size, distribution, and movement patterns. By working within a common framework, there is less chance that comparative analyses are colored by model details rather than general principles.     

Block duplications of a zeta-zeta-alpha-theta gene set in the rabbit alpha-like globin gene cluster

In order to understand the coordinate regulation between the alpha-like and beta-like globins during the developmental switches in hemoglobin synthesis, we have studied the rabbit alpha-like globin gene family. A cluster of six linked genes arranged 5'-zeta 1-alpha 1-theta 1-zeta 2-zeta 3-theta 2-3' has been isolated as a set of overlapping clones from a library of rabbit genomic DNA. Blot-hybridization analysis of genomic DNA not only confirms this linkage arrangement but also reveals the presence of additional zeta and theta genes. We propose that this gene cluster was generated by a block duplication of a set of alpha-like genes; the proposed duplication unit is zeta-zeta-alpha-theta. Further duplications of a zeta-zeta-theta set are also proposed to have occurred. As expected for a duplicated locus, the rabbit alpha-like gene cluster contains long blocks of internal homology. The Z homology block is about 7.2 kilobase pairs long and contains the zeta genes; the T homology block is about 4.7 kilobase pairs long and contains a theta gene. Surprisingly, both Z and T homology blocks are flanked by a common junction sequence (J) which contains a region very similar to the 3'-untranslated sequence of an alpha-globin gene. Analysis of the J sequences suggests a recombination mechanism by which the alpha gene could have been deleted from the second set of genes in the cluster (zeta 2-zeta 3-theta 2). The relationships among the genes in characterized alpha-like gene clusters in mammals are summarized. The rabbit gene cluster differs from those of other mammals principally in the loss of a gene orthologous to the human psi alpha 1 and in the block duplication of the zeta-zeta-alpha-theta gene set.     

Synthesis of Polyamide Nucleic Acids (Pnas), Pna /Dna-Chimeras and Phosphonic Ester Nucleic Acids (Phonas)

Abstract

The synthesis of polyamide nucleic acids (PNAs) and derivatives thereof by different synthetic routes is described. The first strategy makes use of 9-Fluorenylmethoxycarbonyl (Fmoc)/monomethoxytrityl (Mmt) protected building blocks, whereas the second approach involves the use of Mmt/acyl protected monomers, which allows the preparation of PNADNA chimera. Additionally, a block coupling strategy is presented for the synthesis of novel phosphonic ester nucleic acids (PHONAs).     

Block alignment: New representation and comparison method to study evolution of genomes

Genomes are not random sequences because natural selection has injected information in biological sequences for billions of years. Inspired by this idea, we developed a simple method to compare genomes considering nucleotide counts in subsequences (blocks) instead of their exact sequences.We introduce the Block Alignment method for comparing two genomes and based on this comparison method, define a similarity score and a distance. The presented model ignores nucleotide order in the sequence. On the other hand, in this block comparison method, due to exclusion of point mutations and small size variations, there is no need for high coverage sequencing which is responsible for the high costs of data production and storage; moreover, the sequence comparisons could be performed with higher speed.Phylogenetic trees of two sets of bacterial genomes were constructed and the results were in full agreement with their already constructed phylogenetic trees. Furthermore, a weighted and directed similarity network of each set of bacterial genomes was inferred ab initio by this model. Remarkably, the communities of these networks are in agreement with the clades of the corresponding phylogenetic trees which means these similarity networks also contain phylogenetic information about the genomes. Moreover, the block comparison method was used to distinguish rob(15;21)c-associated iAMP21 and sporadic iAMP21 rearrangements in subgroups of chromosome 21 in acute lymphoblastic leukemia. Our results show a meaningful difference between the number of contigs that mapped to chromosomes 15 and 21 in these cases. Furthermore, the presented block alignment model can select the candidate blocks to perform more accurate analysis and it is capable to find conserved blocks on a set of genomes.     

A novel peptide derived from vascular endothelial growth factor prevents amyloid beta aggregation and toxicity

Amyloid-β oligomers (Aβo) are the most pathologically relevant Aβ species in Alzheimer's disease (AD), because they induce early synaptic dysfunction that leads to learning and memory impairments. In contrast, increasing VEGF (Vascular Endothelial Growth Factor) brain levels have been shown to improve learning and memory processes, and to alleviate Aβ-mediated synapse dysfunction. Here, we designed a new peptide, the blocking peptide (BP), which is derived from an Aβo-targeted domain of the VEGF protein, and investigated its effect on Aβ-associated toxicity. Using a combination of biochemical, 3D and ultrastructural imaging, and electrophysiological approaches, we demonstrated that BP strongly interacts with Aβo and blocks Aβ fibrillar aggregation process, leading to the formation of Aβ amorphous aggregates. BP further impedes the formation of structured Aβo and prevents their pathogenic binding to synapses. Importantly, acute BP treatment successfully rescues long-term potentiation (LTP) in the APP/PS1 mouse model of AD, at an age when LTP is highly impaired in hippocampal slices. Moreover, BP is also able to block the interaction between Aβo and VEGF, which suggests a dual mechanism aimed at both trapping Aβo and releasing VEGF to alleviate Aβo-induced synaptic damage. Our findings provide evidence for a neutralizing effect of the BP on Aβ aggregation process and pathogenic action, highlighting a potential new therapeutic strategy.     

Small iminium ions block gramicidin channels in lipid bilayers.

Guanidinium and acetamidinium, when added to the bathing solution in concentrations of approximately 0.1M, cause brief blocks in the single channel potassium currents from channels formed in planar lipid bilayers by gramicidin A. Single channel lifetimes are not affected indicating that the channel structure is not modified by the blockers. Guanidinium block durations and interblock times are approximately exponential in distribution. Block frequencies increase with guanidinium concentration whereas block durations are unaffected. Increases in membrane potential cause an increase in block frequency as expected for a positively charged blocker but a decrease in block duration suggesting that the block is relieved when the blocker passes through the channel. At low pH, urea, formamide, and acetamide cause similar blocks suggesting that the protonated species of these molecules also block. Arginine and several amines do not block. This indicates that only iminium ions which are small enough to enter the channel can cause blocks in gramicidin channels.     

Automated assembly of protein blocks for database searching.

A system is described for finding and assembling the most highly conserved regions of related proteins for database searching. First, an automated version of Smith's algorithm for finding motifs is used for sensitive detection of multiple local alignments. Next, the local alignments are converted to blocks and the best set of non-overlapping blocks is determined. When the automated system was applied successively to all 437 groups of related proteins in the PROSITE catalog, 1764 blocks resulted; these could be used for very sensitive searches of sequence databases. Each block was calibrated by searching the SWISS-PROT database to obtain a measure of the chance distribution of matches, and the calibrated blocks were concatenated into a database that could itself be searched. Examples are provided in which distant relationships are detected either using a set of blocks to search a sequence database or using sequences to search the database of blocks. The practical use of the blocks database is demonstrated by detecting previously unknown relationships between oxidoreductases and by evaluating a proposed relationship between HIV Vif protein and thiol proteases.     

Modelling gastrointestinal bioelectric activity

The development of an anatomically realistic biophysically based model of the human gastrointestinal (GI) tract is presented. A major objective of this work is to develop a modelling framework that can be used to integrate the physiological, anatomical and medical knowledge of the GI system. The anatomical model was developed by fitting derivative continuous meshes to digitised data taken from images of the visible man. Structural information, including fibre distributions of the smooth muscle layers and the arrangement of the networks of interstitial cells of Cajal, were incorporated using published information. A continuum modelling framework was used to simulate electrical activity from the single cell to the whole organ and body. Also computed was the external magnetic field generated from the GI electrical activity.

The set of governing equations were solved using a combination of numerical techniques. Activity at the (continuum) cell level was solved using a high-resolution trilinear finite element procedure that had been defined from the previously fitted C¹ continuous anatomical mesh. Multiple dipolar sources were created from the excitation waves which were embedded within a coupled C¹ continuous torso model to produce both the cutaneous electrical field and the external magnetic field.

Initial simulations were performed using a simplified geometry to test the implementation of the numerical solution procedure. The numerical procedures were shown to rapidly converge with mesh refinement. In the process of this testing, errors in a long standing analytic solution were identified and are corrected in Appendix B.

Results of single cell activity were compared to published results illustrating that the key features of the slow wave activity were successfully replicated. Simulations using a two-dimensional slice through the gastric wall produced slow wave activity that agreed with the known frequency and propagation characteristics. Three-dimensional simulations were also performed using the full stomach mesh and results illustrated the slow wave propagation throughout the stomach musculature.     

Multilevel Nonlinear Mixed-Effect Crown Ratio Models for Individual Trees of Mongolian Oak (Quercus mongolica) in Northeast China

In this study, an individual tree crown ratio (CR) model was developed with a data set from a total of 3134 Mongolian oak (Quercus mongolica) trees within 112 sample plots allocated in Wangqing Forest Bureau of northeast China. Because of high correlation among the observations taken from the same sampling plots, the random effects at levels of both blocks defined as stands that have different site conditions and plots were taken into account to develop a nested two-level nonlinear mixed-effect model. Various stand and tree characteristics were assessed to explore their contributions to improvement of model prediction. Diameter at breast height, plot dominant tree height and plot dominant tree diameter were found to be significant predictors. Exponential model with plot dominant tree height as a predictor had a stronger ability to account for the heteroskedasticity. When random effects were modeled at block level alone, the correlations among the residuals remained significant. These correlations were successfully reduced when random effects were modeled at both block and plot levels. The random effects from the interaction of blocks and sample plots on tree CR were substantially large. The model that took into account both the block effect and the interaction of blocks and sample plots had higher prediction accuracy than the one with the block effect and population average considered alone. Introducing stand density into the model through dummy variables could further improve its prediction. This implied that the developed method for developing tree CR models of Mongolian oak is promising and can be applied to similar studies for other tree species.     

A method for the analysis of domain movements in large biomolecular complexes

A new method for the analysis of domain movements in large, multichain, biomolecular complexes is presented. The method is applicable to any molecule for which two atomic structures are available that represent a conformational change indicating a possible domain movement. The method is blind to atomic bonding and atom type and can, therefore, be applied to biomolecular complexes containing different constituent molecules such as protein, RNA, or DNA. At the heart of the method is the use of blocks located at grid points spanning the whole molecule. The rotation vector for the rotation of atoms from each block between the two conformations is calculated. Treating components of these vectors as coordinates means that each block is associated with a point in a “rotation space” and that blocks with atoms that rotate together, perhaps as part of the same rigid domain, will have colocated points. Thus a domain can be identified from the clustering of points from blocks that span it. Domain pairs are accepted for analysis of their relative movements in terms of screw axes based upon a set of reasonable criteria. Here, we report on the application of the method to biomolecules covering a considerable size range: hemoglobin, liver alcohol dehydrogenase, S‐Adenosylhomocysteine hydrolase, aspartate transcarbamylase, and the 70S ribosome. The results provide a depiction of the conformational change within each molecule that is easily understood, giving a perspective that is expected to lead to new insights. Of particular interest is the allosteric mechanism in some of these molecules. Results indicate that common boundaries between subunits and domains are good regions to focus on as movement in one subunit can be transmitted to another subunit through such interfaces. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Models of neocortical layer 5b pyramidal cells capturing a wide range of dendritic and perisomatic active properties

The thick-tufted layer 5b pyramidal cell extends its dendritic tree to all six layers of the mammalian neocortex and serves as a major building block for the cortical column. L5b pyramidal cells have been the subject of extensive experimental and modeling studies, yet conductance-based models of these cells that faithfully reproduce both their perisomatic Na(+)-spiking behavior as well as key dendritic active properties, including Ca(2+) spikes and back-propagating action potentials, are still lacking. Based on a large body of experimental recordings from both the soma and dendrites of L5b pyramidal cells in adult rats, we characterized key features of the somatic and dendritic firing and quantified their statistics. We used these features to constrain the density of a set of ion channels over the soma and dendritic surface via multi-objective optimization with an evolutionary algorithm, thus generating a set of detailed conductance-based models that faithfully replicate the back-propagating action potential activated Ca(2+) spike firing and the perisomatic firing response to current steps, as well as the experimental variability of the properties. Furthermore, we show a useful way to analyze model parameters with our sets of models, which enabled us to identify some of the mechanisms responsible for the dynamic properties of L5b pyramidal cells as well as mechanisms that are sensitive to morphological changes. This automated framework can be used to develop a database of faithful models for other neuron types. The models we present provide several experimentally-testable predictions and can serve as a powerful tool for theoretical investigations of the contribution of single-cell dynamics to network activity and its computational capabilities.