Structure-based energetics of protein interfaces guides foot-and-mouth disease virus vaccine design

Virus capsids are primed for disassembly, yet capsid integrity is key to generating a protective immune response. Foot-and-mouth disease virus (FMDV) capsids comprise identical pentameric protein subunits held together by tenuous noncovalent interactions and are often unstable. Chemically inactivated or recombinant empty capsids, which could form the basis of future vaccines, are even less stable than live virus. Here we devised a computational method to assess the relative stability of protein-protein interfaces and used it to design improved candidate vaccines for two poorly stable, but globally important, serotypes of FMDV: O and SAT2. We used a restrained molecular dynamics strategy to rank mutations predicted to strengthen the pentamer interfaces and applied the results to produce stabilized capsids. Structural analyses and stability assays confirmed the predictions, and vaccinated animals generated improved neutralizing-antibody responses to stabilized particles compared to parental viruses and wild-type capsids.

Repulsion-induced surface-migration, by ballistics and bounce

The motion of adsorbate molecules across surfaces is fundamental to self-assembly, material growth, and heterogeneous catalysis. Recent Scanning Tunneling Microscopy studies have demonstrated the electron-induced long-range surface-migration of ethylene, benzene, and related molecules, moving tens of Angstroms across Si(100). We present a model of the previously unexplained long-range recoil of chemisorbed ethylene across the surface of silicon. The molecular dynamics reveal two key elements for directed long-range migration: first ‘ballistic’ motion that causes the molecule to leave the ab initio slab of the surface traveling 3–8 Å above it out of range of its roughness, and thereafter skipping-stone ‘bounces’ that transport it further to the observed long distances. Using a previously tested Impulsive Two-State model, we predict comparable long-range recoil of atomic chlorine following electron-induced dissociation of chlorophenyl chemisorbed at Cu(110)

Thermodynamic Stabillity of Hydrogen-Bonded Systems in Polar and Nonpolar Environments

The thermodynamic properties of a selected set of benchmark hydrogen-bonded systems (acetic acid dimer and the complexes of acetic acid with acetamide and methanol) was studied with the goal of obtaining detailed information on solvent effects on the hydrogen-bonded interactions using water, chloroform, and n-heptane as representatives for a wide range in the dielectric constant. Solvent effects were investigated using both explicit and implicit solvation models. For the explicit description of the solvent, molecular dynamics and Monte Carlo simulations in the isothermal isobaric (NpT) ensemble combined with the free energy perturbation technique were performed to determine solvation free energies. Within the implicit solvation approach, the polarizable continuum model and the conductor-like screening model were applied. Combination of gas phase results with the results obtained from the different solvation models through an appropriate thermodynamic cycle allows estimation of complexation free energies, enthalpies, and the respective entropic contributions in solution. Owing to the strong solvation effects of water the cyclic acetic acid dimer is not stable in aqueous solution. In less polar solvents the double hydrogen bond structure of the acetic acid dimer remains stable. This finding is in agreement with previous theoretical and experimental results. A similar trend as for the acetic acid dimer is also observed for the acetamide complex. The methanol complex was found to be thermodynamically unstable in gas phase as well as in any of the three solvents. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 31: 2046-2055, 2010

Electronic properties of a methane-water solution

Electronic properties of a methane-water solution were investigated by a sequential quantum mechanical/molecular dynamics approach. Upon hydration methane acquires an induced dipole moment of similar to 0.5 +/- 0.2 D. This is related to polarisation effects and to weak methane-water hydrogen bond interactions. From gas phase to solution, the first vertical excitation and ionisation energies of methane are red-shifted by 0.45 +/- 0.25 and 0.87 +/- 0.40 eV, respectively. We also report results for the dynamic polarisability of methane in water. In comparison with water, no difference was found for the average monomeric dipole moment of water molecules in close interaction with methane. (c) 2011 Elsevier B.V. All rights reserved.

Classical simulation of deposition of thiophene oligomers on TiO(2)-anatase: Relevance of long-range electrostatic interactions

We performed classical molecular dynamics simulations of the vapor-deposition of alpha-T4 oligomers on the TiO(2)-anatase (101) surface, comparing different sets of charges associated with the atoms of the model. The potential energy surfaces for alpha-T4 and TiO(2) were described by re-parametrizations of the Universal force field with charges given by the charge equilibration (QEq) scheme, or with fixed charges obtained by an ab initio method using the Hirshfeld partition. The two sets of charges lead to completely different results for the interface formation, and for the characteristics of the organic film, with a clearly defined alpha-T4 contact layer in the QEq case, and a more homogeneous molecular distribution when using Hirshfeld charges. The main reason for the discrepancy was found to be the incorrect charge assignment given by QEq to the sulfur and alpha-carbon atoms in thiophenes, and highlight the relevance of long-range interactions in the organization of molecular films. (C) 2009 Elsevier B.V. All rights reserved.

Spin current in the Mobius cyclacene belts

We report on density functional theory studies of the electronic structure and magnetic properties of Mobius-[n]cyclacenes. The geometry of Mobius bands presents a modulation of bond lengths that is needed to accommodate the twist. This modulation takes the form of bond alternation defects analogous to those of solitons in polyacetylene. The ground state of all Mobius bands is a triplet, with a spin density distribution that follows the bond length modulation. A molecular dynamics simulation of the Mobius cyclacene at 300 K shows that the twist travels around the belt inducing a magnetic current. (C) 2009 Elsevier B. V. All rights reserved.

Selective Wrapping and Supramolecular Structures of Polyfluorene-Carbon Nanotube Hybrids

We report on the photophysical properties of single-walled carbon nanotube (SWNT) suspensions In toluene solutions of poly[9,9-dioctylfluorenyl-2,7-diyl](PFO). Steady-state and time-resolved photoluminescence spectroscopy in the near-infrared and visible spectral regions are used to study the interaction of the dispersed SWNTs with the wrapped polymer. Molecular dynamics simulations of the PFO-SWNT hybrids in toluene were carried out to evaluate the energetics of different wrapping geometries. The simulated fluorescence spectra in the visible region were obtained by the quantum chemical ZINDO-CI method, using a sampling of structures obtained from the dynamics trajectories. The tested schemes consider polymer chains aligned along the nanotube axis, where chirality has a minimal effect, or forming helical structures, where a preference for high chiral angles is evidenced. Moreover, toluene affects the polymer structure favoring the helical conformation. Simulations show that the most stable hybrid system is the PFO-wrapped (8,6) nanotube, in agreement with the experimentally observed selectivity.

Computer simulations of copper and gold nanowires and single-wall nanowires

Copper and gold nanowires under tension evolve to form linear atomic chains (LACs), and the study and understanding of this evolution is an important subject for the development of nanocontacts. Here we study the differences and similarities between copper and gold nanowires (NWs) under stress along the [111] crystallographic direction until their rupture using tight-binding molecular dynamics. In both metals, the first significant rearrangement occurs due to one inside atom that goes to the NW` surface. In an attempt to better understand this effect, for both metals we also consider hollow NW`s where the inside atoms were excluded after the initial relaxation to create single-wall NW`s (SWNWs). The dynamical evolution of these SWNWs provides insight on the formation of the constriction that evolves to form LACs. Studying the calculated forces supported by the NW`s we show that SWNWs can sustain larger forces before the first major rearrangement in the copper and gold when compared to the original NW`s.

Adsorption Behavior of 5-Fluorouracil on Au(111): An In Situ STM Study

The biological effects of chemical substitution of DNA bases triggered several investigations of their physicochemical properties This paper studies the adsorption behavior of a halogenated uracil, 5-fluorouracil (5FU). at the electrochemical interface of Au(111) and sulfuric acid solution. Upon modulation of the electric field across the interface, four distinct phases could be inferred by means of cyclic voltammetry (CV) At negative potentials relative to the SCE electrode, limited by the threshold of hydrogen evolution, no molecular species could be detected by scanning tunneling microscopy (STM) at the reconstructed Au(111)-(23 x root 3) surface, indicating that any physisorbed molecules are randomly distributed Incursion into more positive potentials increases the surface population but doer not form any two-dimensional (2D) physisorbed ordered structure Instead, we observed metastable structures that are only detectable. on surfaces with high defect density At sufficiently high positive potentials. limited by gold oxidation, the molecules are chemisorbed in a (3 x 2 root 3) ordered structure. with the aromatic ring perpendicular to the surface We report the densest chemisorbed monolayer for pyrimidine-derivative molecules (area per molecule 0 14 +/- 0 04 nm(2)). A comparison of the adsorption behavior of uracil derivatives has been made based on recent results of chemical substitution and solvent effects. We propose that pi-stacking is enhanced when halogens are incorporated in the uracil structure, in a similar fashion to what is observed in then crystal structure

Differential effects of TR ligands on hormone dissociation rates: Evidence for multiple ligand entry/exit pathways

Some nuclear receptor (NR) ligands promote dissociation of radiolabeled bound hormone from the buried ligand binding cavity (LBC) more rapidly than excess unlabeled hormone itself This result was interpreted to mean that challenger ligands bind allosteric sites on the LBD to induce hormone dissociation, and recent findings indicate that ligands bind weakly to multiple sites on the LBD surface. Here we show, that a large fraction of thyroid hormone receptor (TR) ligands promote rapid dissociation (T(1/2) < 2 h) of , radiolabeled T(3) vs. T(3) (T(1/2), approximate to 5-7 h). We cannot discern relationships between this effect and ligand size, activity or affinity for TR beta. One ligand, GC-24, binds the TR LBC and (weakly) to the TR beta-LBD surface that mediates dimer/heterodimer interaction, but we cannot link this interaction to rapid T(3) dissociation. Instead, several lines of evidence suggest that the challenger ligand must interact with the buried LBC to promote rapid T(3) release. Since previous molecular dynamics simulations suggest that TR ligands leave the LBC by several routes, we propose that a subset of challenger ligands binds and stabilizes a partially unfolded intermediate state of TR that arises during T(3) release and that this effect enhances hormone dissociation. (C) 2009 Elsevier Ltd. All rights reserved.

Interaction of a C-terminal peptide of Bos taurus diacylglycerol acyltransferase 1 with model membranes

Diacylglycerol acyltransferase 1 (DGAT1) catalyzes the final and dedicated step in the synthesis of triacylglycerol, which is believed to involve the lipids oleoyl coenzyme A (OCoA) and dioleoyl-sn-glycerol (DOG) as substrates. In this work we investigated the interaction of a specific peptide, referred to as SIT2, on the C-terminal of DGAT1 (HKWCIRHFYKP) with model membranes made with OCoA and DOG in Langmuir monolayers and liposomes. According to the circular dichroism and fluorescence data, conformational changes on SIT2 were seen only on liposomes containing OCoA and DOG. In Langmuir monolayers, SIT2 causes the isotherms of neat OCoA and DOG monolayers to be expanded, but has negligible effect on mixed monolayers of OCoA and DOG. This synergistic interaction between SIT2 and DOG + OCoA may be rationalized in terms of a molecular model in which SIT2 may serve as a linkage between the two lipids. Our results therefore provide molecular-level evidence for the interaction between this domain and the substrates OCoA and DOG for the synthesis of triacylglycerol. (C) 2009 Elsevier B.V. All rights reserved.

Structural and thermodynamic analysis of thrombin:suramin interaction in solution and crystal phases

Suramin is a hexasulfonated naphthylurea which has been recently characterized as a non-competitive inhibitor of human alpha-thrombin activity over fibrinogen, although its binding site and mode of interaction with the enzyme remain elusive. Here, we determined two X-ray structure of the thrombin: suramin complex, refined at 2.4 angstrom resolution. While a single thrombin: suramin complex was found in the asymmetric unit cell of the crystal, some of the crystallographic contacts with symmetrically related molecules are mediated by both the enzyme and the ligand. Molecular dynamics simulations with the 1:1 complex demonstrate a large rearrangement of suramin in the complex, but with the protein scaffold and the more extensive protein-ligand regions keep unchanged. Small-angle X-ray scattering measurements at high micromolar concentration demonstrate a suramin-induced dimerization of the enzyme. These data indicating a dissimilar binding mode in the monomeric and oligomeric states, with a monomeric, 1:1 complex to be more likely to exist at the thrombin physiological, nanomolar concentration range. Collectively, close understanding on the structural basis for interaction is given which might establish a basis for design of suramin analogues targeting thrombin. Crown Copyright (C) 2009 Published by Elsevier B.V. All rights reserved.

Catalytic mechanism of inulinase from Arthrobacter sp S37

Detailed catalytic roles of the conserved Glu323, Asp460, and Glu519 of Arthrobacter sp. S37 inulinase (EnIA), a member of the glycoside hydrolase family 32, were investigated by site-directed mutagenesis and pH-dependence studies of the enzyme efficiency and homology modeling were carried out for EnIA and for D460E mutant. The enzyme efficiency (k(cat)/K-m) of the E323A and E519A mutants was significantly lower than that of the wild-type due to a substantial decrease in k(cat), but not due to variations in K-m, consistent with their putative roles as nucleophile and acid/base catalyst, respectively. The D460A mutant was totally inactive, whereas the D460E and D460N mutants were active to some extent, revealing Asp460 as a catalytic residue and demonstrating that the presence of a carboxylate group in this position is a prerequisite for catalysis. The pH-dependence studies indicated that the pK(a) of the acid/base catalyst decreased from 9.2 for the wild-type enzyme to 7.0 for the D460E mutant, implicating Asp460 as the residue that interacts with the acid/base catalyst Glu519 and elevates its pK(a). Homology modeling and molecular dynamics simulation of the wild-type enzyme and the D460E mutant shed light on the structural roles of Glu323, Asp460, and Glu519 in the catalytic activity of the enzyme. (C) 2008 Elsevier Inc. All rights reserved.

Analysis of Agonist and Antagonist Effects on Thyroid Hormone Receptor Conformation by Hydrogen/Deuterium Exchange

Thyroid hormone receptors (TRs) are ligand-gated transcription factors with critical roles in development and metabolism. Although x-ray structures of TR ligand-binding domains (LBDs) with agonists are available, comparable structures without ligand (apo-TR) or with antagonists are not. It remains important to understand apo-LBD conformation and the way that it rearranges with ligands to develop better TR pharmaceuticals. In this study, we conducted hydrogen/deuterium exchange on TR LBDs with or without agonist (T(3)) or antagonist (NH(3)). Both ligands reduce deuterium incorporation into LBD amide hydrogens, implying tighter overall folding of the domain. As predicted, mass spectroscopic analysis of individual proteolytic peptides after hydrogen/deuterium exchange reveals that ligand increases the degree of solvent protection of regions close to the buried ligand-binding pocket. However, there is also extensive ligand protection of other regions, including the dimer surface at H10-H11, providing evidence for allosteric communication between the ligand-binding pocket and distant interaction surfaces. Surprisingly, C-terminal activation helix H12, which is known to alter position with ligand, remains relatively protected from solvent in all conditions suggesting that it is packed against the LBD irrespective of the presence or type of ligand. T(3), but not NH(3), increases accessibility of the upper part of H3-H5 to solvent, and we propose that TR H12 interacts with this region in apo-TR and that this interaction is blocked by T(3) but not NH(3.) We present data from site-directed mutagenesis experiments and molecular dynamics simulations that lend support to this structural model of apo-TR and its ligand-dependent conformational changes. (Molecular Endocrinology 25: 15-31, 2011)

On the Denaturation Mechanisms of the Ligand Binding Domain of Thyroid Hormone Receptors

The ligand binding domain (LBD) of nuclear hormone receptors adopts a very compact, mostly alpha-helical structure that binds specific ligands with very high affinity. We use circular dichroism spectroscopy and high-temperature molecular dynamics Simulations to investigate unfolding of the LBDs of thyroid hormone receptors (TRs). A molecular description of the denaturation mechanisms is obtained by molecular dynamics Simulations of the TR alpha and TR beta LBDs in the absence and in the presence of the natural ligand Triac. The Simulations Show that the thermal unfolding of the LBD starts with the loss of native contacts and secondary Structure elements, while the Structure remains essentially compact, resembling a molten globule state. This differs From most protein denaturation simulations reported to date and suggests that the folding mechanism may start with the hydrophobic collapse of the TR LBDs. Our results reveal that the stabilities of the LBDs of the TR alpha and TR beta Subtypes are affected to different degrees by the binding of the isoform selective ligand Triac and that ligand binding confers protection against thermal denaturation and unfolding in a subtype specific manner. Our Simulations indicate two mechanisms by which the ligand stabilizes the LBD: (1) by enhancing the interactions between H8 and H 11, and the interaction of the region between H I and the Omega-loop with the core of the LBD, and (2) by shielding the hydrophobic H6 from hydration.