Sequence-level population simulations over large genomic regions

Simulation is an invaluable tool for investigating the effects of various population genetics modeling assumptions on resulting patterns of genetic diversity, and for assessing the performance of statistical techniques, for example those designed to detect and measure the genomic effects of selection. It is also used to investigate the effectiveness of various design options for genetic association studies. Backward-in-time simulation methods are computationally efficient and have become widely used since their introduction in the 1980s. The forward-in-time approach has substantial advantages in terms of accuracy and modeling flexibility, but at greater computational cost. We have developed flexible and efficient simulation software and a rescaling technique to aid computational efficiency that together allow the simulation of sequence-level data over large genomic regions in entire diploid populations under various scenarios for demography, mutation, selection, and recombination, the latter including hotspots and gene conversion. Our forward evolution of genomic regions (FREGENE) software is freely available from www.ebi.ac.uk/projects/BARGEN together with an ancillary program to generate phenotype labels, either binary or quantitative. In this article we discuss limitations of coalescent-based simulation, introduce the rescaling technique that makes large-scale forward-in-time simulation feasible, and demonstrate the utility of various features of FREGENE, many not previously available.     

Defective ribonucleoside diphosphate reductase impairs replication fork progression in Escherichia coli

The observed lengthening of the C period in the presence of a defective ribonucleoside diphosphate reductase has been assumed to be due solely to the low deoxyribonucleotide supply in the nrdA101 mutant strain. We show here that the nrdA101 mutation induces DNA double-strand breaks at the permissive temperature in a recB-deficient background, suggesting an increase in the number of stalled replication forks that could account for the slowing of replication fork progression observed in the nrdA101 strain in a Rec(+) context. These DNA double-strand breaks require the presence of the Holliday junction resolvase RuvABC, indicating that they have been generated from stalled replication forks that were processed by the specific reaction named "replication fork reversal." Viability results supported the occurrence of this process, as specific lethality was observed in the nrdA101 recB double mutant and was suppressed by the additional inactivation of ruvABC. None of these effects seem to be due to the limitation of the deoxyribonucleotide supply in the nrdA101 strain even at the permissive temperature, as we found the same level of DNA double-strand breaks in the nrdA(+) strain growing under limited (2-microg/ml) or under optimal (5-microg/ml) thymidine concentrations. We propose that the presence of an altered NDP reductase, as a component of the replication machinery, impairs the progression of the replication fork, contributing to the lengthening of the C period in the nrdA101 mutant at the permissive temperature.     

Influence of channel subunit composition on L-type Ca2+ current kinetics and cardiac wave stability

Previous studies have demonstrated that the slope of the function relating the action potential duration (APD) and the diastolic interval, known as the APD restitution curve, plays an important role in the initiation and maintenance of ventricular fibrillation. Since the APD restitution slope critically depends on the kinetics of the L-type Ca(2+) current, we hypothesized that manipulation of the subunit composition of these channels may represent a powerful strategy to control cardiac arrhythmias. We studied the kinetic properties of the human L-type Ca(2+) channel (Ca(v)1.2) coexpressed with the alpha(2)delta-subunit alone (alpha(1C) + alpha(2)delta) or in combination with beta(2a), beta(2b), or beta(3) subunits (alpha(1C) + alpha(2)delta + beta), using Ca(2+) as the charge carrier. We then incorporated the kinetic properties observed experimentally into the L-type Ca(2+) current mathematical model of the cardiac action potential to demonstrate that the APD restitution slope can be selectively controlled by altering the subunit composition of the Ca(2+) channel. Assuming that beta(2b) most closely resembles the native cardiac L-type Ca(2+) current, the absence of beta, as well as the coexpression of beta(2a), was found to flatten restitution slope and stabilize spiral waves. These results imply that subunit modification of L-type Ca(2+) channels can potentially be used as an antifibrillatory strategy.     

Sulfenamides as prodrugs of NH-acidic compounds: a new prodrug option for the amide bond

The objective of this report is to introduce the novel concept of utilizing sulfenamides as prodrugs for compounds containing an NH-acidic functionality, particularly weakly acidic amide-type functionalities (amides, ureas, carbamates, etc.). Included are the syntheses and physicochemical characterizations of some model sulfenamides to illustrate the promise of this new prodrug technology.     

Wnt/beta-catenin signaling in cancer stemness and malignant behavior

Stem cells are defined by their intrinsic capacity to self-renew and differentiate. Cancer stem cells retain both these features but have lost homeostatic mechanisms which maintain normal cell numbers. The canonical Wnt/beta-catenin signaling pathway plays a central role in modulating the delicate balance between stemness and differentiation in several adult stem cell niches such as the hair follicles in the skin, the mammary gland, and the intestinal crypt. Accordingly, constitutive Wnt signaling activation, resulting from mutations in genes encoding its downstream components, underlies tumorigenesis in these tissues. In the majority of sporadic colorectal cancer cases, the rate-limiting event is either loss of APC function or oncogenic beta-catenin mutations. However, although the presence of these initiating mutations would predict nuclear beta-catenin accumulation throughout the tumor mass, heterogeneous intracellular distributions of this key Wnt signaling molecule are observed within primary tumors and their metastases. In particular, tumor cells located at the invasive front and those migrating into the adjacent stromal tissues show nuclear beta-catenin staining. Hence, different levels of Wnt signaling activity reflect tumor heterogeneity and are likely to account for distinct cellular activities such as proliferation and epithelial-mesenchymal transitions, which prompt tumor growth and malignant behavior, respectively. Several intrinsic (cell-autonomous and/or autocrine) and extrinsic (paracrine, derived from the tumor microenvironment) factors may explain this heterogeneity of Wnt/beta-catenin signaling activity within the tumor mass.     

Synthesis and biological evaluation of non-peptide alpha(v)beta(3)/alpha(5)beta(1) integrin dual antagonists containing 5,6-dihydropyridin-2-one scaffolds

Small constrained non-peptidic molecules consisting of a polyfunctionalized rigid core, carrying appendages corresponding to arginine and aspartic acid side chains, have been recently reported to be promising for drug development. In this work, the 5,6-dihydropyridin-2-one was envisaged as a scaffold to turn into potential integrin ligands, introducing a carboxylic acid and a basic appendage. The synthesis and the antiadhesion activity of a small library of peptidomimetics capable to recognize alpha(v)beta(3) and alpha(5)beta(1) integrins has been herein reported.     

Squamate richness in the Brazilian Cerrado and its environmental–climatic associations

We investigate patterns of species richness of squamates (lizards, snakes, and amphisbaenians) in the Brazilian Cerrado, identifying areas of particularly high richness, and testing predictions of large‐scale richness hypotheses by analysing the relationship between species richness and environmental climatic variables. We used point localities from museum collections to produce maps of the predicted distributions for 237 Cerrado squamate species, using niche‐modelling techniques. We superimposed distributions of all species on a composite map, depicting richness across the ecosystem. Then, we performed a multiple regression analysis using eigenvector‐based spatial filtering (Principal Coordinate of Neighbour Matrices) to assess environmental–climatic variables that are best predictors of species richness. We found that the environmental–climatic and spatial filters multiple regression model explained 78% of the variation in Cerrado squamate richness (r² = 0.78; F = 32.66; P < 0.01). Best predictors of species richness were: annual precipitation, precipitation seasonality, altitude, net primary productivity, and precipitation during the driest quarter. A model selection approach revealed that several mechanisms related to the different diversity hypothesis might work together to explain richness variation in the Cerrado. Areas of higher species richness in Cerrado were located mainly in the south‐west, north, extreme east, and scattered areas in the north‐west portions of the biome. Partitioning of energy among species, habitat differentiation, and tolerance to variable environments may be the primary ecological factors determining variation in squamate richness across the Cerrado. High richness areas in northern Cerrado, predicted by our models, are still poorly sampled, and biological surveys are warranted in that region. The south‐western region of the Cerrado exhibits high species richness and is also undergoing high levels of deforestation. Therefore, maintenance of existing reserves, establishment of ecological corridors among reserves, and creation of new reserves are urgently needed to ensure conservation of species in these areas.     

Differences in the degree of inhibition of NDP reductase by chemical inactivation and by the thermosensitive mutation <Emphasis Type="Italic">nrdA101</Emphasis> in <Emphasis Type="Italic">Escherichia coli</Emphasis> suggest an effect on chromosome segregation

NDP reductase activity can be inhibited either by treatment with hydroxyurea or by incubation of an nrdA _ts mutant strain at the non-permissive temperature. Both methods inhibit replication, but experiments on these two types of inhibition yielded very different results. The chemical treatment immediately inhibited DNA synthesis but did not affect the cell and nucleoid appearance, while the incubation of an nrdA101 mutant strain at the non-permissive temperature inhibited DNA synthesis after more than 50 min, and resulted in aberrant chromosome segregation, long filaments, and a high frequency of anucleate cells. These phenotypes are not induced by SOS. In view of these results, we suggest there is an indirect relationship between NDP reductase and the chromosome segregation machinery through the maintenance of the proposed replication hyperstructure.     

Directional transition from initiation to elongation in bacterial translation

The transition of the 30S initiation complex (IC) to the translating 70S ribosome after 50S subunit joining provides an important checkpoint for mRNA selection during translation in bacteria. Here, we study the timing and control of reactions that occur during 70S IC formation by rapid kinetic techniques, using a toolbox of fluorescence-labeled translation components. We present a kinetic model based on global fitting of time courses obtained with eight different reporters at increasing concentrations of 50S subunits. IF1 and IF3 together affect the kinetics of subunit joining, but do not alter the elemental rates of subsequent steps of 70S IC maturation. After 50S subunit joining, IF2-dependent reactions take place independent of the presence of IF1 or IF3. GTP hydrolysis triggers the efficient dissociation of fMet-tRNA^fMet from IF2 and promotes the dissociation of IF2 and IF1 from the 70S IC, but does not affect IF3. The presence of non-hydrolyzable GTP analogs shifts the equilibrium towards a stable 70S–mRNA–IF1–IF2–fMet-tRNA^fMet complex. Our kinetic analysis reveals the molecular choreography of the late stages in translation initiation.