1
|
Mous S, Poitevin F, Hunter MS, Asthagiri DN, Beck TL. Structural biology in the age of X-ray free-electron lasers and exascale computing. Curr Opin Struct Biol 2024; 86:102808. [PMID: 38547555 DOI: 10.1016/j.sbi.2024.102808] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2023] [Revised: 02/07/2024] [Accepted: 03/07/2024] [Indexed: 05/19/2024]
Abstract
Serial femtosecond X-ray crystallography has emerged as a powerful method for investigating biomolecular structure and dynamics. With the new generation of X-ray free-electron lasers, which generate ultrabright X-ray pulses at megahertz repetition rates, we can now rapidly probe ultrafast conformational changes and charge movement in biomolecules. Over the last year, another innovation has been the deployment of Frontier, the world's first exascale supercomputer. Synergizing extremely high repetition rate X-ray light sources and exascale computing has the potential to accelerate discovery in biomolecular sciences. Here we outline our perspective on each of these remarkable innovations individually, and the opportunities and challenges in yoking them within an integrated research infrastructure.
Collapse
Affiliation(s)
- Sandra Mous
- Linac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, 94025, CA, USA
| | - Frédéric Poitevin
- Linac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, 94025, CA, USA
| | - Mark S Hunter
- Linac Coherent Light Source, SLAC National Accelerator Laboratory, Menlo Park, 94025, CA, USA.
| | - Dilipkumar N Asthagiri
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, 37830-6012, TN, USA
| | - Thomas L Beck
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, 37830-6012, TN, USA.
| |
Collapse
|
2
|
Beck TL, Carloni P, Asthagiri DN. All-Atom Biomolecular Simulation in the Exascale Era. J Chem Theory Comput 2024; 20:1777-1782. [PMID: 38382017 DOI: 10.1021/acs.jctc.3c01276] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
Exascale supercomputers have opened the door to dynamic simulations, facilitated by AI/ML techniques, that model biomolecular motions over unprecedented length and time scales. This new capability holds the potential to revolutionize our understanding of fundamental biological processes. Here we report on some of the major advances that were discussed at a recent CECAM workshop in Pisa, Italy, on the topic with a primary focus on atomic-level simulations. First, we highlight examples of current large-scale biomolecular simulations and the future possibilities enabled by crossing the exascale threshold. Next, we discuss challenges to be overcome in optimizing the usage of these powerful resources. Finally, we close by listing several grand challenge problems that could be investigated with this new computer architecture.
Collapse
Affiliation(s)
- Thomas L Beck
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States
| | - Paolo Carloni
- INM-9/IAS-5 Computational Biomedicine, Forschungszentrum Jülich, Wilhelm-Johnen-Straße, D-54245 Jülich, Germany
- Department of Physics, RWTH Aachen University, D-52078 Aachen, Germany
| | - Dilipkumar N Asthagiri
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, United States
| |
Collapse
|
3
|
Clark JA, Prabhu VM, Douglas JF. Molecular Dynamics Simulation of the Influence of Temperature and Salt on the Dynamic Hydration Layer in a Model Polyzwitterionic Polymer PAEDAPS. J Phys Chem B 2023; 127:8185-8198. [PMID: 37668318 PMCID: PMC10578162 DOI: 10.1021/acs.jpcb.3c03654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/06/2023]
Abstract
We investigate the hydration of poly(3-[2-(acrylamido) ethyldimethylammonio] propanesulfonate) over a range of temperatures in pure water and with the inclusion of 0.1 mol/L NaCl using atomistic molecular dynamics simulation. Drawing on concepts drawn from the field of glass-forming liquids, we use the Debye-Waller parameter () for describing the water mobility gradient around the polybetaine backbone extending to an overall distance ≈18 Å. The water mobility in this layer is defined through the mean-square water molecule displacement at a time on the order of water's β-relaxation time. The brushlike topology of polybetaines leads to two regions in the dynamic hydration layer. The inner region of ≈10.5 Å is explored by pendant group conformational motions, and the outer region of ≈7.5 Å represents an extended layer of reduced water mobility relative to bulk water. The dynamic hydration layer extends far beyond the static hydration layer, adjacent to the polymer.
Collapse
Affiliation(s)
- Jennifer A. Clark
- Materials Science and Engineering Division, Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, Maryland 20899, United States
| | - Vivek M. Prabhu
- Materials Science and Engineering Division, Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, Maryland 20899, United States
| | - Jack F. Douglas
- Materials Science and Engineering Division, Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, Maryland 20899, United States
| |
Collapse
|
4
|
Herman CE, Valiya Parambathu A, Asthagiri DN, Lenhoff AM. Polarizability Plays a Decisive Role in Modulating Association between Molecular Cations and Anions. J Phys Chem Lett 2023; 14:7020-7026. [PMID: 37523856 DOI: 10.1021/acs.jpclett.3c01566] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
Electrostatic interactions involving proteins depend on not only the ionic charges involved but also their chemical identities. Here we examine the origins of incompletely understood differences in the strength of association of different pairs of monovalent molecular ions that are relevant to protein-protein and protein-ligand interactions. Cationic analogues of the basic amino acid side chains are simulated, along with oxyanionic analogues of cation-exchange ligands and acidic amino acids. Experimentally observed association trends with respect to the cations, but not anions, are captured by a nonpolarizable model. An effective continuum correction to account for electronic polarizability can capture both trends better but at the expense of fidelity to the underlying free energy landscape for ion-pair association. A polarizable model proves decisive in capturing experimentally suggested trends with respect to both cations and anions; critically, the free energy landscape for ion-pair association is itself altered, thus altering configurational sampling.
Collapse
Affiliation(s)
- Chase E Herman
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, Delaware 19716, United States
| | - Arjun Valiya Parambathu
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, Delaware 19716, United States
| | - Dilipkumar N Asthagiri
- Oak Ridge National Laboratory, 1 Bethel Valley Road, Oak Ridge, Tennessee 37830, United States
| | - Abraham M Lenhoff
- Department of Chemical and Biomolecular Engineering, University of Delaware, 150 Academy Street, Newark, Delaware 19716, United States
| |
Collapse
|
5
|
Adhikari RS, Parambathu AV, Chapman WG, Asthagiri DN. Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory. J Phys Chem B 2022; 126:9607-9616. [PMID: 36354351 DOI: 10.1021/acs.jpcb.2c05725] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)─an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration─physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.
Collapse
Affiliation(s)
- Rohan S Adhikari
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas77005, United States
| | - Arjun Valiya Parambathu
- Department of Chemical and Biomolecular Engineering, University of Delaware, Newark, Delaware19711, United States
| | - Walter G Chapman
- Department of Chemical and Biomolecular Engineering, Rice University, Houston, Texas77005, United States
| | | |
Collapse
|
6
|
Gomez DT, Pratt LR, Asthagiri DN, Rempe SB. Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory. Acc Chem Res 2022; 55:2201-2212. [PMID: 35829622 PMCID: PMC9386901 DOI: 10.1021/acs.accounts.2c00078] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here, we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F- over Cl- in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory-experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.
Collapse
Affiliation(s)
- Diego T. Gomez
- Department
of Chemical & Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States,
| | - Lawrence R. Pratt
- Department
of Chemical & Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States,
| | - Dilipkumar N. Asthagiri
- Department
of Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005, United States,
| | - Susan B. Rempe
- Center
for Integrated Nanotechnologies, Sandia
National Laboratories, Albuquerque, New Mexico 87185, United States,
| |
Collapse
|
7
|
Car-Parrinello molecular dynamics study of CuF, AgF, CuPF6 and AgPF6 in acetonitrile solvent and Cluster-Continuum calculation of the solvation free energy of Cu(I), Ag(I) and Li(I). J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119368] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
|
8
|
Pliego JR. Copper-Catalyzed Aromatic Fluorination of 2-(2-bromophenyl)pyridine via Cu(I)/Cu(III) Mechanism in Acetonitrile Solvent: Cluster-Continuum Free Energy Profile and Microkinetic Analysis. MOLECULAR CATALYSIS 2022. [DOI: 10.1016/j.mcat.2022.112560] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
9
|
Milovanović B, Petković M, Popov I, Etinski M. Water-Mediated Interactions Enhance Alkaline Earth Cation Chelation in Neighboring Cavities of a Cytosine Quartet in the DNA Quadruplex. J Phys Chem B 2021; 125:11996-12005. [PMID: 34694801 DOI: 10.1021/acs.jpcb.1c05598] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Larger Coulombic repulsion between divalent cations compared to the monovalent counterparts dictates the cation-cation distance in the central ion channel of quadruplexes. In this work, density functional theory and a continuum solvation model were employed to study bond energies of alkaline earth cations in adjacent cavities of the central ion channel. Four crystallized tetramolecular quadruplexes with various geometric constraints and structural motifs available in the Protein Data Bank were examined in order to understand how the cation binding affinities could be increased in aqueous solution. A cytosine quartet sandwiched between guanine quartets has a larger bond energy of the second alkaline earth cation in comparison with guanine and uracil quartets. Four highly conserved hydrogen-bonded water molecules in the center of the cytosine quartet are responsible for a higher electrostatic interaction with the cations in comparison with guanines' carbonyl groups. The reported findings are valuable for the design of synthetic quadruplexes templated with divalent cations for optoelectronic applications.
Collapse
Affiliation(s)
- Branislav Milovanović
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16 11000 Belgrade, Serbia
| | - Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16 11000 Belgrade, Serbia
| | - Igor Popov
- Institute for Multidisciplinary Research, University of Belgrade, Kneza Višeslava 1, 11030 Belgrade, Serbia.,Institut of Physics, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16 11000 Belgrade, Serbia
| |
Collapse
|