Cooperative RNA Polymerase Molecules Behavior on a Stochastic Sequence-Dependent Model for Transcription Elongation

The transcription process is crucial to life and the enzyme RNA polymerase (RNAP) is the major component of the transcription machinery. The development of single-molecule techniques, such as magnetic and optical tweezers, atomic-force microscopy and single-molecule fluorescence, increased our understanding of the transcription process and complements traditional biochemical studies. Based on these studies, theoretical models have been proposed to explain and predict the kinetics of the RNAP during the polymerization, highlighting the results achieved by models based on the thermodynamic stability of the transcription elongation complex. However, experiments showed that if more than one RNAP initiates from the same promoter, the transcription behavior slightly changes and new phenomenona are observed. We proposed and implemented a theoretical model that considers collisions between RNAPs and predicts their cooperative behavior during multi-round transcription generalizing the Bai et al. stochastic sequence-dependent model. In our approach, collisions between elongating enzymes modify their transcription rate values. We performed the simulations in Mathematica® and compared the results of the single and the multiple-molecule transcription with experimental results and other theoretical models. Our multi-round approach can recover several expected behaviors, showing that the transcription process for the studied sequences can be accelerated up to 48% when collisions are allowed: the dwell times on pause sites are reduced as well as the distance that the RNAPs backtracked from backtracking sites. © 2013 Costa et al.

Modelling polymers with side chains: MEH-PPV and P3HT

Modelling polymers with side chains is always a challenge once the degrees of freedom are very high. In this study, we present a successful methodology to model poly[2-methoxy-5-(2′-ethyl-hexyloxy)-p-phenylenevinylene] (MEH-PPV) and poly[3-hexylthiophene] (P3HT) in solutions, taking into account the influence of side chains on the polymer conformation. Molecular dynamics and semi-empirical quantum mechanical methods were used for structure optimisation and evaluation of optical properties. The methodology allows to describe structural and optical characteristics of the polymers in a satisfactory way, as well as to evaluate some usual simplifications adopted for modelling these systems. Effective conjugation lengths of 8-14.6 and 21 monomers were obtained for MEH-PPV and P3HT, respectively, in accordance with experimental findings. In addition, anti/syn conformations of these polymers could be predicted based on intrinsic interactions of the lateral branches. © 2013 Copyright Taylor and Francis Group, LLC.

Fetal-Maternal Interactions in the Synepitheliochorial Placenta Using the eGFP Cloned Cattle Model

Background: To investigate mechanisms of fetal-maternal cell interactions in the bovine placenta, we developed a model of transgenic enhanced Green Fluorescent Protein (t-eGFP) expressing bovine embryos produced by nuclear transfer (NT) to assess the distribution of fetal-derived products in the bovine placenta. In addition, we searched for male specific DNA in the blood of females carrying in vitro produced male embryos. Our hypothesis is that the bovine placenta is more permeable to fetal-derived products than described elsewhere. Methodology/Principal Findings: Samples of placentomes, chorion, endometrium, maternal peripheral blood leukocytes and blood plasma were collected during early gestation and processed for nested-PCR for eGFP and testis-specific Y-encoded protein (TSPY), western blotting and immunohistochemistry for eGFP detection, as well as transmission electron microscopy to verify the level of interaction between maternal and fetal cells. TSPY and eGFP DNA were present in the blood of cows carrying male pregnancies at day 60 of pregnancy. Protein and mRNA of eGFP were observed in the trophoblast and uterine tissues. In the placentomes, the protein expression was weak in the syncytial regions, but intense in neighboring cells on both sides of the fetal-maternal interface. Ultrastructurally, our samples from t-eGFP expressing NT pregnancies showed to be normal, such as the presence of interdigitating structures between fetal and maternal cells. In addition, channels-like structures were present in the trophoblast cells. Conclusions/Significance: Data suggested that there is a delivery of fetal contents to the maternal system on both systemic and local levels that involved nuclear acids and proteins. It not clear the mechanisms involved in the transfer of fetal-derived molecules to the maternal system. This delivery may occur through nonclassical protein secretion; throughout transtrophoblastic-like channels and/or by apoptotic processes previously described. In conclusion, the bovine synepitheliochorial placenta displays an intimate fetal-maternal interaction, similar to other placental types for instance human and mouse. © 2013 Pereira et al.

Simulações por dinâmica molecular fine-e coarse-grained das interações intermoleculares entre peptídeos antimicrobianos da família Mastoparano e membranas modelo

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Estudos da interação de peptídeos antimicrobianos com modelo de membrana por simulações de dinâmica molecular

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Computer simulation study of collective dynamics in the glass former Ca(NO3)(2)center dot 4H(2)O

Time correlation functions of current fluctuations were calculated by molecular dynamics (MD) simulations in order to investigate sound waves of high wavevectors in the glass-forming liquid Ca(NO3)(2)center dot 4H(2)O. Dispersion curves, omega(k), were obtained for longitudinal (LA) and transverse acoustic (TA) modes, and also for longitudinal optic (LO) modes. Spectra of LA modes calculated by MD simulations were modeled by a viscoelastic model within the memory function framework. The viscoelastic model is used to rationalize the change of slope taking place at k similar to 0.3 angstrom(-1) in the omega(k) curve of acoustic modes. For still larger wavevectors, mixing of acoustic and optic modes is observed. Partial time correlation functions of longitudinal mass currents were calculated separately for the ions and the water molecules. The wavevector dependence of excitation energies of the corresponding partial LA modes indicates the coexistence of a relatively stiff subsystem made of cations and anions, and a softer subsystem made of water molecules. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4751548]

Molecular and Electronic Structure of the Peptide Subunit of Geobacter sulfurreducens Conductive Phi from First Principles

The respiration of metal oxides by the bacterium Geobacter sulfurreducens requires the assembly of a small peptide (the GS pilin) into conductive filaments termed pili. We gained insights into the contribution of the GS pilin to the pilus conductivity by developing a homology model and performing molecular dynamics simulations of the pilin peptide in vacuo and in solution. The results were consistent with a predominantly helical peptide containing the conserved a-helix region required for pilin assembly but carrying a short carboxy-terminal random-coiled segment rather than the large globular head of other bacterial pilins. The electronic structure of the pain was also explored from first principles and revealed a biphasic charge distribution along the pilin and a low electronic HOMO-LUMO gap, even in a wet environment. The low electronic band gap was the result of strong electrostatic fields generated by the alignment of the peptide bond dipoles in the pilin's alpha-helix and by charges from ions in solution and amino acids in the protein. The electronic structure also revealed some level of orbital delocalization in regions of the pilin containing aromatic amino acids and in spatial regions of high resonance where the HOMO and LUMO states are, which could provide an optimal environment for the hopping of electrons under thermal fluctuations. Hence, the structural and electronic features of the pilin revealed in these studies support the notion of a pilin peptide environment optimized for electron conduction.

Quantum spherical model with competing interactions

We analyse the phase diagram of a quantum mean spherical model in terms of the temperature T, a quantum parameter g, and the ratio p = -J(2)/J(1) where J(1) > 0 refers to ferromagnetic interactions between first-neighbour sites along the d directions of a hypercubic lattice, and J(2) < 0 is associated with competing anti ferromagnetic interactions between second neighbours along m <= d directions. We regain a number of known results for the classical version of this model, including the topology of the critical line in the g = 0 space, with a Lifshitz point at p = 1/4, for d > 2, and closed-form expressions for the decay of the pair correlations in one dimension. In the T = 0 phase diagram, there is a critical border, g(c) = g(c) (p) for d >= 2, with a singularity at the Lifshitz point if d < (m + 4)/2. We also establish upper and lower critical dimensions, and analyse the quantum critical behavior in the neighborhood of p = 1/4. 2012 (C) Elsevier B.V. All rights reserved.

The Role of Molecular Conformation and Polarizable Embedding for One- and Two-Photon Absorption of Disperse Orange 3 in Solution

Solvent effects on the one- and two-photon absorption (IPA and 2PA) of disperse orange 3 (DO3) in dimethyl sulfoxide (DMSO) are studied using a discrete polarizable embedding (PE) response theory. The scheme comprises a quantum region containing the chromophore and an atomically granulated classical region for the solvent accounting for full interactions within and between the two regions. Either classical molecular dynamics (MD) or hybrid Car-Parrinello (CP) quantum/classical (QM/MM) molecular dynamics simulations are employed to describe the solvation of DO3 in DMSO, allowing for an analysis of the effect of the intermolecular short-range repulsion, long-range attraction, and electrostatic interactions on the conformational changes of the chromophore and also the effect of the solute-solvent polarization. PE linear response calculations are performed to verify the character, solvatochromic shift, and overlap of the two lowest energy transitions responsible for the linear absorption spectrum of DO3 in DMSO in the visible spectral region. Results of the PE linear and quadratic response calculations, performed using uncorrelated solute-solvent configurations sampled from either the classical or hybrid CP QM/MM MD simulations, are used to estimate the width of the line shape function of the two electronic lowest energy excited states, which allow a prediction of the 2PA cross-sections without the use of empirical parameters. Appropriate exchange-correlation functionals have been employed in order to describe the charge-transfer process following the electronic transitions of the chromophore in solution.

Computational chemistry studies of UV induced processes in human skin

This thesis presents and uses the techniques of computational chemistry to explore two different processes induced in human skin by ultraviolet light. The first is the transformation of urocanic acid into a immunosuppressing agent, and the other is the enzymatic action of the 8-oxoguanine glycosylase enzyme. The photochemistry of urocanic acid is investigated by time-dependent density functional theory. Vertical absorption spectra of the molecule in different forms and environments is assigned and candidate states for the photochemistry at different wavelengths are identified. Molecular dynamics simulations of urocanic acid in gas phase and aqueous solution reveals considerable flexibility under experimental conditions, particularly for for the cis isomer where competition between intra- and inter-molecular interactions increases flexibility. A model to explain the observed gas phase photochemistry of urocanic acid is developed and it is shown that a reinterpretation in terms of a mixture between isomers significantly enhances the agreement between theory and experiment , and resolves several peculiarities in the spectrum. A model for the photochemistry in the aqueous phase of urocanic acid is then developed, in which two excited states governs the efficiency of photoisomerization. The point of entrance into a conical intersection seam is shown to explain the wavelength dependence of photoisomerization quantum yield. Finally some mechanistic aspects of the DNA repair enzyme 8-oxoguanine glycosylase is investigated with density functional theory. It is found that the critical amino acid of the active site can provide catalytic power in several different manners, and that a recent proposal involving a SN1 type of mechanism seems the most efficient one.

Variance Reduction in Analytical Chemistry : New Numerical Methods in Chemometrics and Molecular Simulation

This thesis is based on five papers addressing variance reduction in different ways. The papers have in common that they all present new numerical methods. Paper I investigates quantitative structure-retention relationships from an image processing perspective, using an artificial neural network to preprocess three-dimensional structural descriptions of the studied steroid molecules. Paper II presents a new method for computing free energies. Free energy is the quantity that determines chemical equilibria and partition coefficients. The proposed method may be used for estimating, e.g., chromatographic retention without performing experiments. Two papers (III and IV) deal with correcting deviations from bilinearity by so-called peak alignment. Bilinearity is a theoretical assumption about the distribution of instrumental data that is often violated by measured data. Deviations from bilinearity lead to increased variance, both in the data and in inferences from the data, unless invariance to the deviations is built into the model, e.g., by the use of the method proposed in paper III and extended in paper IV. Paper V addresses a generic problem in classification; namely, how to measure the goodness of different data representations, so that the best classifier may be constructed. Variance reduction is one of the pillars on which analytical chemistry rests. This thesis considers two aspects on variance reduction: before and after experiments are performed. Before experimenting, theoretical predictions of experimental outcomes may be used to direct which experiments to perform, and how to perform them (papers I and II). After experiments are performed, the variance of inferences from the measured data are affected by the method of data analysis (papers III-V).

Supramolecular order and dynamics of discotic materials studied by solid-state NMR recoupling methods

The topic of this thesis is the investigation of structure,order and dynamics in discotic mesogens by advancedsolid-state NMR spectroscopy. Most of the discotic mesogensunder investigation are hexa-peri-hexabenzocoronene (HBC)derivatives which are of particular interest for potentialdevice applications due to their high one-dimensional chargecarrier mobilities. The supramolecular stacking arrangement of the discoticcores was investigated by 2D 1H-1H double-quantum (DQ)methods, which were modified by incorporating the WATERGATEsuppression technique into the experiments in order toovercome severe phase problems arising from the strongsignal of the long alkyl sidechains. Molecular dynamics and sample orientation was probed throughthe generation of sideband patterns by reconversion rotorencoding in 2D recoupling experiments. These experimentswere extended by new recoupling schemes to enable thedistinction of motion and orientation effects. The solid-state NMR studies presented in this work aim tothe understanding of structure-property relationships in theinvestigated discotic materials, while the experimentsapplied to these materials include new recoupling schemeswhich make the desired information on molecular orientationand dynamics accessible without isotope labelling.

Development of empirical potentials for amorphous silica

Development of empirical potentials for amorphous silica Amorphous silica (SiO2) is of great importance in geoscience and mineralogy as well as a raw material in glass industry. Its structure is characterized as a disordered continuous network of SiO4 tetrahedra. Many efforts have been undertaken to understand the microscopic properties of silica by classical molecular dynamics (MD) simulations. In this method the interatomic interactions are modeled by an effective potential that does not take explicitely into account the electronic degrees of freedom. In this work, we propose a new methodology to parameterize such a potential for silica using ab initio simulations, namely Car-Parrinello (CP) method [Phys. Rev. Lett. 55, 2471 (1985)]. The new potential proposed is compared to the BKS potential [Phys. Rev. Lett. 64, 1955 (1990)] that is considered as the benchmark potential for silica. First, CP simulations have been performed on a liquid silica sample at 3600 K. The structural features so obtained have been compared to the ones predicted by the classical BKS potential. Regarding the bond lengths the BKS tends to underestimate the Si-O bond whereas the Si-Si bond is overestimated. The inter-tetrahedral angular distribution functions are also not well described by the BKS potential. The corresponding mean value of theSiOSi angle is found to be ≃ 147◦, while the CP yields to aSiOSi angle centered around 135◦. Our aim is to fit a classical Born-Mayer/Coulomb pair potential using ab initio calculations. To this end, we use the force-matching method proposed by Ercolessi and Adams [Europhys. Lett. 26, 583 (1994)]. The CP configurations and their corresponding interatomic forces have been considered for a least square fitting procedure. The classical MD simulations with the resulting potential have lead to a structure that is very different from the CP one. Therefore, a different fitting criterion based on the CP partial pair correlation functions was applied. Using this approach the resulting potential shows a better agreement with the CP data than the BKS ones: pair correlation functions, angular distribution functions, structure factors, density of states and pressure/density were improved. At low temperature, the diffusion coefficients appear to be three times higher than those predicted by the BKS model, however showing a similar temperature dependence. Calculations have also been carried out on crystalline samples in order to check the transferability of the potential. The equilibrium geometry as well as the elastic constants of α-quartz at 0 K are well described by our new potential although the crystalline phases have not been considered for the parameterization. We have developed a new potential for silica which represents an improvement over the pair potentials class proposed so far. Furthermore, the fitting methodology that has been developed in this work can be applied to other network forming systems such as germania as well as mixtures of SiO2 with other oxides (e.g. Al2O3, K2O, Na2O).

Molecular organization of n-cyanobiphenyl liquid crystals on a molybdenite surface

The alignement and anchoring of liquid crystals on solid surfaces is a key problem for modern device technology that until now has been treated empirically, but that can now be tackled by atomistic computer simulations. Molecular dynamics (MD) simulations were used in this thesis work to study two films of 7 and 8 n-alkyl-4’cyanobiphenyl (7CB and 8CB) liquid crystals , with a thickness of 15 nm, confined between two (001) surfaces of MoS2 (molybdenite). The isotropic and nematic phases of both liquid crystals were simulated, and the resulting structures characterized structurally. A new force field was designed to model the interactions between the liquid crystal (LC) molecules and the surface of molybdenite, while an accurate force field developed previously was used to model the 7CB and 8CB molecules. The results show that the (001) molybdenite surface induces a planar orientation in both the liquid crystals. For the nematic phase of 8CB, one of the two solid/LC interfaces is composed of a first layer of molecules aligned parallel to the surface, followed by a second layer of molecules aligned perpendicular to the surface (also called, homeotropic). The effect of the surface appears to be local in nature as it is confined to the first 15 Angström of the LC film. Conversely, for the nematic phase of 7CB, a planar ordering is established into the LC film. The LC molecules at the interface with the molybdenite appear to align preferentially their alkyl chains toward the solid substrate. The resulting tilt angle of molecules was found to be in good agreement with experimental measurements available in literature. Despite the fact that the MD simulations spanned a time range of more than 100 ns, the nematic phases of both 7CB and 8CB were found not to be completely formed. In order to confirm the findings presented in this thesis, we propose to extend the current study.

Dynamics, rheology and critical properties of colloidal fluid mixtures

Liquids under the influence of external fields exhibit a wide range of intriguing phenomena that can be markedly different from the behaviour of a quiescent system. This work considers two different systems — a glassforming Yukawa system and a colloid-polymer mixture — by Molecular Dynamics (MD) computer simulations coupled to dissipative particle dynamics. The former consists of a 50-50 binary mixture of differently-sized, like-charged colloids interacting via a screened Coulomb (Yukawa) potential. Near the glass transition the influence of an external shear field is studied. In particular, the transition from elastic response to plastic flow is of interest. At first, this model is characterised in equilibrium. Upon decreasing temperature it exhibits the typical dynamics of glassforming liquids, i.e. the structural relaxation time τα grows strongly in a rather small temperature range. This is discussed with respect to the mode-coupling theory of the glass transition (MCT). For the simulation of bulk systems under shear, Lees-Edwards boundary conditions are applied. At constant shear rates γ˙ ≫ 1/τα the relevant time scale is given by 1/γ˙ and the system shows shear thinning behaviour. In order to understand the pronounced differences between a quiescent system and a system under shear, the response to a suddenly commencing or terminating shear flow is studied. After the switch-on of the shear field the shear stress shows an overshoot, marking the transition from elastic to plastic deformation, which is connected to a super-diffusive increase of the mean squared displacement. Since the average static structure only depends on the value of the shear stress, it does not discriminate between those two regimes. The distribution of local stresses, in contrast, becomes broader as soon as the system starts flowing. After a switch-off of the shear field, these additional fluctuations are responsible for the fast decay of stresses, which occurs on a time scale 1/γ˙ . The stress decay after a switch-off in the elastic regime, on the other hand, happens on the much larger time scale of structural relaxation τα. While stresses decrease to zero after a switch-off for temperatures above the glass transition, they decay to a finite value for lower temperatures. The obtained results are important for advancing new theoretical approaches in the framework of mode-coupling theory. Furthermore, they suggest new experimental investigations on colloidal systems. The colloid-polymer mixture is studied in the context of the behaviour near the critical point of phase separation. For the MD simulations a new effective model with soft interaction potentials is introduced and its phase diagram is presented. Here, mainly the equilibrium properties of this model are characterised. While the self-diffusion constants of colloids and polymers do not change strongly when the critical point is approached, critical slowing down of interdiffusion is observed. The order parameter fluctuations can be determined through the long-wavelength limit of static structure factors. For this strongly asymmetric mixture it is shown how the relevant structure factor can be extracted by a diagonalisation of a matrix that contains the partial static structure factors. By presenting first results of this model under shear it is demonstrated that it is suitable for non-equilibrium simulations as well.