Theoretical Study of Electrical Conduction Through a Molecule Connected to Metallic Nanocontacts

We present a theoretical study of electron transport through a molecule connected to two metallic nanocontacts. The system investigated is 1,4 benzene-dithiolate (BDT) chemically bonded to two Au contacts. The surface chemistry is modeled by representing the tips of the Au contacts as two atomic clusters and treating the molecule-cluster complex as a single entity in an extended Huckel tight binding scheme. We model the tips using several different cluster geometries. An ideal lead is attached to each cluster, and the lead to lead transmission is calculated. The role of the molecule-cluster interaction in transport is analyzed by using single channel leads. We then extend the calculations to multi-channel leads that are a more realistic model of the tip's environment. Using the finite-voltage, finite temperature Landauer formula, we calculate the differential conductance for the different systems studied. The similarities and differences between the predictions of the present class of models and recent experimental work are discussed.Comment: 23 pages, 11 figures, accepted PR

Correlation between sequence hydrophobicity and surface-exposure pattern of database proteins

Hydrophobicity is thought to be one of the primary forces driving the folding of proteins. On average, hydrophobic residues occur preferentially in the core, whereas polar residues tends to occur at the surface of a folded protein. By analyzing the known protein structures, we quantify the degree to which the hydrophobicity sequence of a protein correlates with its pattern of surface exposure. We have assessed the statistical significance of this correlation for several hydrophobicity scales in the literature, and find that the computed correlations are significant but far from optimal. We show that this less than optimal correlation arises primarily from the large degree of mutations that naturally occurring proteins can tolerate. Lesser effects are due in part to forces other than hydrophobicity and we quantify this by analyzing the surface exposure distributions of all amino acids. Lastly we show that our database findings are consistent with those found from an off-lattice hydrophobic-polar model of protein folding.Comment: 16 pages, 2 tables, 8 figure

Non-Equilibrium Polar Localization of Proteins in Bacterial Cells

Many proteins are observed to localize to the poles within bacterial cells. Some bacteria show unipolar localization, yet under different conditions bipolar patterns can emerge. One mechanism for spontaneous polar localization has been shown to involve the combination of protein aggregation and nucleoid occlusion. Whether the different observed patterns represent global energy minima for the cellular system remains to be determined. In this paper we show that for a model consisting only of protein aggregation along with an excluded volume effect due to the DNA polymer, that unipolar patterns are the global energy ground state regardless of protein concentration and DNA density. We extend the model to allow for proteins to be added to the cellular volume at a constant rate and show that bipolar (or multi-foci) patterns emerge as the result of the system being kinetically trapped in a local energy minimum. Lastly we also consider the situation of a growing cell that starts with a pre-existing aggregate at one of the poles and determine conditions under which either unipolar or bipolar patterns can exist at the point when it is ready to divide. This work sheds new interpretations on recently published experimental data and suggests experiments to test whether such a mechanism can drive patterning in bacteria

A Model for Cell Population Size Control Using Asymmetric Division

In multicellular organisms one can find examples where a growing tissue divides up until some final fixed cell number. Asymmetric division plays a prevalent feature in tissue differentiation in these organisms, where the daughters of each asymmetric division inherit unequal amounts of a fate determining molecule and as a result follow different developmental fates. In some tissues the accumulation or decrease of cell cycle regulators acts as an intrinsic timing mechanism governing proliferation. Here we present a minimal model based on asymmetric division and dilution of a cell-cycle regulator that can generate any final population size that might be needed. We show that within the model there are a variety of growth mechanisms from linear to non-linear that can lead to the same final cell count. Interestingly, when we include noise at division we find that there are special final cell population sizes that can be generated with high confidence that are flanked by population sizes that are less robust to division noise. When we include further perturbations in the division process we find that these special populations can remain relatively stable and in some cases even improve in their fidelity

Operational Principles for the Dynamics of the In Vitro ParA-ParB System

In many bacteria the ParA-ParB protein system is responsible for actively segregating DNA during replication. ParB proteins move by interacting with DNA bound ParA-ATP, stimulating their unbinding by catalyzing hydrolysis, that leads to rectified motion due to the creation of a wake of depleted ParA. Recent in vitro experiments have shown that a ParB covered magnetic bead can move with constant speed over a DNA covered substrate that is bound by ParA. It has been suggested that the formation of a gradient in ParA leads to diffusion-ratchet like motion of the ParB bead but how it forms and generates a force is still a matter of exploration. Here we develop a deterministic model for the in vitro ParA-ParB system and show that a ParA gradient can spontaneously form due to any amount of initial spatial noise in bound ParA. The speed of the bead is independent of this noise but depends on the ratio of the range of ParA-ParB force on the bead to that of removal of surface bound ParA by ParB. We find that at a particular ratio the speed attains a maximal value. We also consider ParA rebinding (including cooperativity) and ParA surface diffusion independently as mechanisms for ParA recovery on the surface. Depending on whether the DNA covered surface is undersaturated or saturated with ParA, we find that the bead can accelerate persistently or potentially stall. Our model highlights key requirements of the ParA-ParB driving force that are necessary for directed motion in the in vitro system that may provide insight into the in vivo dynamics of the ParA-ParB system

Conservation of regulatory elements between two species of Drosophila

BACKGROUND: One of the important goals in the post-genomic era is to determine the regulatory elements within the non-coding DNA of a given organism's genome. The identification of functional cis-regulatory modules has proven difficult since the component factor binding sites are small and the rules governing their arrangement are poorly understood. However, the genomes of suitably diverged species help to predict regulatory elements based on the generally accepted assumption that conserved blocks of genomic sequence are likely to be functional. To judge the efficacy of strategies that prefilter by sequence conservation it is important to know to what extent the converse assumption holds, namely that functional elements common to both species will fall within these conserved blocks. The recently completed sequence of a second Drosophila species provides an opportunity to test this assumption for one of the experimentally best studied regulatory networks in multicellular organisms, the body patterning of the fly embryo. RESULTS: We find that 50%–70% of known binding sites reside in conserved sequence blocks, but these percentages are not greatly enriched over what is expected by chance. Finally, a computational genome-wide search in both species for regulatory modules based on clusters of binding sites suggests that genes central to the regulatory network are consistently recovered. CONCLUSIONS: Our results indicate that binding sites remain clustered for these "core modules" while not necessarily residing in conserved blocks. This is an important clue as to how regulatory information is encoded in the genome and how modules evolve

DNA Segregation Under Par Protein Control

The spatial organization of DNA is mediated by the Par protein system in some bacteria. ParB binds specifically to the parS sequence on DNA and orchestrates its motion by interacting with ParA bound to the nucleoid. In the case of plasmids, a single ParB bound plasmid is observed to execute oscillations between cell poles while multiple plasmids eventually settle at equal distances from each other along the cell’s length. While the potential mechanism underlying the ParA-ParB interaction has been discussed, it remains unclear whether ParB-complex oscillations are stable limit cycles or merely decaying transients to a fixed point. How are dynamics affected by substrate length and the number of complexes? We present a deterministic model for ParA-ParB driven DNA segregation where the transition between stable arrangements and oscillatory behaviour depends only on five parameters: ParB-complex number, substrate length, ParA concentration, ParA hydrolysis rate and the ratio of the lengthscale over which the ParB complex stimulates ParA hydrolysis to the lengthscale over which ParA interacts with the ParB complex. When the system is buffered and the ParA rebinding rate is constant we find that ParB-complex dynamics is independent of substrate length and complex number above a minimum system size. Conversely, when ParA resources are limited, we find that changing substrate length and increasing complex number leads to counteracting mechanisms that can both generate or subdue oscillatory dynamics. We argue that cells may be poised near a critical level of ParA so that they can transition from oscillatory to fixed point dynamics as the cell cycle progresses so that they can both measure their size and faithfully partition their genetic material. Lastly, we show that by modifying the availability of ParA or depletion zone size, we can capture some of the observed differences in ParB-complex positioning between replicating chromosomes in B. subtilis cells and low-copy plasmids in E. coli cells

Identifying Proteins of High Designability via Surface-Exposure Patterns

Using an off-lattice model, we fully enumerate folded conformations of polypeptide chains of up to N = 19 monomers. Structures are found to differ markedly in designability, defined as the number of sequences with that structure as a unique lowest-energy conformation. We find that designability is closely correlated with the pattern of surface exposure of the folded structure. For longer chains, complete enumeration of structures is impractical. Instead, structures can be randomly sampled, and relative designability estimated either from designability within the random sample, or directly from surface-exposure pattern. We compare the surface-exposure patterns of those structures identified as highly designable to the patterns of naturally occurring proteins.Comment: 17 pages, 12 figure