Evaluating parameterization protocols for hydration free energy calculations with the AMOEBA polarizable force field

Hydration free energy (HFE) calculations are often used to assess the performance of biomolecular force fields and the quality of assigned parameters. The AMOEBA polarizable force field moves beyond traditional pairwise additive models of electrostatics and may be expected to improve upon predictions of thermodynamic quantities such as HFEs over and above fixed point charge models. The recent SAMPL4 challenge evaluated the AMOEBA polarizable force field in this regard, but showed substantially worse results than those using the fixed point charge GAFF model. Starting with a set of automatically generated AMOEBA parameters for the SAMPL4 dataset, we evaluate the cumulative effects of a series of incremental improvements in parameterization protocol, including both solute and solvent model changes. Ultimately the optimized AMOEBA parameters give a set of results that are not statistically significantly different from those of GAFF in terms of signed and unsigned error metrics. This allows us to propose a number of guidelines for new molecule parameter derivation with AMOEBA, which we expect to have benefits for a range of biomolecular simulation applications such as protein ligand binding studie

Strategies to calculate water binding free energies in protein–ligand complexes

Water molecules are commonplace in protein binding pockets, where they can typically form a complex between the protein and a ligand or become displaced upon ligand binding. As a result, it is often of great interest to establish both the binding free energy and location of such molecules. Several approaches to predicting the location and affinity of water molecules to proteins have been proposed and utilized in the literature, although it is often unclear which method should be used under what circumstances. We report here a comparison between three such methodologies, Just Add Water Molecules (JAWS), Grand Canonical Monte Carlo (GCMC), and double-decoupling, in the hope of understanding the advantages and limitations of each method when applied to enclosed binding sites. As a result, we have adapted the JAWS scoring procedure, allowing the binding free energies of strongly bound water molecules to be calculated to a high degree of accuracy, requiring significantly less computational effort than more rigorous approaches. The combination of JAWS and GCMC offers a route to a rapid scheme capable of both locating and scoring water molecules for rational drug design

Three dimensional modeling via photographs for documentation of a village bath

24th International CIPA Symposium; Strasbourg; France; 2 September 2013 through 6 September 2013The aim of this study is supporting the conceptual discussions of architectural restoration with three dimensional modeling of monuments based on photogrammetric survey. In this study, a 16th century village bath in Ulamiş, Seferihisar, and Izmir is modeled for documentation. Ulamiş is one of the historical villages within which Turkish population first settled in the region of Seferihisar - Urla. The methodology was tested on an antique monument; a bath with a cubical form. Within the limits of this study, only the exterior of the bath was modeled. The presentation scale for the bath was determined as 1 / 50, considering the necessities of designing structural interventions and architectural ones within the scope of a restoration project. The three dimensional model produced is a realistic document presenting the present situation of the ruin. Traditional plan, elevation and perspective drawings may be produced from the model, in addition to the realistic textured renderings and wireframe representations. The model developed in this study provides opportunity for presenting photorealistic details of historical morphologies in scale. Compared to conventional drawings, the renders based on the 3d models provide an opportunity for conceiving architectural details such as color, material and texture. From these documents, relatively more detailed restitution hypothesis can be developed and intervention decisions can be taken. Finally, the principles derived from the case study can be used for 3d documentation of historical structures with irregular surfaces

Conformation and dynamics of human urotensin II and urotensin related peptide in aqueous solution

Conformation and dynamics of the vasoconstrictive peptides human urotensin II (UII) and urotensin related peptide (URP) have been investigated by both unrestrained and enhanced-sampling molecular-dynamics (MD) simulations and NMR spectroscopy. These peptides are natural ligands of the G-protein coupled urotensin II receptor (UTR) and have been linked to mammalian pathophysiology. UII and URP cannot be characterized by a single structure but exist as an equilibrium of two main classes of ring conformations, <i>open</i> and <i>folded</i>, with rapidly interchanging subtypes. The <i>open</i> states are characterized by turns of various types centered at K⁸Y⁹ or F⁶W⁷ predominantly with no or only sparsely populated transannular hydrogen bonds. The <i>folded</i> conformations show multiple turns stabilized by highly populated transannular hydrogen bonds comprising centers F⁶W⁷K⁸ or W⁷K⁸Y⁹. Some of these conformations have not been characterized previously. The equilibrium populations that are experimentally difficult to access were estimated by replica-exchange MD simulations and validated by comparison of experimental NMR data with chemical shifts calculated with density-functional theory. UII exhibits approximately 72% <i>open</i>:28% <i>folded</i> conformations in aqueous solution. URP shows very similar ring conformations as UII but differs in an <i>open:folded</i> equilibrium shifted further toward <i>open</i> conformations (86:14) possibly arising from the absence of folded N-terminal tail-ring interaction. The results suggest that the different biological effects of UII and URP are not caused by differences in ring conformations but rather by different interactions with UTR

CD1b-restricted GEM T cell responses are modulated by Mycobacterium tuberculosis mycolic acid meromycolate chains

Tuberculosis, caused by Mycobacterium tuberculosis, remains a major human pandemic. Germline-encoded mycolyl lipid-reactive (GEM) T cells are donor-unrestricted and recognize CD1b-presented mycobacterial mycolates. However, the molecular requirements governing mycolate antigenicity for the GEM T cell receptor (TCR) remain poorly understood. Here, we demonstrate CD1b expression in tuberculosis granulomas and reveal a central role for meromycolate chains in influencing GEM-TCR activity. Meromycolate fine structure influences T cell responses in TB-exposed individuals, and meromycolate alterations modulate functional responses by GEM-TCRs. Computational simulations suggest that meromycolate chain dynamics regulate mycolate head group movement, thereby modulating GEM-TCR activity. Our findings have significant implications for the design of future vaccines that target GEM T cells

The ELBA Force Field for Coarse-Grain Modeling of Lipid Membranes

A new coarse-grain model for molecular dynamics simulation of lipid membranes is presented. Following a simple and conventional approach, lipid molecules are modeled by spherical sites, each representing a group of several atoms. In contrast to common coarse-grain methods, two original (interdependent) features are here adopted. First, the main electrostatics are modeled explicitly by charges and dipoles, which interact realistically through a relative dielectric constant of unity (). Second, water molecules are represented individually through a new parametrization of the simple Stockmayer potential for polar fluids; each water molecule is therefore described by a single spherical site embedded with a point dipole. The force field is shown to accurately reproduce the main physical properties of single-species phospholipid bilayers comprising dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylethanolamine (DOPE) in the liquid crystal phase, as well as distearoylphosphatidylcholine (DSPC) in the liquid crystal and gel phases. Insights are presented into fundamental properties and phenomena that can be difficult or impossible to study with alternative computational or experimental methods. For example, we investigate the internal pressure distribution, dipole potential, lipid diffusion, and spontaneous self-assembly. Simulations lasting up to 1.5 microseconds were conducted for systems of different sizes (128, 512 and 1058 lipids); this also allowed us to identify size-dependent artifacts that are expected to affect membrane simulations in general. Future extensions and applications are discussed, particularly in relation to the methodology's inherent multiscale capabilities