Sequence-function relationships in E. coli transcriptional regulation

Understanding how transcriptional regulatory sequence maps to regulatory function remains a difficult problem in regulatory biology. Given a particular DNA sequence for a bacterial promoter region, we would like to be able to say which transcription factors bind there, how strongly they bind, and whether they interact with each other and/or RNA polymerase, with the ultimate objective of integrating knowledge of these parameters into a prediction of gene expression levels. The theoretical framework of statistical thermodynamics provides a useful framework for doing so, enabling us to predict how gene expression levels depend on transcription factor binding energies and concentrations. We used thermodynamic models, coupled with models of the sequence-dependent binding energies of transcription factors and RNAP, to construct a genotype to phenotype map for the level of repression exhibited by the lac promoter, and tested it experimentally using a set of promoter variants from E. coli strains isolated from different natural environments. For this work, we sought to ``reverse engineer'' naturally occurring promoter sequences to understand how variations in promoter sequence affects gene expression. The natural inverse of this approach is to ``forward engineer'' promoter sequences to obtain targeted levels of gene expression. We used a high precision model of RNAP-DNA sequence dependent binding energy, coupled with a thermodynamic model relating binding energy to gene expression, to predictively design and verify a suite of synthetic E. coli promoters whose expression varied over nearly three orders of magnitude.

However, although thermodynamic models enable predictions of mean levels of gene expression, it has become evident that cell-to-cell variability or ``noise'' in gene expression can also play a biologically important role. In order to address this aspect of gene regulation, we developed models based on the chemical master equation framework and used them to explore the noise properties of a number of common E. coli regulatory motifs; these properties included the dependence of the noise on parameters such as transcription factor binding strength and copy number. We then performed experiments in which these parameters were systematically varied and measured the level of variability using mRNA FISH. The results showed a clear dependence of the noise on these parameters, in accord with model predictions.

Finally, one shortcoming of the preceding modeling frameworks is that their applicability is largely limited to systems that are already well-characterized, such as the lac promoter. Motivated by this fact, we used a high throughput promoter mutagenesis assay called Sort-Seq to explore the completely uncharacterized transcriptional regulatory DNA of the E. coli mechanosensitive channel of large conductance (MscL). We identified several candidate transcription factor binding sites, and work is continuing to identify the associated proteins.

Wurtzite-to-tetragonal structure phase transformation and size effect in ZnO nanorods

The deformation of [0001]-oriented ZnO nanorods with hexagonal cross sections under uniaxial tensile loading is analyzed through molecular statistical thermodynamics (MST) simulations. The focus is on the size dependence of mechanical behavior in ZnO nanorods with diameters ranging from 1.95 to 17.5 nm. An irreversible phase transformation from the wurtzite (P6(3)mc space group) structure to a tetragonal structure (P4(2)/mnm space group) occurs during the tensile loading process. Young's modulus before the transformation demonstrates a size dependence consistent with what is observed in experiments. A stronger size dependence of response is seen after the transformation and is attributed to the polycrystalline nature of the transformed structure. A comparison of the MST and molecular dynamics (MD) methods shows that MST is 60 times faster than MD and yields results consistent with the results of MD.

Molecular weight dependence of phase Behavior of PEO/P(EO-b-DMS) blends: Application of Sanchez-Lacombe lattice fluid theory

Molecular weight dependence of phase separation behavior of the Poly (ethylene oxide) (PEO)/Poly(ethylene oxide-block-dimethylsiloxane) (P(EO-b-DMS)) blends was investigated by both experimental and theoretical methods. The cloud point curves of PEO/P(EO-b-DMS) blends were obtained by turbidity method. Based on Sanchez-Lacombe lattice fluid theory (SLLFT), the adjustable parameter, epsilon*(12)/k (quantifying the interaction energy between different components), was evaluated by fitting the experimental data in phase diagrams. To calculate the spinodals, binodals, and the volume changes of mixing for these blends, three modified combining rules of the scaling parameters for the block copolymer were introduced.

Effect of chain-length dependence of interaction parameter on spinodals for polydisperse polymer blends

The chain-length dependence of the Flory-Huggins (FH) interaction parameter is introduced into the FH lattice theory for polydisperse polymer-blend systems. The spinodals are calculated for the model polymer blends with different chain lengths and distributions. It is found that all the related variables r(n), r(w), r(z), and chain-length distribution, have effects on the spinodals for polydisperse polymer blends.

What has brought on the effects of number-average molecular weight on the spinodals in polymer mixtures?

On the basis of the thermodynamics of Gibbs, the spinodal for the quasibinary system was derived in the framework of the Sanchez-Lacombe lattice fluid theory. All of the spinodals were calculated based on a model polydisperse polymer mixture, where each polymer contains three different molecular weight subcomponents. According to our calculations, the spinodal depends on both weight-average ((M) over bar (w)) and number-average ((M) over bar (n)) molecular weights, whereas that of the z-average molecular weight is invisible. Moreover, the extreme of the spinodal decreases when the polydispersity index (eta = (M) over bar (w)/(M) over bar (n)) of the polymer increases. The effect of polydispersity on the spinodal decreases when the molecular weight gets larger and can be negligible at a certain large molecular weight. It is well-known that the influence of polydispersity on the phase equilibrium (coexisting curve, cloud point curves) is much more pronounced than on the spinodal. The effect of M, on the spinodal is discussed as it results from the infuluence of composition temperatures, molecular weight, and the latter's distribution on free volume. An approximate expression, which is in the assumptions of v* v(1)* = v(2)* and 1/r --> 0 for both of the polymers, was also derived for simplification. It can be used in high molecular weight, although it failed to make visible the effect of number-average molecular weight on the spinodal.

Comment on "Kinetics of spinodal decomposition in the ising model with vacancy diffusion"

In a recent paper [Phys. Rev. B 50, 3477 (1994)], P. Fratzl and O. Penrose present the results of the Monte Carlo simulation of the spinodal decomposition problem (phase separation) using the vacancy dynamics mechanism. They observe that the t1/3 growth regime is reached faster than when using the standard Kawasaki dynamics. In this Comment we provide a simple explanation for the phenomenon based on the role of interface diffusion, which they claim is irrelevant for the observed behavior.

Vacancy-assisted domain-growth in asymmetric binary alloys: A Monte Carlo Study

A Monte Carlo simulation study of the vacancy-assisted domain growth in asymmetric binary alloys is presented. The system is modeled using a three-state ABV Hamiltonian which includes an asymmetry term. Our simulated system is a stoichiometric two-dimensional binary alloy with a single vacancy which evolves according to the vacancy-atom exchange mechanism. We obtain that, compared to the symmetric case, the ordering process slows down dramatically. Concerning the asymptotic behavior it is algebraic and characterized by the Allen-Cahn growth exponent x51/2. The late stages of the evolution are preceded by a transient regime strongly affected by both the temperature and the degree of asymmetry of the alloy. The results are discussed and compared to those obtained for the symmetric case.

Estudio teórico de propiedades termodinámicas y optoelectronicas de nuevos materiales fotovoltaicos de banda intermedia = Theoretical study of thermodynamic and optoelectronic properties of new intermediate band photovoltaic materials

Los materiales de banda intermedia han atraido la atención de la comunidad científica en el campo de la energía solar fotovoltaica en los últimos años. Sin embargo, con el objetivo de entender los fundamentos de las células solares de banda intermedia, se debe llevar a cabo un estudio profundo de la características de los materiales. Esto se puede hacer mediante un modelo teórico usando Primeros Principios. A partir de este enfoque se pueden obtener resultados tales como la estructura electrónica y propiedades ópticas, entre otras, de los semiconductores fuertemente dopados y sus precursores. Con el fin de desentrañar las estructuras de estos sistemas electrónicos, esta tesis presenta un estudio termodinámico y optoelectrónico de varios materiales fotovoltaicos. Específicamente se caracterizaron los materiales avanzados de banda intermedia y sus precursores. El estudio se hizo en términos de caracterización teórica de la estructura electrónica, la energética del sistema, entre otros. Además la estabilidad se obtuvo usando configuraciones adaptadas a la simetría del sistema y basado en la combinatoria. Las configuraciones de los sitios ocupados por defectos permiten obtener información sobre un espacio de configuraciones donde las posiciones de los dopantes sustituidos se basan en la simetría del sólido cristalino. El resultado puede ser tratado usando elementos de termodinámica estadística y da información de la estabilidad de todo el espacio simétrico. Además se estudiaron otras características importantes de los semiconductores de base. En concreto, el análisis de las interacciones de van der Waals fueron incluidas en el semiconductor en capas SnS2, y el grado de inversión en el caso de las espinelas [M]In2S4. En este trabajo además realizamos una descripción teórica exhaustiva del sistema CdTe:Bi. Este material de banda-intermedia muestra características que son distintas a las de los otros materiales estudiados. También se analizó el Zn como agente modulador de la posición de las sub-bandas prohibidas en el material de banda-intermedia CuGaS2:Ti. Analizándose además la viabilidad termodinámica de la formación de este compuesto. Finalmente, también se describió el GaN:Cr como material de banda intermedia, en la estructura zinc-blenda y en wurtztite, usando configuraciones de sitios ocupados de acuerdo a la simetría del sistema cristalino del semiconductor de base. Todos los resultados, siempre que fue posible, fueron comparados con los resultados experimentales. ABSTRACT The intermediate-band materials have attracted the attention of the scientific community in the field of the photovoltaics in recent years. Nevertheless, in order to understand the intermediate-band solar cell fundamentals, a profound study of the characteristics of the materials is required. This can be done using theoretical modelling from first-principles. The electronic structure and optical properties of heavily doped semiconductors and their precursor semiconductors are, among others, results that can be obtained from this approach. In order to unravel the structures of these crystalline systems, this thesis presents a thermodynamic and optoelectronic study of several photovoltaic materials. Specifically advanced intermediate-band materials and their precursor semiconductors were characterized. The study was made in terms of theoretical characterization of the electronic structure, energetics among others. The stability was obtained using site-occupancy-disorder configurations adapted to the symmetry of the system and based on combinatorics. The site-occupancy-disorder method allows the formation of a configurational space of substitutional dopant positions based on the symmetry of the crystalline solid. The result, that can be treated using statistical thermodynamics, gives information of the stability of the whole space of symmetry of the crystalline lattice. Furthermore, certain other important characteristics of host semiconductors were studied. Specifically, the van der Waal interactions were included in the SnS2 layered semiconductor, and the inversion degree in cases of [M]In2S4 spinels. In this work we also carried out an exhaustive theoretical description of the CdTe:Bi system. This intermediate-band material shows characteristics that are distinct from those of the other studied intermediate-band materials. In addition, Zn was analysed as a modulator of the positions of the sub-band gaps in the CuGaS2:Ti intermediate-band material. The thermodynamic feasibility of the formation of this compound was also carried out. Finally GaN:Cr intermediate-band material was also described both in the zinc-blende and the wurtztite type structures, using the symmetry-adapted-space of configurations. All results, whenever possible, were compared with experimental results.

Disjoining pressure and the film-height-dependent surface tension of thin liquid films: New insight from capillary wave fluctuations

In this paper we review simulation and experimental studies of thermal capillary wave fluctuations as an ideal means for probing the underlying disjoining pressure and surface tensions, and more generally, fine details of the Interfacial Hamiltonian Model. We discuss recent simulation results that reveal a film-height-dependent surface tension not accounted for in the classical Interfacial Hamiltonian Model. We show how this observation may be explained bottom-up from sound principles of statistical thermodynamics and discuss some of its implications

An information theory model of hydrophobic interactions.

A molecular model of poorly understood hydrophobic effects is heuristically developed using the methods of information theory. Because primitive hydrophobic effects can be tied to the probability of observing a molecular-sized cavity in the solvent, the probability distribution of the number of solvent centers in a cavity volume is modeled on the basis of the two moments available from the density and radial distribution of oxygen atoms in liquid water. The modeled distribution then yields the probability that no solvent centers are found in the cavity volume. This model is shown to account quantitatively for the central hydrophobic phenomena of cavity formation and association of inert gas solutes. The connection of information theory to statistical thermodynamics provides a basis for clarification of hydrophobic effects. The simplicity and flexibility of the approach suggest that it should permit applications to conformational equilibria of nonpolar solutes and hydrophobic residues in biopolymers.

Grand canonical Monte Carlo simulation study of methane adsorption at an open graphite surface and in slitlike carbon pores at 273 K

Grand canonical Monte Carlo (GCMC) simulation was used for the systematic investigation of the supercritical methane adsorption at 273 K on an open graphite surface and in slitlike micropores of different sizes. For both considered adsorption systems the calculated excess adsorption isotherms exhibit a maximum. The effect of the pore size on the maximum surface excess and isosteric enthalpy of adsorption for methane storage at 273 K is discussed. The microscopic detailed picture of methane densification near the homogeneous graphite wall and in slitlike pores at 273 K is presented with selected local density profiles and snapshots. Finally, the reliable pore size distributions, obtained in the range of the microporosity, for two pitch-based microporous activated carbon fibers are calculated from the local excess adsorption isotherms obtained via the GCMC simulation. The current systematic study of supercritical methane adsorption both on an open graphite surface and in slitlike micropores performed by the GCMC summarizes recent investigations performed at slightly different temperatures and usually a lower pressure range by advanced methods based on the statistical thermodynamics.

Prediction of vapor-moisture equilibriums for a wood-moisture system using a modified UNIQUAC model

A modified UNIQUAC model has been extended to describe and predict the equilibrium relative humidity and moisture content for wood. The method is validated over a range of moisture content from oven-dried state to fiber saturation point, and over a temperature range of 20-70 degrees C. Adjustable parameters and binary interaction parameters of the UNIQUAC model were estimated from experimental data for Caribbean pine and Hoop pine as well as data available in the literature. The two group-interaction parameters for the wood-moisture system were consistent with using function group contributions for H2O, -OH and -CHO. The result reconfirms that the main contributors to water adsorption in cell walls are the hydroxyl groups of the carbohydrates in cellulose and hemicelluloses. This provides some physical insight into the intermolecular force and energy between bound water and the wood material. (c) 2006 Elsevier Ltd. All rights reserved.

NMR microscopy: A tool for measuring monomer diffusion in supercritical CO2

The diffusion of styrene into linear low density polyethylene in a solution of supercritical CO2 was investigated using NMR microimaging. For both pure styrene and styrene dissolved in supercritical CO2, the diffusion was found to follow Fickian kinetics. Supercritical CO2 was found to enhance the rate and extent of diffusion of styrene into the substrate by up to three times under the conditions of this investigation, compared to pure styrene. NMR imaging was used to measure the concentration profiles of the styrene penetrants in real time, and the results were fitted to a Fickian model for diffusion. At a CO2 pressure of 150 bar and temperature of 40 degrees C, the diffusion coefficient of a 30 wt-% solution of styrene into LLDPE was calculated to be 1 X 10(-11) m(2). s(-1). This is significantly faster than the diffusion coefficient measured for pure styrene diffusion at 40 degrees C (3 x 10(-12) m(2). s(-1)). The diffusion coefficients determined by gravimetric analysis were slightly higher than those determined by the imaging method. This was probably due to residual styrene and/or polystyrene adhering to the surface of the films in the gravimetric technique.