1
|
Abstract
Thole-style mutual induction models for molecular polarization have been adopted by several popular polarizable force fields (FFs) for their simplicity and transferability. The atomic polarizability parameters of these models are typically derived by fitting to ab initio or/and experimental molecular polarizabilities. In this work, we improve upon Thole polarizability parameters by employing both high-level quantum mechanics molecular polarizabilities and electrostatic potential (ESP) responses on three-dimensional grids. Our results indicate that the two approaches to derive atomic polarizability parameters are both effective, while the ESP approaches can also capture the polarization for the atoms with lone pair electrons. The resulting polarizability parameters have been validated on a set of over 7200 molecules covering the most common elements found in organic molecules (C, H, O, N, P, S, F, Cl, Br, and I). These parameters have also been tested on the experimentally measured molecular polarizabilities of 422 molecules. The final set of parameters derived in this work show notable improvement over the current AMOEBA set. The result is a highly transferable, expanded set of atomic polarizabilities defined by the local chemical environment in the form of SMARTS patterns. These parameters can be used directly in molecular mechanics polarizable potential energy functions such as AMOEBA, AMOEBA+, and other Thole-style models.
Collapse
Affiliation(s)
| | | | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| |
Collapse
|
2
|
Tollefson MR, Litman JM, Qi G, O'Connell CE, Wipfler MJ, Marini RJ, Bernabe HV, Tollefson WTA, Braun TA, Casavant TL, Smith RJH, Schnieders MJ. Structural Insights into Hearing Loss Genetics from Polarizable Protein Repacking. Biophys J 2019; 117:602-612. [PMID: 31327459 DOI: 10.1016/j.bpj.2019.06.030] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2019] [Revised: 06/10/2019] [Accepted: 06/25/2019] [Indexed: 12/21/2022] Open
Abstract
Hearing loss is associated with ∼8100 mutations in 152 genes, and within the coding regions of these genes are over 60,000 missense variants. The majority of these variants are classified as "variants of uncertain significance" to reflect our inability to ascribe a phenotypic effect to the observed amino acid change. A promising source of pathogenicity information is biophysical simulation, although input protein structures often contain defects because of limitations in experimental data and/or only distant homology to a template. Here, we combine the polarizable atomic multipole optimized energetics for biomolecular applications force field, many-body optimization theory, and graphical processing unit acceleration to repack all deafness-associated proteins and thereby improve average structure MolProbity score from 2.2 to 1.0. We then used these optimized wild-type models to create over 60,000 structures for missense variants in the Deafness Variation Database, which are being incorporated into the Deafness Variation Database to inform deafness pathogenicity prediction. Finally, this work demonstrates that advanced polarizable atomic multipole force fields are efficient enough to repack the entire human proteome.
Collapse
Affiliation(s)
- Mallory R Tollefson
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa; Molecular Otolaryngology & Renal Research Laboratories, Department of Otolaryngology-Head and Neck Surgery, University of Iowa Hospitals and Clinics, Iowa City, Iowa
| | - Jacob M Litman
- Department of Biochemistry, University of Iowa, Iowa City, Iowa
| | - Guowei Qi
- Department of Biochemistry, University of Iowa, Iowa City, Iowa
| | - Claire E O'Connell
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa
| | - Matthew J Wipfler
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa
| | - Robert J Marini
- Molecular Otolaryngology & Renal Research Laboratories, Department of Otolaryngology-Head and Neck Surgery, University of Iowa Hospitals and Clinics, Iowa City, Iowa
| | - Hernan V Bernabe
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa; Molecular Otolaryngology & Renal Research Laboratories, Department of Otolaryngology-Head and Neck Surgery, University of Iowa Hospitals and Clinics, Iowa City, Iowa
| | | | - Terry A Braun
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa
| | - Thomas L Casavant
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa
| | - Richard J H Smith
- Molecular Otolaryngology & Renal Research Laboratories, Department of Otolaryngology-Head and Neck Surgery, University of Iowa Hospitals and Clinics, Iowa City, Iowa.
| | - Michael J Schnieders
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa; Department of Biochemistry, University of Iowa, Iowa City, Iowa.
| |
Collapse
|
3
|
Kondratick CM, Litman JM, Shaffer KV, Washington MT, Dieckman LM. Crystal structures of PCNA mutant proteins defective in gene silencing suggest a novel interaction site on the front face of the PCNA ring. PLoS One 2018; 13:e0193333. [PMID: 29499038 PMCID: PMC5834165 DOI: 10.1371/journal.pone.0193333] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2017] [Accepted: 02/08/2018] [Indexed: 11/19/2022] Open
Abstract
Proliferating cell nuclear antigen (PCNA), a homotrimeric protein, is the eukaryotic sliding clamp that functions as a processivity factor for polymerases during DNA replication. Chromatin association factor 1 (CAF-1) is a heterotrimeric histone chaperone protein that is required for coupling chromatin assembly with DNA replication in eukaryotes. CAF-1 association with replicating DNA, and the targeting of newly synthesized histones to sites of DNA replication and repair requires its interaction with PCNA. Genetic studies have identified three mutant forms of PCNA in yeast that cause defects in gene silencing and exhibit altered association of CAF-1 to chromatin in vivo, as well as inhibit binding to CAF-1 in vitro. Three of these mutant forms of PCNA, encoded by the pol30-6, pol30-8, and the pol30-79 alleles, direct the synthesis of PCNA proteins with the amino acid substitutions D41A/D42A, R61A/D63A, and L126A/I128A, respectively. Interestingly, these double alanine substitutions are located far away from each other within the PCNA protein. To understand the structural basis of the interaction between PCNA and CAF-1 and how disruption of this interaction leads to reduced gene silencing, we determined the X-ray crystal structures of each of these mutant PCNA proteins. All three of the substitutions caused disruptions of a surface cavity on the front face of the PCNA ring, which is formed in part by three loops comprised of residues 21–24, 41–44, and 118–134. We suggest that this cavity is a novel binding pocket required for the interaction between PCNA and CAF-1, and that this region in PCNA also represents a potential binding site for other PCNA-binding proteins.
Collapse
Affiliation(s)
- Christine M. Kondratick
- Department of Biochemistry, University of Iowa College of Medicine, Iowa City, IA, United States of America
| | - Jacob M. Litman
- Department of Biochemistry, University of Iowa College of Medicine, Iowa City, IA, United States of America
| | - Kurt V. Shaffer
- Department of Chemistry, Creighton University, Omaha, NE, United States of America
| | - M. Todd Washington
- Department of Biochemistry, University of Iowa College of Medicine, Iowa City, IA, United States of America
| | - Lynne M. Dieckman
- Department of Chemistry, Creighton University, Omaha, NE, United States of America
- * E-mail:
| |
Collapse
|
4
|
Nessler IJ, Litman JM, Schnieders MJ. Toward polarizable AMOEBA thermodynamics at fixed charge efficiency using a dual force field approach: application to organic crystals. Phys Chem Chem Phys 2016; 18:30313-30322. [PMID: 27524378 PMCID: PMC5102770 DOI: 10.1039/c6cp02595a] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
First principles prediction of the structure, thermodynamics and solubility of organic molecular crystals, which play a central role in chemical, material, pharmaceutical and engineering sciences, challenges both potential energy functions and sampling methodologies. Here we calculate absolute crystal deposition thermodynamics using a novel dual force field approach whose goal is to maintain the accuracy of advanced multipole force fields (e.g. the polarizable AMOEBA model) while performing more than 95% of the sampling in an inexpensive fixed charge (FC) force field (e.g. OPLS-AA). Absolute crystal sublimation/deposition phase transition free energies were determined using an alchemical path that grows the crystalline state from a vapor reference state based on sampling with the OPLS-AA force field, followed by dual force field thermodynamic corrections to change between FC and AMOEBA resolutions at both end states (we denote the three step path as AMOEBA/FC). Importantly, whereas the phase transition requires on the order of 200 ns of sampling per compound, only 5 ns of sampling was needed for the dual force field thermodynamic corrections to reach a mean statistical uncertainty of 0.05 kcal mol-1. For five organic compounds, the mean unsigned error between direct use of AMOEBA and the AMOEBA/FC dual force field path was only 0.2 kcal mol-1 and not statistically significant. Compared to experimental deposition thermodynamics, the mean unsigned error for AMOEBA/FC (1.4 kcal mol-1) was more than a factor of two smaller than uncorrected OPLS-AA (3.2 kcal mol-1). Overall, the dual force field thermodynamic corrections reduced condensed phase sampling in the expensive force field by a factor of 40, and may prove useful for protein stability or binding thermodynamics in the future.
Collapse
Affiliation(s)
- Ian J Nessler
- Department of Chemical Engineering, University of Iowa, Iowa City, IA 52242, USA
| | - Jacob M Litman
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242, USA
| | - Michael J Schnieders
- Department of Biochemistry, University of Iowa, Iowa City, IA 52242, USA and Department of Biomedical Engineering, University of Iowa, Iowa City, IA 52242, USA.
| |
Collapse
|
5
|
LuCore SD, Litman JM, Powers KT, Gao S, Lynn AM, Tollefson WTA, Fenn TD, Washington MT, Schnieders MJ. Dead-End Elimination with a Polarizable Force Field Repacks PCNA Structures. Biophys J 2015; 109:816-26. [PMID: 26287633 PMCID: PMC4547145 DOI: 10.1016/j.bpj.2015.06.062] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2015] [Revised: 06/07/2015] [Accepted: 06/29/2015] [Indexed: 11/15/2022] Open
Abstract
A balance of van der Waals, electrostatic, and hydrophobic forces drive the folding and packing of protein side chains. Although such interactions between residues are often approximated as being pairwise additive, in reality, higher-order many-body contributions that depend on environment drive hydrophobic collapse and cooperative electrostatics. Beginning from dead-end elimination, we derive the first algorithm, to our knowledge, capable of deterministic global repacking of side chains compatible with many-body energy functions. The approach is applied to seven PCNA x-ray crystallographic data sets with resolutions 2.5-3.8 Å (mean 3.0 Å) using an open-source software. While PDB_REDO models average an Rfree value of 29.5% and MOLPROBITY score of 2.71 Å (77th percentile), dead-end elimination with the polarizable AMOEBA force field lowered Rfree by 2.8-26.7% and improved mean MOLPROBITY score to atomic resolution at 1.25 Å (100th percentile). For structural biology applications that depend on side-chain repacking, including x-ray refinement, homology modeling, and protein design, the accuracy limitations of pairwise additivity can now be eliminated via polarizable or quantum mechanical potentials.
Collapse
Affiliation(s)
- Stephen D LuCore
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa
| | - Jacob M Litman
- Department of Biochemistry, University of Iowa, Iowa City, Iowa
| | - Kyle T Powers
- Department of Biochemistry, University of Iowa, Iowa City, Iowa
| | - Shibo Gao
- Department of Biochemistry, University of Iowa, Iowa City, Iowa
| | - Ava M Lynn
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa
| | | | | | | | - Michael J Schnieders
- Department of Biomedical Engineering, University of Iowa, Iowa City, Iowa; Department of Biochemistry, University of Iowa, Iowa City, Iowa.
| |
Collapse
|