Glycogen phosphorylase inhibitors: a free energy perturbation analysis of glucopyranose spirohydantoin analogues

GP catalyzes the phosphorylation of glycogen to Glc-1-P. Because of its fundamental role in the metabolism of glycogen, GP has been the target for a systematic structure-assisted design of inhibitory compounds, which could be of value in the therapeutic treatment of type 2 diabetes mellitus. The most potent catalytic-site inhibitor of GP identified to date is spirohydantoin of glucopyranose (hydan). In this work, we employ MD free energy simulations to calculate the relative binding affinities for GP of hydan and two spirohydantoin analogues, methyl-hydan and n-hydan, in which a hydrogen atom is replaced by a methyl- or amino group, respectively. The results are compared with the experimental relative affinities of these ligands, estimated by kinetic measurements of the ligand inhibition constants. The calculated binding affinity for methyl-hydan (relative to hydan) is 3.75 +/- 1.4 kcal/mol, in excellent agreement with the experimental value (3.6 +/- 0.2 kcal/mol). For n-hydan, the calculated value is 1.0 +/- 1.1 kcal/mol, somewhat smaller than the experimental result (2.3 +/- 0.1 kcal/mol). A free energy decomposition analysis shows that hydan makes optimum interactions with protein residues and specific water molecules in the catalytic site. In the other two ligands, structural perturbations of the active site by the additional methyl- or amino group reduce the corresponding binding affinities. The computed binding free energies are sensitive to the preference of a specific water molecule for two well-defined positions in the catalytic site. The behavior of this water is analyzed in detail, and the free energy profile for the translocation of the water between the two positions is evaluated. The results provide insights into the role of water molecules in modulating ligand binding affinities. A comparison of the interactions between a set of ligands and their surrounding groups in X-ray structures is often used in the interpretation of binding free energy differences and in guiding the design of new ligands. For the systems in this work, such an approach fails to estimate the order of relative binding strengths, in contrast to the rigorous free energy treatment.

The Landau free energy of hard ellipses obtained from microscopic simulations

Systems of two-dimensional hard ellipses of varying aspect ratios and packing fractions are studied by Monte Carlo simulations in the generalised canonical ensemble. From this microscopic model, we extract a coarse-grained macroscopic Landau-de Gennes free energy as a function of packing fraction and orientational order parameter. We separate the free energy into the ideal orientational entropy of non-interacting two-dimensional spins and an excess free energy associated with excluded volume interactions. We further explore the isotropic-nematic phase transition using our empirical expression for the free energy and find that the nature of the phase transition is continuous for the aspect ratios we studied.

Resolving the energy paradox of chaperone/usher-mediated fibre assembly

Periplasmic chaperone/usher machineries are used for assembly of filamentous adhesion organelles of Gram-negative pathogens in a process that has been suggested to be driven by folding energy. Structures of mutant chaperone-subunit complexes revealed a final folding transition (condensation of the subunit hydrophobic core) on the release of organelle subunit from the chaperone-subunit pre-assembly complex and incorporation into the final fibre structure. However, in view of the large interface between chaperone and subunit in the pre-assembly complex and the reported stability of this complex, it is difficult to understand how final folding could release sufficient energy to drive assembly. In the present paper, we show the X-ray structure for a native chaperone-fibre complex that, together with thermodynamic data, shows that the final folding step is indeed an essential component of the assembly process. We show that completion of the hydrophobic core and incorporation into the fibre results in an exceptionally stable module, whereas the chaperone-subunit preassembly complex is greatly destabilized by the high-energy conformation of the bound subunit. This difference in stabilities creates a free energy potential that drives fibre formation.

Finite-stretching corrections to the Milner-Witten-Cates theory for polymer brushes

This paper investigates finite-stretching corrections to the classical Milner-Witten-Cates theory for semi-dilute polymer brushes in a good solvent. The dominant correction to the free energy originates from an entropic repulsion caused by the impenetrability of the grafting surface, which produces a depletion of segments extending a distance $\mu \propto L^{-1}$ from the substrate, where $L$ is the classical brush height. The next most important correction is associated with the translational entropy of the chain ends, which creates the well-known tail where a small population of chains extend beyond the classical brush height by a distance $\xi \propto L^{-1/3}$. The validity of these corrections is confirmed by quantitative comparison with numerical self-consistent field theory.

Selectivity of bis(calix[4]diquinone) ionophores towards metal ions in solvent dimethylsulfoxide: a molecular mechanics and molecular dynamics study

Molecular modelling studies have been carried out on two bis(calix[4]diqu(inone) ionophores, each created from two (calix[4]diquinone)arenes bridged at their bottom rims via alkyl chains (CH2)(n), 1: n = 3, 2; n = 4, in order to understand the reported selectivity of these ligands towards different sized metal ions such as Na+, K+, Rb+, and Cs+ in dmso solution. Conformational. analyses have been carried out which show that in the lowest energy conformations of the two macrocycles, the individual calix[4]diquinones exhibit a combination of partial cone, 1,3-alternate and cone conformations. The interactions of these alkali metals with the macrocycles have been studied in the gas phase and in a periodic box of solvent dmso by molecular mechanics and molecular dynamics calculations. Molecular mechanics calculations have been carried out on the mode of entry of the ions into the macrocycles and suggest that this is likely to occur from the side of the central cavity, rather than through the main axis of the calix[4]diquinones. There are energy barriers of ca. 19 kcal mol(-1) for this entry path in the gas phase, but in solution no energy barrier is found. Molecular dynamics simulations show that in both 1 and 2, though particularly in the latter macrocycle, one or two solvent molecules are bonded to the metal throughout the course of the simulation, often to the exclusion, of one or more of the ether oxygen atoms. By contrast the carbonyl oxygen atoms remain bonded to the metal atoms throughout with bond lengths that remain significantly less than those to the ether oxygen atoms. Free energy perturbation studies have been carried out in dmso and indicate that for 1, the selectivity follows the order Rb+ approximate to K+ > Cs+ >> Na+, which is partially in agreement with the experimental results. The energy differences are small and indeed the ratio between stability constants found for Cs+ and K+ complexes is only 0.60, showing that 1 has only a slight preference for K+. For the larger receptor 2, which is better suited to metal complexation, the binding affinity follows the pattern Cs+ >> Rb+ >> K+ >> Na+, with energy differences of 5.75, 2.61, 2.78 kcal mol(-1) which is perfectly consistent with experimental results.

Isothermal titration calorimetry study of epicatechin binding to serum albumin

The interaction of epicatechin with bovine serum albumin (BSA) was studied by isothermal titration calorimetry. The binding constant (K) and associated thermodynamic binding parameters (n, Delta H) were determined for the interaction at three solution concentrations of BSA using a binding model assuming independent binding sites. These data show weak non-covalent binding of epicatechin to BSA. The interaction energetics varied with BSA concentration in the calorimeter cell, suggesting that the binding of epicatechin induced BSA aggregation. The free energy (Delta G) remained constant within a range of 2 kJ mol(-1) and negative entropy was observed, indicating an enthalpy driven exothermic interaction. It is concluded that the non-covalent epicatechin-BSA complex is formed by hydrogen bonding. (c) 2006 Elsevier B.V. All rights reserved.

Fast and accurate SCFT calculations for periodic block-copolymer morphologies using the spectral method with Anderson mixing

We study the numerical efficiency of solving the self-consistent field theory (SCFT) for periodic block-copolymer morphologies by combining the spectral method with Anderson mixing. Using AB diblock-copolymer melts as an example, we demonstrate that this approach can be orders of magnitude faster than competing methods, permitting precise calculations with relatively little computational cost. Moreover, our results raise significant doubts that the gyroid (G) phase extends to infinite $\chi N$. With the increased precision, we are also able to resolve subtle free-energy differences, allowing us to investigate the layer stacking in the perforated-lamellar (PL) phase and the lattice arrangement of the close-packed spherical (S$_{cp}$) phase. Furthermore, our study sheds light on the existence of the newly discovered Fddd (O$^{70}$) morphology, showing that conformational asymmetry has a significant effect on its stability.

Computer simulations of structures, energetics and dynamics of myoglobin center dot center dot center dot ligand complexes

Myoglobin has been studied in considerable detail using different experimental and computational techniques over the past decades. Recent developments in time-resolved spectroscopy have provided experimental data amenable to detailed atomistic simulations. The main theme of the present review are results on the structures, energetics and dynamics of ligands ( CO, NO) interacting with myoglobin from computer simulations. Modern computational methods including free energy simulations, mixed quantum mechanics/molecular mechanics simulations, and reactive molecular dynamics simulations provide insight into the dynamics of ligand dynamics in confined spaces complementary to experiment. Application of these methods to calculate and understand experimental observations for myoglobin interacting with CO and NO are presented and discussed.

Computational chemistry for elucidating protein function: energetics and dynamics of myoglobin-ligand systems

State-of-the-art computational methodologies are used to investigate the energetics and dynamics of photodissociated CO and NO in myoglobin (Mb···CO and Mb···NO). This includes the combination of molecular dynamics, ab initio MD, free energy sampling, and effective dynamics methods to compare the results with studies using X-ray crystallography and ultrafast spectroscopy metho ds. It is shown that modern simulation techniques along with careful description of the intermolecular interactions can give quantitative agreement with experiments on complex molecular systems. Based on this agreement predictions for as yet uncharacterized species can be made.

Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

Emotion and anticipation in an enactive framework for cognition (response to Andy Clark)

Undeniably, anticipation plays a crucial role in cognition. By what means, to what extent, and what it achieves remain open questions. In a recent BBS target article, Clark (in press) depicts an integrative model of the brain that builds on hierarchical Bayesian models of neural processing (Rao and Ballard, 1999; Friston, 2005; Brown et al., 2011), and their most recent formulation using the free-energy principle borrowed from thermodynamics (Feldman and Friston, 2010; Friston, 2010; Friston et al., 2010). Hierarchical generative models of cognition, such as those described by Clark, presuppose the manipulation of representations and internal models of the world, in as much detail as is perceptually available. Perhaps surprisingly, Clark acknowledges the existence of a “virtual version of the sensory data” (p. 4), but with no reference to some of the historical debates that shaped cognitive science, related to the storage, manipulation, and retrieval of representations in a cognitive system (Shanahan, 1997), or accounting for the emergence of intentionality within such a system (Searle, 1980; Preston and Bishop, 2002). Instead of demonstrating how this Bayesian framework responds to these foundational questions, Clark describes the structure and the functional properties of an action-oriented, multi-level system that is meant to combine perception, learning, and experience (Niedenthal, 2007).

Oppositely charged polyelectrolytes: complex formation and effects of chain asymmetry

The formation of complexes in solutions of oppositely charged polyions has been studied by Monte Carlo simulations. The amount as well as the length, and thus, the absolute charge of one of the polyions have been varied. There is an increasing tendency to form large clusters as the excess of one kind of polyion decreases. When all polyions have the same length, this tendency reaches a maximum near, but off, equivalent amounts of the two types of polyions. When one kind of polyion is made shorter, the propensity to form large clusters decreases and the fluctuations in cluster charge increases. Simple free-energy expressions have been formulated on the basis of a set of simple rules that help rationalize the observations. By calculating cluster distributions in both grand canonical and canonical ensembles, it has been possible to show the extent of finite-size effects in the simulations.

Homoleptic copper(II) and copper(I) complexes of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline. Six-coordinate copper(I) in solution

Reaction of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline (L) with Cu(ClO(4))(2)center dot 6H(2)O in methanol in 3:1 M ratio at room temperature yields light green [CuL(3)](ClO(4))(2)center dot H(2)O (1). The X-ray crystal structure of the hemi acetonitrile solvate [CuL(3)](ClO(4))(2)center dot 0.5CH(3)CN has been determined which shows Jahn-Teller distortion in the CuN(6) core present in the cation [CuL(3)](2+). Complex 1 gives an axial EPR spectrum in acetonitrile-toluene glass with g(parallel to) = 2.262 (A(parallel to) = 169 x 10 (4) cm (1)) and g(perpendicular to) = 2.069. The Cu(II/I) potential in 1 in CH(2)Cl(2) at a glassy carbon electrode is 0.32 V versus NHE. This potential does not change with the addition of extra L in the medium implicating generation of a six-coordinate copper(I) species [CuL(3)](+) in solution. B3LYP/LanL2DZ calculations show that the six Cu-N bond distances in [CuL(3)](+) are 2.33, 2.25, 2.32, 2.25, 2.28 and 2.25 angstrom while the ideal Cu(I)-N bond length in a symmetric Cu(I)N(6) moiety is estimated as 2.25 angstrom. Reaction of L with Cu(CH(3)CN)(4)ClO(4) in dehydrated methanol at room temperature even in 4:1 M proportion yields [CuL(2)]ClO(4) (2). Its (1)H NMR spectrum indicates that the metal in [CuL(2)](+) is tetrahedral. The Cu(II/I) potential in 2 is found to be 0.68 V versus NHE in CH(2)Cl(2) at a glassy carbon electrode. In presence of excess L, 2 yields the cyclic voltammogram of 1. From (1)H NMR titration, the free energy of binding of L to [CuL(2)](+) to produce [CuL(3)](+) in CD(2)Cl(2) at 298 K is estimated as -11.7 (+/-0.2) kJ mol (1).

Band gap control via tuning of inversion degree in CdIn2S4 spinel

Based on theoretical arguments we propose a possible route for controlling the band-gap in the promising photovoltaic material CdIn2S4. Our ab initio calculations show that the experimental degree of inversion in this spinel (fraction of tetrahedral sites occupied by In) corresponds approximately to the equilibrium value given by the minimum of the theoretical inversion free energy at a typical synthesis temperature. Modification of this temperature, or of the cooling rate after synthesis, is then expected to change the inversion degree, which in turn sensitively tunes the electronic band-gap of the solid, as shown here by Heyd-Scuseria-Ernzerhof screened hybrid functional calculations.

Reversible state transition in nanoconfined aqueous solutions

Using molecular dynamics simulations, we find a reversible transition between the dispersion and aggregation states of solute molecules in aqueous solutions confined in nanoscale geometry, which is not observed in macroscopic systems. The nanoscale confinement also leads to a significant increase of the critical aggregation concentration (CAC). A theoretical model based on Gibbs free energy calculation is developed to describe the simulation results. It indicates that the reversible state transition is attributed to the low free energy barrier (of order kBT) in between two energy minima corresponding to the dispersion and aggregation states, and the enhancement of the CAC results from the fact that at lower concentrations the number of solute molecules is not large enough to allow the formation of a stable cluster in the confined systems.