Free energy and molecular dynamics calculations for the cubic-tetragonal phase transition in zirconia

The high-temperature cubic-tetragonal phase transition of pure stoichiometric zirconia is studied by molecular dynamics (MD) simulations and within the framework of the Landau theory of phase transformations. The interatomic forces are calculated using an empirical, self-consistent, orthogonal tight-binding model, which includes atomic polarizabilities up to the quadrupolar level. A first set of standard MD calculations shows that, on increasing temperature, one particular vibrational frequency softens. The temperature evolution of the free-energy surfaces around the phase transition is then studied with a second set of calculations. These combine the thermodynamic integration technique with constrained MD simulations. The results seem to support the thesis of a second-order phase transition but with unusual, very anharmonic behavior above the transition temperature.

Time-dependent tight binding

Starting from a Lagrangian mean-field theory, a set of time-dependent tight-binding equations is derived to describe dynamically and self-consistently an interacting system of quantum electrons and classical nuclei. These equations conserve norm, total energy and total momentum. A comparison with other tight-binding models is made. A previous tight-binding result for forces on atoms in the presence of electrical current flow is generalized to the time-dependent domain and is taken beyond the limit of local charge neutrality.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

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.

Locating binding poses in protein-ligand systems using reconnaissance metadynamics

A molecular dynamics-based protocol is proposed for finding and scoring protein-ligand binding poses. This protocol uses the recently developed reconnaissance metadynamics method, which employs a self-learning algorithm to construct a bias that pushes the system away from the kinetic traps where it would otherwise remain. The exploration of phase space with this algorithm is shown to be roughly six to eight times faster than unbiased molecular dynamics and is only limited by the time taken to diffuse about the surface of the protein. We apply this method to the well-studied trypsin-benzamidine system and show that we are able to refind all the poses obtained from a reference EADock blind docking calculation. These poses can be scored based on the length of time the system remains trapped in the pose. Alternatively, one can perform dimensionality reduction on the output trajectory and obtain a map of phase space that can be used in more expensive free-energy calculations.

Binding of Calcium and Carbonate to Polyacrylates

Polyacrylate molecules can be used to slow the growth of calcium carbonate. However, little is known about the mechanism by which the molecules impede the growth rate. A recent computational study (Bulo et al. Macromolecules 2007, 40, 3437) used metadynamics to investigate the binding of calcium to polyacrylate chains and has thrown some light on the coiling and precipitation of these polymers. We extend these simulations to examine the binding of calcium and carbonate to polyacrylate chains. We show that calcium complexed with both carbonate and polyacrylate is a very stable species. The free energies of calcium-carbonate-polyacrylate complexes, with different polymer configurations, are calculated, and differences in the free energy of the binding of carbonate are shown to be due to differences in the amount of steric hindrance about the calcium, which prevents the approach of the carbonate ion.

Organic synthesis in the interstellar medium by low-energy carbon irradiation

We present a first principles molecular dynamics (FPMD) study of the interaction of low-energy neutral carbon projectiles with amorphous solid water clusters at 30 K. Reactions involving the carbon atom at an initial energy of 11 and 1.7 eV with 30-molecule clusters have been investigated. Simulations indicate that the formation of hydroxymethylene, an intermediate in formaldehyde production, dominates at the higher energy. The reaction proceeds by fragmenting a water molecule, binding the carbon to the OH radical, and saturating the C valence with a hydrogen atom that can arise from the originally dissociated water molecule, or through a chain of proton transfer events. We identified several possible pathways for the formation of HCOH. When the initial collision occurs at the periphery of the cluster, we observe the formation of CO and the evaporation of water molecules. At the lower energy water fragmentation is not favorable, thus leading to the formation of weakly bound carbon-water complexes. © 2013 American Chemical Society.

Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice-leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point. 

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

An Energy Descriptor To Quantify Methane Selectivity in Fischer-Tropsch Synthesis:A Density Functional Theory Study

Selectivity is a fundamental issue in heterogeneous catalysis. In this study, the CH(4) selectivity in Fischer-Tropsch synthesis is chosen to be investigated: CH4 selectivity on Rh, Co, Ru, Fe, and Re surfaces is computed by first-principles methods. In conjunction with kinetic analyses, we are able to derive the effective barrier difference between methane formation and chain growth (Delta E(eff)) to quantify the CH(4) selectivity. By using this energy descriptor, the ranking of methane selectivity predicted from density functional theory (DFT) calculations is consistent with experimental work. Moreover, a linear correlation between Delta E(eff) and the chemisorption energy of C + 4H (Delta H) is found. This fundamental finding possesses the following significance: (i) it shows that the selectivity, which appears to have kinetic characteristics, is largely determined by thermodynamic properties; and (ii) it suggests that an increase of the binding strength of C + 4H will suppress methane selectivity.

Ab initio diffusional potential energy surface for CO chemisorption on Pd{110} at high coverage:Coupled translation and rotation

The ground state potential energy surface for CO chemisorption across Pd{110} has been calculated using density functional theory with gradient corrections at monolayer coverage. The most stable site corresponds well with the experimental adsorption heat, and it is found that the strength of binding to sites is in the following order: pseudo-short-bridge>atop>long-bridge>hollow. Pathways and transition states for CO surface diffusion, involving a correlation between translation and orientation, are proposed and discussed. (C) 1997 American Institute of Physics.

Free fatty acid receptors:structural models and elucidation of ligand binding interactions

BACKGROUND: The free fatty acid receptors (FFAs), including FFA1 (orphan name: GPR40), FFA2 (GPR43) and FFA3 (GPR41) are G protein-coupled receptors (GPCRs) involved in energy and metabolic homeostasis. Understanding the structural basis of ligand binding at FFAs is an essential step toward designing potent and selective small molecule modulators.

RESULTS: We analyse earlier homology models of FFAs in light of the newly published FFA1 crystal structure co-crystallized with TAK-875, an ago-allosteric ligand, focusing on the architecture of the extracellular binding cavity and agonist-receptor interactions. The previous low-resolution homology models of FFAs were helpful in highlighting the location of the ligand binding site and the key residues for ligand anchoring. However, homology models were not accurate in establishing the nature of all ligand-receptor contacts and the precise ligand-binding mode. From analysis of structural models and mutagenesis, it appears that the position of helices 3, 4 and 5 is crucial in ligand docking. The FFA1-based homology models of FFA2 and FFA3 were constructed and used to compare the FFA subtypes. From docking studies we propose an alternative binding mode for orthosteric agonists at FFA1 and FFA2, involving the interhelical space between helices 4 and 5. This binding mode can explain mutagenesis results for residues at positions 4.56 and 5.42. The novel FFAs structural models highlight higher aromaticity of the FFA2 binding cavity and higher hydrophilicity of the FFA3 binding cavity. The role of the residues at the second extracellular loop used in mutagenesis is reanalysed. The third positively-charged residue in the binding cavity of FFAs, located in helix 2, is identified and predicted to coordinate allosteric modulators.

CONCLUSIONS: The novel structural models of FFAs provide information on specific modes of ligand binding at FFA subtypes and new suggestions for mutagenesis and ligand modification, guiding the development of novel orthosteric and allosteric chemical probes to validate the importance of FFAs in metabolic and inflammatory conditions. Using our FFA homology modelling experience, a strategy to model a GPCR, which is phylogenetically distant from GPCRs with the available crystal structures, is discussed.