1
|
Sampson JM, Cannon DA, Duan J, Epstein JCK, Sergeeva AP, Katsamba PS, Mannepalli SM, Bahna FA, Adihou H, Guéret SM, Gopalakrishnan R, Geschwindner S, Rees DG, Sigurdardottir A, Wilkinson T, Dodd RB, De Maria L, Mobarec JC, Shapiro L, Honig B, Buchanan A, Friesner RA, Wang L. Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations. bioRxiv 2024:2024.04.22.590325. [PMID: 38712280 PMCID: PMC11071377 DOI: 10.1101/2024.04.22.590325] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2024]
Abstract
Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.
Collapse
Affiliation(s)
| | | | - Jianxin Duan
- Schrödinger, GmbH, Life Sciences Software, Mannheim, Germany
| | | | - Alina P. Sergeeva
- Columbia University, Department of Systems Biology, New York, NY, USA
| | | | - Seetha M. Mannepalli
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA, 10027
| | - Fabiana A. Bahna
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA, 10027
| | - Hélène Adihou
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, Gothenburg, Sweden
- Max Planck Institute of Molecular Physiology, AstraZeneca-MPI Satellite Unit, Dortmund, Germany
| | - Stéphanie M. Guéret
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, Gothenburg, Sweden
- Max Planck Institute of Molecular Physiology, AstraZeneca-MPI Satellite Unit, Dortmund, Germany
| | - Ranganath Gopalakrishnan
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism, BioPharmaceuticals R&D, Gothenburg, Sweden
- Max Planck Institute of Molecular Physiology, AstraZeneca-MPI Satellite Unit, Dortmund, Germany
| | - Stefan Geschwindner
- AstraZeneca, Mechanistic and Structural Biology, Discovery Sciences, R&D, Cambridge, UK
| | | | | | | | - Roger B. Dodd
- AstraZeneca, Biologics Engineering, R&D, Cambridge, UK
| | - Leonardo De Maria
- AstraZeneca, Medicinal Chemistry, Research and Early Development, Respiratory and Immunology, BioPharmaceuticals R&D, Gothenburg, Sweden
| | - Juan Carlos Mobarec
- AstraZeneca, Mechanistic and Structural Biology, Discovery Sciences, R&D, Cambridge, UK
| | - Lawrence Shapiro
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA, 10027
- Columbia University, Department of Biochemistry and Molecular Biophysics, New York, NY, USA
| | - Barry Honig
- Columbia University, Department of Systems Biology, New York, NY, USA
- Columbia University, Zuckerman Mind Brain Behavior Institute, New York, NY, USA, 10027
- Columbia University, Department of Biochemistry and Molecular Biophysics, New York, NY, USA
- Columbia University, Department of Medicine, New York, NY, USA
| | | | | | - Lingle Wang
- Schrödinger, Inc., Life Sciences Software, New York, NY, USA
| |
Collapse
|
2
|
Yi X, Zhang L, Friesner RA, McDermott A. Predicted and Experimental NMR Chemical Shifts at Variable Temperatures: The Effect of Protein Conformational Dynamics. J Phys Chem Lett 2024; 15:2270-2278. [PMID: 38381862 DOI: 10.1021/acs.jpclett.3c02589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [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] [Indexed: 02/23/2024]
Abstract
NMR chemical shifts provide a sensitive probe of protein structure and dynamics but remain challenging to predict and interpret. We examine the effect of protein conformational distributions on 15N chemical shifts for dihydrofolate reductase (DHFR), comparing QM/MM predicted shifts with experimental shifts in solution as well as frozen distributions. Representative snapshots from MD trajectories exhibit variation in predicted 15N chemical shifts of up to 25 ppm. The average over the fluctuations is in significantly better agreement with room temperature solution experimental values than the prediction for any single optimal conformations. Meanwhile, solid-state NMR (SSNMR) measurements of frozen solutions at 105 K exhibit broad lines whose widths agree well with the widths of distributions of predicted shifts for samples from the trajectory. The backbone torsion angle ψi-1 varies over 60° on the picosecond time scale, compensated by φi. These fluctuations can explain much of the shift variation.
Collapse
Affiliation(s)
- Xu Yi
- Department of Chemistry, Columbia University, New York, New York 10025, United States
| | - Lichirui Zhang
- Department of Chemistry, Columbia University, New York, New York 10025, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10025, United States
| | - Ann McDermott
- Department of Chemistry, Columbia University, New York, New York 10025, United States
| |
Collapse
|
3
|
Cao Y, Halls MD, Friesner RA. Highly efficient implementation of analytic nonadiabatic derivative couplings within the pseudospectral method. J Chem Phys 2024; 160:084106. [PMID: 38385510 DOI: 10.1063/5.0188277] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Accepted: 01/31/2024] [Indexed: 02/23/2024] Open
Abstract
A pseudospectral implementation of nonadiabatic derivative couplings in the Tamm-Dancoff approximation is reported, and the accuracy and efficiency of the pseudospectral nonadiabatic derivative couplings are studied. Our results demonstrate that the pseudospectral method provides mean absolute errors of 0.2%-1.9%, while providing a significant speedup. Benchmark calculations on fullerenes (Cn, n up to 100) using B3LYP achieved 10- to 15-fold, 8- to 17-fold, and 43- to 75-fold speedups for 6-31G**, 6-31++G**, and cc-pVTZ basis sets, respectively, when compared to the conventional spectral method.
Collapse
Affiliation(s)
- Yixiang Cao
- Schrödinger Inc., 1540 Broadway, 24th Floor, New York, New York 10036, USA
| | - Mathew D Halls
- Schrödinger Inc., 9868 Scranton, Suite 3200, San Diego, California 92121, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
4
|
Xu Y, Chen J, Aydt AP, Zhang L, Sergeyev I, Keeler EG, Choi B, He S, Reichman DR, Friesner RA, Nuckolls C, Steigerwald ML, Roy X, McDermott AE. Electron and Spin Delocalization in [Co 6 Se 8 (PEt 3 ) 6 ] 0/+1 Superatoms. Chemphyschem 2024; 25:e202300064. [PMID: 38057144 DOI: 10.1002/cphc.202300064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2023] [Revised: 11/01/2023] [Indexed: 12/08/2023]
Abstract
Molecular clusters can function as nanoscale atoms/superatoms, assembling into superatomic solids, a new class of solid-state materials with designable properties through modifications on superatoms. To explore possibilities on diversifying building blocks, here we thoroughly studied one representative superatom, Co6 Se8 (PEt3 )6 . We probed its structural, electronic, and magnetic properties and revealed its detailed electronic structure as valence electrons delocalize over inorganic [Co6 Se8 ] core while ligands function as an insulated shell. 59 Co SSNMR measurements on the core and 31 P, 13 C on the ligands show that the neutral Co6 Se8 (PEt3 )6 is diamagnetic and symmetric, with all ligands magnetically equivalent. Quantum computations cross-validate NMR results and reveal degenerate delocalized HOMO orbitals, indicating aromaticity. Ligand substitution keeps the inorganic core nearly intact. After losing one electron, the unpaired electron in [Co6 Se8 (PEt3 )6 ]+1 is delocalized, causing paramagnetism and a delocalized electron spin. Notably, this feature of electron/spin delocalization over a large cluster is attractive for special single-electron devices.
Collapse
Affiliation(s)
- Yunyao Xu
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Jia Chen
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Alexander P Aydt
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Lichirui Zhang
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Ivan Sergeyev
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Eric G Keeler
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Bonnie Choi
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Shoushou He
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - David R Reichman
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Colin Nuckolls
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | | | - Xavier Roy
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| | - Ann E McDermott
- Department of Chemistry, Columbia University New York, New York, 10027, USA
| |
Collapse
|
5
|
Coskun D, Lihan M, Rodrigues JPGLM, Vass M, Robinson D, Friesner RA, Miller EB. Using AlphaFold and Experimental Structures for the Prediction of the Structure and Binding Affinities of GPCR Complexes via Induced Fit Docking and Free Energy Perturbation. J Chem Theory Comput 2024; 20:477-489. [PMID: 38100422 DOI: 10.1021/acs.jctc.3c00839] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2023]
Abstract
Free energy perturbation (FEP) remains an indispensable method for computationally assaying prospective compounds in advance of synthesis. However, before FEP can be deployed prospectively, it must demonstrate retrospective recapitulation of known experimental data where the subtle details of the atomic ligand-receptor model are consequential. An open question is whether AlphaFold models can serve as useful initial models for FEP in the regime where there exists a congeneric series of known chemical matter but where no experimental structures are available either of the target or of close homologues. As AlphaFold structures are provided without a bound ligand, we employ induced fit docking to refine the AlphaFold models in the presence of one or more congeneric ligands. In this work, we first validate the performance of our latest induced fit docking technology, IFD-MD, on a retrospective set of public experimental GPCR structures with 95% of cross-docks producing a pose with a ligand RMSD ≤ 2.5 Å in the top two predictions. We then apply IFD-MD and FEP on AlphaFold models of the somatostatin receptor family of GPCRs. We use AlphaFold models produced prior to the availability of any experimental structure from this family. We arrive at FEP-validated models for SSTR2, SSTR4, and SSTR5, with RMSE around 1 kcal/mol, and explore the challenges of model validation under scenarios of limited ligand data, ample ligand data, and categorical data.
Collapse
Affiliation(s)
- Dilek Coskun
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Muyun Lihan
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | | | - Márton Vass
- Schrödinger Technologies Limited, Davidson House, First Floor, Reading RG1 3 EU, U.K
| | - Daniel Robinson
- Schrödinger Technologies Limited, Davidson House, First Floor, Reading RG1 3 EU, U.K
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3110, New York, New York 10036, United States
| | - Edward B Miller
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| |
Collapse
|
6
|
Weber JL, Vuong H, Friesner RA, Reichman DR. Expanding the Design Space of Constraints in Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2023; 19:7567-7576. [PMID: 37889331 DOI: 10.1021/acs.jctc.3c00654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2023]
Abstract
We formulate and characterize a new constraint for auxiliary-field quantum Monte Carlo (AFQMC) applicable for general fermionic systems, which allows for the accumulation of phase in the random walk but disallows walkers with a magnitude of phase greater than π with respect to the trial wave function. For short imaginary times, before walkers accumulate sizable phase values, this approach is equivalent to exact free projection, allowing one to observe the accumulation of bias associated with the constraint and thus estimate its magnitude a priori. We demonstrate the stability of this constraint over arbitrary imaginary times and system sizes, highlighting the removal of noise due to the fermionic sign problem. Benchmark total energies for a variety of weakly and strongly correlated molecular systems reveal a distinct bias with respect to standard phaseless AFQMC, with a comparative increase in accuracy given sufficient quality of the trial wave function for the set of studied cases. We then take this constraint, termed linecut AFQMC (lc-AFQMC), and systematically release it (lcR-AFQMC), providing a route to obtain a smooth bridge between constrained AFQMC and the exact free projection results.
Collapse
Affiliation(s)
- John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Hung Vuong
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
7
|
Debnath S, Neufeld VA, Jacobson LD, Rudshteyn B, Weber JL, Berkelbach TC, Friesner RA. Accurate Quantum Chemical Reaction Energies for Lithium-Mediated Electrolyte Decomposition and Evaluation of Density Functional Approximations. J Phys Chem A 2023; 127:9178-9184. [PMID: 37878768 PMCID: PMC10795021 DOI: 10.1021/acs.jpca.3c04369] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2023]
Abstract
An important concern related to the performance of Li-ion batteries is the formation of a solid electrolyte interphase on the surface of the anode. This film is formed from the decomposition of electrolytes and can have important effects on the stability and performance. Here, we evaluate the decomposition pathway of ethylene carbonate and related organic electrolyte molecules using a series of density functional approximations and correlated wave function (WF) methods, including the coupled-cluster theory with single, double, and perturbative triple excitations [CCSD(T)] and auxiliary-field quantum Monte Carlo (AFQMC). We find that the transition state barrier associated with ring opening varies widely across different functionals, ranging from 3.01 to 17.15 kcal/mol, which can be compared to the value of 12.84 kcal/mol predicted by CCSD(T). This large variation underscores the importance of benchmarking against accurate WF methods. A performance comparison of all of the density functionals used in this study reveals that the M06-2X-D3 (a meta-hybrid GGA), CAM-B3LYP-D3 (a range-separated hybrid), and B2GP-PLYP-D3 (a double hybrid) perform the best, with average errors of about 1.50-1.60 kcal/mol compared to CCSD(T). We also compared the performance of the WF methods that are more scalable than CCSD(T), finding that DLPNO-CCSD(T) and phaseless AFQMC with a DFT trial wave function exhibit average errors of 1.38 and 1.74 kcal/mol, respectively.
Collapse
Affiliation(s)
- Sibali Debnath
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Verena A Neufeld
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | | | - Benjamin Rudshteyn
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Schrödinger, Inc., New York, New York 10036, United States
| | - John L Weber
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Schrödinger, Inc., New York, New York 10036, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Schrödinger, Inc., New York, New York 10036, United States
| |
Collapse
|
8
|
Neugebauer H, Vuong HT, Weber JL, Friesner RA, Shee J, Hansen A. Toward Benchmark-Quality Ab Initio Predictions for 3d Transition Metal Electrocatalysts: A Comparison of CCSD(T) and ph-AFQMC. J Chem Theory Comput 2023; 19:6208-6225. [PMID: 37655473 DOI: 10.1021/acs.jctc.3c00617] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Generating accurate ab initio ionization energies for transition metal complexes is an important step toward the accurate computational description of their electrocatalytic reactions. Benchmark-quality data is required for testing existing theoretical methods and developing new ones but is complicated to obtain for many transition metal compounds due to the potential presence of both strong dynamical and static electron correlation. In this regime, it is questionable whether the so-called gold standard, coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), provides the desired level of accuracy─roughly 1-3 kcal/mol. In this work, we compiled a test set of 28 3d metal-containing molecules relevant to homogeneous electrocatalysis (termed 3dTMV) and computed their vertical ionization energies (ionization potentials) with CCSD(T) and phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) in the def2-SVP basis set. A substantial effort has been made to converge away the phaseless bias in the ph-AFQMC reference values. We assess a wide variety of multireference diagnostics and find that spin-symmetry breaking of the CCSD wave function and the PBE0 density functional correlate well with our analysis of multiconfigurational wave functions. We propose quantitative criteria based on symmetry breaking to delineate correlation regimes inside of which appropriately performed CCSD(T) can produce mean absolute deviations from the ph-AFQMC reference values of roughly 2 kcal/mol or less and outside of which CCSD(T) is expected to fail. We also present a preliminary assessment of density functional theory (DFT) functionals on the 3dTMV set.
Collapse
Affiliation(s)
- Hagen Neugebauer
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Hung T Vuong
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
9
|
Sergeeva AP, Katsamba PS, Liao J, Sampson JM, Bahna F, Mannepalli S, Morano NC, Shapiro L, Friesner RA, Honig B. Free Energy Perturbation Calculations of Mutation Effects on SARS-CoV-2 RBD::ACE2 Binding Affinity. J Mol Biol 2023:168187. [PMID: 37355034 DOI: 10.1016/j.jmb.2023.168187] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2023] [Revised: 06/13/2023] [Accepted: 06/15/2023] [Indexed: 06/26/2023]
Abstract
The strength of binding between human angiotensin converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of viral spike protein plays a role in the transmissibility of the SARS-CoV-2 virus. In this study we focus on a subset of RBD mutations that have been frequently observed in infected individuals and probe binding affinity changes to ACE2 using surface plasmon resonance (SPR) measurements and free energy perturbation (FEP) calculations. Our SPR results are largely in accord with previous studies but discrepancies do arise due to differences in experimental methods and to protocol differences even when a single method is used. Overall, we find that FEP performance is superior to that of other computational approaches examined as determined by agreement with experiment and, in particular, by its ability to identify stabilizing mutations. Moreover, the calculations successfully predict the observed cooperative stabilization of binding by the Q498R N501Y double mutant present in Omicron variants and offer a physical explanation for the underlying mechanism. Overall, our results suggest that despite the significant computational cost, FEP calculations may offer an effective strategy to understand the effects of interfacial mutations on protein-protein binding affinities and, hence, in a variety of practical applications such as the optimization of neutralizing antibodies.
Collapse
Affiliation(s)
- Alina P Sergeeva
- Department of Systems Biology, Columbia University Medical Center, New York, NY 10032
| | - Phinikoula S Katsamba
- Zuckerman Mind Brain and Behavior Institute, Columbia University, New York, NY 10027, USA
| | - Junzhuo Liao
- Department of Chemistry, Columbia University, New York, NY 10027, USA
| | - Jared M Sampson
- Department of Chemistry, Columbia University, New York, NY 10027, USA; Schrödinger, Inc., New York, NY 10036, USA
| | - Fabiana Bahna
- Zuckerman Mind Brain and Behavior Institute, Columbia University, New York, NY 10027, USA
| | - Seetha Mannepalli
- Zuckerman Mind Brain and Behavior Institute, Columbia University, New York, NY 10027, USA
| | - Nicholas C Morano
- Zuckerman Mind Brain and Behavior Institute, Columbia University, New York, NY 10027, USA
| | - Lawrence Shapiro
- Zuckerman Mind Brain and Behavior Institute, Columbia University, New York, NY 10027, USA; Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA.
| | | | - Barry Honig
- Department of Systems Biology, Columbia University Medical Center, New York, NY 10032; Zuckerman Mind Brain and Behavior Institute, Columbia University, New York, NY 10027, USA; Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA; Department of Medicine, Columbia University, New York, NY 10032.
| |
Collapse
|
10
|
Chen W, Cui D, Jerome SV, Michino M, Lenselink EB, Huggins DJ, Beautrait A, Vendome J, Abel R, Friesner RA, Wang L. Enhancing Hit Discovery in Virtual Screening through Absolute Protein-Ligand Binding Free-Energy Calculations. J Chem Inf Model 2023; 63:3171-3185. [PMID: 37167486 DOI: 10.1021/acs.jcim.3c00013] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
In the hit identification stage of drug discovery, a diverse chemical space needs to be explored to identify initial hits. Contrary to empirical scoring functions, absolute protein-ligand binding free-energy perturbation (ABFEP) provides a theoretically more rigorous and accurate description of protein-ligand binding thermodynamics and could, in principle, greatly improve the hit rates in virtual screening. In this work, we describe an implementation of an accurate and reliable ABFEP method in FEP+. We validated the ABFEP method on eight congeneric compound series binding to eight protein receptors including both neutral and charged ligands. For ligands with net charges, the alchemical ion approach is adopted to avoid artifacts in electrostatic potential energy calculations. The calculated binding free energies correlate with experimental results with a weighted average of R2 = 0.55 for the entire dataset. We also observe an overall root-mean-square error (RMSE) of 1.1 kcal/mol after shifting the zero-point of the simulation data to match the average experimental values. Through ABFEP calculations using apo versus holo protein structures, we demonstrated that the protein conformational and protonation state changes between the apo and holo proteins are the main physical factors contributing to the protein reorganization free energy manifested by the overestimation of raw ABFEP calculated binding free energies using the holo structures of the proteins. Furthermore, we performed ABFEP calculations in three virtual screening applications for hit enrichment. ABFEP greatly improves the hit rates as compared to docking scores or other methods like metadynamics. The good performance of ABFEP in rank ordering compounds demonstrated in this work confirms it as a useful tool to improve the hit rates in virtual screening, thus facilitating hit discovery.
Collapse
Affiliation(s)
- Wei Chen
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Di Cui
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Steven V Jerome
- Schrödinger, Inc., 10201 Wateridge Circle, Suite 220, San Diego, California 92121, United States
| | - Mayako Michino
- Tri-Institutional Therapeutics Discovery Institute, 413 E. 69th Street, New York, New York 10065, United States
| | | | - David J Huggins
- Tri-Institutional Therapeutics Discovery Institute, 413 E. 69th Street, New York, New York 10065, United States
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York 10065, United States
| | - Alexandre Beautrait
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Jeremie Vendome
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Robert Abel
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Lingle Wang
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| |
Collapse
|
11
|
Shee J, Weber JL, Reichman DR, Friesner RA, Zhang S. On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond. J Chem Phys 2023; 158:140901. [PMID: 37061483 PMCID: PMC10089686 DOI: 10.1063/5.0134009] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within "chemical accuracy" (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob's ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 - N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.
Collapse
Affiliation(s)
- James Shee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, USA
| |
Collapse
|
12
|
Yi X, Zhang L, Friesner RA, McDermott A. Predicted and Experimental NMR Chemical Shifts at Variable Temperatures: The Effect of Protein Conformational Dynamics. bioRxiv 2023:2023.01.25.525502. [PMID: 36747635 PMCID: PMC9900828 DOI: 10.1101/2023.01.25.525502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
NMR chemical shifts provide a sensitive probe of protein structure and dynamics. Prediction of shifts, and therefore interpretation of shifts, particularly for the frequently measured amidic 15 N sites, remains a tall challenge. We demonstrate that protein 15 N chemical shift prediction from QM/MM predictions can be improved if conformational variation is included via MD sampling, focusing on the antibiotic target, E. coli Dihydrofolate reductase (DHFR). Variations of up to 25 ppm in predicted 15 N chemical shifts are observed over the trajectory. For solution shifts the average of fluctuations on the low picosecond timescale results in a superior prediction to a single optimal conformation. For low temperature solid state measurements, the histogram of predicted shifts for locally minimized snapshots with specific solvent arrangements sampled from the trajectory explains the heterogeneous linewidths; in other words, the conformations and associated solvent are 'frozen out' at low temperatures and result in inhomogeneously broadened NMR peaks. We identified conformational degrees of freedom that contribute to chemical shift variation. Backbone torsion angles show high amplitude fluctuations during the trajectory on the low picosecond timescale. For a number of residues, including I60, ψ varies by up to 60º within a conformational basin during the MD simulations, despite the fact that I60 (and other sites studied) are in a secondary structure element and remain well folded during the trajectory. Fluctuations in ψ appear to be compensated by other degrees of freedom in the protein, including φ of the succeeding residue, resulting in "rocking" of the amide plane with changes in hydrogen bonding interactions. Good agreement for both room temperature and low temperature NMR spectra provides strong support for the specific approach to conformational averaging of computed chemical shifts.
Collapse
|
13
|
Coskun D, Chen W, Clark AJ, Lu C, Harder ED, Wang L, Friesner RA, Miller EB. Reliable and Accurate Prediction of Single-Residue p Ka Values through Free Energy Perturbation Calculations. J Chem Theory Comput 2022; 18:7193-7204. [PMID: 36384001 DOI: 10.1021/acs.jctc.2c00954] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Accurate prediction of the pKa's of protein residues is crucial to many applications in biological simulation and drug discovery. Here, we present the use of free energy perturbation (FEP) calculations for the prediction of single-protein residue pKa values. We begin with an initial set of 191 residues with experimentally determined pKa values. To isolate sampling limitations from force field inaccuracies, we develop an algorithm to classify residues whose environments are significantly affected by crystal packing effects. We then report an approach to identify buried histidines that require significant sampling beyond what is achieved in typical FEP calculations. We therefore define a clean data set not requiring algorithms capable of predicting major conformational changes on which other pKa prediction methods can be tested. On this data set, we report an RMSE of 0.76 pKa units for 35 ASP residues, 0.51 pKa units for 44 GLU residues, and 0.67 pKa units for 76 HIS residues.
Collapse
Affiliation(s)
- Dilek Coskun
- Schrödinger, Inc., 1540 Broadway, New York, New York10036, United States
| | - Wei Chen
- Schrödinger, Inc., 1540 Broadway, New York, New York10036, United States
| | - Anthony J Clark
- Schrödinger, Inc., 1540 Broadway, New York, New York10036, United States
| | - Chao Lu
- Schrödinger, Inc., 1540 Broadway, New York, New York10036, United States
| | - Edward D Harder
- Schrödinger, Inc., 1540 Broadway, New York, New York10036, United States
| | - Lingle Wang
- Schrödinger, Inc., 1540 Broadway, New York, New York10036, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3110, New York, New York10036, United States
| | - Edward B Miller
- Schrödinger, Inc., 1540 Broadway, New York, New York10036, United States
| |
Collapse
|
14
|
Dajnowicz S, Agarwal G, Stevenson JM, Jacobson LD, Ramezanghorbani F, Leswing K, Friesner RA, Halls MD, Abel R. High-Dimensional Neural Network Potential for Liquid Electrolyte Simulations. J Phys Chem B 2022; 126:6271-6280. [PMID: 35972463 DOI: 10.1021/acs.jpcb.2c03746] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Liquid electrolytes are one of the most important components of Li-ion batteries, which are a critical technology of the modern world. However, we still lack the computational tools required to accurately calculate key properties of these materials (viscosity and ionic diffusivity) from first principles necessary to support improved designs. In this work, we report a machine learning-based force field for liquid electrolyte simulations, which bridges the gap between the accuracy of range-separated hybrid density functional theory and the efficiency of classical force fields. Predictions of material properties made with this force field are quantitatively accurate compared to experimental data. Our model uses the QRNN deep neural network architecture, which includes both long-range interactions and global charge equilibration. The training data set is composed solely of non-periodic density functional theory (DFT), allowing the practical use of an accurate theory (here, ωB97X-D3BJ/def2-TZVPD), which would be prohibitively expensive for generating large data sets with periodic DFT. In this report, we focus on seven common carbonates and LiPF6, but this methodology has very few assumptions and can be readily applied to any liquid electrolyte system. This provides a promising path forward for large-scale atomistic modeling of many important battery chemistries.
Collapse
Affiliation(s)
| | - Garvit Agarwal
- Schrödinger, Inc., New York, New York 10036, United States
| | | | | | | | - Karl Leswing
- Schrödinger, Inc., New York, New York 10036, United States
| | - Richard A Friesner
- Schrödinger, Inc., New York, New York 10036, United States.,Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Mathew D Halls
- Schrödinger, Inc., San Diego, California 92121, United States
| | - Robert Abel
- Schrödinger, Inc., New York, New York 10036, United States
| |
Collapse
|
15
|
Xu T, Zhu K, Beautrait A, Vendome J, Borrelli KW, Abel R, Friesner RA, Miller EB. Induced-Fit Docking Enables Accurate Free Energy Perturbation Calculations in Homology Models. J Chem Theory Comput 2022; 18:5710-5724. [PMID: 35972903 DOI: 10.1021/acs.jctc.2c00371] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Homology models have been used for virtual screening and to understand the binding mode of a known active compound; however, rarely have the models been shown to be of sufficient accuracy, comparable to crystal structures, to support free-energy perturbation (FEP) calculations. We demonstrate here that the use of an advanced induced-fit docking methodology reliably enables predictive FEP calculations on congeneric series across homology models ≥30% sequence identity. Furthermore, we show that retrospective FEP calculations on a congeneric series of drug-like ligands are sufficient to discriminate between predicted binding modes. Results are presented for a total of 29 homology models for 14 protein targets, showing FEP results comparable to those obtained using experimentally determined crystal structures for 86% of homology models with template structure sequence identities ranging from 30 to 50%. Implications for the use and validation of homology models in drug discovery projects are discussed.
Collapse
Affiliation(s)
- Tianchuan Xu
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| | - Kai Zhu
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| | | | - Jeremie Vendome
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| | - Kenneth W Borrelli
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| | - Robert Abel
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3110, New York, New York 10036, United States
| | - Edward B Miller
- Schrödinger, Inc., 1540 Broadway, New York, New York 10036, United States
| |
Collapse
|
16
|
Weber JL, Vuong H, Devlaminck PA, Shee J, Lee J, Reichman DR, Friesner RA. A Localized-Orbital Energy Evaluation for Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2022; 18:3447-3459. [PMID: 35507769 DOI: 10.1021/acs.jctc.2c00111] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) has recently emerged as a promising method for the production of benchmark-level simulations of medium- to large-sized molecules because of its accuracy and favorable polynomial scaling with system size. Unfortunately, the memory footprints of standard energy evaluation algorithms are nontrivial, which can significantly impact timings on graphical processing units (GPUs) where memory is limited. Previous attempts to reduce scaling by taking advantage of the low-rank structure of the Coulombic integrals have been successful but exhibit high prefactors, making their utility limited to very large systems. Here we present a complementary cubic-scaling route to reduce memory and computational scaling based on the low rank of the Coulombic interactions between localized orbitals, focusing on the application to ph-AFQMC. We show that the error due to this approximation, which we term localized-orbital AFQMC (LO-AFQMC), is systematic and controllable via a single variable and that the method is computationally favorable even for small systems. We present results demonstrating robust retention of accuracy versus both experiment and full ph-AFQMC for a variety of test cases chosen for their potential difficulty for localized-orbital-based methods, including the singlet-triplet gaps of the polyacenes benzene through pentacene, the heats of formation for a set of Platonic hydrocarbon cages, and the total energy of ferrocene, Fe(Cp)2. Finally, we reproduce our previous result for the gas-phase ionization energy of Ni(Cp)2, agreeing with full ph-AFQMC to within statistical error while using less than 1/15th of the computer time.
Collapse
Affiliation(s)
- John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Hung Vuong
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Pierre A Devlaminck
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
17
|
Rudshteyn B, Weber JL, Coskun D, Devlaminck PA, Zhang S, Reichman DR, Shee J, Friesner RA. Calculation of Metallocene Ionization Potentials via Auxiliary Field Quantum Monte Carlo: Toward Benchmark Quantum Chemistry for Transition Metals. J Chem Theory Comput 2022; 18:2845-2862. [PMID: 35377642 PMCID: PMC9123894 DOI: 10.1021/acs.jctc.1c01071] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
The accurate ab initio prediction of ionization energies is essential to understanding the electrochemistry of transition metal complexes in both materials science and biological applications. However, such predictions have been complicated by the scarcity of gas phase experimental data, the relatively large size of the relevant molecules, and the presence of strong electron correlation effects. In this work, we apply all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multideterminant trial wave functions to six metallocene complexes to compare the computed adiabatic and vertical ionization energies with experimental results. We find that ph-AFQMC yields mean absolute errors (MAEs) of 1.69 ± 1.02 kcal/mol for the adiabatic energies and 2.85 ± 1.13 kcal/mol for the vertical energies. We also carry out density functional theory (DFT) calculations using a variety of functionals, which yields MAEs of 3.62-6.98 kcal/mol and 3.31-9.88 kcal/mol, as well as one variant of localized coupled cluster calculations (DLPNO-CCSD(T0) with moderate PNO cutoffs), which has MAEs of 4.96 and 6.08 kcal/mol, respectively. We also test the reliability of DLPNO-CCSD(T0) and DFT on acetylacetonate (acac) complexes for adiabatic energies measured in the same manner experimentally, and we find higher MAEs, ranging from 4.56 to 10.99 kcal/mol (with a different ordering) for DFT and 6.97 kcal/mol for DLPNO-CCSD(T0). Finally, by utilizing experimental solvation energies, we show that accurate reduction potentials in solution for the metallocene series can be obtained from the AFQMC gas phase results.
Collapse
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Pierre A Devlaminck
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United States.,Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
18
|
Cao Y, Halls MD, Vadicherla TR, Friesner RA. Pseudospectral implementations of long-range corrected density functional theory. J Comput Chem 2021; 42:2089-2102. [PMID: 34415620 DOI: 10.1002/jcc.26739] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2021] [Revised: 07/11/2021] [Accepted: 08/03/2021] [Indexed: 11/06/2022]
Abstract
We have implemented pseudospectral density-functional theory (DFT) with long-range corrected DFT functionals (PS-LRC) in quantum mechanics package Jaguar, and applied it in the calculations of geometry optimizations, dimmer interaction energies, polarizabilities and first-order hyperpolarizabilities, harmonic vibrational frequencies, S1 and T1 excitation energies, singlet-triplet gaps, charge transfer numbers, oscillator strengths, reaction barrier heights, electron-transfer couplings, and charge-transfer excitation energies. From our accuracy benchmark analysis, PS grids, PS dealiasing functions, PS atomic corrections, PS multigrid strategy, PS length scales, and PS cutoff scheme perform well in PS DFT with LRC density functionals with very small and ignorable deviations when compared to the conventional spectral (CS) method. The timing benchmark study of S1 excitation energy calculations of fullerenes (Cn , n up to 540) demonstrates that PS-LRC achieves 1.4-8.4-fold speedups in SCF, 22-32-fold speedups in Tamm-Dancoff approximation, and 6-15-fold speedups in total wall clock time with an average error 0.004 eV of excitation energies compared to the CS method.
Collapse
|
19
|
Cao Y, Halls MD, Friesner RA. Highly efficient implementation of the analytical gradients of pseudospectral time-dependent density functional theory. J Chem Phys 2021; 155:024115. [PMID: 34266272 DOI: 10.1063/5.0055379] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The accuracy and efficiency of time-dependent density functional theory (TDDFT) excited state gradient calculations using the pseudospectral method are presented. TDDFT excited state geometry optimizations of the G2 test set molecules, the organic fluorophores with large Stokes shifts, and the Pt-complexes show that the pseudospectral method gives average errors of 0.01-0.1 kcal/mol for the TDDFT excited state energy, 0.02-0.06 pm for the bond length, and 0.02-0.12° for the bond angle when compared to the results from conventional TDDFT. TDDFT gradient calculations of fullerenes (Cn, n up to 540) with the B3LYP functional and 6-31G** basis set show that the pseudospectral method provides 8- to 14-fold speedups in the total wall clock time over the conventional methods. The pseudospectral TDDFT gradient calculations with a diffuse basis set give higher speedups than the calculations for the same basis set without diffuse functions included.
Collapse
Affiliation(s)
- Yixiang Cao
- Schrödinger Inc., 120 West 45th Street, Tower 45, 17th Floor, New York, New York 10036, USA
| | - Mathew D Halls
- Schrödinger Inc., 10201 Wateridge Circle, Suite 220, San Diego, California 92121, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| |
Collapse
|
20
|
Lu C, Wu C, Ghoreishi D, Chen W, Wang L, Damm W, Ross GA, Dahlgren MK, Russell E, Von Bargen CD, Abel R, Friesner RA, Harder ED. OPLS4: Improving Force Field Accuracy on Challenging Regimes of Chemical Space. J Chem Theory Comput 2021; 17:4291-4300. [PMID: 34096718 DOI: 10.1021/acs.jctc.1c00302] [Citation(s) in RCA: 495] [Impact Index Per Article: 165.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Chao Lu
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | - Chuanjie Wu
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | - Delaram Ghoreishi
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | - Wei Chen
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | - Lingle Wang
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | - Wolfgang Damm
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | - Gregory A. Ross
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | - Markus K. Dahlgren
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | - Ellery Russell
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | | | - Robert Abel
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Edward D. Harder
- Schrodinger, Incorporated, 120 West 45th Street, New York, New York 10036, United States
| |
Collapse
|
21
|
Miller EB, Murphy RB, Sindhikara D, Borrelli KW, Grisewood MJ, Ranalli F, Dixon SL, Jerome S, Boyles NA, Day T, Ghanakota P, Mondal S, Rafi SB, Troast DM, Abel R, Friesner RA. Reliable and Accurate Solution to the Induced Fit Docking Problem for Protein-Ligand Binding. J Chem Theory Comput 2021; 17:2630-2639. [PMID: 33779166 DOI: 10.1021/acs.jctc.1c00136] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present a reliable and accurate solution to the induced fit docking problem for protein-ligand binding by combining ligand-based pharmacophore docking, rigid receptor docking, and protein structure prediction with explicit solvent molecular dynamics simulations. This novel methodology in detailed retrospective and prospective testing succeeded to determine protein-ligand binding modes with a root-mean-square deviation within 2.5 Å in over 90% of cross-docking cases. We further demonstrate these predicted ligand-receptor structures were sufficiently accurate to prospectively enable predictive structure-based drug discovery for challenging targets, substantially expanding the domain of applicability for such methods.
Collapse
Affiliation(s)
- Edward B Miller
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Robert B Murphy
- Schrödinger, Inc., 10201 Wateridge Circle, San Diego, California 92121, United States
| | - Daniel Sindhikara
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Kenneth W Borrelli
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Matthew J Grisewood
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Fabio Ranalli
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Steven L Dixon
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Steven Jerome
- Schrödinger, Inc., 10201 Wateridge Circle, San Diego, California 92121, United States
| | - Nicholas A Boyles
- Schrödinger, Inc., 101 SW Main Street, Portland, Oregon 97204, United States
| | - Tyler Day
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Phani Ghanakota
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Sayan Mondal
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Salma B Rafi
- Schrödinger, Inc., 101 SW Main Street, Portland, Oregon 97204, United States
| | - Dawn M Troast
- Morphic Therapeutic, 35 Gatehouse Drive, Waltham, Massachusetts 02451, United States
| | - Robert Abel
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3110, New York, New York 10036, United States
| |
Collapse
|
22
|
Zhou T, Tsybovsky Y, Gorman J, Rapp M, Cerutti G, Chuang GY, Katsamba PS, Sampson JM, Schön A, Bimela J, Boyington JC, Nazzari A, Olia AS, Shi W, Sastry M, Stephens T, Stuckey J, Teng IT, Wang P, Wang S, Zhang B, Friesner RA, Ho DD, Mascola JR, Shapiro L, Kwong PD. Cryo-EM Structures of SARS-CoV-2 Spike without and with ACE2 Reveal a pH-Dependent Switch to Mediate Endosomal Positioning of Receptor-Binding Domains. Cell Host Microbe 2020; 28:867-879.e5. [PMID: 33271067 PMCID: PMC7670890 DOI: 10.1016/j.chom.2020.11.004] [Citation(s) in RCA: 239] [Impact Index Per Article: 59.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2020] [Revised: 10/20/2020] [Accepted: 11/09/2020] [Indexed: 01/21/2023]
Abstract
The SARS-CoV-2 spike employs mobile receptor-binding domains (RBDs) to engage the human ACE2 receptor and to facilitate virus entry, which can occur through low-pH-endosomal pathways. To understand how ACE2 binding and low pH affect spike conformation, we determined cryo-electron microscopy structures-at serological and endosomal pH-delineating spike recognition of up to three ACE2 molecules. RBDs freely adopted "up" conformations required for ACE2 interaction, primarily through RBD movement combined with smaller alterations in neighboring domains. In the absence of ACE2, single-RBD-up conformations dominated at pH 5.5, resolving into a solitary all-down conformation at lower pH. Notably, a pH-dependent refolding region (residues 824-858) at the spike-interdomain interface displayed dramatic structural rearrangements and mediated RBD positioning through coordinated movements of the entire trimer apex. These structures provide a foundation for understanding prefusion-spike mechanics governing endosomal entry; we suggest that the low pH all-down conformation potentially facilitates immune evasion from RBD-up binding antibody.
Collapse
Affiliation(s)
- Tongqing Zhou
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Yaroslav Tsybovsky
- Electron Microscopy Laboratory, Cancer Research Technology Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA
| | - Jason Gorman
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Micah Rapp
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA
| | - Gabriele Cerutti
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA
| | - Gwo-Yu Chuang
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Phinikoula S Katsamba
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA
| | - Jared M Sampson
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA; Department of Chemistry, Columbia University, New York, NY 10027, USA
| | - Arne Schön
- Department of Biology, Johns Hopkins University, Baltimore, MD, 21218, USA
| | - Jude Bimela
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA
| | - Jeffrey C Boyington
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Alexandra Nazzari
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Adam S Olia
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Wei Shi
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Mallika Sastry
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Tyler Stephens
- Electron Microscopy Laboratory, Cancer Research Technology Program, Leidos Biomedical Research, Inc., Frederick National Laboratory for Cancer Research, Frederick, MD 21702, USA
| | - Jonathan Stuckey
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - I-Ting Teng
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Pengfei Wang
- Aaron Diamond AIDS Research Center, Columbia University Vagelos College of Physicians and Surgeons, New York, NY 10032, USA
| | - Shuishu Wang
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Baoshan Zhang
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | | | - David D Ho
- Aaron Diamond AIDS Research Center, Columbia University Vagelos College of Physicians and Surgeons, New York, NY 10032, USA
| | - John R Mascola
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA
| | - Lawrence Shapiro
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA; Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA; Aaron Diamond AIDS Research Center, Columbia University Vagelos College of Physicians and Surgeons, New York, NY 10032, USA.
| | - Peter D Kwong
- Vaccine Research Center, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Bethesda, MD 20892, USA; Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA.
| |
Collapse
|
23
|
Weber JL, Churchill EM, Jockusch S, Arthur EJ, Pun AB, Zhang S, Friesner RA, Campos LM, Reichman DR, Shee J. In silico prediction of annihilators for triplet-triplet annihilation upconversion via auxiliary-field quantum Monte Carlo. Chem Sci 2020; 12:1068-1079. [PMID: 34163873 PMCID: PMC8179011 DOI: 10.1039/d0sc03381b] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
The energy of the lowest-lying triplet state (T1) relative to the ground and first-excited singlet states (S0, S1) plays a critical role in optical multiexcitonic processes of organic chromophores. Focusing on triplet–triplet annihilation (TTA) upconversion, the S0 to T1 energy gap, known as the triplet energy, is difficult to measure experimentally for most molecules of interest. Ab initio predictions can provide a useful alternative, however low-scaling electronic structure methods such as the Kohn–Sham and time-dependent variants of Density Functional Theory (DFT) rely heavily on the fraction of exact exchange chosen for a given functional, and tend to be unreliable when strong electronic correlation is present. Here, we use auxiliary-field quantum Monte Carlo (AFQMC), a scalable electronic structure method capable of accurately describing even strongly correlated molecules, to predict the triplet energies for a series of candidate annihilators for TTA upconversion, including 9,10 substituted anthracenes and substituted benzothiadiazole (BTD) and benzoselenodiazole (BSeD) compounds. We compare our results to predictions from a number of commonly used DFT functionals, as well as DLPNO-CCSD(T0), a localized approximation to coupled cluster with singles, doubles, and perturbative triples. Together with S1 estimates from absorption/emission spectra, which are well-reproduced by TD-DFT calculations employing the range-corrected hybrid functional CAM-B3LYP, we provide predictions regarding the thermodynamic feasibility of upconversion by requiring (a) the measured T1 of the sensitizer exceeds that of the calculated T1 of the candidate annihilator, and (b) twice the T1 of the annihilator exceeds its S1 energetic value. We demonstrate a successful example of in silico discovery of a novel annihilator, phenyl-substituted BTD, and present experimental validation via low temperature phosphorescence and the presence of upconverted blue light emission when coupled to a platinum octaethylporphyrin (PtOEP) sensitizer. The BTD framework thus represents a new class of annihilators for TTA upconversion. Its chemical functionalization, guided by the computational tools utilized herein, provides a promising route towards high energy (violet to near-UV) emission. Electronic structure theories such as AFQMC can accurately predict the low-lying excited state energetics of organic chromophores involved in triplet–triplet annihilation upconversion. A novel class of benzothiadiazole annihilators is discovered.![]()
Collapse
Affiliation(s)
- John L Weber
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Emily M Churchill
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Steffen Jockusch
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Evan J Arthur
- Schrodinger Inc 120 West 45th Street New York NY 1003 USA
| | - Andrew B Pun
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute 162 5th Avenue New York NY 10010 USA.,Department of Physics, College of William and Mary Williamsburg VA 23187 USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Luis M Campos
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - David R Reichman
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - James Shee
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| |
Collapse
|
24
|
Zhou T, Tsybovsky Y, Olia AS, Gorman J, Rapp MA, Cerutti G, Chuang GY, Katsamba PS, Nazzari A, Sampson JM, Schon A, Wang PD, Bimela J, Shi W, Teng IT, Zhang B, Boyington JC, Sastry M, Stephens T, Stuckey J, Wang S, Friesner RA, Ho DD, Mascola JR, Shapiro L, Kwong PD. Cryo-EM Structures Delineate a pH-Dependent Switch that Mediates Endosomal Positioning of SARS-CoV-2 Spike Receptor-Binding Domains. bioRxiv 2020. [PMID: 32637958 DOI: 10.1101/2020.07.04.187989] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The SARS-CoV-2 spike employs mobile receptor-binding domains (RBDs) to engage the ACE2 receptor and to facilitate virus entry. Antibodies can engage RBD but some, such as CR3022, fail to inhibit entry despite nanomolar spike affinity. Here we show the SARS-CoV-2 spike to have low unfolding enthalpy at serological pH and up to 10-times more unfolding enthalpy at endosomal pH, where we observe significantly reduced CR3022 affinity. Cryo-EM structures -at serological and endosomal pH- delineated spike recognition of up to three ACE2 molecules, revealing RBD to freely adopt the 'up' conformation. In the absence of ACE2, single-RBD-up conformations dominated at pH 5.5, resolving into a locked all-down conformation at lower pH. Notably, a pH-dependent refolding region (residues 824-858) at the spike-interdomain interface displayed dramatic structural rearrangements and mediated RBD positioning and spike shedding of antibodies like CR3022. An endosomal mechanism involving spike-conformational change can thus facilitate immune evasion from RBD-'up'-recognizing antibody.
Collapse
|
25
|
Kumar M, Shee J, Rudshteyn B, Reichman DR, Friesner RA, Miller CE, Francisco JS. Multiple Stable Isoprene-Ozone Complexes Reveal Complex Entrance Channel Dynamics in the Isoprene + Ozone Reaction. J Am Chem Soc 2020; 142:10806-10813. [PMID: 32431151 DOI: 10.1021/jacs.0c02360] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Accurately characterizing isoprene ozonolysis continues to challenge atmospheric chemists. The reaction is believed to be a spontaneous, concerted cycloaddition. However, little information is available about the entrance channel and isoprene-ozone complexes thought to define the long-range portion of the reaction coordinate. Our coupled cluster and auxiliary field quantum Monte Carlo calculations predict multiple stable isoprene-ozone van der Waals complexes for trans-isoprene in the gas phase with moderate association energies. These results indicate that long-range dynamics in the isoprene-ozone entrance channel can impact the overall reaction in the troposphere and provide the spectroscopic information necessary to extend the microwave characterization of isoprene ozonolysis to prereactive complexes. At the air-water interface, Born-Oppenheimer molecular dynamics simulations indicate that the cycloaddition reaction between ozone and trans-isoprene follows a stepwise mechanism, which is quite distinct from our proposed gas-phase mechanism and occurs on a femtosecond time scale. The stepwise nature of isoprene ozonolysis on the aqueous surface is more consistent with the DeMore mechanism than with the Criegee mechanism suggested by the gas-phase calculations, suggesting that the reaction media may play an important role in the reaction. Overall, these predictions aim to provide a missing fundamental piece of molecular insight into isoprene ozonolysis, which has broad tropospheric implications due to its critical role as a nighttime source of hydroxyl radicals.
Collapse
Affiliation(s)
- Manoj Kumar
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - James Shee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Benjamin Rudshteyn
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Charles E Miller
- Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California 91109, United States
| | - Joseph S Francisco
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
26
|
Rudshteyn B, Coskun D, Weber JL, Arthur EJ, Zhang S, Reichman DR, Friesner RA, Shee J. Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2020; 16:3041-3054. [DOI: 10.1021/acs.jctc.0c00070] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L. Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrodinger Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
27
|
Li G, Rudshteyn B, Shee J, Weber JL, Coskun D, Bochevarov AD, Friesner RA. Accurate Quantum Chemical Calculation of Ionization Potentials: Validation of the DFT-LOC Approach via a Large Data Set Obtained from Experiments and Benchmark Quantum Chemical Calculations. J Chem Theory Comput 2020; 16:2109-2123. [PMID: 32150400 DOI: 10.1021/acs.jctc.9b00875] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Density functional theory (DFT) is known to often fail when calculating thermodynamic values, such as ionization potentials (IPs), due to nondynamical error (i.e., the self-interaction term). Localized orbital corrections (LOCs), derived from assigning corresponding corrections for the atomic orbitals, bonds, and paired and unpaired electrons, are utilized to correct the IPs calculated from DFT. Some of the assigned parameters, which are physically due to the contraction of and change of the environment around a bond, depend on identifying the location in the molecule from which the electron is removed using differences in the charge density between neutral and oxidized species. In our training set, various small organic and inorganic molecules from the literature with the reported experimental IP were collected using the NIST database. For certain molecules with uncertain or no experimental measurements, we obtain the IP using coupled cluster theory and auxiliary field quantum Monte Carlo. After applying these corrections, as generated by least-squares regression, LOC reduces the mean absolute deviation (MAD) of the training set from 0.143 to 0.046 eV (R2 = 0.895), and LOC reduces the MAD of the test set from 0.192 to 0.097 eV (R2 = 0.833).
Collapse
Affiliation(s)
- Guangqi Li
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Benjamin Rudshteyn
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | | | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
28
|
Shee J, Arthur EJ, Zhang S, Reichman DR, Friesner RA. Singlet–Triplet Energy Gaps of Organic Biradicals and Polyacenes with Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2019; 15:4924-4932. [DOI: 10.1021/acs.jctc.9b00534] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrodinger Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
29
|
Clark AJ, Negron C, Hauser K, Sun M, Wang L, Abel R, Friesner RA. Relative Binding Affinity Prediction of Charge-Changing Sequence Mutations with FEP in Protein-Protein Interfaces. J Mol Biol 2019; 431:1481-1493. [PMID: 30776430 PMCID: PMC6453258 DOI: 10.1016/j.jmb.2019.02.003] [Citation(s) in RCA: 55] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2018] [Revised: 01/17/2019] [Accepted: 02/04/2019] [Indexed: 01/16/2023]
Abstract
Building on the substantial progress that has been made in using free energy perturbation (FEP) methods to predict the relative binding affinities of small molecule ligands to proteins, we have previously shown that results of similar quality can be obtained in predicting the effect of mutations on the binding affinity of protein–protein complexes. However, these results were restricted to mutations which did not change the net charge of the side chains due to known difficulties with modeling perturbations involving a change in charge in FEP. Various methods have been proposed to address this problem. Here we apply the co-alchemical water approach to study the efficacy of FEP calculations of charge changing mutations at the protein–protein interface for the antibody–gp120 system investigated previously and three additional complexes. We achieve an overall root mean square error of 1.2 kcal/mol on a set of 106 cases involving a change in net charge selected by a simple suitability filter using side-chain predictions and solvent accessible surface area to be relevant to a biologic optimization project. Reasonable, although less precise, results are also obtained for the 44 more challenging mutations that involve buried residues, which may in some cases require substantial reorganization of the local protein structure, which can extend beyond the scope of a typical FEP simulation. We believe that the proposed prediction protocol will be of sufficient efficiency and accuracy to guide protein engineering projects for which optimization and/or maintenance of a high degree of binding affinity is a key objective.
Collapse
Affiliation(s)
- Anthony J Clark
- Schrodinger Inc., 120 W 45th Street, New York, NY 10036, USA.
| | | | - Kevin Hauser
- Schrodinger Inc., 120 W 45th Street, New York, NY 10036, USA
| | - Mengzhen Sun
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3178, New York, NY 10027, USA
| | - Lingle Wang
- Schrodinger Inc., 120 W 45th Street, New York, NY 10036, USA
| | - Robert Abel
- Schrodinger Inc., 120 W 45th Street, New York, NY 10036, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3178, New York, NY 10027, USA
| |
Collapse
|
30
|
Roos K, Wu C, Damm W, Reboul M, Stevenson JM, Lu C, Dahlgren MK, Mondal S, Chen W, Wang L, Abel R, Friesner RA, Harder ED. OPLS3e: Extending Force Field Coverage for Drug-Like Small Molecules. J Chem Theory Comput 2019; 15:1863-1874. [PMID: 30768902 DOI: 10.1021/acs.jctc.8b01026] [Citation(s) in RCA: 584] [Impact Index Per Article: 116.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Katarina Roos
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
- Department of Cell and Molecular Biology, Uppsala University, Biomedical Centre, Box 596, SE-751 24 Uppsala, Sweden
| | - Chuanjie Wu
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Wolfgang Damm
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Mark Reboul
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - James M. Stevenson
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Chao Lu
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Markus K. Dahlgren
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Sayan Mondal
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Wei Chen
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Lingle Wang
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Robert Abel
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Edward D. Harder
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| |
Collapse
|
31
|
Abel R, Manas ES, Friesner RA, Farid RS, Wang L. Modeling the value of predictive affinity scoring in preclinical drug discovery. Curr Opin Struct Biol 2018; 52:103-110. [PMID: 30321805 DOI: 10.1016/j.sbi.2018.09.002] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2018] [Revised: 09/02/2018] [Accepted: 09/07/2018] [Indexed: 12/31/2022]
Abstract
Drug discovery is widely recognized to be a difficult and costly activity in large part due to the challenge of identifying chemical matter which simultaneously optimizes multiple properties, one of which is affinity for the primary biological target. Further, many of these properties are difficult to predict ahead of expensive and time-consuming compound synthesis and experimental testing. Here we highlight recent work to develop compound affinity prediction models, and extensively investigate the value such models may provide to preclinical drug discovery. We demonstrate that the ability of these models to improve the overall probability of success is crucially dependent on the shape of the error distribution, not just the root-mean-square error. In particular, while scoring more molecule ideas generally improves the probability of project success when the error distribution is Gaussian, fat-tail distributions such as a Cauchy distribution, can lead to a situation where scoring more ideas actually decreases the overall probability of success.
Collapse
Affiliation(s)
- Robert Abel
- Schrodinger, Inc., 120 West 45th Street, New York, NY 10036, United States.
| | - Eric S Manas
- GlaxoSmithKline, 1250 South Collegeville Road, Collegeville, PA 19426, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027, United States
| | - Ramy S Farid
- Schrodinger, Inc., 120 West 45th Street, New York, NY 10036, United States
| | - Lingle Wang
- Schrodinger, Inc., 120 West 45th Street, New York, NY 10036, United States
| |
Collapse
|
32
|
Shee J, Arthur EJ, Zhang S, Reichman DR, Friesner RA. Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units. J Chem Theory Comput 2018; 14:4109-4121. [DOI: 10.1021/acs.jctc.8b00342] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrödinger
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
33
|
Jacobson LD, Bochevarov AD, Watson MA, Hughes TF, Rinaldo D, Ehrlich S, Steinbrecher TB, Vaitheeswaran S, Philipp DM, Halls MD, Friesner RA. Automated Transition State Search and Its Application to Diverse Types of Organic Reactions. J Chem Theory Comput 2017; 13:5780-5797. [PMID: 28957627 DOI: 10.1021/acs.jctc.7b00764] [Citation(s) in RCA: 98] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Transition state search is at the center of multiple types of computational chemical predictions related to mechanistic investigations, reactivity and regioselectivity predictions, and catalyst design. The process of finding transition states in practice is, however, a laborious multistep operation that requires significant user involvement. Here, we report a highly automated workflow designed to locate transition states for a given elementary reaction with minimal setup overhead. The only essential inputs required from the user are the structures of the separated reactants and products. The seamless workflow combining computational technologies from the fields of cheminformatics, molecular mechanics, and quantum chemistry automatically finds the most probable correspondence between the atoms in the reactants and the products, generates a transition state guess, launches a transition state search through a combined approach involving the relaxing string method and the quadratic synchronous transit, and finally validates the transition state via the analysis of the reactive chemical bonds and imaginary vibrational frequencies as well as by the intrinsic reaction coordinate method. Our approach does not target any specific reaction type, nor does it depend on training data; instead, it is meant to be of general applicability for a wide variety of reaction types. The workflow is highly flexible, permitting modifications such as a choice of accuracy, level of theory, basis set, or solvation treatment. Successfully located transition states can be used for setting up transition state guesses in related reactions, saving computational time and increasing the probability of success. The utility and performance of the method are demonstrated in applications to transition state searches in reactions typical for organic chemistry, medicinal chemistry, and homogeneous catalysis research. In particular, applications of our code to Michael additions, hydrogen abstractions, Diels-Alder cycloadditions, carbene insertions, and an enzyme reaction model involving a molybdenum complex are shown and discussed.
Collapse
Affiliation(s)
- Leif D Jacobson
- Schrödinger, Inc. , 120 West 45th St., New York, New York 10036, United States
| | - Art D Bochevarov
- Schrödinger, Inc. , 120 West 45th St., New York, New York 10036, United States
| | - Mark A Watson
- Schrödinger, Inc. , 120 West 45th St., New York, New York 10036, United States
| | - Thomas F Hughes
- Schrödinger, Inc. , 120 West 45th St., New York, New York 10036, United States
| | - David Rinaldo
- Schrödinger GmbH , Dynamostrasse 13, D-68165 Mannheim, Germany
| | - Stephan Ehrlich
- Schrödinger GmbH , Dynamostrasse 13, D-68165 Mannheim, Germany
| | | | - S Vaitheeswaran
- Schrödinger, Inc. , 222 Third St., Suite 2230, Cambridge, Massachusetts 02142, United States
| | - Dean M Philipp
- Schrödinger, Inc. , 101 SW Main St., Suite 1300, Portland, Oregon 97204, United States
| | - Mathew D Halls
- Schrödinger, Inc. , 5820 Oberlin Dr., Suite 203, San Diego, California 92121, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University , 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
34
|
Abel R, Wang L, Mobley DL, Friesner RA. A Critical Review of Validation, Blind Testing, and Real- World Use of Alchemical Protein-Ligand Binding Free Energy Calculations. Curr Top Med Chem 2017; 17:2577-2585. [DOI: 10.2174/1568026617666170414142131] [Citation(s) in RCA: 62] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2016] [Revised: 01/16/2017] [Accepted: 01/17/2017] [Indexed: 11/22/2022]
Abstract
Protein-ligand binding is among the most fundamental phenomena underlying all molecular
biology, and a greater ability to more accurately and robustly predict the binding free energy of a small
molecule ligand for its cognate protein is expected to have vast consequences for improving the efficiency
of pharmaceutical drug discovery. We briefly reviewed a number of scientific and technical advances
that have enabled alchemical free energy calculations to recently emerge as a preferred approach,
and critically considered proper validation and effective use of these techniques. In particular, we characterized
a selection bias effect which may be important in prospective free energy calculations, and introduced
a strategy to improve the accuracy of the free energy predictions.
Collapse
Affiliation(s)
- Robert Abel
- Schrödinger, Inc., 120 West 45th Street, New York, NY 10036, United States
| | - Lingle Wang
- Schrödinger, Inc., 120 West 45th Street, New York, NY 10036, United States
| | - David L. Mobley
- Departments of Pharmaceutical Sciences and Chemistry, University of California, Irvine, CA 92697, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027, United States
| |
Collapse
|
35
|
Abstract
A principal goal of drug discovery project is to design molecules that can tightly and selectively bind to the target protein receptor. Accurate prediction of protein-ligand binding free energies is therefore of central importance in computational chemistry and computer aided drug design. Multiple recent improvements in computing power, classical force field accuracy, enhanced sampling methods, and simulation setup have enabled accurate and reliable calculations of protein-ligands binding free energies, and position free energy calculations to play a guiding role in small molecule drug discovery. In this Account, we outline the relevant methodological advances, including the REST2 (Replica Exchange with Solute Temperting) enhanced sampling, the incorporation of REST2 sampling with convential FEP (Free Energy Perturbation) through FEP/REST, the OPLS3 force field, and the advanced simulation setup that constitute our FEP+ approach, followed by the presentation of extensive comparisons with experiment, demonstrating sufficient accuracy in potency prediction (better than 1 kcal/mol) to substantially impact lead optimization campaigns. The limitations of the current FEP+ implementation and best practices in drug discovery applications are also discussed followed by the future methodology development plans to address those limitations. We then report results from a recent drug discovery project, in which several thousand FEP+ calculations were successfully deployed to simultaneously optimize potency, selectivity, and solubility, illustrating the power of the approach to solve challenging drug design problems. The capabilities of free energy calculations to accurately predict potency and selectivity have led to the advance of ongoing drug discovery projects, in challenging situations where alternative approaches would have great difficulties. The ability to effectively carry out projects evaluating tens of thousands, or hundreds of thousands, of proposed drug candidates, is potentially transformative in enabling hard to drug targets to be attacked, and in facilitating the development of superior compounds, in various dimensions, for a wide range of targets. More effective integration of FEP+ calculations into the drug discovery process will ensure that the results are deployed in an optimal fashion for yielding the best possible compounds entering the clinic; this is where the greatest payoff is in the exploitation of computer driven design capabilities. A key conclusion from the work described is the surprisingly robust and accurate results that are attainable within the conventional classical simulation, fixed charge paradigm. No doubt there are individual cases that would benefit from a more sophisticated energy model or dynamical treatment, and properties other than protein-ligand binding energies may be more sensitive to these approximations. We conclude that an inflection point in the ability of MD simulations to impact drug discovery has now been attained, due to the confluence of hardware and software development along with the formulation of "good enough" theoretical methods and models.
Collapse
Affiliation(s)
- Robert Abel
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Lingle Wang
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Edward D. Harder
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - B. J. Berne
- Department
of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department
of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
36
|
Shee J, Zhang S, Reichman DR, Friesner RA. Chemical Transformations Approaching Chemical Accuracy via Correlated Sampling in Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2017; 13:2667-2680. [DOI: 10.1021/acs.jctc.7b00224] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- James Shee
- Department
of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Shiwei Zhang
- Department
of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - David R. Reichman
- Department
of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department
of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
37
|
Abel R, Mondal S, Masse C, Greenwood J, Harriman G, Ashwell MA, Bhat S, Wester R, Frye L, Kapeller R, Friesner RA. Accelerating drug discovery through tight integration of expert molecular design and predictive scoring. Curr Opin Struct Biol 2017; 43:38-44. [DOI: 10.1016/j.sbi.2016.10.007] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2016] [Accepted: 10/07/2016] [Indexed: 01/08/2023]
|
38
|
Wang L, Deng Y, Wu Y, Kim B, LeBard DN, Wandschneider D, Beachy M, Friesner RA, Abel R. Accurate Modeling of Scaffold Hopping Transformations in Drug Discovery. J Chem Theory Comput 2016; 13:42-54. [PMID: 27933808 DOI: 10.1021/acs.jctc.6b00991] [Citation(s) in RCA: 88] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Lingle Wang
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Yuqing Deng
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Yujie Wu
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Byungchan Kim
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - David N. LeBard
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Dan Wandschneider
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Mike Beachy
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Richard A. Friesner
- Department
of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Robert Abel
- Schrödinger,
Inc., 120 West 45th Street, New York, New York 10036, United States
| |
Collapse
|
39
|
Clark AJ, Gindin T, Zhang B, Wang L, Abel R, Murret CS, Xu F, Bao A, Lu NJ, Zhou T, Kwong PD, Shapiro L, Honig B, Friesner RA. Free Energy Perturbation Calculation of Relative Binding Free Energy between Broadly Neutralizing Antibodies and the gp120 Glycoprotein of HIV-1. J Mol Biol 2016; 429:930-947. [PMID: 27908641 PMCID: PMC5383735 DOI: 10.1016/j.jmb.2016.11.021] [Citation(s) in RCA: 71] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2016] [Revised: 11/07/2016] [Accepted: 11/23/2016] [Indexed: 12/01/2022]
Abstract
Direct calculation of relative binding affinities between antibodies and antigens is a long-sought goal. However, despite substantial efforts, no generally applicable computational method has been described. Here, we describe a systematic free energy perturbation (FEP) protocol and calculate the binding affinities between the gp120 envelope glycoprotein of HIV-1 and three broadly neutralizing antibodies (bNAbs) of the VRC01 class. The protocol has been adapted from successful studies of small molecules to address the challenges associated with modeling protein–protein interactions. Specifically, we built homology models of the three antibody–gp120 complexes, extended the sampling times for large bulky residues, incorporated the modeling of glycans on the surface of gp120, and utilized continuum solvent-based loop prediction protocols to improve sampling. We present three experimental surface plasmon resonance data sets, in which antibody residues in the antibody/gp120 interface were systematically mutated to alanine. The RMS error in the large set (55 total cases) of FEP tests as compared to these experiments, 0.68 kcal/mol, is near experimental accuracy, and it compares favorably with the results obtained from a simpler, empirical methodology. The correlation coefficient for the combined data set including residues with glycan contacts, R2 = 0.49, should be sufficient to guide the choice of residues for antibody optimization projects, assuming that this level of accuracy can be realized in prospective prediction. More generally, these results are encouraging with regard to the possibility of using an FEP approach to calculate the magnitude of protein–protein binding affinities. Natural bNAbs against HIV-1 are a promising starting point for therapies. FEP provides a promising technique to computationally predict antibody structure optimizations. We show that FEP applied to bNAb mutations predicts binding affinity changes. FEP can provide an efficient method to identify potency-enhancing mutations to bNAbs.
Collapse
Affiliation(s)
- Anthony J Clark
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3178, New York, NY 10027, USA
| | - Tatyana Gindin
- Department of Pathology, Columbia University Medical Center, 630 W. 168th St, New York, NY 10032, USA
| | - Baoshan Zhang
- Vaccine Research Center, NIAID, National Institutes of Health, 40 Convent Drive, Bethesda, MD 20892, USA
| | - Lingle Wang
- Schrodinger Inc., 120 W 45th Street, New York, NY 10036, USA
| | - Robert Abel
- Schrodinger Inc., 120 W 45th Street, New York, NY 10036, USA
| | - Colleen S Murret
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3178, New York, NY 10027, USA
| | - Fang Xu
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3178, New York, NY 10027, USA
| | - Amy Bao
- Vaccine Research Center, NIAID, National Institutes of Health, 40 Convent Drive, Bethesda, MD 20892, USA
| | - Nina J Lu
- Vaccine Research Center, NIAID, National Institutes of Health, 40 Convent Drive, Bethesda, MD 20892, USA
| | - Tongqing Zhou
- Vaccine Research Center, NIAID, National Institutes of Health, 40 Convent Drive, Bethesda, MD 20892, USA
| | - Peter D Kwong
- Department of Biochemistry and Biophysics, Columbia University Medical Center, 701 West 168th Street, New York, NY 10032, USA; Vaccine Research Center, NIAID, National Institutes of Health, 40 Convent Drive, Bethesda, MD 20892, USA
| | - Lawrence Shapiro
- Department of Biochemistry and Biophysics, Columbia University Medical Center, 701 West 168th Street, New York, NY 10032, USA; Vaccine Research Center, NIAID, National Institutes of Health, 40 Convent Drive, Bethesda, MD 20892, USA
| | - Barry Honig
- Department of Biochemistry and Molecular Biophysics, Center for Computational Biology and Bioinformatics, Department of Systems Biology, Department of Medicine, Howard Hughes Medical Institute, Columbia University, 1130 Street Nicholas Avenue, Room 815, New York, NY 10032, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, MC 3178, New York, NY 10027, USA.
| |
Collapse
|
40
|
Abstract
The organizational structure is described for a new program that permits computa tions on a variety of quantum mechanical problems in chemical dynamics and spec troscopy. Particular attention was devoted to developing and using algorithms that exploit the capabilities of current vector supercomputers. A key component in this procedure is the recursive transformation of the large sparse Hamiltonian matrix into a much smaller tridiagonal matrix. An ap plication to time-dependent laser-mole cule energy transfer is presented. Rate of energy deposition in the multimode mole cule for systematic variations in the mo lecular intermode coupling parameters is emphasized.
Collapse
|
41
|
Clark AJ, Tiwary P, Borrelli K, Feng S, Miller EB, Abel R, Friesner RA, Berne BJ. Prediction of Protein-Ligand Binding Poses via a Combination of Induced Fit Docking and Metadynamics Simulations. J Chem Theory Comput 2016; 12:2990-8. [PMID: 27145262 DOI: 10.1021/acs.jctc.6b00201] [Citation(s) in RCA: 146] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ligand docking is a widely used tool for lead discovery and binding mode prediction based drug discovery. The greatest challenges in docking occur when the receptor significantly reorganizes upon small molecule binding, thereby requiring an induced fit docking (IFD) approach in which the receptor is allowed to move in order to bind to the ligand optimally. IFD methods have had some success but suffer from a lack of reliability. Complementing IFD with all-atom molecular dynamics (MD) is a straightforward solution in principle but not in practice due to the severe time scale limitations of MD. Here we introduce a metadynamics plus IFD strategy for accurate and reliable prediction of the structures of protein-ligand complexes at a practically useful computational cost. Our strategy allows treating this problem in full atomistic detail and in a computationally efficient manner and enhances the predictive power of IFD methods. We significantly increase the accuracy of the underlying IFD protocol across a large data set comprising 42 different ligand-receptor systems. We expect this approach to be of significant value in computationally driven drug design.
Collapse
Affiliation(s)
- Anthony J Clark
- Department of Chemistry, Columbia University , New York, New York 10027, United States
| | - Pratyush Tiwary
- Department of Chemistry, Columbia University , New York, New York 10027, United States
| | - Ken Borrelli
- Schrödinger, Inc. , 120 West 45th Street, New York, New York 10036, United States
| | - Shulu Feng
- Schrödinger, Inc. , 120 West 45th Street, New York, New York 10036, United States
| | - Edward B Miller
- Schrödinger, Inc. , 120 West 45th Street, New York, New York 10036, United States
| | - Robert Abel
- Schrödinger, Inc. , 120 West 45th Street, New York, New York 10036, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University , New York, New York 10027, United States
| | - B J Berne
- Department of Chemistry, Columbia University , New York, New York 10027, United States
| |
Collapse
|
42
|
Murphy RB, Repasky MP, Greenwood JR, Tubert-Brohman I, Jerome S, Annabhimoju R, Boyles NA, Schmitz CD, Abel R, Farid R, Friesner RA. WScore: A Flexible and Accurate Treatment of Explicit Water Molecules in Ligand–Receptor Docking. J Med Chem 2016; 59:4364-84. [DOI: 10.1021/acs.jmedchem.6b00131] [Citation(s) in RCA: 58] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Robert B. Murphy
- Schrödinger, Inc., 101 SW Main Street, Portland Oregon 97204, United States
| | - Matthew P. Repasky
- Schrödinger, Inc., 101 SW Main Street, Portland Oregon 97204, United States
| | - Jeremy R. Greenwood
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Ivan Tubert-Brohman
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Steven Jerome
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | | | - Nicholas A. Boyles
- Schrödinger, Inc., 101 SW Main Street, Portland Oregon 97204, United States
| | | | - Robert Abel
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Ramy Farid
- Schrödinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Richard A. Friesner
- Department of
Chemistry, Columbia University, New York, 3000 Broadway,
MC 3110, New York 10036, United States
| |
Collapse
|
43
|
Cao Y, Hughes T, Giesen D, Halls MD, Goldberg A, Vadicherla TR, Sastry M, Patel B, Sherman W, Weisman AL, Friesner RA. Highly efficient implementation of pseudospectral time-dependent density-functional theory for the calculation of excitation energies of large molecules. J Comput Chem 2016; 37:1425-41. [DOI: 10.1002/jcc.24350] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2015] [Revised: 02/02/2016] [Accepted: 02/03/2016] [Indexed: 12/25/2022]
Affiliation(s)
- Yixiang Cao
- Schrödinger Inc; 120 West 45th Street, Tower 45, 17th Floor New York New York 10036
| | - Thomas Hughes
- Schrödinger Inc; 120 West 45th Street, Tower 45, 17th Floor New York New York 10036
| | - Dave Giesen
- Schrödinger Inc; 120 West 45th Street, Tower 45, 17th Floor New York New York 10036
| | - Mathew D. Halls
- Schrödinger Inc; 120 West 45th Street, Tower 45, 17th Floor New York New York 10036
| | - Alexander Goldberg
- Schrödinger Inc; 8910 University Lane, Suite 270 San Diego California 92122
| | - Tati Reddy Vadicherla
- Schrödinger; Plot No. 573, B & C, Road No. 1 Jubilee Hills, Hyderabad Telangana 5000096 India
| | - Madhavi Sastry
- Schrödinger; Sanali Infopark, 8-2-120/113 Banjara Hills, Hyderabad Andhra Pradesh 500034 India
| | - Bhargav Patel
- Schrödinger; Sanali Infopark, 8-2-120/113 Banjara Hills, Hyderabad Andhra Pradesh 500034 India
| | - Woody Sherman
- Schrödinger Inc; 245 First St. Riverview II, 18th Floor Cambridge Massachusetts 02142
| | - Andrew L. Weisman
- Department of Applied Physics and Applied Mathematics; Columbia University; New York New York 10027
| | | |
Collapse
|
44
|
Coskun D, Jerome SV, Friesner RA. Evaluation of the Performance of the B3LYP, PBE0, and M06 DFT Functionals, and DBLOC-Corrected Versions, in the Calculation of Redox Potentials and Spin Splittings for Transition Metal Containing Systems. J Chem Theory Comput 2016; 12:1121-8. [DOI: 10.1021/acs.jctc.5b00782] [Citation(s) in RCA: 51] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Dilek Coskun
- Department
of Chemistry, Columbia University, New York, New York 10027, United States
| | | | - Richard A. Friesner
- Department
of Chemistry, Columbia University, New York, New York 10027, United States
| |
Collapse
|
45
|
Jerome SV, Hughes TF, Friesner RA. Successful application of the DBLOC method to the hydroxylation of camphor by cytochrome p450. Protein Sci 2016; 25:277-85. [PMID: 26441133 PMCID: PMC4815313 DOI: 10.1002/pro.2819] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2015] [Revised: 10/01/2015] [Accepted: 10/02/2015] [Indexed: 01/19/2023]
Abstract
The activation barrier for the hydroxylation of camphor by cytochrome P450 was computed using a mixed quantum mechanics/molecular mechanics (QM/MM) model of the full protein-ligand system and a fully QM calculation using a cluster model of the active site at the B3LYP/LACVP*/LACV3P** level of theory, which consisted of B3LYP/LACV3P** single point energies computed at B3LYP/LACVP* optimized geometries. From the QM/MM calculation, a barrier height of 17.5 kcal/mol was obtained, while the experimental value was known to be less than or equal to 10 kcal/mol. This process was repeated using the D3 correction for hybrid DFT in order to investigate whether the inadequate treatment of dispersion interaction was responsible for the overestimation of the barrier. While the D3 correction does reduce the computed barrier to 13.3 kcal/mol, it was still in disagreement with experiment. After application of a series of transition metal optimized localized orbital corrections (DBLOC) and without any refitting of parameters, the barrier was further reduced to 10.0 kcal/mol, which was consistent with the experimental results. The DBLOC method to CH bond activation in methane monooxygenase (MMO) was also applied, as a second, independent test. The barrier in MMO was known, by experiment, to be 15.4 kcal/mol. After application of the DBLOC corrections to the MMO barrier compute by B3LYP, in a previous study, and accounting for dispersion with Grimme's D3 method, the unsigned deviation from experiment was improved from 3.2 to 2.3 kcal/mol. These results suggested that the combination of dispersion plus localized orbital corrections could yield significant quantitative improvements in modeling the catalytic chemistry of transition-metal containing enzymes, within the limitations of the statistical errors of the model, which appear to be on the order of approximately 2 kcal/mole.
Collapse
Affiliation(s)
- Steven V Jerome
- Department of Chemistry, Columbia University, New York, New York, 10027
| | | | | |
Collapse
|
46
|
Knight JL, Krilov G, Borrelli KW, Williams J, Gunn JR, Clowes A, Cheng L, Friesner RA, Abel R. Leveraging Data Fusion Strategies in Multireceptor Lead Optimization MM/GBSA End-Point Methods. J Chem Theory Comput 2015; 10:3207-20. [PMID: 26588291 DOI: 10.1021/ct500189s] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Accurate and efficient affinity calculations are critical to enhancing the contribution of in silico modeling during the lead optimization phase of a drug discovery campaign. Here, we present a large-scale study of the efficacy of data fusion strategies to leverage results from end-point MM/GBSA calculations in multiple receptors to identify potent inhibitors among an ensemble of congeneric ligands. The retrospective analysis of 13 congeneric ligand series curated from publicly available data across seven biological targets demonstrates that in 90% of the individual receptor structures MM/GBSA scores successfully identify subsets of inhibitors that are more potent than a random selection, and data fusion strategies that combine MM/GBSA scores from each of the receptors significantly increase the robustness of the predictions. Among nine different data fusion metrics based on consensus scores or receptor rankings, the SumZScore (i.e., converting MM/GBSA scores into standardized Z-Scores within a receptor and computing the sum of the Z-Scores for a given ligand across the ensemble of receptors) is found to be a robust and physically meaningful metric for combining results across multiple receptors. Perhaps most surprisingly, even with relatively low to modest overall correlations between SumZScore and experimental binding affinities, SumZScore tends to reliably prioritize subsets of inhibitors that are at least as potent as those that are prioritized from a "best" single receptor identified from known compounds within the congeneric series.
Collapse
Affiliation(s)
- Jennifer L Knight
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Goran Krilov
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Kenneth W Borrelli
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Joshua Williams
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - John R Gunn
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Alec Clowes
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Luciano Cheng
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| | - Richard A Friesner
- Columbia University , Department of Chemistry, 3000 Broadway, MC 3110, New York, New York 10027, United States
| | - Robert Abel
- Schrödinger, 120 West 45th Street, 17th Floor, Tower 45, New York, New York 10036-4041, United States
| |
Collapse
|
47
|
Zhu K, Shirts MR, Friesner RA, Jacobson MP. Multiscale Optimization of a Truncated Newton Minimization Algorithm and Application to Proteins and Protein-Ligand Complexes. J Chem Theory Comput 2015; 3:640-8. [PMID: 26637042 DOI: 10.1021/ct600129f] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We optimize a truncated Newton (TN) minimization algorithm and computer package, TNPACK, developed for macromolecular minimizations by applying multiscale methods, analogous to those used in molecular dynamics (e.g., r-RESPA). The molecular mechanics forces are divided into short- and long-range components, with the long-range forces updated only intermittently in the iterative evaluations. This algorithm, which we refer to as MSTN, is implemented as a modification to the TNPACK package and is tested on energy minimizations of protein loops, entire proteins, and protein-ligand complexes and compared with the unmodified truncated Newton algorithm, a quasi-Newton algorithm (LBFGS), and a conjugate gradient algorithm (CG+). In vacuum minimizations, the speedup of MSTN relative to the unmodified TN algorithm (TNPACK) depends on system size and the distance cutoffs used for defining the short- and long-range interactions and the long-range force updating frequency, but it is 4 to 5 times greater in the work reported here. This algorithm works best for the minimization of small portions of a protein and shows some degradation (speedup factor of 2-3) for the minimization of entire proteins. The MSTN algorithm is faster than the quasi-Newton and conjugate gradient algorithms by approximately 1 order of magnitude. We also present a modification of the algorithm which permits minimizations with a generalized Born implicit solvent model, using a self-consistent procedure that increases the computational expense, relative to a vacuum, by only a small factor (∼3-4).
Collapse
Affiliation(s)
- Kai Zhu
- Department of Chemistry and Center for Biomolecular Simulation, Columbia University, New York, New York 10027, and Department of Pharmaceutical Chemistry, University of California, San Francisco, California 94158-2517
| | - Michael R Shirts
- Department of Chemistry and Center for Biomolecular Simulation, Columbia University, New York, New York 10027, and Department of Pharmaceutical Chemistry, University of California, San Francisco, California 94158-2517
| | - Richard A Friesner
- Department of Chemistry and Center for Biomolecular Simulation, Columbia University, New York, New York 10027, and Department of Pharmaceutical Chemistry, University of California, San Francisco, California 94158-2517
| | - Matthew P Jacobson
- Department of Chemistry and Center for Biomolecular Simulation, Columbia University, New York, New York 10027, and Department of Pharmaceutical Chemistry, University of California, San Francisco, California 94158-2517
| |
Collapse
|
48
|
Harder E, Damm W, Maple J, Wu C, Reboul M, Xiang JY, Wang L, Lupyan D, Dahlgren MK, Knight JL, Kaus JW, Cerutti DS, Krilov G, Jorgensen WL, Abel R, Friesner RA. OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins. J Chem Theory Comput 2015; 12:281-96. [PMID: 26584231 DOI: 10.1021/acs.jctc.5b00864] [Citation(s) in RCA: 1992] [Impact Index Per Article: 221.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The parametrization and validation of the OPLS3 force field for small molecules and proteins are reported. Enhancements with respect to the previous version (OPLS2.1) include the addition of off-atom charge sites to represent halogen bonding and aryl nitrogen lone pairs as well as a complete refit of peptide dihedral parameters to better model the native structure of proteins. To adequately cover medicinal chemical space, OPLS3 employs over an order of magnitude more reference data and associated parameter types relative to other commonly used small molecule force fields (e.g., MMFF and OPLS_2005). As a consequence, OPLS3 achieves a high level of accuracy across performance benchmarks that assess small molecule conformational propensities and solvation. The newly fitted peptide dihedrals lead to significant improvements in the representation of secondary structure elements in simulated peptides and native structure stability over a number of proteins. Together, the improvements made to both the small molecule and protein force field lead to a high level of accuracy in predicting protein-ligand binding measured over a wide range of targets and ligands (less than 1 kcal/mol RMS error) representing a 30% improvement over earlier variants of the OPLS force field.
Collapse
Affiliation(s)
- Edward Harder
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Wolfgang Damm
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Jon Maple
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Chuanjie Wu
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Mark Reboul
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Jin Yu Xiang
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Lingle Wang
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Dmitry Lupyan
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Markus K Dahlgren
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Jennifer L Knight
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Joseph W Kaus
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - David S Cerutti
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Goran Krilov
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - William L Jorgensen
- Department of Chemistry, Yale University , New Haven, Connecticut 06520, United States
| | - Robert Abel
- Schrodinger, Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University , 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
49
|
Wang L, Wu Y, Deng Y, Kim B, Pierce L, Krilov G, Lupyan D, Robinson S, Dahlgren MK, Greenwood J, Romero DL, Masse C, Knight JL, Steinbrecher T, Beuming T, Damm W, Harder E, Sherman W, Brewer M, Wester R, Murcko M, Frye L, Farid R, Lin T, Mobley DL, Jorgensen WL, Berne BJ, Friesner RA, Abel R. Accurate and reliable prediction of relative ligand binding potency in prospective drug discovery by way of a modern free-energy calculation protocol and force field. J Am Chem Soc 2015; 137:2695-703. [PMID: 25625324 DOI: 10.1021/ja512751q] [Citation(s) in RCA: 770] [Impact Index Per Article: 85.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Designing tight-binding ligands is a primary objective of small-molecule drug discovery. Over the past few decades, free-energy calculations have benefited from improved force fields and sampling algorithms, as well as the advent of low-cost parallel computing. However, it has proven to be challenging to reliably achieve the level of accuracy that would be needed to guide lead optimization (∼5× in binding affinity) for a wide range of ligands and protein targets. Not surprisingly, widespread commercial application of free-energy simulations has been limited due to the lack of large-scale validation coupled with the technical challenges traditionally associated with running these types of calculations. Here, we report an approach that achieves an unprecedented level of accuracy across a broad range of target classes and ligands, with retrospective results encompassing 200 ligands and a wide variety of chemical perturbations, many of which involve significant changes in ligand chemical structures. In addition, we have applied the method in prospective drug discovery projects and found a significant improvement in the quality of the compounds synthesized that have been predicted to be potent. Compounds predicted to be potent by this approach have a substantial reduction in false positives relative to compounds synthesized on the basis of other computational or medicinal chemistry approaches. Furthermore, the results are consistent with those obtained from our retrospective studies, demonstrating the robustness and broad range of applicability of this approach, which can be used to drive decisions in lead optimization.
Collapse
Affiliation(s)
- Lingle Wang
- Schrödinger, Inc. , 120 West 45th Street, New York, New York 10036, United States
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
50
|
Stafford KA, Trbovic N, Butterwick JA, Abel R, Friesner RA, Palmer AG. Conformational preferences underlying reduced activity of a thermophilic ribonuclease H. J Mol Biol 2014; 427:853-866. [PMID: 25550198 DOI: 10.1016/j.jmb.2014.11.023] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2014] [Accepted: 11/14/2014] [Indexed: 11/30/2022]
Abstract
The conformational basis for reduced activity of the thermophilic ribonuclease HI enzyme from Thermus thermophilus, compared to its mesophilic homolog from Escherichia coli, is elucidated using a combination of NMR spectroscopy and molecular dynamics (MD) simulations. Explicit-solvent all-atom MD simulations of the two wild-type proteins and an E. coli mutant in which a glycine residue is inserted after position 80 to mimic the T. thermophilus protein reproduce the differences in conformational dynamics determined from (15)N spin-relaxation NMR spectroscopy of three loop regions that surround the active site and contain functionally important residues: the glycine-rich region, the handle region, and the β5/αE loop. Examination of the MD trajectories indicates that the thermophilic protein samples conformations productive for substrate binding and activity less frequently than the mesophilic enzyme, although these differences may manifest as either increased or decreased relative flexibility of the different regions. Additional MD simulations indicate that mutations increasing activity of the T. thermophilus enzyme at mesophilic temperatures do so by reconfiguring the local environments of the mutated sites to more closely resemble active conformations. Taken together, the results show that both locally increased and decreased flexibility contribute to an overall reduction in activity of T. thermophilus ribonuclease H compared to its mesophilic E. coli homolog.
Collapse
Affiliation(s)
- Kate A Stafford
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA
| | - Nikola Trbovic
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA
| | - Joel A Butterwick
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA
| | - Robert Abel
- Department of Chemistry, Columbia University, New York, NY 10027, USA
| | | | - Arthur G Palmer
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York, NY 10032, USA.
| |
Collapse
|