Beyond Ehrenfest: correlated non-adiabatic molecular dynamics

A method for introducing correlations between electrons and ions that is computationally affordable is described. The central assumption is that the ionic wavefunctions are narrow, which makes possible a moment expansion for the full density matrix. To make the problem tractable we reduce the remaining many-electron problem to a single-electron problem by performing a trace over all electronic degrees of freedom except one. This introduces both one- and two-electron quantities into the equations of motion. Quantities depending on more than one electron are removed by making a Hartree-Fock approximation. Using the first-moment approximation, we perform a number of tight binding simulations of the effect of an electric current on a mobile atom. The classical contribution to the ionic kinetic energy exhibits cooling and is independent of the bias. The quantum contribution exhibits strong heating, with the heating rate proportional to the bias. However, increased scattering of electrons with increasing ionic kinetic energy is not observed. This effect requires the introduction of the second moment.

Ultrafast coherent control and quantum encoding of molecular vibrations in D2+ using intense laser pulses

Intense, few-femtosecond pulse technology has enabled studies of the fastest vibrational relaxation processes. The hydrogen group vibrations can be imaged and manipulated using intense infrared pulses. Through numerical simulation, we demonstrate an example of ultrafast coherent control that could be effected with current experimental facilities, and observed using high-resolution time-of-flight spectroscopy. The proposal is a pump-probe-type technique to manipulate the D2+ ion with ultrashort pulse sequences. The simulations presented show that vibrational selection can be achieved through pulse delay. We find that the vibrational system can be purified to a two-level system thus realizing a vibrational qubit. A novel scheme for the selective transfer of population between these two levels, based on a Raman process and conditioned upon the delay time of a second control-pulse is outlined, and may enable quantum encoding with this system.

Intense-field dissociation dynamics of D2+ molecular ions using ultrafast laser pulses

The dynamics of dissociation of pre-ionized D2+ molecules using intense (10^12–10^15 W cm-2), ultrashort (50 fs), infrared (? = 790 nm) laser pulses are examined. Use of an intensity selective scan technique has allowed the deuterium energy spectrum to be measured over a broad range of intensity. It is found that the dominant emission shifts to lower energies as intensity is increased, in good agreement with corresponding wavepacket simulations. The results are consistent with an interpretation in terms of bond softening, which at high intensity (approximately >3 × 10^14 W cm-2) becomes dominated by dissociative ionization. Angular distribution measurements reveal the presence of slow molecular dissociation, an indication that vibrational trapping mechanisms occur in this molecule.

Quantum revivals in ultrashort intense field dissociation of molecular ions

Ultrashort (<15 fs) high intensity (1014-1016 W cm-2) laser pulses have provided novel methods for investigation of the dynamics of simple molecular ions such as H2+ and D2+. In this paper we report on simulations carried out for the D2+ molecular ion, within the Born- Oppenheimer and two-state approximations. These simulations allow one to investigate the dissociation dynamics of the D2+ molecular ion when subjected to such ultrashort, intense laser pulses. In particular, these simulations are compared to the results from recent pump-probe experiments, in which, the nuclear vibrational motion of D2+ has been imaged. Simulations suggest that the nature of the dissociation process, be it 1- or 2-photon, may be influenced by the tuning of the pump-probe delay time.

Pathways to state-selective control of vibrational wavepackets in the deuterium molecular ion

The capability of intense ultrashort laser pulses to initiate, control and image vibrational wavepacket dynamics in the deuterium molecular ion has been simulated with a view to inform and direct future femtosecond pump-control-probe experiments. The intense-field coherent control of the vibrational superposition has been studied as a function of pulse intensity and delay time, to provide an indication of key constraints for experimental studies. For selected cases of the control mechanism, probing of the subsequent vibrational wavepacket dynamics has been simulated via the photodissociation (PD) channel. Such PD probing is shown to elucidate the modified wavepacket dynamics where the position of the quantum revival is sensitive to the control process. Through Fourier transform analysis the PD yield is also shown to provide a characterisation of the vibrational distribution. It has been shown that a simple 'critical R cut-off' approximation can be used to reproduce the effect of a probe pulse interaction, providing a convenient and efficient alternative to intensive computer simulations of the PD mechanism in the deuterium molecular ion.

SO and CS observations of molecular clouds - II. Analysis and modelling of the abundance ratios - probing O-2/CO with SO/CS?

We here analyse the observational SO and CS data presented in Nilsson ct al. (2000). The SO/CS integrated intensity ratio maps are presented for 19 molecular clouds, together with tables of relevant ratios at strategic positions, where we have also observed (SO)-S-34 and/or (CS)-S-34. The SO/CS abundance ratio as calculated from an LTE analysis is highly varying within and between the sources. Our isotopomer observations and Monte Carlo simulations verify that this is not an artifact due to optical depth problems. The variation of the maximum SO/CS abundance ratio between the clouds is 0.2-7. The largest variations within a cloud are found for the most nearby objects, possibly indicating resolution effects. We have also performed time dependent chemical simulations. We compare the simulations with our observed SO/CS abundance ratios and suggest a varying oxygen to carbon initial abundance, differing temporal evolution, density differences and X-ray sources associated with young stellar objects as possible explanations to the variations. In particular, the observed variation of the maximum SO/CS abundance ratio between the clouds can be explained by using initial O/C+ abundance ratios in the range 1.3-2.5. We finally derive a relationship between the SO/CS and O-2/CO abundance ratios, which may be used as a guide to find the most promising interstellar O-2 search targets.

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.

Neutron diffraction, NMR and molecular dynamics study of glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate

beta-D-glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate in a 6 : 1 molar ratio (ionic liquid : glucose) has been studied by neutron scattering, NMR and molecular dynamics simulations. Good agreement was found between simulated neutron scattering profiles generated for isotopically substituted liquid systems and those experimentally determined as well as between simulated and experimental diffusion coefficients obtained by Pulsed Field Gradient NMR spectroscopy. The overriding glucose-ionic liquid interactions in the liquid are hydrogen-bonding between acetate oxygens and sugar hydroxyl groups. The ionic liquid cation was found to play only a minor role in the solvation of the sugar and does not participate in hydrogen-bonding with the sugar to any significant degree. NOESY experiments lend further evidence that there is no direct interaction between sugar hydroxyl groups and acidic hydrogens on the ionic liquid cation.

Binding modes of diketo-acid inhibitors of HIV-1 integrase: A comparative molecular dynamics simulation study

HIV-1 integrase (IN) has become an attractive target since drug resistance against HIV-1 reverse transcriptase (RT) and protease (PR) has appeared. Diketo acid (DKA) inhibitors are potent and selective inhibitors of HIV-1 IN: however the action mechanism is not well understood. Here, to study the inhibition mechanism of DKAs we performed 10 ns comparative molecular dynamics simulations on HIV-1 IN bound with three most representative DMA inhibitors: Shionogi inhibitor, S-1360 and two Merck inhibitors L-731,988 and L-708,906. Our simulations show that the acidic part of S-1360 formed salt bridge and cation-pi interactions with Lys159. In addition, the catalytic Glu152 in S-1360 was pushed away from the active site to form an ion-pair interaction with Arg199. The Merck inhibitors can maintain either one or both of these ion-pair interaction features. The difference in potencies of the DMA inhibitors is thus attributed to the different binding modes at the catalytic site. Such structural information at atomic level, not only demonstrates the action modes of DMA inhibitors but also provides a novel starting point for structural-based design of HIV-1 IN inhibitors.

The Plasmodium falciparum Malaria M1 Alanyl Aminopeptidase (PfA-M1): Insights of Catalytic Mechanism and Function from MD Simulations

Malaria caused by several species of Plasmodium is major parasitic disease of humans, causing 1-3 million deaths worldwide annually. The widespread resistance of the human parasite to current drug therapies is of major concern making the identification of new drug targets urgent. While the parasite grows and multiplies inside the host erythrocyte it degrades the host cell hemoglobin and utilizes the released amino acids to synthesize its own proteins. The P. falciparum malarial M1 alanyl-aminopeptidase (PfA-M1) is an enzyme involved in the terminal stages of hemoglobin digestion and the generation of an amino acid pool within the parasite. The enzyme has been validated as a potential drug target since inhibitors of the enzyme block parasite growth in vitro and in vivo. In order to gain further understanding of this enzyme, molecular dynamics simulations using data from a recent crystal structure of PfA-M1 were performed. The results elucidate the pentahedral coordination of the catalytic Zn in these metallo-proteases and provide new insights into the roles of this cation and important active site residues in ligand binding and in the hydrolysis of the peptide bond. Based on the data, we propose a two-step catalytic mechanism, in which the conformation of the active site is altered between the Michaelis complex and the transition state. In addition, the simulations identify global changes in the protein in which conformational transitions in the catalytic domain are transmitted at the opening of the N-terminal 8 angstrom-long channel and at the opening of the 30 angstrom-long C-terminal internal chamber that facilitates entry of peptides to the active site and exit of released amino acids. The possible implications of these global changes with regard to enzyme function are discussed.

Molecular dynamics studies of the STAT3 homodimer:DNA complex: relationships between STAT3 mutations and protein-DNA recognition.

Signal Transducers and Activators of Transcription (STAT) proteins are a group of latent cytoplasmic transcription factors involved in cytokine signaling. STAT3 is a member of the STAT family and is expressed at elevated levels in a large number of diverse human cancers and is now a validated target for anticancer drug discovery.. Understanding the dynamics of the STAT3 dimer interface, accounting for both protein-DNA and protein-protein interactions, with respect to the dynamics of the latent unphosphorylated STAT3 monomer, is important for designing potential small-molecule inhibitors of the activated dimer. Molecular dynamics (MD) simulations have been used to study the activated STAT3 homodimer:DNA complex and the latent unphosphorylated STAT3 monomer in an explicit water environment. Analysis of the data obtained from MD simulations over a 50 ns time frame has suggested how the transcription factor interacts with DNA, the nature of the conformational changes, and ways in which function may be affected. Examination of the dimer interface, focusing on the protein-DNA interactions, including involvement of water molecules, has revealed the key residues contributing to the recognition events involved in STAT3 protein-DNA interactions. This has shown that the majority of mutations in the DNA-binding domain are found at the protein-DNA interface. These mutations have been mapped in detail and related to specific protein-DNA contacts. Their structural stability is described, together with an analysis of the model as a starting-point for the discovery of novel small-molecule STAT3 inhibitors.

The molecular dynamics of Trypanosoma brucei UDP-galactose 4'-epimerase:a drug target for African sleeping sickness

During the past century, several epidemics of human African trypanosomiasis, a deadly disease caused by the protist Trypanosoma brucei, have afflicted sub-Saharan Africa. Over 10 000 new victims are reported each year, with hundreds of thousands more at risk. As current drug treatments are either highly toxic or ineffective, novel trypanocides are urgently needed. The T. brucei galactose synthesis pathway is one potential therapeutic target. Although galactose is essential for T. brucei survival, the parasite lacks the transporters required to intake galactose from the environment. UDP-galactose 4'-epimerase (TbGalE) is responsible for the epimerization of UDP-glucose to UDP-galactose and is therefore of great interest to medicinal chemists. Using molecular dynamics simulations, we investigate the atomistic motions of TbGalE in both the apo and holo states. The sampled conformations and protein dynamics depend not only on the presence of a UDP-sugar ligand, but also on the chirality of the UDP-sugar C4 atom. This dependence provides important insights into TbGalE function and may help guide future computer-aided drug discovery efforts targeting this protein.

A Neutron Diffraction and Molecular Dynamics Investigation of Acetate-based Ionic Liquids as Solvents for Glucose

The liquid state structure of the ionic liquid, 1-ethyl-3-methylimidazolium acetate, and the solute/solvent structure of glucose dissolved in the ionic liquid at a 1: 6 molar ratio have been investigated at 323 K by molecular dynamics simulations and neutron diffraction experiments using H/D isotopically substituted materials. Interactions between hydrogen-bond donating cation sites and polar, directional hydrogen-bond accepting acetate anions are examined. Ion-ion radial distribution functions for the neat ionic liquid, calculated from both MD and derived from the empirical potential structure refinement model to the experimental data, show the alternating shell-structure of anions around the cation, as anticipated. Spatial probability distributions reveal the main anion-to-cation features as in-plane interactions of anions with imidazolium ring hydrogens and cation-cation planar stacking. Interestingly, the presence of the polarised hydrogen-bond acceptor anion leads to increased anion-anion tail-tail structuring within each anion shell, indicating the onset of hydrophobic regions within the anion regions of the liquid.

Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.