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.

Impact of energy-using product policies

Energy-using products (EuPs), such as domestic appliances, audio-visual and ICT equipment contribute significantly to CO2 emissions, both in the domestic and non-domestic sectors. Policies that encourage the use of more energy efficient products can therefore generate significant reductions in overall energy consumption and hence, CO2 emissions. To the extent that these policies cause an increase the average production cost of EuPs, they may impose economic costs on producers, or on consumers, or on both. In this theoretical paper, an adaptation of a simple vertical product differentiation model – in which products are characterised in terms of their quality and their energy consumption – is used to analyse the impact of the different EuP polices on product innovation and to assess the resultant economic impacts on producers and consumers. It is shown that whereas the imposition of a binding product standard for energy efficiency unambiguously reduces aggregate profit and increases the average market price in the absence of any learning effects, the introduction or strengthening of demand-side measures (such as energy labelling) may reduce, or increase, aggregate profit. Even in the case where the overall impact is unambiguously negative, the effects of product innovation and learning can be in either direction.

Synechococcus sp. PCC 7002 CpcB lyase null mutations alter phycocyanin chromophore function and, consequently, affect the redistribution of excitation energy via the light state transition /

Phycobilisomes are the major light harvesting complexes for cyanobacteria and phycocyanin is the primary phycobiliprotein of the phycobilisome rod. The phycocyanobilin lyases responsible for chromophorylating the phycocyanin p subunit (CpcB) have been recently identified in the cyanobacterium Synechococcus sp. PCC 7002. Surprisingly, mutants missing the CpcB lyases were nevertheless capable of producing pigmented phycocyanin. 10K absorbance measurements revealed that the energy states of the p phycocyanin chromophores were only subtly shifted; however, 77K steady state fluorescence emission spectroscopy showed excitation energy transfer involving the targeted chromophores to be highly disrupted. Such evidence suggests that phycobilin orientation within the binding domain is specifically modified. We hypothesized that alternate, less specific lyases are able to act on the p binding sites. A phycocyanin linker-polypeptide deficient mutant was similarly characterized. The light state transition, a short term adaptation of the photosynthetic light harvesting apparatus resulting in the redistribution of excitation energy among the photo systems, was shown to be dominated by the reallocation of phycocyanin-absorbed excitation energy. Treatment with a high M phosphate buffer effectively prevented the redistribution of both chlorophyll a- and phycobilisome- absorbed excitation energy, suggesting that the two effects are not strictly independent. The mutant strains required a larger redistribution of excitation energy between light states, perhaps to compensate for their loss in phycobilisome antenna function.

Orientation of Non-Native 2-Monosubstituted and 2,3-Disubstituted 1,4-Naphthoquinones in the A1 binding site of PSI and the effect on the rate of electron transfer : an electron paramagnetic resonance spectroscopy study

The proce-ss ofoxygenic photosynthesis is vital to life on Earth. the central event in photosynthesis is light induced electron transfer that converts light into energy for growth. Ofparticular significance is the membrane bound multisubunit protein known as Photosystem I (PSI). PSI is a reaction centre that is responsible for the transfer of electrons across the membrane to reduce NADP+ to NADPH. The recent publication ofa high resolution X-ray structure of PSI has shown new information about the structure, in particular the electron transfer cofactors, which allows us to study it in more detail. In PSI, the secondary acceptor is crucial for forward electron transfer. In this thesis, the effect of removing the native acceptor phylloquinone and replacing it with a series of structurally related quinones was investigated via transient electron paramagnetic resonance (EPR) experiments. The orientation of non native quinones in the binding site and their ability to function in the electron transfer process was determined. It was found that PSI will readily accept alkyl naphthoquinones and anthraquinone. Q band EPR experiments revealed that the non-native quinones are incorporated into the binding site with the same orientation of the headgroup as in the native system. X band EPR spectra and deuteration experiments indicate that monosubstituted naphthoquinones are bound to the Al site with their side group in the position occupied by the methyl group in native PSI (meta to the hydrogen bonded carbonyl oxygen). X band EPR experiments show that 2, 3- disubstituted methyl naphthoquinones are also incorporated into the Al site in the same orientation as phylloquinone, even with the presence of a halogen- or sulfur-containing side chain in the position normally occupied by the phytyl tail ofphylloquinone. The exception to this is 2-bromo-3-methyl --.- _. -. - -- - - 4 _._ _ _ - _ _ naphthoquinone which has a poorly resolved spectrum, making determination of the orientation difficuh. All of the non-native quinones studied act as efficient electron acceptors. However, forward electron transfer past the quinone could only be demonstrated for anthraquinone, which has a more negative midpoint potential than phylloquinone. In the case of anthraquinone, an increased rate of forward electron transfer compared to native PSI was found. From these results we can conclude that the rate ofelectron transfer from Al to Fx in native PSI lies in the normal region ofthe Marcus Curve.