Ian R. Gould

Lipid14: The Amber Lipid Force Field

The AMBER lipid force field has been updated to create Lipid14, allowing tensionless simulation of a number of lipid types with the AMBER MD package. The modular nature of this force field allows numerous combinations of head and tail groups to create different lipid types, enabling the easy insertion of new lipid species. The Lennard-Jones and torsion parameters of both the head and tail groups have been revised and updated partial charges calculated. The force field has been validated by simulating bilayers of six different lipid types for a total of 0.5 μs each without applying a surface tension; with favorable comparison to experiment for properties such as area per lipid, volume per lipid, bilayer thickness, NMR order parameters, scattering data, and lipid lateral diffusion. As the derivation of this force field is consistent with the AMBER development philosophy, Lipid14 is compatible with the AMBER protein, nucleic acid, carbohydrate, and small molecule force fields.

The Structure of Imidazolium-Based Ionic Liquids: Insights From Ion-Pair Interactions

A large number of ab-initio (B3LYP/6–31++G(d,p)) computed ion-pair structures have been examined in order to determine if such calculations are capable of offering insight into the physical properties of the liquid state, particularly viscosity and melting point. Ion pairings based around the 1-butyl-3-methylimidazolium (C4C1im) cations and a range of anions (Cl, BF4, and N(Tf)2 where N(Tf)2 is bis(trifluoromethylsulfonly)imide) were chosen because of the range of viscosities exhibited by the corresponding ionic liquids. We have used these results to build up a ‘picture’ of the ionic liquid structure which is consistent with molecular dynamics simulations and experimental evidence. However, further work is required to established if such an analysis could be predictive. Nevertheless, we establish clear relationships relating ion-pair association energy, a derived ‘connectivity index’, and the diversity of structures with viscosity and melting point. Our calculations indicate that ions in C4C1imCl form a strong, highly connected and regular array thus rationalizing the high viscosity and melting point. In contrast the ion-pairs of C4C1imN(Tf)2 form a weakly interacting, highly disordered, and low connectivity network consistent with the low viscosity and melting point. C4C1imBF4 lies midway between these two extremes.

GAFFlipid: a General Amber Force Field for the accurate molecular dynamics simulation of phospholipid

Previous attempts to simulate phospholipid bilayers using the General Amber Force Field (GAFF) yielded many bilayer characteristics in agreement with experiment, however when using a tensionless NPT ensemble the bilayer is seen to compress to an undesirable extent resulting in low areas per lipid and high order parameters in comparison to experiment. In this work, the GAFF Lennard-Jones parameters for the simulation of acyl chains are corrected to allow the accurate and stable simulation of pure lipid bilayers. Lipid bilayers comprised of six phospholipid types were simulated for timescales approaching a quarter of a microsecond under tensionless constant pressure conditions using Graphics Processing Units. Structural properties including area per lipid, volume per lipid, bilayer thickness, order parameter and headgroup hydration show favourable agreement with available experimental values. Expanding the system size from 72 to 288 lipids and a more experimentally realistic 2 × 288 lipid bilayer stack induces little change in the observed properties. This preliminary work is intended for combination with the new AMBER Lipid11 modular force field as part of on-going attempts to create a modular phospholipid AMBER force field allowing tensionless NPT simulations of complex lipid bilayers.

Assessment of standard force field models against high-quality ab initio potential curves for prototypes of pi-pi, CH/pi, and SH/pi interactions.

Several popular force fields, namely, CHARMM, AMBER, OPLS‐AA, and MM3, have been tested for their ability to reproduce highly accurate quantum mechanical potential energy curves for noncovalent interactions in the benzene dimer, the benzene‐CH4 complex, and the benzene‐H2S complex. All of the force fields are semi‐quantitatively correct, but none of them is consistently reliable quantitatively. Re‐optimization of Lennard‐Jones parameters and symmetry‐adapted perturbation theory analysis for the benzene dimer suggests that better agreement cannot be expected unless more flexible functional forms (particularly for the electrostatic contributions) are employed for the empirical force fields.

Structure and dynamics of phospholipid bilayers using recently developed general all-atom force fields

Two fully hydrated pure‐species phospholipids bilayers, 1,2‐dimyristoyl‐sn‐glycero‐3‐phosphocholine (DMPC) and 1,2‐dioleoyl‐sn‐glycero‐3‐phosphorylcholine (DOPC), in the fluid phase and explicit solvent have been studied using molecular dynamics simulation. Atom interactions were modeled using recently developed force fields based on AMBER with full atomistic details. Several representative liquid phase properties for the structure and dynamics of lipids with different length of hydrocarbon chains and different level of saturation have been reproduced without artificially biasing the system in order to match experimental data. In particular, as the new GAFF (General Amber Force Field) has not been explicitly developed to reproduce lipid characteristics and is naturally compatible with standard AMBER nucleic acids and proteins parameters, it is here proven a promising tool to study mixed lipid–protein processes as protein activity dependence on membrane composition, permeation of solute across membranes, and other cellular processes.

Glyoxylate as a backbone linkage for a prebiotic ancestor of RNA

The origin of the first RNA polymers is central to most current theories for the origin of life. Difficulties associated with the prebiotic formation of RNA have lead to the general consensus that a simpler polymer preceded RNA. However, polymers proposed as possible ancestors to RNA are not much easier to synthesize than RNA itself. One particular problem with the prebiotic synthesis of RNA is the formation of phosphoester bonds in the absence of chemical activation. Here we demonstrate that glyoxylate (the ionized form of glyoxylic acid), a plausible prebiotic molecule, represents a possible ancestor of the phosphate group in modern RNA. Although in low yields (∼ 1%), acetals are formed from glyoxylate and nucleosides under neutral conditions, provided that metal ions are present (e.g., Mg2+), and provided that water is removed by evaporation at moderate temperatures (e.g., 65 ∘C), i.e. under “drying conditions”. Such acetals are termed ga-dinucleotides and possess a linkage that is analogous to the backbone in RNA in both structure and electrostatic charge. Additionally, an energy-minimized model of a gaRNA duplex predicts a helical structure similar to that of A-form RNA. We propose that glyoxylate-acetal linkages would have had certain advantages over phosphate linkages for early self-replicating polymers, but that the distinct functional properties of phosphoester and phosphodiester bonds would have eventually lead to the replacement of glyoxylate by phosphate.

Direct identification and decongestion of Fermi resonances by control of pulse time ordering in two-dimensional IR spectroscopy

We show that it is possible to both directly measure and directly calculate Fermi resonance couplings in benzene. The measurement method used was a particular form of two-dimensional infrared spectroscopy (2D-IR) known as doubly vibrationally enhanced four wave mixing. By using different pulse orderings, vibrational cross peaks could be measured either purely at the frequencies of the base vibrational states or split by the coupling energy. This capability is a feature currently unique to this particular form of 2D-IR and can be helpful in the decongestion of complex spectra. Five cross peaks of the ring breathing mode nu13 with a range of combination bands were observed spanning a region of 1500-4550 cm(-1). The coupling energy was measured for two dominant states of the nu13+nu16 Fermi resonance tetrad. Dephasing rates were measured in the time domain for nu13 and the two (nu13+nu16) Fermi resonance states. The electronic and mechanical vibrational anharmonic coefficients were calculated to second and third orders, respectively, giving information on relative intensities of the cross peaks and enabling the Fermi resonance states of the combination band nu13+nu16 at 3050-3100 cm(-1) to be calculated. The excellent agreement between calculated and measured spectral intensities and line shapes suggests that assignment of spectral features from ab initio calculations is both viable and practicable for this form of spectroscopy.

Imaging phase separation in model lipid membranes through the use of BODIPY based molecular rotors

In order to fully understand the dynamics of processes within biological lipid membranes, it is necessary to possess an intimate knowledge of the physical state and ordering of lipids within the membrane. Here we report the use of three molecular rotors based on meso-substituted boron-dipyrrin (BODIPY) in combination with fluorescence lifetime spectroscopy to investigate the viscosity and phase behaviour of model lipid bilayers. In phase-separated giant unilamellar vesicles, we visualise both liquid-ordered (Lo) and liquid-disordered (Ld) phases using fluorescence lifetime imaging microscopy (FLIM), determining their associated viscosity values, and investigate the effect of composition on the viscosity of these phases. Additionally, we use molecular dynamics simulations to investigate the orientation of the BODIPY probes within the bilayer, as well as using molecular dynamics simulations and fluorescence correlation spectroscopy (FCS) to compare diffusion coefficients with those predicted from the fluorescence lifetimes of the probes.

Low-lying excited states of 4-dimethylaminobenzonitrile: equilibrium geometries, vibrational frequencies and charge transfer character

Abstract Equilibrium geometries and vibrational frequencies of the ground and low-lying excited states of 4-(dimethyl-amino)benzonitrile (DMABN) are reported (6–31G ∗∗ basis set, ab initio RHF and CI-singles methods). The results suggest that: (i) the ground state structure is planar (C s symmetry); (ii) the (relaxed) locally excited state has the dimethylamino group twisted somewhat about the NAr bond (C 1 symmetry); and (iii) the TICT state is actually a transition state between two mirror image locally excited state structures in the gas phase, where the MeN bonds are perpendicular to the NAr bond.

Functional asymmetry in the lysyl-tRNA synthetase explored by molecular dynamics, free energy calculations and experiment

BackgroundCharging of transfer-RNA with cognate amino acid is accomplished by the aminoacyl-tRNA synthetases, and proceeds through an aminoacyl adenylate intermediate. The lysyl-tRNA synthetase has evolved an active site that specifically binds lysine and ATP. Previous molecular dynamics simulations of the heat-inducible Escherichia coli lysyl-tRNA synthetase, LysU, have revealed differences in the binding of ATP and aspects of asymmetry between the nominally equivalent active sites of this dimeric enzyme. The possibility that this asymmetry results in different binding affinities for the ligands is addressed here by a parallel computational and biochemical study.ResultsBiochemical experiments employing isothermal calorimetry, steady-state fluorescence and circular dichroism are used to determine the order and stoichiometries of the lysine and nucleotide binding events, and the associated thermodynamic parameters. An ordered mechanism of substrate addition is found, with lysine having to bind prior to the nucleotide in a magnesium dependent process. Two lysines are found to bind per dimer, and trigger a large conformational change. Subsequent nucleotide binding causes little structural rearrangement and crucially only occurs at a single catalytic site, in accord with the simulations.Molecular dynamics based free energy calculations of the ATP binding process are used to determine the binding affinities of each site. Significant differences in ATP binding affinities are observed, with only one active site capable of realizing the experimental binding free energy. Half-of-the-sites models in which the nucleotide is only present at one active site achieve their full binding potential irrespective of the subunit choice. This strongly suggests the involvement of an anti-cooperative mechanism. Pathways for relaying information between the two active sites are proposed.ConclusionsThe asymmetry uncovered here appears to be a common feature of oligomeric aminoacyl-tRNA synthetases, and may play an important functional role. We suggest a manner in which catalytic efficiency could be improved by LysU operating in an alternating sites mechanism.