Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

The effect of halogen atom on the molecular quadratic nonlinearity of half sandwich complexes in the presence of an acceptor

Half sandwich complexes of the type [CpM(CO)(n)X] {X=Cl, Br, I; If, M=Fe, Ru; n=2 and if M=Mo; n=3} and [CpNiPPh3X] {X=Cl, Br, I} have been synthesized and their second order molecular nonlinearity (beta) measured at 1064 nm in CHCl3 by the hyper-Rayleigh scattering technique. Iron complexes consistently display larger beta values than ruthenium complexes while nickel complexes have marginally larger beta values than iron complexes. In the presence of an acceptor ligand such as CO or PPh3, the role of the halogen atom is that of a pi donor. The better overlap of Cl orbitals with Fe and Ni metal centres make Cl a better pi donor than Br or I in the respective complexes. Consequently, M-pi interaction is stronger in Fe/Ni-Cl complexes. The value of beta decreases as one goes down the halogen group. For the complexes of 4d metal ions where the metal-ligand distance is larger, the influence of pi orbital overlap appears to be less important, resulting in moderate changes in beta as a function of halogen substitution. (C) 2006 Elsevier B.V. All rights reserved.

Three-dimensional structure of heat shock protein 90 from Plasmodium falciparum: molecular modelling approach to rational drug design against malaria

We have recently implicated heat shock protein 90 from Plasmodium falciparum (PfHsp90) as a potential drug target against malaria. Using inhibitors specific to the nucleotide binding domain of Hsp90, we have shown potent growth inhibitory effects on development of malarial parasite in human erythrocytes. To gain better understanding of the vital role played by PfHsp90 in parasite growth, we have modeled its three dimensional structure using recently described full length structure of yeast Hsp90. Sequence similarity found between PfHsp90 and yeast Hsp90 allowed us to model the core structure with high confidence. The superimposition of the predicted structure with that of the template yeast Hsp90 structure reveals an RMSD of 3.31 angstrom. The N-terminal and middle domains showed the least RMSD (1.76 angstrom) while the more divergent C-terminus showed a greater RMSD (2.84 angstrom) with respect to the template. The structure shows overall conservation of domains involved in nucleotide binding, ATPase activity, co-chaperone binding as well as inter-subunit interactions. Important co-chaperones known to modulate Hsp90 function in other eukaryotes are conserved in malarial parasite as well. An acidic stretch of amino acids found in the linker region, which is uniquely extended in PfHsp90 could not be modeled in this structure suggesting a flexible conformation. Our results provide a basis to compare the overall structure and functional pathways dependent on PfHsp90 in malarial parasite. Further analysis of differences found between human and parasite Hsp90 may make it possible to design inhibitors targeted specifically against malaria.

Effect of substitution on molecular conformation and packing features in a series of aryl substituted ethyl-6-methyl-4-phenyl-2-thioxo-1,2,3,4-tetrahydropyrimidine-5-carboxylates

The supramolecular structures of eight aryl protected ethyl-6-methyl-4-phenyl-2-thioxo-1,2,3,4 tetrahydropyrimidine-5-carboxyl ates were analyzed in order to understand the effect of variations in functional groups on molecular geometry, conformation and packing of molecules in the crystalline lattice. It is observed that the existence of a short intra-molecular C-H center dot center dot center dot pi interaction between the aromatic hydrogen of the aryl ring with the isolated double bond of the six-membered tetrahydropyrimidine ring is a key feature which imparts additional stability to the molecular conformation in the solid state. The compounds pack via the cooperative involvement of both N-H center dot center dot center dot S=C and N-H center dot center dot center dot O=C intermolecular dimers forming a sheet like structure. In addition, weak C-H center dot center dot center dot O and C-H center dot center dot center dot pi intermolecular interactions provide additional stability to the crystal packing.

Diffusion of Pure CH4 and Its Binary Mixture with CO2 in Faujasite NaY: A Combination of Neutron Scattering Experiments and Molecular Dynamics Simulations

The self-diffusion properties of pure CH4 and its binary mixture with CO2 within MY zeolite have been investigated by combining an experimental quasi-elastic neutron scattering (QENS) technique and classical Molecular dynamics simulations. The QENS measurements carried out at 200 K led to an unexpected self-diffusivity profile for Pure CH4 with the presence of a maximum for a loading of 32 CH4/unit cell, which was never observed before for the diffusion of apolar species in azeolite system With large windows. Molecular dynamics simulations were performed using two distinct microscopic models for representing the CH4/NaY interactions. Depending on the model, we are able to fairly reproduce either the magnitude or the profile of the self-diffusivity.Further analysis allowed LIS to provide some molecular insight into the diffusion mechanism in play. The QENS measurements report only a slight decrease of the self-diffusivity of CH4 in the presence of CO2 when the CO2 loading increases. Molecular dynamics simulations successfully capture this experimental trend and suggest a plausible microscopic diffusion mechanism in the case of this binary mixture.

Clinician perspective on molecular profiling of non-small-cell lung cancer

The article by Meric-Bernstam et al1 that was recently published in Journal of Clinical Oncology raises important questions about the clinical application of large-scale genomic testing. We congratulate the authors for this ambitious study, which successfully profiled 2,000 consecutive patients with advanced cancer. The next-generation sequencing (NGS) platform was used for 1,749 of 2,000 patients (87.5%). Of 789 patients with potentially actionable mutations, 83 (11%, or 4% of screened population) were enrolled in a genomically matched clinical study. As the editorial2 accompanying the article by Meric-Bernstam et al1 pointed out, the 4% figure, albeit disappointing, may be an underestimate because cancers such as lung adenocarcinoma and melanoma, for which ≥ 50% of patients have actionable mutations, were under-represented. ...

Molecular line and continuum studies of the early stages of star formation

New stars form in dense interstellar clouds of gas and dust called molecular clouds. The actual sites where the process of star formation takes place are the dense clumps and cores deeply embedded in molecular clouds. The details of the star formation process are complex and not completely understood. Thus, determining the physical and chemical properties of molecular cloud cores is necessary for a better understanding of how stars are formed. Some of the main features of the origin of low-mass stars, like the Sun, are already relatively well-known, though many details of the process are still under debate. The mechanism through which high-mass stars form, on the other hand, is poorly understood. Although it is likely that the formation of high-mass stars shares many properties similar to those of low-mass stars, the very first steps of the evolutionary sequence are unclear. Observational studies of star formation are carried out particularly at infrared, submillimetre, millimetre, and radio wavelengths. Much of our knowledge about the early stages of star formation in our Milky Way galaxy is obtained through molecular spectral line and dust continuum observations. The continuum emission of cold dust is one of the best tracers of the column density of molecular hydrogen, the main constituent of molecular clouds. Consequently, dust continuum observations provide a powerful tool to map large portions across molecular clouds, and to identify the dense star-forming sites within them. Molecular line observations, on the other hand, provide information on the gas kinematics and temperature. Together, these two observational tools provide an efficient way to study the dense interstellar gas and the associated dust that form new stars. The properties of highly obscured young stars can be further examined through radio continuum observations at centimetre wavelengths. For example, radio continuum emission carries useful information on conditions in the protostar+disk interaction region where protostellar jets are launched. In this PhD thesis, we study the physical and chemical properties of dense clumps and cores in both low- and high-mass star-forming regions. The sources are mainly studied in a statistical sense, but also in more detail. In this way, we are able to examine the general characteristics of the early stages of star formation, cloud properties on large scales (such as fragmentation), and some of the initial conditions of the collapse process that leads to the formation of a star. The studies presented in this thesis are mainly based on molecular line and dust continuum observations. These are combined with archival observations at infrared wavelengths in order to study the protostellar content of the cloud cores. In addition, centimetre radio continuum emission from young stellar objects (YSOs; i.e., protostars and pre-main sequence stars) is studied in this thesis to determine their evolutionary stages. The main results of this thesis are as follows: i) filamentary and sheet-like molecular cloud structures, such as infrared dark clouds (IRDCs), are likely to be caused by supersonic turbulence but their fragmentation at the scale of cores could be due to gravo-thermal instability; ii) the core evolution in the Orion B9 star-forming region appears to be dynamic and the role played by slow ambipolar diffusion in the formation and collapse of the cores may not be significant; iii) the study of the R CrA star-forming region suggests that the centimetre radio emission properties of a YSO are likely to change with its evolutionary stage; iv) the IRDC G304.74+01.32 contains candidate high-mass starless cores which may represent the very first steps of high-mass star and star cluster formation; v) SiO outflow signatures are seen in several high-mass star-forming regions which suggest that high-mass stars form in a similar way as their low-mass counterparts, i.e., via disk accretion. The results presented in this thesis provide constraints on the initial conditions and early stages of both low- and high-mass star formation. In particular, this thesis presents several observational results on the early stages of clustered star formation, which is the dominant mode of star formation in our Galaxy.

Understanding Membranes through the Molecular Design of Lipids

Lipids are amphiphilic molecules that are composed of hydrophilic and hydrophobic regions. A typical membranous aggregate (vesicles, water-filled lipid nanospheres) is formed upon the self-organization of lipids in water from a diverse collection of amphiphiles producing a dynamic supramolecular structure that shows phase behavior and ordering as required for specific biological functions. The determination of various physical properties of lipid aggregates is the key to determining structure-function relationships. Over the years, we have designed and synthesized a wide variety of lipid molecular systems for the investigation of their membrane-forming properties and have used them for purposes such as gene delivery and enzyme activation. In this feature article, we focus on our work on various types of lipids including ion-paired amphiphiles, cholesterol-based lipids, aromatic lipids, macrocyclic lipids containing disulfide tethers; cationic dimeric lipids, and so forth. The emphasis is oil experimental design and bottom-line conclusions.

Molecular dynamics studies of nanoparticle impacts

This thesis concerns the dynamics of nanoparticle impacts on solid surfaces. These impacts occur, for instance, in space, where micro- and nanometeoroids hit surfaces of planets, moons, and spacecraft. On Earth, materials are bombarded with nanoparticles in cluster ion beam devices, in order to clean or smooth their surfaces, or to analyse their elemental composition. In both cases, the result depends on the combined effects of countless single impacts. However, the dynamics of single impacts must be understood before the overall effects of nanoparticle radiation can be modelled. In addition to applications, nanoparticle impacts are also important to basic research in the nanoscience field, because the impacts provide an excellent case to test the applicability of atomic-level interaction models to very dynamic conditions. In this thesis, the stopping of nanoparticles in matter is explored using classical molecular dynamics computer simulations. The materials investigated are gold, silicon, and silica. Impacts on silicon through a native oxide layer and formation of complex craters are also simulated. Nanoparticles up to a diameter of 20 nm (315000 atoms) were used as projectiles. The molecular dynamics method and interatomic potentials for silicon and gold are examined in this thesis. It is shown that the displacement cascade expansionmechanism and crater crown formation are very sensitive to the choice of atomic interaction model. However, the best of the current interatomic models can be utilized in nanoparticle impact simulation, if caution is exercised. The stopping of monatomic ions in matter is understood very well nowadays. However, interactions become very complex when several atoms impact on a surface simultaneously and within a short distance, as happens in a nanoparticle impact. A high energy density is deposited in a relatively small volume, which induces ejection of material and formation of a crater. Very high yields of excavated material are observed experimentally. In addition, the yields scale nonlinearly with the cluster size and impact energy at small cluster sizes, whereas in macroscopic hypervelocity impacts, the scaling 2 is linear. The aim of this thesis is to explore the atomistic mechanisms behind the nonlinear scaling at small cluster sizes. It is shown here that the nonlinear scaling of ejected material yield disappears at large impactor sizes because the stopping mechanism of nanoparticles gradually changes to the same mechanism as in macroscopic hypervelocity impacts. The high yields at small impactor size are due to the early escape of energetic atoms from the hot region. In addition, the sputtering yield is shown to depend very much on the spatial initial energy and momentum distributions that the nanoparticle induces in the material in the first phase of the impact. At the later phases, the ejection of material occurs by several mechanisms. The most important mechanism at high energies or at large cluster sizes is atomic cluster ejection from the transient liquid crown that surrounds the crater. The cluster impact dynamics detected in the simulations are in agreement with several recent experimental results. In addition, it is shown that relatively weak impacts can induce modifications on the surface of an amorphous target over a larger area than was previously expected. This is a probable explanation for the formation of the complex crater shapes observed on these surfaces with atomic force microscopy. Clusters that consist of hundreds of thousands of atoms induce long-range modifications in crystalline gold.

Purification, characterization and molecular cloning of a monocot mannose-binding lectin from Remusatia vivipara with nematicidal activity

A mannose-binding lectin (RVL) was purified from the tubers of Remusatia vivipara, a monocot plant by single-step affinity chromatography on asialofetuin-Sepharose 4B. RVL agglutinated only rabbit erythrocytes and was inhibited by mucin, asialomucin, asialofetuin and thyroglobulin. Lectin activity was stable up to 80A degrees C and under wide range of pH (2.0-9.3). SDS-PAGE and gel filtration results showed the lectin is a homotetramer of Mr 49.5 kDa, but MALDI analysis showed two distinct peaks corresponding to subunit mass of 12 kDa and 12.7 kDa. Also the N-terminal sequencing gave two different sequences indicating presence of two polypeptide chains. Cloning of RVL gene indicated posttranslational cleavage of RVL precursor into two mature polypeptides of 116 and 117 amino-acid residues. Dynamic light scattering (DLS) and gel filtration studies together confirmed the homogeneity of the purified lectin and supported RVL as a dimer with Mr 49.5 kDa derived from single polypeptide precursor of 233 amino acids. Purified RVL exerts potent nematicidal activity on Meloidogyne incognita, a root knot nematode. Fluorescent confocal microscopic studies demonstrated the binding of RVL to specific regions of the alimentary-tract and exhibited a potent toxic effect on M. incognita. RVL-mucin complex failed to interact with the gut confirming the receptor mediated lectin interaction. Very high mortality (88%) rate was observed at lectin concentration as low as 30 A mu g/ml, suggesting its potential application in the development of nematode resistant transgenic-crops.

Ion channel formation by zervamicin-IIB. A molecular modelling study.

Zervamicin-IIB (Zrv-IIB) is a 16 residue peptaibol which forms voltage-activated, multiple conductance level channels in planar lipid bilayers. A molecular model of Zrv-IIB channels is presented. The structure of monomeric Zrv-IIB is based upon the crystal structure of Zervamicin-Leu. The helical backbone is kinked by a hydroxyproline residue at position 10. Zrv-IIB channels are modelled as helix bundles of from 4 to 8 parallel helices surrounding a central pore. The monomers are packed with their C-terminal helical segments in close contact, and the bundles are stabilized by hydrogen bonds between glutamine 11 and hydroxyproline 10 of adjacent helices. Interaction energy profiles for movement of three different probes species (K+, Cl- and water) through the central pore are analyzed. The conformations of: (a) the sidechain of glutamine 3; (b) the hydroxyl group of hydroxyproline 10; and (c) the C-terminal hydroxyl group are "optimized" in order to maximize favourable interactions between the channel and the probes, resulting in favourable interaction energy profiles for all three. This suggests that conformational flexibility of polar sidechains enables the channel lining to mimic an aqueous environment.

Charge density based classification of intermolecular interactions in molecular crystals

Evaluation of intermolecular interactions in terms of both experimental and theoretical charge density analyses has produced a unified picture with which to classify strong and weak hydrogen bonds, along with van der Waals interactions, into three regions.

Molecular markers for discriminating Streptococcus pyogenes and S.dysgalactiae subspecies equisimilis

Given the increasing aetiological importance of Streptococcus dysgalactiae subspecies equisimilis in diseases which are primarily attributed to S. pyogenes, molecular markers are essential to distinguish these species and delineate their epidemiology more precisely. Many clinical microbiology laboratories rely on agglutination reactivity and biochemical tests to distinguish them. These methods have limitations which are particularly exacerbated when isolates with mixed properties are encountered. In order to provide additional distinguishing parameters that could be used to unequivocally discriminate these two common pathogens, we assess here three molecular targets: the speB gene, intergenic region upstream of the scpG gene (IRSG) and virPCR. Of these, the former two respectively gave positive and negative results for S. pyogenes, and negative and positive results for S. dysgalactiae subsp. equisimilis. Thus,a concerted use of these nucleic acid-based methods is particularly helpful in epidemiological surveillance to accurately assess the relative contribution of these species to streptococcal infections and diseases.