Ab initio molecular orbital studies of closed shell flavins

Ab initio calculations with a (7 s 3 p) basis set are performed on uracil, lumazine, alloxazine and various isoalloxazines. The results as total energies and charge distributions are discussed in relation to the biochemical behaviour of the flavins. The calculations correctly predict equilibrium situations in the alloxazine-isoalloxazine system and explain the high affinity for nucleophilic addition at N₅ in the flavins. The reduction of flavins and their reoxidation by oxygen are discussed.     

Ab initio simulations of Cu binding sites on the N-terminal region of prion protein

The human prion protein binds Cu²⁺ ions in the octarepeat domain of the N-terminal tail up to full occupancy at pH 7.4. Recent experiments have shown that the HGGG octarepeat subdomain is responsible for holding the metal bound in a square-planar configuration. By using first principle ab initio molecular dynamics simulations of the Car–Parrinello type, the coordination of copper to the binding sites of the prion protein octarepeat region is investigated. Simulations are carried out for a number of structured binding sites. Results for the complexes Cu(HGGGW)(wat), Cu(HGGG), and [Cu(HGGG)]₂ are presented. While the presence of a Trp residue and a water molecule does not seem to affect the nature of the copper coordination, high stability of the bond between copper and the amide nitrogen of deprotonated Gly residues is confirmed in all cases. For the more interesting [Cu(HGGG)]₂ complex, a dynamically entangled arrangement of the two domains with exchange of amide nitrogen bonds between the two copper centers emerges, which is consistent with the short Cu–Cu distance observed in experiments at full copper occupancy.     

Ab initio calculations on the thermodynamic properties of azaborospiropentanes

Following our recent studies of the thermodynamic properties of azaspiropentane and borospiropentane, in consideration of their usefulness as new potential high energy materials, we follow up with ab initio calculations on the thermodynamic properties of azaborospiropentanes. Properties reported in this study include optimized structural parameters, vibrational frequencies, enthalpies of formation, specific enthalpies of combustion, proton affinities, and hydride affinities. Our results indicate that azatriborospiropentane gives off most energy when combusted, as evidenced by its specific enthalpy of combustion of about −52 kJ per gram. Figure Optimized geometry for R-azatriborospiropentane (10) Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ab initio and dielectric studies of succinic acid and maleic acid in 1,4-dioxane

In a continuing effort to understand the hydrogen bond through the study of dielectric and computational conformal studies of dilute solutions, succinic acid and maleic acid are studied in solutions of 1,4-dioxane solvent. Dielectric studies give an account of the net dipole moment of the system under study, which is then compared with the values obtained from conformal analysis. The dielectric measurements were made at 303 K at a frequency of 9.83 GHz using a X-band microwave test bench in order to determine the relaxation times and the dipole moments. The static dielectric permittivity and the high frequency dielectric permittivity were measured using a LCR meter and an Abbe's refractometer, respectively. The results are inspected in comparison with the dipole moment results of ab initio calculations of some of the conformers in gas phase and in liquid phase. Gaussian-03 software package with 6–31G(d) basis set optimisation was used for this purpose. Onsager's reaction field model is used to study the solvation of the dicarboxylic acids in 1,4-dioxane. The results are interpreted in terms of the intermolecular and intramolecular hydrogen bond interactions in the dilute systems.     

Interaction of copper organometallic precursors with barrier layers of Ti, Ta and W and their nitrides: a first-principles molecular dynamics study

Processes for the deposition of copper films on transition metal barrier layers by means CVD using organometallic precursors are often found to lead to poor adhesion characteristics of the grown film. By means of first-principles molecular dynamics simulations, we show that the source of the problem is the strong reactivity of the surfaces toward the precursors, which decompose spontaneously upon contact with the surface leading to contamination of the interface. Our simulations consider Ti, Ta, and W as barrier layers, and Cu(hfac)-(tmvs) as precursor. In contrast, we show that surfaces of these metals properly passivated with nitrogen, in such a way that only N atoms are exposed on the surface, are much less active and do not lead to decomposition of the precursor. We propose this passivation procedure as a practical solution to the adhesion problem. Figure CupraSelect on the WN (100) surface     

Ab initio RNA folding by discrete molecular dynamics: from structure prediction to folding mechanisms

RNA molecules with novel functions have revived interest in the accurate prediction of RNA three-dimensional (3D) structure and folding dynamics. However, existing methods are inefficient in automated 3D structure prediction. Here, we report a robust computational approach for rapid folding of RNA molecules. We develop a simplified RNA model for discrete molecular dynamics (DMD) simulations, incorporating base-pairing and base-stacking interactions. We demonstrate correct folding of 150 structurally diverse RNA sequences. The majority of DMD-predicted 3D structures have <4 A deviations from experimental structures. The secondary structures corresponding to the predicted 3D structures consist of 94% native base-pair interactions. Folding thermodynamics and kinetics of tRNA(Phe), pseudoknots, and mRNA fragments in DMD simulations are in agreement with previous experimental findings. Folding of RNA molecules features transient, non-native conformations, suggesting non-hierarchical RNA folding. Our method allows rapid conformational sampling of RNA folding, with computational time increasing linearly with RNA length. We envision this approach as a promising tool for RNA structural and functional analyses.     

Pair Potentials from ab initio Calculations for use in MD Simulations of Molten Alkali Carbonates

Ab initio calculations at the Hartree-Fock SCF level have been carried out to determine the pair interaction between the alkali ions and the carbonate ion. A distinction has been made between terms in the metal ion - carbonate ion interaction which have different physical origins, such as static coulomb interaction, short-range repulsion and electronic polarization. The additivity of the pair interaction is investigated in 3-body calculations. It is shown that for these 3-body systems pairwise addition of 2-body interactions from which polarization effects have been omitted is superior to pairwise addition of the full Hartree-Fock interactions. A model potential based on these modified interactions has been constructed. Results of MD simulations show that both structural and dynamical properties are well described by these pair potentials.     

The role of flexibility and hydration on the sequence-specific DNA recognition by the Tn916 integrase protein: a molecular dynamics analysis

The N-terminal domain of the Tn916 integrase protein (INT-DBD) is responsible for DNA binding in the process of strand cleavage and joining reactions required for transposition of the Tn916 conjugative transposon. Site-specific association is facilitated by numerous protein-DNA contacts from the face of a three-stranded beta-sheet inserted into the major groove. The protein undergoes a subtle conformational transition and is slightly unfolded in the protein-DNA complex. The conformation of many charged residues is poorly defined by NMR data but mutational studies have indicated that removal of polar side chains decreases binding affinity, while non-polar contacts are malleable. Based on analysis of the binding enthalpy and binding heat capacity, we have reasoned that dehydration of the protein-DNA interface is incomplete. This study presents results from a molecular dynamics investigation of the INT-DBD-DNA complex aimed at a more detailed understanding of the role of conformational dynamics and hydration in site-specific binding. Comparison of simulations (total of 13 ns) of the free protein and of the bound protein conformation (in isolation or DNA-bound) reveals intrinsic flexibility in certain parts of the molecule. Conformational adaptation linked to partial unfolding appears to be induced by protein-DNA contacts. The protein-DNA hydrogen-bonding network is highly dynamic. The simulation identifies protein-DNA interactions that are poorly resolved or only surmised from the NMR ensemble. Single water molecules and water clusters dynamically optimize the complementarity of polar interactions at the 'wet' protein-DNA interface. The simulation results are useful to establish a qualitative link between experimental data on individual residue's contribution to binding affinity and thermodynamic properties of INT-DBD alone and in complex with DNA.     

Ab initio computational study of beta-cellobiose conformers using B3LYP/6-311++G**

The molecular structure of 27 conformers of beta-cellobiose were studied in vacuo through gradient geometry optimization using B3LYP density functionals and the 6-311++G** basis set. The conformationally dependent geometry changes and energies were explored as well as the hydrogen-bonding network. The lowest electronic energy structures found were not those suggested from available crystallographic and NMR solution data, where the glycosidic dihedral angles fall in the region (phi, psi) approximately (40 degrees, -20 degrees ). Rather, 'flipped' conformations in which the dihedral angles are in the range (phi, psi) approximately (180 degrees, 0 degrees ) are energetically more stable by approximately 2.5 kcal/mol over the 'experimentally accepted' structure. Further, when the vibrational free energy, deltaG, obtained from the calculated frequencies, is compared throughout the series, structures with (phi, psi) in the experimentally observed range still have higher free energy ( approximately 2.0 kcal/mol) than 'flipped' forms. The range of bridging dihedral angles of the 'normal' conformers, resulting from the variance in the phi dihedral is larger than that found in the 'flipped' forms. Due to this large flat energy surface for the normal conformations, we surmise that the summation of populations of these conformations will favor the 'normal' conformations, although evidence suggests that polar solvent effects may play the dominant role in providing stability for the 'normal' forms. Even though some empirical studies previously found the 'flipped' conformations to be lowest in energy, these studies have been generally discredited because they were in disagreement with experimental results. Most of the DFT/ab initio conformations reported here have not been reported previously in the ab initio literature, in part because the use of less rigorous theoretical methods, i.e. smaller basis sets, have given results in general agreement with experimental data, that is, they energetically favored the 'normal' forms. These are the first DFT/ab initio calculations at this level of theory, apparently because of the length and difficulty of carrying out optimizations at these high levels.     

Ab initio investigations on the geometric and electronic structures of a diblock molecular diode under the influence of an external bias

Theoretical investigations on the diblock molecular diode, thiophene–thiazole compound, have been carried out at the Hartree–Fock (HF) level by considering the interaction under the external bias. They demonstrate that the electronic structures of this kind of diode molecule are essentially different from those based on the Aviram and Ratner model, in which donor and acceptor π-conjugated segments are separated by an insulating σ-bonded segment, in terms of the energy levels of the frontier molecular orbitals as well as their spatial distributions. The introduction of the external bias modifies both the geometric and electronic structures. In particular, the spatial distributions of the frontier molecular orbitals are also shifted under the external bias. Moreover, all these features show a strong dependence on the polarity of the applied bias due to the build in intrinsic molecular asymmetric structures, which could be used to intuitively interpret the asymmetrical current–voltage behaviours of molecules.     

Introducing a novel method based on the imperialistic competitive algorithm to determine fluorine intermolecular potential from ab initio calculations and calculation of some properties via MD simulations

In this work, the possibility of obtaining an accurate site-site potential model suitable for use in molecular dynamics (MD) simulations of fluorine from ab initio calculations has been explored. The exploration was made on ab initio calculations. To reduce the ab initio pair potentials into a site-site potential, a higher significance was assigned to the configuration which is more stable. For this purpose, the imperialistic competitive algorithm (ICA) was implemented as a powerful optimisation tool. The calculated second virial coefficients were compared to the experimental values to test the quality of the presented intermolecular potential. The relative error for the calculated second virial coefficient ranged from 0.1 to 5.6%. MD simulations were used to evaluate the ability of the proposed intermolecular potential function. The relative error for the MD simulations ranged from 0.5 to 5.2%. The results are in good agreement with experimental data.     

Quantum simulation for muoniated and deuterated methyl radicals in implicit water solvent: combined ab initio path integral molecular dynamics and the polarizable continuum model simulation study

We have combined ab initio path integral molecular dynamics (PIMD) simulation and the polarizable continuum model (PCM) method to efficiently incorporate solvent effects into nuclear quantum fluctuation of molecular systems. Our combined ab initio PIMD–PCM simulation was applied to muoniated and deuterated methyl radical immersed in implicit water solvent to gain information on solvent and isotope effects from one simulation run. We found that solvent effects lead to the bond elongation and a decrease in the magnitude of isotropic hyperfine coupling constants. These are consistent with the trends in conventional static calculations and experiments. In addition, the performance of cavity models (universal force field, united atom specified for Kohn–Sham and these hybrid models) and the conservation of the PIMD–PCM Hamiltonian were accessed. We confirmed that solvent effects on nuclear quantum fluctuation are efficiently computed using our combined simulation of quantum solute in implicit solvent.     

Mechanical features of various silkworm crystalline considering hydration effect via molecular dynamics simulations

Silk materials are receiving significant attention as base materials for various functional nanomaterials and nanodevices, due to its exceptionally high mechanical properties, biocompatibility, and degradable characteristics. Although crystalline silk regions are composed of various repetitive motifs with differing amino acid sequences, how the effect of humidity works differently on each of the motifs and their structural characteristics remains unclear. We report molecular dynamics (MD) simulations on various silkworm fibroins composed of major motifs (i.e. (GAGAGS)_n, (GAGAGA)_n, and (GAGAGY)_n) at varying degrees of hydration, and reveal how each major motifs of silk fibroins change at each degrees of hydration using MD simulations and their structural properties in mechanical perspective via steered molecular dynamics simulations. Our results explain what effects humidity can have on nanoscale materials and devices consisting of crystalline silk materials.     

Ion-specific diffusion rates through transmembrane protein channels. A molecular dynamics study

The migration of different alkali metal cations through a transmembrane model channel is simulated by means of the molecular dynamics technique. The parameters of the model are chosen in close relation to the gramicidin A channel. Coulomb- and van der Waals-type potentials between the ions and flexible carbonyl groups of the pore-forming molecule are used to describe the ion channel interaction. The diffusion properties of the ions are obtained from three-dimensional trajectory calculations. The diffusion rates for the different ions Li+, Na+, K+ and Rb+ are affected not only by the mass of the particles but also very strongly by their size. The latter effect is more pronounced for rigid channels, i.e., for binding vibrational frequencies of the CO groups with v greater than 400 cm-1. In this range the selectivity sequence for the diffusion rates is the inverse of that expected from normal rate theory but agrees with that found in experiments for gramicidin A.     

Ab initio study of phenyl benzoate: structure,conformational analysis,dipole moment,IR and Raman vibrational spectra

The molecular structure (bond distances and angles), conformational properties, dipole moment and vibrational spectroscopic data (vibrational frequencies, IR and Raman intensities) of phenyl benzoate were calculated using Hartree–Fock (HF), density functional (DFT), and second order Møller–Plesset perturbation theory (MP2) with basis sets ranging from 6-31G* to 6-311++G**. The theoretical results are discussed mainly in terms of comparisons with available experimental data. For geometric data, good agreement between theory and experiment is obtained for the MP2, B3LYP and B3PW91 levels with basis sets including diffuse functions. The B3LYP/6-31+G* theory level estimates the shape of the experimental functions for phenyl torsion around the Ph–O and Ph–C bonds well, but reproduces the height of the rotational barriers poorly. The B3LYP/6-31+G* harmonic force constants were scaled by applying the scaled quantum mechanical force field (SQM) technique. The calculated vibrational spectra were interpreted and band assignments were reported. They are in excellent agreement with experimental IR and Raman spectra.Figure Calculated and experimental (GED) potential energy functions for torsional motion of phenyl benzoate relative to the minimum value. a The potential function for torsion about the O₃–C₄ bond. b The potential function for torsion about the C₂–C₁₀ bond.     

Effect of chirality,length and diameter of carbon nanotubes on the adsorption of 20 amino acids: a molecular dynamics simulation study

We carried out molecular dynamics simulations to study the adsorption of all the 20 amino acids (AAs; aromatic, polar and non-polar) on the surface of chiral, zigzag and armchair single-walled carbon nanotubes. The adsorption was occurring in all systems. In the aromatic AAs, the π–π stacking and the semi-hydrogen bond formation cause a strong interaction with the carbon nanotubes (CNTs). We also investigated the chirality, length and diameter dependencies on adsorption energies. We found that all AAs have more tendency to adsorption on the chiral and zigzag CNTs over the armchair. The results show that increasing both the diameter and the length causes the enhancement of the adsorption energy. But, the effect of the length is more than of the diameter. For example, the adsorption energy of Trp on the surface of CNT (4,1), with 2 nm length, is 20.4 kcal/mol. When the length of CNT becomes twice, the adsorption energy increases by 24 ± 0.3%. But by doubling the diameter, the adsorption energy increased only by 9.8 ± 0.25%.     

Molecular dynamics study of the dissociation of an antigen-antibody complex in solution

Preliminary results are reported of a molecular dynamics calculation of free energy variations during the dissociation of an antigen-antibody complex, hen egg-white lysozyme — Fab D1.3, using atomic coordinates determined by the group of R. J. Poljak, and explicit handling of solvent molecules. After equilibration of the complex in solution at 300 K, a dissociation path was generated by a directed dynamics protocol. Then the thermodynamic perturbation method was used for computing the derivative of the free energy of the system with respect to dissociation coordinate, both for the undissociated complex and in two points along the path. 200-ps molecular dynamics simulations were carried out at each of these points. The results obtained are discussed, with special emphasis on the role of interstitial water in the appearance of a hydrophobic activation free energy.Université Paul Sabatier, and URA 505 of C.N.R.S.