The diamond pyramid structure in electroless copper deposit, its atomic model and molecular dynamics simulation

In this seminar, I will talk about the discovery of the diamond pyramid structures in the electroless copper deposits on both epoxy and stainless steel substrates. The surface morphology of the structure was characterized with scanning electron microscopy (SEM). According to the morphological feature of the structure, an atom model was brought forward in order to describe the possible mechanism of forming such structure. Molecular dynamics simulations were then carried out to investigate the growing process of the diamond pyramid structure. The final structures of the simulation were compared with the SEM images and the atomic model. The radial distribution function of the final structures of the simulation was compared with that calculated from the X-ray diffraction pattern of the electroless copper deposit sample.

Molecular dynamics simulation investigation on the plastic flow behaviour of silicon during nanometric cutting

Molecular dynamics (MD) simulation was carried out to acquire an in-depth understanding of the flow behaviour of single crystal silicon during nanometric cutting on three principal crystallographic planes and at different cutting temperatures. The key findings were that (i) the substrate material underneath the cutting tool was observed for the first time to experience a rotational flow akin to fluids at all the tested temperatures up to 1200 K. (ii) The degree of flow in terms of vorticity was found higher on the (1 1 1) crystal plane signifying better machinability on this orientation in accord with the current pool of knowledge (iii) an increase in the machining temperature reduces the springback effect and thereby the elastic recovery and (iv) the cutting orientation and the cutting temperature showed significant dependence on the location of the stagnation region in the cutting zone of the substrate.

Self-healing of cracks during ductile regime machining of silicon: Insights from molecular dynamics simulation

During nanoindentation and ductile-regime machining of silicon, a phenomenon known as “self-healing” takes place in that the microcracks, microfractures, and small spallings generated during the machining are filled by the plastically flowing ductile phase of silicon. However, this phenomenon has not been observed in simulation studies. In this work, using a long-range potential function, molecular dynamics simulation was used to provide an improved explanation of this mechanism. A unique phenomenon of brittle cracking was discovered, typically inclined at an angle of 45° to 55° to the cut surface, leading to the formation of periodic arrays of nanogrooves being filled by plastically flowing silicon during cutting. This observation is supported by the direct imaging. The simulated X-ray diffraction analysis proves that in contrast to experiments, Si-I to Si-II (beta tin) transformation during ductile-regime cutting is highly unlikely and solid-state amorphisation of silicon caused solely by the machining stress rather than the cutting temperature is the key to its brittle-ductile transition observed during the MD simulations

Excess Electron Localization in Solvated DNA Bases

We present a first-principles molecular dynamics study of an excess electron in condensed phase models of solvated DNA bases. Calculations on increasingly large microsolvated clusters taken from liquid phase simulations show that adiabatic electron affinities increase systematically upon solvation, as for optimized gas-phase geometries. Dynamical simulations after vertical attachment indicate that the excess electron, which is initially found delocalized, localizes around the nucleobases within a 15 fs time scale. This transition requires small rearrangements in the geometry of the bases.

Ion Association in [bmim][PF6]/Naphthalene Mixtures: An Experimental and Computational Study

Mixtures of room temperature ionic liquids (IL) with neutral organic molecules provide a valuable testing ground to investigate the interplay of the ionic and molecular-dipolar state in dense Coulomb systems at near ambient conditions. In the present study, the viscosity eta and the ionic conductivity a of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6])/naphthalene mixtures at T = 80 degrees C have been measured at 10 stoichiometries spanning the composition range from pure naphthalene to pure [bmim][PF6]. The viscosity grows nearly monotonically with increasing IL mole fraction (x), whereas the conductivity per ion displays a clear peak at x approximate to 15%. The origin of this maximum has been investigated using molecular dynamics simulations based on a classical force field. Snapshots of the simulated samples show that the conductivity maximum is due to the gradual transition in the IL component from an ionic state at high x to a dipolar fluid made of neutral ion pairs at low x. At concentrations x <0.20 the ion pairs condense into molecular-thin filaments bound by dipolar forces and extending in between nanometric droplets of IL. These results are confirmed and complemented by the computation of dynamic and transport properties in [bmim][PF6]/naphthalene mixtures at low IL concentration.

A crystallographic and modelling study of a human telomeric RNA (TERRA) quadruplex

DNA telomeric repeats in mammalian cells are transcribed to guanine-rich RNA sequences, which adopt parallel-stranded G-quadruplexes with a propeller-like fold. The successful crystallization and structure analysis of a bimolecular human telomeric RNA G-quadruplex, folded into the same crystalline environment as an equivalent DNA oligonucleotide sequence, is reported here. The structural basis of the increased stability of RNA telomeric quadruplexes over DNA ones and their preference for parallel topologies is described here. Our findings suggest that the 2'-OH hydroxyl groups in the RNA quadruplex play a significant role in redefining hydration structure in the grooves and the hydrogen bonding networks. The preference for specific nucleotides to populate the C3'-endo sugar pucker domain is accommodated by alterations in the phosphate backbone, which leads to greater stability through enhanced hydrogen bonding networks. Molecular dynamics simulations on the DNA and RNA quadruplexes are consistent with these findings. The computations, based on the native crystal structure, provide an explanation for RNA G-quadruplex ligand binding selectivity for a group of naphthalene diimide ligands as compared to the DNA G-quadruplex.

Molecular Modeling on Inhibitor Complexes and Active-Site Dynamics of Cytochrome P450 C17, a Target for Prostate Cancer Therapy

A molecular model for the P450 enzyme cytochrome P450 C17 (CYP17) is presented based on sequence alignments of multiple template structures and homology modeling. This enzyme plays a central role in the biosynthesis of testosterone and is emerging as a major target in prostate cancer, with the recently developed inhibitor abiraterone currently in advanced clinical trials. The model is described in detail, together with its validation, by providing structural explanations to available site-directed mutagenesis data. The CYP17 molecule in this model is in the form of a triangular prism, with an edge of similar to 55 angstrom and a thickness of similar to 37 angstrom. It is predominantly helical, comprising 13 alpha helices interspersed by six 3(10) helices and 11 beta-sheets. Multinanosecond molecular dynamics simulations in explicit solvent have been carried out, and principal components analysis has been used to reveal the details of dynamics around the active site. Coarse-grained methods have also been used to verify low-frequency motions, which have been correlated with active-site gating. The work also describes the results of docking synthetic inhibitors, including the drug abiraterone and the natural substrate pregnenolone, in the CYP17 active site together with molecular dynamics simulations on the complexes. (C) 2010 Elsevier Ltd. All rights reserved.

A molecular model for drug binding to tandem repeats of telomeric G-quadruplexes.

The extreme 3'-ends of human telomeres consist of 150–250 nucleotides of single-stranded DNA sequence together with associated proteins. Small-molecule ligands can compete with these proteins and induce a conformational change in the DNA to a four-stranded quadruplex arrangement, which is also no longer a substrate for the telomerase enzyme. The modified telomere ends provide signals to the DNA-damage-response system and trigger senescence and apoptosis. Experimental structural data are available on such quadruplex complexes comprising up to four telomeric DNA repeats, but not on longer systems that are more directly relevant to the single-stranded overhang in human cells. The present paper reports on a molecular modelling study that uses Molecular Dynamics simulation methods to build dimer and tetramer quadruplex repeats. These incorporate ligand-binding sites and are models for overhang–ligand complexes.

Comparison of Dynamics of Extracellular Accesses to the β(1) and β(2) Adrenoceptors Binding Sites Uncovers the Potential of Kinetic Basis of Antagonist Selectivity.

From the molecular mechanism of antagonist unbinding in the ß(1) and ß(2) adrenoceptors investigated by steered molecular dynamics, we attempt to provide further possibilities of ligand subtype and subspecies selectivity. We have simulated unbinding of ß(1) -selective Esmolol and ß(2) -selective ICI-118551 from both receptors to the extracellular environment and found distinct molecular features of unbinding. By calculating work profiles, we show different preference in antagonist unbinding pathways between the receptors, in particular, perpendicular to the membrane pathway is favourable in the ß(1) adrenoceptor, whereas the lateral pathway involving helices 5, 6 and 7 is preferable in the ß(2) adrenoceptor. The estimated free energy change of unbinding based on the preferable pathway correlates with the experimental ligand selectivity. We then show that the non-conserved K347 (6.58) appears to facilitate in guiding Esmolol to the extracellular surface via hydrogen bonds in the ß(1) adrenoceptor. In contrast, hydrophobic and aromatic interactions dominate in driving ICI-118551 through the easiest pathway in the ß(2) adrenoceptor. We show how our study can stimulate design of selective antagonists and discuss other possible molecular reasons of ligand selectivity, involving sequential binding of agonists and glycosylation of the receptor extracellular surface. © 2012 John Wiley & Sons A/S.

Conformational Dynamics of the Human Propeller Telomeric DNA Quadruplex on a Microsecond Time Scale

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. While G-quartet stems have been well characterised, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 µs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.

Density Functional study of the photoactive yellow protein's chromophore

We have investigated the structural and electronic properties of p-coumaric acid, the chromophore of the photoactive yellow protein (PYP), by means of first-principles molecular dynamics based on density functional theory (DFT). We have studied the chromophore both in the vacuum and in an extended model which includes the nearest residues in the binding pocket of PYP, as derived from crystallographic data. We have characterized the ground state of the isolated chromophore in its protonated and deprotonated forms and computed the energy barrier involved in the trans to cis isomerization process around the carbon-carbon double bond. A comparison of the optimized structures of the chromophore in the vacuum and in the extended protein model, both in the trans (ground state of PYP in the dark) and cis (first light-activated intermediate) configuration, shows how the protein environment affects the chromophore in the first step of the photocycle. Our model gives an energy storage of 25 kcal/mol associated with the trans-to-cia photoisomerization. Finally, we have elucidated the nature of the electronic excitation relevant for the photochemistry of PYP by means of time-dependent DFT calculations.

Mechanistic Study of 1,3-Butadiene Formation in Acetylene Hydrogenation over the Pd-Based Catalysts Using Density Functional Calculations

Green oil, which leads to the deactivation of the catalysts used for the selective hydrogenation of acetylene, has long been observed but its formation mechanism is not fully understood. In this work, the formation of 1,3-butadiene, known to be the precursor of green oil, on both Pd(111) and Pd(211) surfaces is examined using density functional theory calculations. The pathways containing C-2 + C-2 coupling reactions as well as the corresponding hydrogenation reactions are studied in detail. Three pathways for 1,3-butadiene production, namely coupling plus hydrogenation and further hydrogenation, hydrogenation plus coupling plus hydrogenation, and a two step hydrogenation followed by coupling, are determined. By comparing the effective barriers, we identify the favored pathway on both surfaces. A general understanding toward the deactivation process of the industrial catalysts is also provided. In addition, the effects of the formation of subsurface carbon atoms as well as the Ag alloying on the 1,3-butadiene formation on Pd-based catalysts are also investigated and compared with experimental results.

Influence of surface structures, subsurface carbon and hydrogen, and surface alloying on the activity and selectivity of acetylene hydrogenation on Pd surfaces: A density functional theory study

The selective hydrogenation of acetylene to ethylene on several Pd surfaces (Pd(111), Pd(100), Pd(211), and Pd(211)-defect) and Pd surfaces with subsurface species (carbon and hydrogen) as well as a number of Pd-based alloys (Pd-M/Pd(111) and Pd-M/Pd(211) (M = Cu, Ag and Au)) are investigated using density functional theory calculations to understand both the acetylene hydrogenation activity and the selectivity of ethylene formation. All the hydrogenation barriers are calculated, and the reaction rates on these surfaces are obtained using a two-step model. Pd(211) is found to have the highest activity for acetylene hydrogenation while Pd(100) gives rise to the lowest activity. In addition, more open surfaces result in over-hydrogenation to form ethane, while the close-packed surface (Pd(111)) is the most selective. However, we also find that the presence of subsurface carbon and hydrogen significantly changes the reactivity and selectivity of acetylene toward hydrogenation on Pd surfaces. On forming surface alloys of Pd with Cu, Ag and Au, the selectivity for ethylene is also found to be changed. A new energy decomposition method is used to quantitatively analyze the factors in determining the changes in selectivity. These surface modifiers are found to block low coordination unselective sites, leading to a decreased ethane production. (C) 2013 The Authors. Published by Elsevier Inc. All rights reserved.

Selective Hydrogenation of Acetylene over Pd-Boron Catalysts: A Density Functional Theory Study

Boron-modified Pd catalysts have shown excellent performance for the selective hydrogenation of alkynes experimentally. In the current work, we investigated the hydrogenation of acetylene on boron-modified Pd(111) and Pd(211) surfaces, utilizing density functional theory calculations. The activity of acetylene hydrogenation has been studied by estimating the effective barrier of the whole process. The selectivity of ethylene formation is investigated from a comparison between the desorption and the hydrogenation of ethylene as well as comparison between the ethylene and the 1,3-butadiene formation. Formation of subsurface carbon and hydrogen on both boron-modified Pd(111) and Pd(211) surfaces has also been evaluated, since these have been reported to affect both the activity and the selectivity of acetylene hydrogenation to produce ethylene on Pd surfaces. Our results provide some important insights into the Pd B catalysts for selective hydrogenation of acetylene and also for more complex hydrogenation systems, such as stereoselective hydrogenation of longer chain alkynes and selective hydrogenation of vegetable oil.

Mn ion dissolution from MnS:a density functional theory study

The dissolution of MnS inclusions could induce pitting corrosion in stainless steels, but its dissolution mechanism is poorly understood at the atomic scale. With the help of ab initio molecular dynamics calculations, one inevitable step in the dissolution of MnS is studied by simulating the process of one Mn ion leaving the surface. The reaction mechanism is determined to contain three steps with two large barriers and a small one, leading to two slow steps in the Mn ion dissolution. Comparing to the Na ion dissolution from NaCl, the barriers of the Mn ion dissolution are much larger, which is a reflection of their different electronic structures.