Charge transfer in DNA: hole charge is confined to a single base pair due to solvation effects

We include solvation effects in tight-binding Hamiltonians for hole states in DNA. The corresponding linear-response parameters are derived from accurate estimates of solvation energy calculated for several hole charge distributions in DNA stacks. Two models are considered: (A) the correction to a diagonal Hamiltonian matrix element depends only on the charge localized on the corresponding site and (B) in addition to this term, the reaction field due to adjacent base pairs is accounted for. We show that both schemes give very similar results. The effects of the polar medium on the hole distribution in DNA are studied. We conclude that the effects of polar surroundings essentially suppress charge delocalization in DNA, and hole states in (GC)n sequences are localized on individual guanines

Electronic and vibrational contributions to first hyperpolarizability of donor-acceptor-substituted azobenzene

In this study we report on the electronic and vibrational (hyper)polarizabilities of donor���acceptorsubstituted azobenzene. It is observed that both electronic and vibrational contributions to the electricdipole first hyperpolarizability of investigated photoactive molecule substantially depend on the conformation. The contributions to the nuclear relaxation first hyperpolarizability are found to be quite important in the case of two considered isomers (cis and trans). Although the double-harmonic term is found to be the largest in terms of magnitude, it is shown that the total value of the nuclear relaxation contribution to vibrational first hyperpolarizability is a result of subtle interplay of higher-order contributions. As a part of the study, we also assess the performance of long-range-corrected densityfunctional theory in determining vibrational contributions to electric dipole (hyper)polarizabilities. In most cases, the applied long-range-corrected exchange correlation potentials amend the drawbacks of their conventional counterparts

Fronts from complex two-dimensional dispersal kernels: theory and application to Reid's paradox

Bimodal dispersal probability distributions with characteristic distances differing by several orders of magnitude have been derived and favorably compared to observations by Nathan [Nature (London) 418, 409 (2002)]. For such bimodal kernels, we show that two-dimensional molecular dynamics computer simulations are unable to yield accurate front speeds. Analytically, the usual continuous-space random walks (CSRWs) are applied to two dimensions. We also introduce discrete-space random walks and use them to check the CSRW results (because of the inefficiency of the numerical simulations). The physical results reported are shown to predict front speeds high enough to possibly explain Reid's paradox of rapid tree migration. We also show that, for a time-ordered evolution equation, fronts are always slower in two dimensions than in one dimension and that this difference is important both for unimodal and for bimodal kernels

Functional analysis of a conserved DNA methyltransferase in caulobacter crescentus

CcrM is a DNA methyltransferase that methylates the adenine in GANTC motifs in the chromo-some of the bacterial model Caulobacter crescentus. The loss of the CcrM homolog is lethal in C. crescentus and in several other species of Alphaproteobacteria. In this research, we used different experimental and bioinformatic approaches to determine why CcrM is so critical to the physiology of C. crescentus. We first showed that CcrM is a resident orphan DNA methyltransferase in non-Rickettsiales Alphaproteobacteria and that its gene is strictly conserved in this clade (with only one ex��ception among the genomes sequenced so far). In C. crescentus, cells depleted in CcrM in rich medium quickly lose viability and present an elongated phenotype characteristic of an im��pairment in cell division. Using minimal medium instead of rich medium as selective and main��tenance substrate, we could generate a AccrM mutant that presents a viability comparable to the wild type strain and only mild morphological defects. On the basis of a transcriptomic ap��proach, we determined that several genes essential for cell division were downregulated in the AccrM strain in minimal medium. We offered decisive arguments to support that the efficient transcription of two of these genes, ftsZ and mipZ, coding respectively for the Z-ring forming GTPase FtsZ and an inhibitor of FtsZ polymerization needed for the correct positioning of the Z- ring at mid-cell, requires the methylation of an adenine in a conserved GANTC motif located in their core promoter region. We propose a model, according to which the genome of C. crescentus encodes a transcriptional activator that requires a methylated adenine in a GANTC context to bind to DNA and suggest that this transcriptional regulator might be the global cell-cycle regulator GcrA. In addition, combining a classic genetic approach and in vitro evolution experiments, we showed that the mortality and cell division defects of the AccrM strain in rich medium are mainly due to limiting intracellular levels of the FtsZ protein. We also studied the dynamics of GANTC methylation in C. crescentus using the SMRT technol��ogy developed by Pacific Biosciences. Our findings support the commonly accepted model, accord��ing to which the methylation state of GANTC motifs varies during the cell cycle of C. crescentus: before the initiation of DNA replication, the GANTC motifs are fully-methylated (methylated on both strands); when the DNA gets replicated, the GANTC motifs become hemi-methylated (methyl��ated on one strand only) and this occurs at different times during replication for different loci along the chromosome depending on their position relative to the origin of replication; the GANTC mo��tifs are only remethylated after DNA replication has finished as a consequence of the massive and short-lived expression of CcrM in predivisional cells. About 30 GANTC motifs in the C. crescentus chromosome were found to be undermethylated in most of the bacterial population; these might be protected from CcrM activity by DNA binding proteins and some of them could be involved in methylation-based bistable transcriptional switches. - CcrM est une ADN m��thyltransf��rase qui m��thyle les ad��nines dans le contexte GANTC dans le g��nome de la bact��rie mod��le Caulobacter crescentus. La perte de l'homologue de CcrM chez C. crescentus et chez plusieurs autres esp��ces d'Alphaproteobact��ries est l��tale. Dans le courant de cette recherche, nous tentons de d��terminer pourquoi la prot��ine CcrM est cruciale pour la survie de C. crescentus. Nous d��montrons d'abord que CcrM est une ad��nine m��thyltransf��rase orpheline r��sidente, dont le g��ne fait partie du g��nome minimal partag�� par les Alphaprot��obact��ries non-Rickettsiales (�� une exception pr��s). Lorsqu'une souche de C. crescentus est priv��e de CcrM, sa viabilit�� d��cro��t rapi��dement et ses cellules pr��sentent une morphologie allong��e qui sugg��re que la division cellulaire est inhib��e. Nous sommes parvenus �� cr��er une souche AccrM en utilisant un milieu minimum, au lieu du milieu riche classiquement employ��, comme milieu de s��lection et de maintenance pour la souche. Lorsque nous avons ��tudi�� le transcriptome de cette souche de C. crescentus priv��e de CcrM, nous avons pu constater que plusieurs g��nes essentiels pour le bon d��roulement de la division cellulaire bact��rienne ��taient r��prim��s. En particulier, l'expression ad��quate des g��nes ftsZ et mipZ - qui codent, respectivement, pour FtsZ, la prot��ine qui constitue, au milieu de la cellule, un anneau prot��ique qui initie le processus de division et pour MipZ, un inhibiteur de la polym��risation de FtsZ qui est indispensable pour le bon positionnement de l'anneau FtsZ - est d��pendante de la pr��sence d'une ad��nine m��thyl��e dans un motif GANTC conserv�� situ�� dans leur r��gion promotrice. Nous pr��sentons un mod��le selon lequel le g��nome de C. crescentus code pour un facteur de transcription qui exige la pr��sence d'une ad��nine m��thyl��e dans un contexte GANTC pour s'attacher �� l'ADN et nous sugg��rons qu'il pourrait s'agir du r��gulateur global du cycle cellulaire GcrA. En outre, nous montrons, en combinant la g��n��tique classique et une approche bas��e sur l'��volution exp��rimentale, que la mortalit�� et l'inhibition de la division cellulaire caract��ristiques de la souche ��ccrMeη milieu riche sont dues �� des niveaux excessivement bas de prot��ine FtsZ. Nous avons aussi ��tudi�� la dynamique de la m��thylation du chromosome de C. crescentus sur la base de la technologie SMRT d��velopp��e par Pacific Biosciences. Nous confirmons le mod��le commun��ment accept��, qui affirme que l'��tat de m��thylation des motifs GANTC change durant le cycle cellulaire de C. crescentus: les motifs GANTC sont compl��tement m��thyl��s (m��thyl��s sur les deux brins) avant de d��but de la r��plication de l'ADN; ils deviennent h��mi-m��thyl��s (m��thyl��s sur un brin seulement) une fois r��pliqu��s, ce qui arrive �� diff��rents moments durant la r��plication pour diff��rents sites le long du chromosome en fonction de leur position par rapport �� l'origine de r��pli-cation; finalement, les motifs GANTC sont rem��thyl��s apr��s la fin de la r��plication du chromosome lorsque la prot��ine CcrM est massivement, mais tr��s transitoirement, produite. Par ailleurs, nous identifions dans le chromosome de C. crescentus environ 30 motifs GANTC qui restent en perma-nence non-m��thyl��s dans une grande partie de la population bact��rienne; ces motifs sont probable-ment prot��g��s de l'action de CcrM par des prot��ines qui s'attachent �� l'ADN et certains d'entre eux pourraient ��tre impliqu��s dans des m��canismes de r��gulation g��n��rant une transcription bistable.

Distributed all-atom simulations of intrinsically unstructured proteins

The Computational Biophysics Group at the Universitat Pompeu Fabra (GRIB-UPF) hosts two unique computational resources dedicated to the execution of large scale molecular dynamics (MD) simulations: (a) the ACMD molecular-dynamics software, used on standard personal computers with graphical processing units (GPUs); and (b) the GPUGRID. net computing network, supported by users distributed worldwide that volunteer GPUs for biomedical research. We leveraged these resources and developed studies, protocols and open-source software to elucidate energetics and pathways of a number of biomolecular systems, with a special focus on flexible proteins with many degrees of freedom. First, we characterized ion permeation through the bactericidal model protein Gramicidin A conducting one of the largest studies to date with the steered MD biasing methodology. Next, we addressed an open problem in structural biology, the determination of drug-protein association kinetics; we reconstructed the binding free energy, association, and dissaciociation rates of a drug like model system through a spatial decomposition and a Makov-chain analysis. The work was published in the Proceedings of the National Academy of Sciences and become one of the few landmark papers elucidating a ligand-binding pathway. Furthermore, we investigated the unstructured Kinase Inducible Domain (KID), a 28-peptide central to signalling and transcriptional response; the kinetics of this challenging system was modelled with a Markovian approach in collaboration with Frank Noe���s group at the Freie University of Berlin. The impact of the funding includes three peer-reviewed publication on high-impact journals; three more papers under review; four MD analysis components, released as open-source software; MD protocols; didactic material, and code for the hosting group.

T-Cell Receptors Binding Orientation over Peptide/MHC Class I Is Driven by Long-Range Interactions.

Crystallographic data about T-Cell Receptor - peptide - major histocompatibility complex class I (TCRpMHC) interaction have revealed extremely diverse TCR binding modes triggering antigen recognition. Understanding the molecular basis that governs TCR orientation over pMHC is still a considerable challenge. We present a simplified rigid approach applied on all non-redundant TCRpMHC crystal structures available. The CHARMM force field in combination with the FACTS implicit solvation model is used to study the role of long-distance interactions between the TCR and pMHC. We demonstrate that the sum of the coulomb interactions and the electrostatic solvation energies is sufficient to identify two orientations corresponding to energetic minima at 0�� and 180�� from the native orientation. Interestingly, these results are shown to be robust upon small structural variations of the TCR such as changes induced by Molecular Dynamics simulations, suggesting that shape complementarity is not required to obtain a reliable signal. Accurate energy minima are also identified by confronting unbound TCR crystal structures to pMHC. Furthermore, we decompose the electrostatic energy into residue contributions to estimate their role in the overall orientation. Results show that most of the driving force leading to the formation of the complex is defined by CDR1,2/MHC interactions. This long-distance contribution appears to be independent from the binding process itself, since it is reliably identified without considering neither short-range energy terms nor CDR induced fit upon binding. Ultimately, we present an attempt to predict the TCR/pMHC binding mode for a TCR structure obtained by homology modeling. The simplicity of the approach and the absence of any fitted parameters make it also easily applicable to other types of macromolecular protein complexes.

Constitutively active mutants of the alpha 1B-adrenergic receptor: role of highly conserved polar amino acids in receptor activation.

Site-directed mutagenesis and molecular dynamics simulations of the alpha 1B-adrenergic receptor (AR) were combined to explore the potential molecular changes correlated with the transition from R (inactive state) to R (active state). Using molecular dynamics analysis we compared the structural/dynamic features of constitutively active mutants with those of the wild type and of an inactive alpha 1B-AR to build a theoretical model which defines the essential features of R and R. The results of site-directed mutagenesis were in striking agreement with the predictions of the model supporting the following hypothesis. (i) The equilibrium between R and R depends on the equilibrium between the deprotonated and protonated forms, respectively, of D142 of the DRY motif. In fact, replacement of D142 with alanine confers high constitutive activity to the alpha 1B-AR. (ii) The shift of R143 of the DRY sequence out of a conserved 'polar pocket' formed by N63, D91, N344 and Y348 is a feature common to all the active structures, suggesting that the role of R143 is fundamental for mediating receptor activation. Disruption of these intramolecular interactions by replacing N63 with alanine constitutively activates the alpha 1B-AR. Our findings might provide interesting generalities about the activation process of G protein-coupled receptors.

Systematic study of the static electrical properties of the CO molecule: influence of the basis set size and correlation energy

The influence of the basis set size and the correlation energy in the static electrical properties of the CO molecule is assessed. In particular, we have studied both the nuclear relaxation and the vibrational contributions to the static molecular electrical properties, the vibrational Stark effect (VSE) and the vibrational intensity effect (VIE). From a mathematical point of view, when a static and uniform electric field is applied to a molecule, the energy of this system can be expressed in terms of a double power series with respect to the bond length and to the field strength. From the power series expansion of the potential energy, field-dependent expressions for the equilibrium geometry, for the potential energy and for the force constant are obtained. The nuclear relaxation and vibrational contributions to the molecular electrical properties are analyzed in terms of the derivatives of the electronic molecular properties. In general, the results presented show that accurate inclusion of the correlation energy and large basis sets are needed to calculate the molecular electrical properties and their derivatives with respect to either nuclear displacements or/and field strength. With respect to experimental data, the calculated power series coefficients are overestimated by the SCF, CISD, and QCISD methods. On the contrary, perturbation methods (MP2 and MP4) tend to underestimate them. In average and using the 6-311 + G(3df) basis set and for the CO molecule, the nuclear relaxation and the vibrational contributions to the molecular electrical properties amount to 11.7%, 3.3%, and 69.7% of the purely electronic ��, ��, and �� values, respectively

An unusual food plant for Cydia pomonella (Linnaeus) (Lepidoptera, Tortricidae) in Mexico

An unusual food plant for Cydia pomonella (Linnaeus) (Lepidoptera, Tortricidae) in Mexico. Larvae of Cydia pomonella (Linnaeus, 1758) were discovered on floral cones of Magnolia schiedeana (Schltdl, 1864) near the natural reserve of La Martinica, Veracruz, M��xico. Magnolia represents an unusual host for this moth species, which is known throughout the world as the "codling moth", a serious pest of fruits of Rosaceae, especially apples. The larvae were identified using taxonomic keys, and identification was corroborated using molecular markers. Further sampling resulted in no additional larvae, hence, the observation was probably that of an ovipositional error by the female, and M. schiedeana is not at risk of attack by this important moth pest.

Enhanced EGFP Fluorescence Emission in Presence of PEG Aqueous Solutions and PIB1000-PEG6000-PIB1000 Copolymer Vesicles.

An EGFP construct interacting with the PIB1000-PEG6000-PIB1000 vesicles surface reported a ~2-fold fluorescence emission enhancement. Because of the constructs nature with the amphiphilic peptide inserted into the PIB core, EGFP is expected to experience a "pure" PEG environment. To unravel this phenomenon PEG/water solutions at different molecular weights and concentrations were used. Already at ~1 : 10 protein/PEG molar ratio the increase in fluorescence emission is observed reaching a plateau correlating with the PEG molecular weight. Parallel experiments in presence of glycerol aqueous solutions did show a slight fluorescence enhancement however starting at much higher concentrations. Molecular dynamics simulations of EGFP in neat water, glycerol, and PEG aqueous solutions were performed showing that PEG molecules tend to "wrap" the protein creating a microenvironment where the local PEG concentration is higher compared to its bulk concentration. Because the fluorescent emission can be perturbed by the refractive index surrounding the protein, the clustering of PEG molecules induces an enhanced fluorescence emission already at extremely low concentrations. These findings can be important when related to the use of EGFP as reported in molecular biology experiments.

Lack of Enantioselectivity in the SULT1A3-catalyzed Sulfoconjugation of Normetanephrine Enantiomers: An In Vitro and Computational Study.

(1R)-Normetanephrine is the natural stereoisomeric substrate for sulfotransferase 1A3 (SULT1A3)-catalyzed sulfonation. Nothing appears known on the enantioselectivity of the reaction despite its potential significance in the metabolism of adrenergic amines and in clinical biochemistry. We confronted the kinetic parameters of the sulfoconjugation of synthetic (1R)-normetanephrine and (1S)-normetanephrine by recombinant human SULT1A3 to a docking model of each normetanephrine enantiomer with SULT1A3 and the 3'-phosphoadenosine-5'-phosphosulfate cofactor on the basis of molecular modeling and molecular dynamics simulations of the stability of the complexes. The K(M) , V(max) , and k(cat) values for the sulfonation of (1R)-normetanephrine, (1S)-normetanephrine, and racemic normetanephrine were similar. In silico models were consistent with these findings as they showed that the binding modes of the two enantiomers were almost identical. In conclusion, SULT1A3 is not substrate-enantioselective toward normetanephrine, an unexpected finding explainable by a mutual adaptability between the ligands and SULT1A3 through an "induced-fit model" in the catalytic pocket. Chirality, 00:000-000, 2012.�� 2012 Wiley Periodicals, Inc.

The role of Cysteine 6.47 in class A GPCRs

Abstract Background: The CWxP motif of transmembrane helix 6 (x: any residue) is highly conserved in class A GPCRs. Within this motif, W6.48 is a big star in the theory of the global ���toggle switch��� because of its key role in the activation mechanism of GPCRs upon ligand binding. With all footlights focused on W6.48, the reason why the preceding residue, C6.47, is largely conserved is still unknown. The present study is aimed to fill up this lack of knowledge by characterizing the role of C6.47 of the CWxP motif. Results: A complete analysis of available crystal structures has been made alongside with molecular dynamics simulations of model peptides to explore a possible structural role for C6.47. Conclusions: We conclude that C6.47 does not modulate the conformation of the TM6 proline kink and propose that C6.47 participates in the rearrangement of the TM6 and TM7 interface accompanying activation.

Unexpected genomic variability in clinical and environmental strains of the pathogenic yeast Candida parapsilosis

Invasive candidiasis is the most commonly reported invasive fungal infection worldwide. Although Candida albicans remains the main cause, the incidence of emerging Candida species, such as C. parapsilosis is increasing. It has been postulated that C. parapsilosis clinical isolates result from a recent global expansion of a virulent clone. However, the availability of a single genome for this species has so far prevented testing this hypothesis at genomic scales. We present here the sequence of three additional strains from clinical and environmental samples. Our analyses reveal unexpected patterns of genomic variation, shared among distant strains, that argue against the clonal expansion hypothesis. All strains carry independent expansions involving an arsenite transporter homolog, pointing to the existence of directional selection in the environment, and independent origins of the two clinical isolates. Furthermore, we report the first evidence for the existence of recombination in this species. Altogether, our results shed new light onto the dynamics of genome evolution in C. parapsilosis.

Analysis of Industrial Granular Flow Applications by Using Advanced Collision Models

Granular flow phenomena are frequently encountered in the design of process and industrial plants in the traditional fields of the chemical, nuclear and oil industries as well as in other activities such as food and materials handling. Multi-phase flow is one important branch of the granular flow. Granular materials have unusual kinds of behavior compared to normal materials, either solids or fluids. Although some of the characteristics are still not well-known yet, one thing is confirmed: the particle-particle interaction plays a key role in the dynamics of granular materials, especially for dense granular materials. At the beginning of this thesis, detailed illustration of developing two models for describing the interaction based on the results of finite-element simulation, dimension analysis and numerical simulation is presented. The first model is used to describing the normal collision of viscoelastic particles. Based on some existent models, more parameters are added to this model, which make the model predict the experimental results more accurately. The second model is used for oblique collision, which include the effects from tangential velocity, angular velocity and surface friction based on Coulomb's law. The theoretical predictions of this model are in agreement with those by finite-element simulation. I n the latter chapters of this thesis, the models are used to predict industrial granular flow and the agreement between the simulations and experiments also shows the validation of the new model. The first case presents the simulation of granular flow passing over a circular obstacle. The simulations successfully predict the existence of a parabolic steady layer and show how the characteristics of the particles, such as coefficients of restitution and surface friction affect the separation results. The second case is a spinning container filled with granular material. Employing the previous models, the simulation could also reproduce experimentally observed phenomena, such as a depression in the center of a high frequency rotation. The third application is about gas-solid mixed flow in a vertically vibrated device. Gas phase motion is added to coherence with the particle motion. The governing equations of the gas phase are solved by using the Large eddy simulation (LES) and particle motion is predicted by using the Lagrangian method. The simulation predicted some pattern formation reported by experiment.

0 AND 02 ADSORPTION ON Ag/Cu[100]

In this work the adsorption mechanisms of atomic and molecular oxygen on Cu(100) surface are studied using ab initio simulation methods. Through the atomistic scale under-standing of the elementary oxidation processes we can further understand the large-scale oxidation. Copper is a material widely used in industry which makes it an interesting subject, and also understanding the oxidation of copper helps us understand the oxidation mechanism of other metals. First we have a look on some theory on surface alloys in general and behaviour of Ag on Cu(100) surface. After that the physical background there is behind the methods of density functional calculations are discussed, and some methods, namely potential energy surfaces and molecular dynamics, are introduced. Then there is a brief look on the numerical details used in the calculations, and after that, the results of the simulations are exhibited.