Monte carlo simulation for molecular gas dynamics

The dynamics of low-density flows is governed by the Boltzmann equation of the kinetic theory of gases. This is a nonlinear integro-differential equation and, in general, numerical methods must be used to obtain its solution. The present paper, after a brief review of Direct Simulation Monte Carlo (DSMC) methods due to Bird, and Belotserkovskii and Yanitskii, studies the details of theDSMC method of Deshpande for mono as well as multicomponent gases. The present method is a statistical particle-in-cell method and is based upon the Kac-Prigogine master equation which reduces to the Boltzmann equation under the hypothesis of molecular chaos. The proposed Markoff model simulating the collisions uses a Poisson distribution for the number of collisions allowed in cells into which the physical space is divided. The model is then extended to a binary mixture of gases and it is shown that it is necessary to perform the collisions in a certain sequence to obtain unbiased simulation.

XPS evidence for molecular charge-transfer doping of graphene

By employing X-ray photoelectron spectroscopy (XPS), we have been able to establish the occurrence of charge-transfer doping in few-layer graphene covered with electron acceptor (TCNE) and donor (TTF) molecules. We have performed quantitative estimates of the extent of charge transfer in these complexes and elucidated the origin of unusual shifts of their Raman G-bands and explained the differences in the dependence of conductivity on n- and p-doping. The study unravels the cause of the apparent difference between the charge-transfer doping and electrochemical doping. (C) 2010 Elsevier B.V. All rights reserved.

Molecular orbital studies on nicotinic acid, nicotinamide and related compounds

Electronic structures of nicotinic, isonicotinic and 2-picolinic acids and their amides have been investigated, using the variable-? Pariser-Parr-Pople (PPP), iterative extended Hückel and MINDO/2 methods. In addition, PPP and MINDO/2 treatments have also been applied to 3-acetylpyridine and protonated nicotinamide. Based on these calculations, dipole moments, electronic transitions, chemical and biological activity are discussed. Comparison is made with experimental results where available.

Molecular features of colorectal cancer predisposition and progression

Both inherited genetic variations and somatically acquired mutations drive cancer development. The aim of this thesis was to gain insight into the molecular mechanisms underlying colorectal cancer (CRC) predisposition and tumor progression. Whereas one-third of CRC may develop in the context of hereditary predisposition, the known highly penetrant syndromes only explain a small fraction of all cases. Genome-wide association studies have shown that ten common single nucleotide polymorphisms (SNPs) modestly predispose to CRC. Our population-based sample series of around thousand CRC cases and healthy controls was genotyped for these SNPs. Tumors of heterozygous patients were analyzed for allelic imbalance, in an attempt to reveal the role of these SNPs in somatic tumor progression. The risk allele of rs6983267 at 8q24 was favored in the tumors significantly more often than the neutral allele, indicating that this germline variant is somatically selected for. No imbalance targeting the risk allele was observed in the remaining loci, suggesting that most of the low-penetrance CRC SNPs mainly play a role in the early stages of the neoplastic process. The ten SNPs were further analyzed in 788 CRC cases, 97 of which had a family history of CRC, to evaluate their combined contribution. A significant association appeared between the overall number of risk alleles and familial CRC and these ten SNPs seem to explain around 9% of the familial clustering of CRC. Finding more CRC susceptibility alleles may facilitate individualized risk prediction and cancer prevention in the future. Microsatellite instability (MSI), resulting from defective mismatch repair function, is a hallmark of Lynch syndrome and observed in a subset of all CRCs. Our aim was to identify microsatellite frameshift mutations that inactivate tumor suppressor genes in MSI CRCs. By sequencing microsatellite repeats of underexpressed genes we found six novel MSI target genes that were frequently mutated in 100 MSI CRCs: 51% in GLYR1, 47% in ABCC5, 43% in WDTC1, 33% in ROCK1, 30% in OR51E2, and 28% in TCEB3. Immunohistochemical staining of GLYR1 revealed defective protein expression in homozygously mutated tumors, providing further support for the loss of function hypothesis. Another mutation screening effort sought to identify MSI target genes with putative oncogenic functions. Microsatellites were similarly sequenced in genes that were overexpressed and, upon mutation, predicted to avoid nonsense-mediated mRNA decay. The mitotic checkpoint kinase TTK harbored protein-elongating mutations in 59% of MSI CRCs and the mutant protein was detected in heterozygous MSI CRC cells. No checkpoint dysregulation or defective protein localization was observable however, and the biological relevance of this mutation may hence be related to other mechanisms. In conclusion, these two large-scale and unbiased efforts identified frequently mutated genes that are likely to contribute to the development of this cancer type and may be utilized in developing diagnostic and therapeutic applications.

Conformations of three heterocyclic perhydropyrrolobenzofurans and polymeric assembly via co-operative inter-molecular C-H center dot center dot center dot O hydrogen bonds

In 1-cyclo-hexyl-6,6,8a-trimethyl-3a,6,7,8a-tetra-hydro-1H-1-benzofuro2, 3-b]pyrrole-2,4(3H,5H)-dione, C19H27NO3, (I), and the isomorphous compounds 6,6,8a-trimethyl-1-phenyl-3a,6,7,8a-tetra-hydro-1H-1-benzofuro2,3-b]p yrrole-2,4(3H,5H)-dione, C19H21NO3, (II), and 6,6,8a-trimethyl-1-(3-pyridyl)-3a,6,7,8a-tetra-hydro-1H-1-benzofuro2, 3-b]pyrrole-2,4(3H,5H)-dione, C18H20N2O3, (III), the tetra-hydro-benzo-dihydro-furo-pyrrolidine ring systems are folded at the cis junction of the five-membered rings, giving rise to a non-planar shape of the tricyclic cores. The dihydro-furan and pyrrolidine rings in (I) are puckered and adopt an envelope conformation. The cyclo-hexene rings adopt a half-chair conformation in all the mol-ecules, while the substituent N-cyclo-hexyl ring in (I) assumes a chair form. Short intra-molecular C-HcO contacts form S(5) and S(6) motifs. The isomorphous compounds (II) and (III) are effectively isostructural, and aggregate into chains via inter-molecular C-HcO hydrogen bonds.

Thermal and photochemical cycloaddition reactions of thiocarbonyls: a qualitative molecular orbital analysis

A comprehensive analysis of thermal and photochemical reactions of thiocarbonyls has been undertaken within the PMO framework employing MINDO/3 orbital energies and wavefunctions. The model is generally successful in rationalizing the observed regiochemistry of such reactions. In particular, the indicated regiochemistry for [4 + 2] thermal cycloadditions of saturated thiones to 2-substituted dienes, for the dimerization of α,β-unsaturated thiones, and for the photochemical cycloadditions of thioketones and thioenones are all in agreement with experimental observations. Interesting predictions are also made concerning cycloadditions of saturated, conjugated, and arylalkyl thiones which have not yet been studied experimentally. The analysis reveals the decisive role played by secondary orbital interactions in determining the observed product selectivity in the photochemical reactions between thioenone and olefins.

X-ray studies on crystalline complexes involving amino acids and peptides. XV. Crystal structures of L-lysine D-glutamate and L-lysine D-aspartate monohydrate and the effect of chirality on molecular aggregation

L-Lysine D-glutamate crystallizes in the monoclinic space group P2(1) with a = 4.902, b = 30.719, c = 9.679 A, beta = 90 degrees and Z = 4. The crystals of L-lysine D-aspartate monohydrate belong to the orthorhombic space group P2(1)2(1)2(1) with a = 5.458, b = 7.152, c = 36.022 A and Z = 4. The structures were solved by the direct methods and refined to R values of 0.125 and 0.040 respectively for 1412 and 1503 observed reflections. The glutamate complex is highly pseudosymmetric. The lysine molecules in it assume a conformation with the side chain staggered between the alpha-amino and the alpha-carboxylate groups. The interactions of the side chain amino groups of lysine in the two complexes are such that they form infinite sequences containing alternating amino and carboxylate groups. The molecular aggregation in the glutamate complex is very similar to that observed in L-arginine D-aspartate and L-arginine D-glutamate trihydrate, with the formation of double layers consisting of both types of molecules. In contrast to the situation in the other three LD complexes, the unlike molecules in L-lysine D-aspartate monohydrate aggregate into alternating layers as in the case of most LL complexes. The arrangement of molecules in the lysine layer is nearly the same as in L-lysine L-aspartate, with head-to-tail sequences as the central feature. The arrangement of aspartate ions in the layers containing them is, however, somewhat unusual. Thus the comparison between the LL and the LD complexes analyzed so far indicates that the reversal of chirality of one of the components in a complex leads to profound changes in molecular aggregation, but these changes could be of more than one type.

Molecular phylogeny of Hemidactylus geckos (Squamata: Gekkonidae) of the Indian subcontinent reveals a unique Indian radiation and an Indian origin of Asian house geckos

Represented by approximately 85 species, Hemidactylus is one of the most diverse and widely distributed genera of reptiles in the world. In the Indian subcontinent, this genus is represented by 28 species out of which at least 13 are endemic to this region. Here, we report the phylogeny of the Indian Hemidactylus geckos based on mitochondrial and nuclear DNA markers sequenced from multiple individuals of widely distributed as well as endemic congeners of India. Results indicate that a majority of the species distributed in India form a distinct clade whose members are largely confined to the Indian subcontinent thus representing a unique Indian radiation. The remaining Hemidactylus geckos of India belong to two other geographical clades representing the Southeast Asian and West-Asian arid zone species. Additionally, the three widely distributed, commensal species (H. brookii, H. frenatus and H. flaviviridis) are nested within the Indian radiation suggesting their Indian origin. Dispersal-vicariance analysis also supports their Indian origin and subsequent dispersal out-of-India into West-Asian arid zone and Southeast Asia. Thus, Indian subcontinent has served as an important arena for diversification amongst the Hemidactylus geckos and in the evolution and spread of its commensal geckos. (C) 2010 Elsevier Inc. All rights reserved.

Droplet Epitaxy of InN Quantum Dots on Si(111) by RF Plasma-Assisted Molecular Beam Epitaxy

InN quantum dots (QDs) were fabricated on Si(111) substrate by droplet epitaxy using an RF plasma-assisted MBE system. Variation of the growth parameters, such as growth temperature and deposition time, allowed us to control the characteristic size and density of the QDs. As the growth temperature was increased from 100 C to 300 degrees C, an enlargement of QD size and a drop in dot density were observed, which was led by the limitation of surface diffusion of adatoms with the limited thermal energy. Atomic force microscopy (AFM) and scanning electron microscopy (SEM) were used to assess the QDs size and density. The chemical bonding configurations of InN QDs were examined by X-ray photo-electron spectroscopy (XPS). Fourier transform infrared (FTIR) spectrum of the deposited InN QDs shows the presence of In-N bond. Temperature-dependent photoluminescence (PL) measurements showed that the emission peak energies of the InN QDs are sensitive to temperature and show a strong peak emission at 0.79 eV.

Association between molecular variation in the phosphoglucose isomerase (Pgi) locus and migration propensity in the Glanville fritillary butterfly

Habitat fragmentation produces patches of suitable habitat surrounded by unfavourable matrix habitat. A species may persist in such a fragmented landscape in an equilibrium between the extinctions and recolonizations of local populations, thus forming a metapopulation. Migration between local populations is necessary for the long-term persistence of a metapopulation. The Glanville fritillary butterfly (Melitaea cinxia) forms a metapopulation in the Åland islands in Finland. There is migration between the populations, the extent of which is affected by several environmental factors and variation in the phenotype of individual butterflies. Different allelic forms of the glycolytic enzyme phosphoglucose isomerase (Pgi) has been identified as a possible genetic factor influencing flight performance and migration rate in this species. The frequency of a certain Pgi allele, Pgi-f, follows the same pattern in relation to population age and connectivity as migration propensity. Furthermore, variation in flight metabolic performance, which is likely to affect migration propensity, has been linked to genetic variation in Pgi or a closely linked locus. The aim of this study was to investigate the association between Pgi genotype and the migration propensity in the Glanville fritillary both at the individual and population levels using a statistical modelling approach. A mark-release-recapture (MRR) study was conducted in a habitat patch network of M. cinxia in Åland to collect data on the movements of individual butterflies. Larval samples from the study area were also collected for population level examinations. Each butterfly and larva was genotyped at the Pgi locus. The MRR data was parameterised with two mathematical models of migration: the Virtual Migration Model (VM) and the spatially explicit diffusion model. VM model predicted and observed numbers of emigrants from populations with high and low frequencies of Pgi-f were compared. Posterior predictive data sets were simulated based on the parameters of the diffusion model. Lack-of-fit of observed values to the model predicted values of several descriptors of movements were detected, and the effect of Pgi genotype on the deviations was assessed by randomizations including the genotype information. This study revealed a possible difference in the effect of Pgi genotype on migration propensity between the two sexes in the Glanville fritillary. The females with and males without the Pgi-f allele moved more between habitat patches, which is probably related to differences in the function of flight in the two sexes. Females may use their high flight capacity to migrate between habitat patches to find suitable oviposition sites, whereas males may use it to acquire mates by keeping a territory and fighting off other intruding males, possibly causing them to emigrate. The results were consistent across different movement descriptors and at the individual and population levels. The effect of Pgi is likely to be dependent on the structure of the landscape and the prevailing environmental conditions.

Complexation of siRNA with Dendrimer: A Molecular Modeling Approach

This paper reports the structural behavior and thermodynamics of the complexation of siRNA with poly(amidoamine) (PAMAM) dendrimers of generation 3 (G3) and 4 (G4) through fully atomistic molecular dynamics (MD) simulations accompanied by free energy calculations and inherent structure determination. We have also done simulation with one siRNA and two dendrimers (2 x G3 or 2xG4) to get the microscopic picture of various binding modes. Our simulation results reveal the formation of stable siRNA-dendrimer complex over nanosecond time scale. With the increase in dendrimcr generation, the charge ratio increases and hence the binding energy between siRNA and dendrimer also increases in accordance with available experimental measurements. Calculated radial distribution functions of amines groups of various subgenerations in a given generation of dendrimer and phosphate in backbone of siRNA reveals that one dendrimer of generation 4 shows better binding with siRNA almost wrapping the dendrimer when compared to the binding with lower generation dendrimer like G3. In contrast, two dendrimers of generation 4 show binding without siRNA wrapping the den-rimer because of repulsion between two dendrimers. The counterion distribution around the complex and the water molecules in the hydration shell of siRNA give microscopic picture of the binding dynamics. We see a clear correlation between water. counterions motions and the complexation i.e. the water molecules and counterions which condensed around siRNA are moved away from the siRNA backbone when dendrimer start binding to the siRNA back hone. As siRNA wraps/bind to the dendrimer counterions originally condensed onto siRNA (Na-1) and dendrimer (Cl-) get released. We give a quantitative estimate of the entropy of counterions and show that there is gain in entropy due to counterions release during the complexation. Furthermore, the free energy of complexation of IG3 and IG4 at two different salt concentrations shows that increase in salt concentration leads to the weakening of the binding affinity of siRNA and dendrimer.