1
|
Degünther J, van der Meer J, Seifert U. General theory for localizing the where and when of entropy production meets single-molecule experiments. Proc Natl Acad Sci U S A 2024; 121:e2405371121. [PMID: 39121164 PMCID: PMC11331124 DOI: 10.1073/pnas.2405371121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2024] [Accepted: 07/11/2024] [Indexed: 08/11/2024] Open
Abstract
The laws of thermodynamics apply to biophysical systems on the nanoscale as described by the framework of stochastic thermodynamics. This theory provides universal, exact relations for quantities like work, which have been verified in experiments where a fully resolved description allows direct access to such quantities. Complementary studies consider partially hidden, coarse-grained descriptions, in which the mean entropy production typically is not directly accessible but can be bounded in terms of observable quantities. Going beyond the mean, we introduce a fluctuating entropy production that applies to individual trajectories in a coarse-grained description under time-dependent driving. Thus, this concept is applicable to the broad and experimentally significant class of driven systems in which not all relevant states can be resolved. We provide a paradigmatic example by studying an experimentally verified protein unfolding process. As a consequence, the entire distribution of the coarse-grained entropy production rather than merely its mean retains spatial and temporal information about the microscopic process. In particular, we obtain a bound on the distribution of the physical entropy production of individual unfolding events.
Collapse
Affiliation(s)
- Julius Degünther
- II. Institut für Theoretische Physik, Universität Stuttgart, Stuttgart70550, Germany
| | - Jann van der Meer
- II. Institut für Theoretische Physik, Universität Stuttgart, Stuttgart70550, Germany
| | - Udo Seifert
- II. Institut für Theoretische Physik, Universität Stuttgart, Stuttgart70550, Germany
| |
Collapse
|
2
|
Okamoto R, Akiyama H, Nakae R, Tanaka Y, Washizu H. Tribo-Chemical Reactions in DLC-Zirconia Sliding under Ethanol Gas Environment. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2024; 40:14953-14963. [PMID: 39004826 DOI: 10.1021/acs.langmuir.4c01114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/16/2024]
Abstract
A reactive force field (ReaxFF) molecular dynamics simulation is performed for the sliding of diamond-like carbon (DLC) and yttria-stabilized zirconia (YSZ) under an ethanol gas environment, motivated by the previous experiment of ultralow friction phenomenon (friction fade-out). We observe (i) dissociation of ethanol into ethoxy and hydrogen, both of which simultaneously adsorb on the YSZ surface, and (ii) dissociation of ethanol into ethyl and hydroxy, the former of which forms a bond with another ethanol molecule and the latter of which adsorbs on the DLC surface. Reaction (i) is enhanced by the sliding motion, but occurs even without it, while reaction (ii) only occurs during sliding with a sufficiently high load pressure. The potentials of mean force for the two reactions are also calculated combining the steered MD and Jarzynski equality. It is shown that the activation energies of reactions (i) and (ii) are significantly lowered by the YSZ and DLC surfaces, respectively, as compared to those in a vacuum. The resultant activation energy is higher for reaction (ii) than for reaction (i).
Collapse
Affiliation(s)
- Ryuichi Okamoto
- Graduate School of Information Science, University of Hyogo, Kobe 650-0047, Japan
| | - Hirotoshi Akiyama
- Graduate School of Information Science, University of Hyogo, Kobe 650-0047, Japan
| | - Rio Nakae
- Graduate School of Information Science, University of Hyogo, Kobe 650-0047, Japan
| | - Yudai Tanaka
- Graduate School of Information Science, University of Hyogo, Kobe 650-0047, Japan
| | - Hitoshi Washizu
- Graduate School of Information Science, University of Hyogo, Kobe 650-0047, Japan
| |
Collapse
|
3
|
Cuong NM, Khanh PN, Nhung LTH, Ha NX, Huong TT, Bauerova K, Kim YH, Tung DD, Thuy TT, Anh NTH. Acetylcholinesterase inhibitory activities of some flavonoids from the root bark of Pinus krempfii Lecomte: in vitro and in silico study. J Biomol Struct Dyn 2024; 42:4888-4901. [PMID: 37325850 DOI: 10.1080/07391102.2023.2223664] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2023] [Accepted: 06/05/2023] [Indexed: 06/17/2023]
Abstract
From the root bark of Pinus krempfii Lecomte, four flavonoids were isolated and evaluated for their inhibitory activities against AChE and BChE enzymes in vitro and in silico. Tectochrysin (1) was found to inhibit AChE with an IC50 value of 33.69 ± 2.80 μM. The docking study results also showed agreement with the in vitro test results. All four compounds also showed the best binding affinity for the AChE enzyme, characterised by binding energy (ΔG) values as low as -8.1 to -9.3 kcal/mol, in which, the compound tectochrysin had the best binding affinity for the AChE protein with a ΔG value of -9.329 kcal/mol. Tectochrysin (1) was also bound to the amino acid Phe295 of AChE with a bond length of 2.8 Å, similar to the control dihydrotanshinone-I. Galangin (2) also showed its in vitro inhibitory activity against BChE with an IC50 value of 82.21 ± 2.70 μM. In silico, it also had the best binding energy value of -9.072 kcal/mol with BChE and formed hydrogen bonds with the His438 (2.85 Å) residues of BChE like the positive control (tacrine). The steered molecular dynamics (SMD) simulation results of these two complexes revealed a mechanistic insight that the protein-ligand complexes showed stable trajectories throughout the 20 and 150 ns simulations. Moreover, the drug likeliness suggested that both flavonoids (1 and 2) were expected to be drug-like and have an LD50 toxicity level of 5. This study has contributed new results for drug discovery and the development of substances with neuroprotective effects, especially for the treatment of Alzheimer's disease.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Nguyen Manh Cuong
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Graduate University of Sciences and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Pham Ngoc Khanh
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Graduate University of Sciences and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Le Thi Hong Nhung
- Faculty of Chemical Technology, Hanoi University of Industry, Bac Tu Liem District, Hanoi, Vietnam
| | - Nguyen Xuan Ha
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Tran Thu Huong
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Katarina Bauerova
- Centre of Experimental Medicine, Institute of Experimental Pharmacology and Toxicology, Slovak Academy of Sciences, Karlova Ves, Slovakia
| | - Young Ho Kim
- College of Pharmacy, Chungnam National University, Daejeon, South Korea
| | | | - Trinh Thi Thuy
- Graduate University of Sciences and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Institute of Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Nguyen Thi Hoang Anh
- Graduate University of Sciences and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Institute of Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| |
Collapse
|
4
|
Truong DT, Ho K, Nhi HTY, Nguyen VH, Dang TT, Nguyen MT. Imidazole[1,5-a]pyridine derivatives as EGFR tyrosine kinase inhibitors unraveled by umbrella sampling and steered molecular dynamics simulations. Sci Rep 2024; 14:12218. [PMID: 38806555 PMCID: PMC11133355 DOI: 10.1038/s41598-024-62743-3] [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: 03/07/2024] [Accepted: 05/21/2024] [Indexed: 05/30/2024] Open
Abstract
Although the use of the tyrosine kinase inhibitors (TKIs) has been proved that it can save live in a cancer treatment, the currently used drugs bring in many undesirable side-effects. Therefore, the search for new drugs and an evaluation of their efficiency are intensively carried out. Recently, a series of eighteen imidazole[1,5-a]pyridine derivatives were synthetized by us, and preliminary analyses pointed out their potential to be an important platform for pharmaceutical development owing to their promising actions as anticancer agents and enzyme (kinase, HIV-protease,…) inhibitors. In the present theoretical study, we further analyzed their efficiency in using a realistic scenario of computational drug design. Our protocol has been developed to not only observe the atomistic interaction between the EGFR protein and our 18 novel compounds using both umbrella sampling and steered molecular dynamics simulations, but also determine their absolute binding free energies. Calculated properties of the 18 novel compounds were in detail compared with those of two known drugs, erlotinib and osimertinib, currently used in cancer treatment. Inspiringly the simulation results promote three imidazole[1,5-a]pyridine derivatives as promising inhibitors into a further step of clinical trials.
Collapse
Affiliation(s)
- Duc Toan Truong
- Laboratory for Chemical Computation and Modeling, Institute for Computational Science and Artificial Intelligence, Van Lang University, Ho Chi Minh City, 70000, Vietnam
- Faculty of Applied Technology, School of Technology, Van Lang University, Ho Chi Minh City, 70000, Vietnam
| | - Kiet Ho
- Institute for Computational Science and Technology (ICST), Quang Trung Software City, Ho Chi Minh City, 70000, Vietnam
| | - Huynh Thi Yen Nhi
- Laboratory for Chemical Computation and Modeling, Institute for Computational Science and Artificial Intelligence, Van Lang University, Ho Chi Minh City, 70000, Vietnam
| | - Van Ha Nguyen
- Faculty of Chemistry, VNU University of Science, Vietnam National University, Hanoi, 19 Le Thanh Tong, Hoan Kiem, Hanoi, 11021, Vietnam
| | - Tuan Thanh Dang
- Faculty of Chemistry, VNU University of Science, Vietnam National University, Hanoi, 19 Le Thanh Tong, Hoan Kiem, Hanoi, 11021, Vietnam
| | - Minh Tho Nguyen
- Laboratory for Chemical Computation and Modeling, Institute for Computational Science and Artificial Intelligence, Van Lang University, Ho Chi Minh City, 70000, Vietnam.
- Faculty of Applied Technology, School of Technology, Van Lang University, Ho Chi Minh City, 70000, Vietnam.
| |
Collapse
|
5
|
Heidari M, Gaichies T, Leibler L, Labousse M. Polymer time crystal: Mechanical activation of reversible bonds by low-amplitude high frequency excitations. SCIENCE ADVANCES 2024; 10:eadn6107. [PMID: 38781335 PMCID: PMC11114238 DOI: 10.1126/sciadv.adn6107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Accepted: 04/17/2024] [Indexed: 05/25/2024]
Abstract
Reversible supramolecular bonds play an important role in materials science and in biological systems. The equilibrium between open and closed bonds and the association rate can be controlled thermally, chemically, by mechanical pulling, by ultrasound, or by catalysts. In practice, these intrinsic equilibrium methods either suffer from a limited range of tunability or may damage the material. Here, we present a nonequilibrium strategy that exploits the dissipative properties of the system to control and change the dynamic properties of sacrificial and reversible networks. We show theoretically and numerically how high-frequency mechanical oscillations of very low amplitude can open or close bonds. This mechanism indicates how reversible bonds could alleviate mechanical fatigue of materials especially at low temperatures where they are fragile. In another area, it suggests that the system can be actively modified by the application of ultrasound to induce gel-fluid transitions and to activate or deactivate adhesion properties.
Collapse
Affiliation(s)
- Maziar Heidari
- Gulliver, CNRS, ESPCI Paris, Université PSL, 75005 Paris, France
| | - Théophile Gaichies
- Gulliver, CNRS, ESPCI Paris, Université PSL, 75005 Paris, France
- Département de chimie, École normale supérieure, Université PSL, 75005 Paris, France
| | - Ludwik Leibler
- Gulliver, CNRS, ESPCI Paris, Université PSL, 75005 Paris, France
| | | |
Collapse
|
6
|
Liu Y, Ghosh TK, Lin G, Chen M. Unbiasing Enhanced Sampling on a High-Dimensional Free Energy Surface with a Deep Generative Model. J Phys Chem Lett 2024; 15:3938-3945. [PMID: 38568182 DOI: 10.1021/acs.jpclett.3c03515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2024]
Abstract
Biased enhanced sampling methods that utilize collective variables (CVs) are powerful tools for sampling conformational ensembles. Due to their large intrinsic dimensions, efficiently generating conformational ensembles for complex systems requires enhanced sampling on high-dimensional free energy surfaces. While temperature-accelerated molecular dynamics (TAMD) can trivially adopt many CVs in a simulation, unbiasing the simulation to generate unbiased conformational ensembles requires accurate modeling of a high-dimensional CV probability distribution, which is challenging for traditional density estimation techniques. Here we propose an unbiasing method based on the score-based diffusion model, a deep generative learning method that excels in density estimation across complex data landscapes. We demonstrate that this unbiasing approach, tested on multiple TAMD simulations, significantly outperforms traditional unbiasing methods and can generate accurate unbiased conformational ensembles. With the proposed approach, TAMD can adopt CVs that focus on improving sampling efficiency and the proposed unbiasing method enables accurate evaluation of ensemble averages of important chemical features.
Collapse
Affiliation(s)
- Yikai Liu
- Department of Mechanical Engineering, Purdue University, West Lafayette, Indiana 47906, United States
| | - Tushar K Ghosh
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47906, United States
| | - Guang Lin
- Department of Mechanical Engineering, Purdue University, West Lafayette, Indiana 47906, United States
| | - Ming Chen
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47906, United States
| |
Collapse
|
7
|
Fei Z, Ma YH. Temperature fluctuations in mesoscopic systems. Phys Rev E 2024; 109:044101. [PMID: 38755872 DOI: 10.1103/physreve.109.044101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Accepted: 02/26/2024] [Indexed: 05/18/2024]
Abstract
Temperature is a fundamental concept in thermodynamics. In macroscopic thermodynamics, systems possess their own intrinsic temperature which equals the reservoir temperature when they equilibrate. In stochastic thermodynamics for simple systems at the microscopic level, thermodynamic quantities other than temperature (a deterministic parameter of the reservoir) are stochastic. To bridge the disparity in the perspectives about temperature between the micro- and macroregimes, we assign a generic mesoscopic N-body system an intrinsic fluctuating temperature T in this work. We simplify the complicated dynamics of numerous particles to one stochastic differential equation with respect to T, where the noise term accounts for finite-size effects arising from random energy transfer between the system and the reservoir. Our analysis reveals that these fluctuations make the extensive quantities (in the thermodynamic limit) deviate from being extensive. Moreover, we derive finite-size corrections, characterized by heat capacity of the system, to the Jarzynski equality. A possible violation of the principle of maximum work that scales with N^{-1} is also discussed. Additionally, we examine the impact of temperature fluctuations in a finite-size Carnot engine. We show that irreversible entropy production resulting from the temperature fluctuations of the working substance diminishes the average efficiency of the cycle as η_{C}-〈η〉∼N^{-1}, highlighting the unattainability of the Carnot efficiency η_{C} for mesoscopic heat engines even under the quasistatic limit. Our general framework paves the way for further exploration of nonequilibrium thermodynamics and the corresponding finite-size effects in a mesoscopic regime.
Collapse
Affiliation(s)
- Zhaoyu Fei
- Department of Physics and Key Laboratory of Optical Field Manipulation of Zhejiang Province, Zhejiang Sci-Tech University, Hangzhou 310018, China
- Graduate School of China Academy of Engineering Physics, No. 10 Xibeiwang East Road, Haidian District, Beijing 100193, China
| | - Yu-Han Ma
- Graduate School of China Academy of Engineering Physics, No. 10 Xibeiwang East Road, Haidian District, Beijing 100193, China
- Department of Physics, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
8
|
Post M, Wolf S, Stock G. Investigation of Rare Protein Conformational Transitions via Dissipation-Corrected Targeted Molecular Dynamics. J Chem Theory Comput 2023; 19:8978-8986. [PMID: 38011829 DOI: 10.1021/acs.jctc.3c01017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
To sample rare events, dissipation-corrected targeted molecular dynamics (dcTMD) applies a constant velocity constraint along a one-dimensional reaction coordinate s, which drives an atomistic system from an initial state into a target state. Employing a cumulant approximation of Jarzynski's identity, the free energy ΔG(s) is calculated from the mean external work and dissipated work of the process. By calculating the friction coefficient Γ(s) from the dissipated work, in a second step, the equilibrium dynamics of the process can be studied by propagating a Langevin equation. While so far dcTMD has been mostly applied to study the unbinding of protein-ligand complexes, here its applicability to rare conformational transitions within a protein and the prediction of their kinetics are investigated. As this typically requires the introduction of multiple collective variables {xj} = x, a theoretical framework is outlined to calculate the associated free energy ΔG(x) and friction Γ(x) from dcTMD simulations along coordinate s. Adopting the α-β transition of alanine dipeptide as well as the open-closed transition of T4 lysozyme as representative examples, the virtues and shortcomings of dcTMD to predict protein conformational transitions and the related kinetics are studied.
Collapse
Affiliation(s)
- Matthias Post
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, Freiburg 79104, Germany
| | - Steffen Wolf
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, Freiburg 79104, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, Freiburg 79104, Germany
| |
Collapse
|
9
|
Cheng CH, Lai PY. Theory of nonequilibrium asymptotic state thermodynamics: Interacting Ehrenfest urn ring as an example. Phys Rev E 2023; 108:064114. [PMID: 38243422 DOI: 10.1103/physreve.108.064114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2023] [Accepted: 11/10/2023] [Indexed: 01/21/2024]
Abstract
A generalized class of the nonequilibrium state, called the nonequilibrium asymptotic state (NEAS), is proposed. The NEAS is constructed within the framework of Fokker-Planck equations in thermodynamic limit. Besides the usual equilibrium state and nonequilibrium steady state (NESS), the class of NEAS could also cover the nonequilibrium periodic state (NEPS) in which its dynamics shows periodicity, the nonequilibrium quasiperiodic state (NEQPS), and nonequilibrium chaotic state (NECS) in which its dynamics becomes chaotic. Based on the theory of NEAS thermodynamics, the corresponding thermodynamics of different NEASs could also be determined. Finally the interacting Ehrenfest urn ring model is used as an example to illustrate how different kinds of NEAS (equilibrium state, uniform NESS, nonuniform NESS, NEPS) in the three-urn case are identified in our framework. In particular, the thermodynamics of NEPS and its phase transitions to other types of NEAS are studied.
Collapse
Affiliation(s)
- Chi-Ho Cheng
- Department of Physics, National Changhua University of Education, Taiwan, Republic of China
| | - Pik-Yin Lai
- Department of Physics and Center for Complex Systems, National Central University, Taiwan, Republic of China
- Physics Division, National Center for Theoretical Sciences, Taipei 10617, Taiwan, Republic of China
| |
Collapse
|
10
|
Liu W, Niu Z, Cheng W, Li X, Duan CK, Yin Z, Rong X, Du J. Experimental Test of the Jarzynski Equality in a Single Spin-1 System Using High-Fidelity Single-Shot Readouts. PHYSICAL REVIEW LETTERS 2023; 131:220401. [PMID: 38101345 DOI: 10.1103/physrevlett.131.220401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2022] [Revised: 10/22/2023] [Accepted: 10/24/2023] [Indexed: 12/17/2023]
Abstract
The Jarzynski equality (JE), which connects the equilibrium free energy with nonequilibrium work statistics, plays a crucial role in quantum thermodynamics. Although practical quantum systems are usually multilevel systems, most tests of the JE were executed in two-level systems. A rigorous test of the JE by directly measuring the work distribution of a physical process in a high-dimensional quantum system remains elusive. Here, we report an experimental test of the JE in a single spin-1 system. We realized nondemolition projective measurement of this three-level system via cascading high-fidelity single-shot readouts and directly measured the work distribution utilizing the two-point measurement protocol. The validity of the JE was verified from the nonadiabatic to adiabatic zone and under different effective temperatures. Our work puts the JE on a solid experimental foundation and makes the nitrogen-vacancy (NV) center system a mature toolbox to perform advanced experiments of stochastic quantum thermodynamics.
Collapse
Affiliation(s)
- Wenquan Liu
- CAS Key Laboratory of Microscale Magnetic Resonance and School of Physical Sciences, University of Science and Technology of China, Hefei 230026, China
- CAS Center for Excellence in Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, China
- School of Science, Beijing University of Posts and Telecommunications, Beijing, 100876, China
| | - Zhibo Niu
- CAS Key Laboratory of Microscale Magnetic Resonance and School of Physical Sciences, University of Science and Technology of China, Hefei 230026, China
- CAS Center for Excellence in Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, China
| | - Wei Cheng
- CAS Key Laboratory of Microscale Magnetic Resonance and School of Physical Sciences, University of Science and Technology of China, Hefei 230026, China
- CAS Center for Excellence in Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, China
| | - Xin Li
- Center for Quantum Technology Research and Key Laboratory of Advanced Optoelectronic Quantum Architecture and Measurements (MOE), School of Physics, Beijing Institute of Technology, Beijing 100081, China
| | - Chang-Kui Duan
- CAS Key Laboratory of Microscale Magnetic Resonance and School of Physical Sciences, University of Science and Technology of China, Hefei 230026, China
- CAS Center for Excellence in Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, China
| | - Zhangqi Yin
- Center for Quantum Technology Research and Key Laboratory of Advanced Optoelectronic Quantum Architecture and Measurements (MOE), School of Physics, Beijing Institute of Technology, Beijing 100081, China
| | - Xing Rong
- CAS Key Laboratory of Microscale Magnetic Resonance and School of Physical Sciences, University of Science and Technology of China, Hefei 230026, China
- CAS Center for Excellence in Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei 230088, China
| | - Jiangfeng Du
- CAS Key Laboratory of Microscale Magnetic Resonance and School of Physical Sciences, University of Science and Technology of China, Hefei 230026, China
- CAS Center for Excellence in Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei 230088, China
| |
Collapse
|
11
|
Moreira SV, Samuelsson P, Potts PP. Stochastic Thermodynamics of a Quantum Dot Coupled to a Finite-Size Reservoir. PHYSICAL REVIEW LETTERS 2023; 131:220405. [PMID: 38101369 DOI: 10.1103/physrevlett.131.220405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 11/07/2023] [Indexed: 12/17/2023]
Abstract
In nanoscale systems coupled to finite-size reservoirs, the reservoir temperature may fluctuate due to heat exchange between the system and the reservoirs. To date, a stochastic thermodynamic analysis of heat, work, and entropy production in such systems is, however, missing. Here we fill this gap by analyzing a single-level quantum dot tunnel coupled to a finite-size electronic reservoir. The system dynamics is described by a Markovian master equation, depending on the fluctuating temperature of the reservoir. Based on a fluctuation theorem, we identify the appropriate entropy production that results in a thermodynamically consistent statistical description. We illustrate our results by analyzing the work production for a finite-size reservoir Szilard engine.
Collapse
Affiliation(s)
- Saulo V Moreira
- Department of Physics and NanoLund, Lund University, Box 118, 22100 Lund, Sweden
- School of Physics, Trinity College Dublin, Dublin 2, Ireland
| | - Peter Samuelsson
- Department of Physics and NanoLund, Lund University, Box 118, 22100 Lund, Sweden
| | - Patrick P Potts
- Department of Physics and Swiss Nanoscience Institute, University of Basel, Klingelbergstrasse 82, 4056 Basel, Switzerland
| |
Collapse
|
12
|
Botti V, Lavecchia di Tocco F, Cannistraro S, Bizzarri AR. Hybridization Kinetics of miR-155 on Gold Surfaces as Investigated by Surface Plasmon Resonance and Atomic Force Spectroscopy. ACS OMEGA 2023; 8:38941-38949. [PMID: 37901511 PMCID: PMC10601050 DOI: 10.1021/acsomega.3c03318] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/12/2023] [Accepted: 08/15/2023] [Indexed: 10/31/2023]
Abstract
miRNAs are short noncoding RNA single strands, with a crucial role in several biological processes. miRNAs are dysregulated in several human diseases, and their detection is an important goal for diagnosis and screening. Innovative biosensors for miRNAs are commonly based on the hybridization process between a miRNA and its corresponding complementary strand (or suitable aptamers) immobilized onto an electrode surface forming a duplex. A detailed description of the hybridization kinetics in working conditions deserves a great deal of interest for the optimization of the biosensing process. Surface plasmon resonance (SPR) and atomic force spectroscopy (AFS) were applied to investigate the hybridization process between miR-155, a multifunctional miRNA that constitutes an important marker overexpressed in several diseases, and its complementary strand (antimiR-155), immobilized on the gold-coated surface of a commercial electrode. Under well-adjusted pH, ionic strength, surface coverage, and concentration, we found that miR-155 has a high affinity for antimiR-155 with kinetics well described by the 1:1 Langmuir model. Both techniques provided an association rate of about 104 M-1 s-1, while a dissociation rate of 10-5 and 10-4 s-1 was assessed by SPR and AFS, respectively. These results allowed us to establish optimized measurement running times for applications in biosensing. An analysis of AFS data also led us to evaluate the binding free energy for the duplex, which was found to be close to that of free molecules in solution. These results could guide in the implementation of fine-tuned working conditions of a biosensor for detecting miRNAs based on correspondent complementary strands.
Collapse
Affiliation(s)
- Valentina Botti
- Biophysics and Nanoscience Centre,
DEB, Università della Tuscia, Largo dell’Università, 01100 Viterbo, Italy
| | | | - Salvatore Cannistraro
- Biophysics and Nanoscience Centre,
DEB, Università della Tuscia, Largo dell’Università, 01100 Viterbo, Italy
| | - Anna Rita Bizzarri
- Biophysics and Nanoscience Centre,
DEB, Università della Tuscia, Largo dell’Università, 01100 Viterbo, Italy
| |
Collapse
|
13
|
Gavryushov S, Bashilov A, Cherashev-Tumanov KV, Kuzmich NN, Burykina TI, Izotov BN. Interaction of Synthetic Cannabinoid Receptor Agonists with Cannabinoid Receptor I: Insights into Activation Molecular Mechanism. Int J Mol Sci 2023; 24:14874. [PMID: 37834323 PMCID: PMC10574015 DOI: 10.3390/ijms241914874] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2023] [Revised: 09/22/2023] [Accepted: 09/27/2023] [Indexed: 10/15/2023] Open
Abstract
Synthetic cannabinoid receptor agonists (SCRAs) have become a wide group of new psychoactive substances since the 2010s. For the last few years, the X-ray structures of the complexes of cannabinoid receptor I (CB1) with SCRAs as well as the complexes of CB1 with its antagonist have been published. Based on those data, SCRA-CB1 interactions are analyzed in detail, using molecular modeling and molecular dynamics simulations. The molecular mechanism of the conformational transformation of the transmembrane domain of CB1 caused by its interaction with SCRA is studied. These conformational changes allosterically modulate the CB1-Gi complex, providing activation of the Gi protein. Based on the X-ray-determined structures of the CB1-ligand complexes, a stable apo conformation of inactive CB1 with a relatively low potential barrier of receptor activation was modeled. For that model, molecular dynamic simulations of SCRA binding to CB1 led to the active state of CB1, which allowed us to explore the key features of this activation and the molecular mechanism of the receptor's structural transformation. The simulated CB1 activation is in accordance with the previously published experimental data for the activation at protein mutations or structural changes of ligands. The key feature of the suggested activation mechanism is the determination of the stiff core of the CB1 transmembrane domain and the statement that the entire conformational transformation of the receptor to the active state is caused by a shift of alpha helix TM7 relative to this core. The shift itself is caused by protein-ligand interactions. It was verified via steered molecular dynamics simulations of the X-ray-determined structures of the inactive receptor, which resulted in the active conformation of CB1 irrespective of the placement of agonist ligand in the receptor's active site.
Collapse
Affiliation(s)
- Sergei Gavryushov
- Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, Vavilova Str. 32, Moscow 119991, Russia
- Institute for Translational Medicine and Biotechnology, Sechenov First Moscow State Medical University, 2-4 Bolshaya Pirogovskaya Str., Moscow 119991, Russia; (A.B.); (K.V.C.-T.); (T.I.B.); (B.N.I.)
| | - Anton Bashilov
- Institute for Translational Medicine and Biotechnology, Sechenov First Moscow State Medical University, 2-4 Bolshaya Pirogovskaya Str., Moscow 119991, Russia; (A.B.); (K.V.C.-T.); (T.I.B.); (B.N.I.)
- Skolkovo Institute of Science and Technology, Bolshoy Boulevard 30, Bld. 1, Moscow 121205, Russia
| | - Konstantin V. Cherashev-Tumanov
- Institute for Translational Medicine and Biotechnology, Sechenov First Moscow State Medical University, 2-4 Bolshaya Pirogovskaya Str., Moscow 119991, Russia; (A.B.); (K.V.C.-T.); (T.I.B.); (B.N.I.)
| | - Nikolay N. Kuzmich
- The Maurice and Vivienne Wohl Institute for Drug Discovery, Weizmann Institute of Science, Rehovot 7610001, Israel;
| | - Tatyana I. Burykina
- Institute for Translational Medicine and Biotechnology, Sechenov First Moscow State Medical University, 2-4 Bolshaya Pirogovskaya Str., Moscow 119991, Russia; (A.B.); (K.V.C.-T.); (T.I.B.); (B.N.I.)
| | - Boris N. Izotov
- Institute for Translational Medicine and Biotechnology, Sechenov First Moscow State Medical University, 2-4 Bolshaya Pirogovskaya Str., Moscow 119991, Russia; (A.B.); (K.V.C.-T.); (T.I.B.); (B.N.I.)
| |
Collapse
|
14
|
Fuller JT, Barnes S, Sadun LA, Ajmera P, Alexandrova AN, Sadun AA. Coenzyme Q10 trapping in mitochondrial complex I underlies Leber's hereditary optic neuropathy. Proc Natl Acad Sci U S A 2023; 120:e2304884120. [PMID: 37733737 PMCID: PMC10523484 DOI: 10.1073/pnas.2304884120] [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: 03/27/2023] [Accepted: 06/21/2023] [Indexed: 09/23/2023] Open
Abstract
How does a single amino acid mutation occurring in the blinding disease, Leber's hereditary optic neuropathy (LHON), impair electron shuttling in mitochondria? We investigated changes induced by the m.3460 G>A mutation in mitochondrial protein ND1 using the tools of Molecular Dynamics and Free Energy Perturbation simulations, with the goal of determining the mechanism by which this mutation affects mitochondrial function. A recent analysis suggested that the mutation's replacement of alanine A52 with a threonine perturbs the stability of a region where binding of the electron shuttling protein, Coenzyme Q10, occurs. We found two functionally opposing changes involving the role of Coenzyme Q10. The first showed that quantum electron transfer from the terminal Fe/S complex, N2, to the Coenzyme Q10 headgroup, docked in its binding pocket, is enhanced. However, this positive adjustment is overshadowed by our finding that the mobility of Coenzyme Q10 in its oxidized and reduced states, entering and exiting its binding pocket, is disrupted by the mutation in a manner that leads to conditions promoting the generation of reactive oxygen species. An increase in reactive oxygen species caused by the LHON mutation has been proposed to be responsible for this optic neuropathy.
Collapse
Affiliation(s)
- Jack T. Fuller
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA90095
| | - Steven Barnes
- Department of Ophthalmology, David Geffen School of Medicine, University of California, Los Angeles, CA90095
- Doheny Eye Institute, Pasadena, CA91103
- Department of Neurobiology, David Geffen School of Medicine, University of California, Los Angeles, CA90095
| | - Lorenzo A. Sadun
- Department of Mathematics, University of Texas at Austin, Austin, TX78712
| | - Pujan Ajmera
- Department of Chemistry and Biochemistry, University of California, Los Angeles, CA90095
| | | | - Alfredo A. Sadun
- Department of Ophthalmology, David Geffen School of Medicine, University of California, Los Angeles, CA90095
- Doheny Eye Institute, Pasadena, CA91103
| |
Collapse
|
15
|
Nguyen HL, Thai NQ, Li MS. Identifying inhibitors of NSP16-NSP10 of SARS-CoV-2 from large databases. J Biomol Struct Dyn 2023; 41:7045-7054. [PMID: 36002258 DOI: 10.1080/07391102.2022.2114941] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Accepted: 08/14/2022] [Indexed: 12/15/2022]
Abstract
The COVID-19 pandemic, which has already claimed millions of lives, continues to pose a serious threat to human health, requiring the development of new effective drugs. Non-structural proteins of SARS-CoV-2 play an important role in viral replication and infection. Among them, NSP16 (non-structured protein 16) and its cofactor NSP10 (non-structured protein 10) perform C2'-O methylation at the 5' end of the viral RNA, which promotes efficient virus replication. Therefore, the NSP16-NSP10 complex becomes an attractive target for drug development. Using a multi-step virtual screening protocol which includes Lipinski's rule, docking, steered molecular dynamics and umbrella sampling, we searched for potential inhibitors from the PubChem and anti-HIV databases. It has been shown that CID 135566620 compound from PubChem is the best candidate with an inhibition constant in the sub-μM range. The Van der Waals interaction was found to be more important than the electrostatic interaction in the binding affinity of this compound to NSP16-NSP10. Further in vitro and in vivo studies are needed to test the activity of the identified compound against COVID-19.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Hoang Linh Nguyen
- Life Science Lab, Institute for Computational Science and Technology, Quang Trung, Software City, Ho Chi Minh City, Vietnam
- Ho Chi Minh City University of Technology (HCMUT), Ho Chi Minh City, Vietnam
- Vietnam National University, Ho Chi Minh City, Vietnam
| | | | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland
| |
Collapse
|
16
|
Zhang ZZ, Tan QS, Wu W. Heat distribution in quantum Brownian motion. Phys Rev E 2023; 108:014138. [PMID: 37583192 DOI: 10.1103/physreve.108.014138] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2023] [Accepted: 07/11/2023] [Indexed: 08/17/2023]
Abstract
We investigate the heat statistics in a relaxation process of quantum Brownian motion described by the Caldeira-Leggett model. By employing the normal mode transformation and the phase-space formulation approach, we can analyze the quantum heat distribution within an exactly dynamical framework beyond the traditional paradigm of Born-Markovian and weak-coupling approximations. It is revealed that the exchange fluctuation theorem for quantum heat generally breaks down in the strongly non-Markovian regime. Our results may improve the understanding about the nonequilibrium thermodynamics of open quantum systems when the usual Markovian treatment is no longer appropriate.
Collapse
Affiliation(s)
- Ze-Zhou Zhang
- Key Laboratory of Quantum Theory and Applications of Ministry of Education, Lanzhou University, Lanzhou 730000, China
- Lanzhou Center for Theoretical Physics and Key Laboratory of Theoretical Physics of Gansu Province, Lanzhou University, Lanzhou 730000, China
| | - Qing-Shou Tan
- Key Laboratory of Hunan Province on Information Photonics and Freespace Optical Communication, College of Physics and Electronics, Hunan Institute of Science and Technology, Yueyang 414000, China
| | - Wei Wu
- Key Laboratory of Quantum Theory and Applications of Ministry of Education, Lanzhou University, Lanzhou 730000, China
- Lanzhou Center for Theoretical Physics and Key Laboratory of Theoretical Physics of Gansu Province, Lanzhou University, Lanzhou 730000, China
| |
Collapse
|
17
|
Wu YX, Chen JF, Pei JH, Zhang F, Quan HT. Heat statistics in the relaxation process of the Edwards-Wilkinson elastic manifold. Phys Rev E 2023; 107:064115. [PMID: 37464632 DOI: 10.1103/physreve.107.064115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2023] [Accepted: 05/23/2023] [Indexed: 07/20/2023]
Abstract
The stochastic thermodynamics of systems with a few degrees of freedom has been studied extensively so far. We would like to extend the study to systems with more degrees of freedom and even further-continuous fields with infinite degrees of freedom. The simplest case for a continuous stochastic field is the Edwards-Wilkinson elastic manifold. It is an exactly solvable model of which the heat statistics in the relaxation process can be calculated analytically. The cumulants require a cutoff spacing to avoid ultraviolet divergence. The scaling behavior of the heat cumulants with time and the system size as well as the large deviation rate function of the heat statistics in the large size limit is obtained.
Collapse
Affiliation(s)
- Yu-Xin Wu
- School of Physics, Peking University, Beijing 100871, China
| | - Jin-Fu Chen
- School of Physics, Peking University, Beijing 100871, China
| | - Ji-Hui Pei
- School of Physics, Peking University, Beijing 100871, China
| | - Fan Zhang
- School of Physics, Peking University, Beijing 100871, China
| | - H T Quan
- School of Physics, Peking University, Beijing 100871, China
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
- Frontiers Science Center for Nano-optoelectronics, Peking University, Beijing 100871, China
| |
Collapse
|
18
|
Kuznets-Speck B, Limmer DT. Inferring equilibrium transition rates from nonequilibrium protocols. Biophys J 2023; 122:1659-1664. [PMID: 36964656 PMCID: PMC10183322 DOI: 10.1016/j.bpj.2023.03.031] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2022] [Revised: 01/08/2023] [Accepted: 03/21/2023] [Indexed: 03/26/2023] Open
Abstract
We develop a theory for inferring equilibrium transition rates from trajectories driven by a time-dependent force using results from stochastic thermodynamics. Applying the Kawasaki relation to approximate the nonequilibrium distribution function in terms of the equilibrium distribution function and the excess dissipation, we formulate a nonequilibrium transition state theory to estimate the rate enhancement over the equilibrium rate due to the nonequilibrium protocol. We demonstrate the utility of our theory in examples of pulling of harmonically trapped particles in one and two dimensions, as well as a semiflexible polymer with a reactive linker in three dimensions. We expect our purely thermodynamic approach will find use in both molecular simulation and force spectroscopy experiments.
Collapse
Affiliation(s)
| | - David T Limmer
- Chemistry Department, University of California, Berkeley, Berkeley, California; Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California; Material Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California; Kavli Energy NanoSciences Institute, University of California, Berkeley, Berkeley, California.
| |
Collapse
|
19
|
Goswami Y, Sastry S. Kinetic reconstruction of free energies as a function of multiple order parameters. J Chem Phys 2023; 158:144502. [PMID: 37061464 DOI: 10.1063/5.0144338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
A vast array of phenomena, ranging from chemical reactions to phase transformations, are analyzed in terms of a free energy surface defined with respect to a single or multiple order parameters. Enhanced sampling methods are typically used, especially in the presence of large free energy barriers, to estimate free energies using biasing protocols and sampling of transition paths. Kinetic reconstructions of free energy barriers of intermediate height have been performed, with respect to a single order parameter, employing the steady state properties of unconstrained simulation trajectories when barrier crossing is achievable with reasonable computational effort. Considering such cases, we describe a method to estimate free energy surfaces with respect to multiple order parameters from a steady state ensemble of trajectories. The approach applies to cases where the transition rates between pairs of order parameter values considered is not affected by the presence of an absorbing boundary, whereas the macroscopic fluxes and sampling probabilities are. We demonstrate the applicability of our prescription on different test cases of random walkers executing Brownian motion in order parameter space with an underlying (free) energy landscape and discuss strategies to improve numerical estimates of the fluxes and sampling. We next use this approach to reconstruct the free energy surface for supercooled liquid silicon with respect to the degree of crystallinity and density, from unconstrained molecular dynamics simulations, and obtain results quantitatively consistent with earlier results from umbrella sampling.
Collapse
Affiliation(s)
- Yagyik Goswami
- Theoretical Sciences Unit and School of Advanced Materials, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur Campus, Bengaluru 560064, India
| | - Srikanth Sastry
- Theoretical Sciences Unit and School of Advanced Materials, Jawaharlal Nehru Centre for Advanced Scientific Research, Jakkur Campus, Bengaluru 560064, India
| |
Collapse
|
20
|
Chen JF, Quan HT. Hierarchical structure of fluctuation theorems for a driven system in contact with multiple heat reservoirs. Phys Rev E 2023; 107:024135. [PMID: 36932622 DOI: 10.1103/physreve.107.024135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Accepted: 02/09/2023] [Indexed: 06/18/2023]
Abstract
For driven open systems in contact with multiple heat reservoirs, we find the marginal distributions of work or heat do not satisfy any fluctuation theorem, but only the joint distribution of work and heat satisfies a family of fluctuation theorems. A hierarchical structure of these fluctuation theorems is discovered from microreversibility of the dynamics by adopting a step-by-step coarse-graining procedure in both classical and quantum regimes. Thus, we put all fluctuation theorems concerning work and heat into a unified framework. We also propose a general method to calculate the joint statistics of work and heat in the situation of multiple heat reservoirs via the Feynman-Kac equation. For a classical Brownian particle in contact with multiple heat reservoirs, we verify the validity of the fluctuation theorems for the joint distribution of work and heat.
Collapse
Affiliation(s)
- Jin-Fu Chen
- School of Physics, Peking University, Beijing 100871, China
| | - H T Quan
- School of Physics, Peking University, Beijing 100871, China
- Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
- Frontiers Science Center for Nano-optoelectronics, Peking University, Beijing 100871, China
| |
Collapse
|
21
|
Chandra H, Arya MK, Verma A. Biodiversity of butterflies (Lepidoptera: Rhopalocera) in the protected landscape of Nandhour, Uttarakhand, India. JOURNAL OF THREATENED TAXA 2023. [DOI: 10.11609/jott.7519.15.1.22448-22470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
An appraisal of butterfly species composition in terms of comparative diversity richness in various habitat types was conducted in and around the Nandhour Wildlife Sanctuary of Terai Arc landscape in Uttarakhand. During the two years of survey period (March 2018–February 2020), a total of 89 species of butterflies belonging to families Nymphalidae 43 species, Pieridae 15 species, Lycaenidae 13 species, Hesperiidae eight species, Papilionidae seven species, and Riodinidae three species were recorded, of which 46 species represented new records for the Nandhour Landscape. Butterfly diversity and richness were highest in dense moist & open dry riverine forests and lowest in human settlements & agricultural land. No significant differences in the number of species were found in moist mixed deciduous forest, subtropical Chir Pine forest, moist Bhabar Sal forest, moist Shiwalik Sal forest, and plantation forest. Eight species are endemic to the Indian Himalayan Region.
Collapse
|
22
|
Guéry-Odelin D, Jarzynski C, Plata CA, Prados A, Trizac E. Driving rapidly while remaining in control: classical shortcuts from Hamiltonian to stochastic dynamics. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2023; 86:035902. [PMID: 36535018 DOI: 10.1088/1361-6633/acacad] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/12/2022] [Accepted: 12/19/2022] [Indexed: 06/17/2023]
Abstract
Stochastic thermodynamics lays down a broad framework to revisit the venerable concepts of heat, work and entropy production for individual stochastic trajectories of mesoscopic systems. Remarkably, this approach, relying on stochastic equations of motion, introduces time into the description of thermodynamic processes-which opens the way to fine control them. As a result, the field of finite-time thermodynamics of mesoscopic systems has blossomed. In this article, after introducing a few concepts of control for isolated mechanical systems evolving according to deterministic equations of motion, we review the different strategies that have been developed to realize finite-time state-to-state transformations in both over and underdamped regimes, by the proper design of time-dependent control parameters/driving. The systems under study are stochastic, epitomized by a Brownian object immersed in a fluid; they are thus strongly coupled to their environment playing the role of a reservoir. Interestingly, a few of those methods (inverse engineering, counterdiabatic driving, fast-forward) are directly inspired by their counterpart in quantum control. The review also analyzes the control through reservoir engineering. Besides the reachability of a given target state from a known initial state, the question of the optimal path is discussed. Optimality is here defined with respect to a cost function, a subject intimately related to the field of information thermodynamics and the question of speed limit. Another natural extension discussed deals with the connection between arbitrary states or non-equilibrium steady states. This field of control in stochastic thermodynamics enjoys a wealth of applications, ranging from optimal mesoscopic heat engines to population control in biological systems.
Collapse
Affiliation(s)
- David Guéry-Odelin
- Laboratoire Collisions, Agrégats, Réactivité, IRSAMC, Université de Toulouse, CNRS, Toulouse, France
| | - Christopher Jarzynski
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD, United States of America
- Institute for Physical Science and Technology, University of Maryland, College Park, MD, United States of America
- Department of Physics, University of Maryland, College Park, MD, United States of America
| | - Carlos A Plata
- Física Teórica, Universidad de Sevilla, Apartado de Correos 1065, E-41080 Sevilla, Spain
| | - Antonio Prados
- Física Teórica, Universidad de Sevilla, Apartado de Correos 1065, E-41080 Sevilla, Spain
| | | |
Collapse
|
23
|
Hepatoprotective Effect of Millettia dielsiana: In Vitro and In Silico Study. MOLECULES (BASEL, SWITZERLAND) 2022; 27:molecules27248978. [PMID: 36558112 PMCID: PMC9787921 DOI: 10.3390/molecules27248978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Revised: 12/07/2022] [Accepted: 12/08/2022] [Indexed: 12/24/2022]
Abstract
In silico docking studies of 50 selected compounds from Millettia dielsiana Harms ex Diels (family Leguminosae) were docked into the binding pocket of the PI3K/mTOR protein. In there, compounds trans-3-O-p-hydroxycinnamoyl ursolic acid (1) and 5,7,4'-trihydroxyisoflavone 7-O-β-D-apiofuranosyl-(1→6)-β-D-glucopyranoside (2) are predicted to be very promising inhibitors against PI3K/mTOR. They direct their cytotoxic activity against Hepatocellular carcinoma with binding affinity (BA) values, the pulling work spent to the co-crystallized ligand from the binding site of PI3K/mTOR (W and Fmax), and the non-equilibrium binding free energy (∆GneqJar) as BA values = -9.237 and -9.083 kcal/mol, W = 83.5 ± 10.6 kcal/mol with Fmax = 336.2 ± 45.3 pN and 126.6 ± 21.7 kcal/mol with Fmax = 430.3 ± 84.0 pN, and ∆GneqJar = -69.86074 and -101.2317 kcal/mol, respectively. In molecular dynamic simulation, the RMSD value of the PI3K/mTOR complex with compounds (1 and 2) was in the range of 0.3 nm to the end of the simulation. Therefore, the compounds (1 and 2) are predicted to be very promising inhibitors against PI3K/mTOR. The crude extract, ethyl acetate fraction and compounds (1 and 2) from Millettia dielsiana exhibited moderate to potent in vitro cytotoxicity on Hepatocellular carcinoma cell line with IC50 values of 81.2 µg/mL, 60.4 µg/mL, 23.1 μM, and 16.3 μM, respectively, and showed relatively potent to potent in vitro antioxidant activity on mouse hepatocytes with ED50 values of 24.4 µg/mL, 19.3 µg/mL, 30.7 μM, and 20.5 μM, respectively. In conclusion, Millettia dielsiana and compounds (1 and 2) are predicted to have very promising cytotoxic activity against Hepatocellular carcinoma and have a hepatoprotective effect.
Collapse
|
24
|
Louwerse MD, Sivak DA. Connections between efficient control and spontaneous transitions in an Ising model. Phys Rev E 2022; 106:064124. [PMID: 36671088 DOI: 10.1103/physreve.106.064124] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Accepted: 12/02/2022] [Indexed: 06/17/2023]
Abstract
A system can be driven between metastable configurations by a time-dependent driving protocol, which uses external control parameters to change the potential energy of the system. Here we investigate the correspondence between driving protocols that are designed to minimize work and the spontaneous transition paths of the system in the absence of driving. We study the spin-inversion reaction in a 2D Ising model, quantifying the timing of each spin flip and heat flow to the system during both a minimum-work protocol and a spontaneous transition. The general order of spin flips during the transition mechanism is preserved between the processes, despite the coarseness of control parameters that are unable to reproduce more detailed features of the spontaneous mechanism. Additionally, external control parameters provide energy to each system component to compensate changes in internal energy, showing how control parameters are tuned during a minimum-work protocol to counteract underlying energetic features. This paper supports a correspondence between minimum-work protocols and spontaneous transition mechanisms.
Collapse
Affiliation(s)
- Miranda D Louwerse
- Department of Chemistry, Simon Fraser University, Burnaby, British Columbia, Canada V5A1S6
| | - David A Sivak
- Department of Physics, Simon Fraser University, Burnaby, British Columbia, Canada V5A1S6
| |
Collapse
|
25
|
Zhong A, DeWeese MR. Limited-control optimal protocols arbitrarily far from equilibrium. Phys Rev E 2022; 106:044135. [PMID: 36397571 DOI: 10.1103/physreve.106.044135] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Accepted: 09/16/2022] [Indexed: 06/16/2023]
Abstract
Recent studies have explored finite-time dissipation-minimizing protocols for stochastic thermodynamic systems driven arbitrarily far from equilibrium, when granted full external control to drive the system. However, in both simulation and experimental contexts, systems often may only be controlled with a limited set of degrees of freedom. Here, going beyond slow- and fast-driving approximations employed in previous studies, we obtain exact finite-time optimal protocols for this limited-control setting. By working with deterministic Fokker-Planck probability density time evolution, we can frame the work-minimizing protocol problem in the standard form of an optimal control theory problem. We demonstrate that finding the exact optimal protocol is equivalent to solving a system of Hamiltonian partial differential equations, which in many cases admit efficiently calculable numerical solutions. Within this framework, we reproduce analytical results for the optimal control of harmonic potentials and numerically devise optimal protocols for two anharmonic examples: varying the stiffness of a quartic potential and linearly biasing a double-well potential. We confirm that these optimal protocols outperform other protocols produced through previous methods, in some cases by a substantial amount. We find that for the linearly biased double-well problem, the mean position under the optimal protocol travels at a near-constant velocity. Surprisingly, for a certain timescale and barrier height regime, the optimal protocol is also nonmonotonic in time.
Collapse
Affiliation(s)
- Adrianne Zhong
- Department of Physics, University of California, Berkeley, Berkeley, California 94720, USA
| | - Michael R DeWeese
- Department of Physics, University of California, Berkeley, Berkeley, California 94720, USA and Redwood Center For Theoretical Neuroscience and Helen Wills Neuroscience Institute, University of California, Berkeley, Berkeley, California 94720, USA
| |
Collapse
|
26
|
Semaan MT, Crutchfield JP. Homeostatic and adaptive energetics: Nonequilibrium fluctuations beyond detailed balance in voltage-gated ion channels. Phys Rev E 2022; 106:044410. [PMID: 36397574 DOI: 10.1103/physreve.106.044410] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2022] [Accepted: 09/15/2022] [Indexed: 06/16/2023]
Abstract
Stochastic thermodynamics has largely succeeded in characterizing both equilibrium and far-from-equilibrium phenomena. Yet many opportunities remain for application to mesoscopic complex systems-especially biological ones-whose effective dynamics often violate detailed balance and whose microscopic degrees of freedom are often unknown or intractable. After reviewing excess and housekeeping energetics-the adaptive and homeostatic components of a system's dissipation-we extend stochastic thermodynamics with a trajectory class fluctuation theorem for nonequilibrium steady-state, nondetailed-balanced complex systems. We then take up the neurobiological examples of voltage-gated sodium and potassium ion channels to apply and illustrate the theory, elucidating their nonequilibrium behavior under a biophysically plausible action potential drive. These results uncover challenges for future experiments and highlight the progress possible understanding the thermodynamics of complex systems-without exhaustive knowledge of every underlying degree of freedom.
Collapse
Affiliation(s)
- Mikhael T Semaan
- Complexity Sciences Center and Department of Physics and Astronomy, University of California, Davis, One Shields Avenue, Davis, California 95616, USA
| | - James P Crutchfield
- Complexity Sciences Center and Department of Physics and Astronomy, University of California, Davis, One Shields Avenue, Davis, California 95616, USA
| |
Collapse
|
27
|
Bassman Oftelie L, Klymko K, Liu D, Tubman NM, de Jong WA. Computing Free Energies with Fluctuation Relations on Quantum Computers. PHYSICAL REVIEW LETTERS 2022; 129:130603. [PMID: 36206437 DOI: 10.1103/physrevlett.129.130603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Revised: 01/07/2022] [Accepted: 08/05/2022] [Indexed: 06/16/2023]
Abstract
As a central thermodynamic property, free energy enables the calculation of virtually any equilibrium property of a physical system, allowing for the construction of phase diagrams and predictions about transport, chemical reactions, and biological processes. Thus, methods for efficiently computing free energies, which in general is a difficult problem, are of great interest to broad areas of physics and the natural sciences. The majority of techniques for computing free energies target classical systems, leaving the computation of free energies in quantum systems less explored. Recently developed fluctuation relations enable the computation of free energy differences in quantum systems from an ensemble of dynamic simulations. While performing such simulations is exponentially hard on classical computers, quantum computers can efficiently simulate the dynamics of quantum systems. Here, we present an algorithm utilizing a fluctuation relation known as the Jarzynski equality to approximate free energy differences of quantum systems on a quantum computer. We discuss under which conditions our approximation becomes exact, and under which conditions it serves as a strict upper bound. Furthermore, we successfully demonstrate a proof of concept of our algorithm using the transverse field Ising model on a real quantum processor. As quantum hardware continues to improve, we anticipate that our algorithm will enable computation of free energy differences for a wide range of quantum systems useful across the natural sciences.
Collapse
Affiliation(s)
| | | | - Diyi Liu
- School of Mathematics, University of Minnesota, Minnesota 55455, USA
| | - Norm M Tubman
- NASA Ames Research Center, Mountain View, California 94035, USA
| | - Wibe A de Jong
- Lawrence Berkeley National Lab, Berkeley, California 94720
| |
Collapse
|
28
|
Materzok T, De Boer D, Gorb S, Müller-Plathe F. Gecko Adhesion on Flat and Rough Surfaces: Simulations with a Multi-Scale Molecular Model. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2022; 18:e2201674. [PMID: 35927024 DOI: 10.1002/smll.202201674] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/16/2022] [Revised: 07/07/2022] [Indexed: 06/15/2023]
Abstract
A multiscale modeling approach is used to develop a particle-based mesoscale gecko spatula model that is able to link atomistic simulations and mesoscale (0.44 µm) simulations. It is used to study the detachment of spatulae from flat as well as nanostructured surfaces. The spatula model is based on microscopical information about spatulae structure and on atomistic molecular simulation results. Target properties for the coarse-graining result from a united-atom model of gecko keratin in periodic boundary conditions (PBC), previously developed by the authors. Pull-off forces necessary to detach gecko keratin under 2D PBC parallel to the surface are previously overestimated when only a small region of a spatula is examined. It is shown here that this is due to the restricted geometry (i.e., missing peel-off mode) and not model parameters. The spatula model peels off when pulled away from a surface, both in the molecular picture of the pull-off process and in the force-extension curve of non-equilibrium simulations mimicking single-spatula detachment studied with atomic force microscopy equipment. The force field and spatula model can reproduce experimental pull-off forces. Inspired by experimental results, the underlying mechanism that causes pull-off forces to be at a minimum on surfaces of varying roughnesses is also investigated. A clear sigmoidal increase in the pull-off force of spatulae with surface roughness shows that adhesion is determined by the ratio between spatula pad area and the area between surface peaks. Experiments showed a correlation with root-mean-square roughness of the surface, but the results of this work indicate that this is not a causality but depends on the area accessible.
Collapse
Affiliation(s)
- Tobias Materzok
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie and Profile Area Thermofluids and Interfaces, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, D-64287, Darmstadt, Germany
| | - Danna De Boer
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie and Profile Area Thermofluids and Interfaces, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, D-64287, Darmstadt, Germany
| | - Stanislav Gorb
- Zoological Institute Functional Morphology and Biomechanics, Kiel University, Am Botanischen Garten 1-9, D-24118, Kiel, Germany
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie and Profile Area Thermofluids and Interfaces, Technische Universität Darmstadt, Alarich-Weiss-Str. 8, D-64287, Darmstadt, Germany
| |
Collapse
|
29
|
The Development of Single Molecule Force Spectroscopy: From Polymer Biophysics to Molecular Machines. Q Rev Biophys 2022; 55:e9. [PMID: 35916314 DOI: 10.1017/s0033583522000087] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
30
|
Thu TTM, Li MS. Protein aggregation rate depends on mechanical stability of fibrillar structure . J Chem Phys 2022; 157:055101. [DOI: 10.1063/5.0088689] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The formation of the fibrillar structure of amyloid proteins/peptides is believed to be associated with neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, etc. Since the rate of aggregation can influence neurotoxicity, finding the key factors that control this rate is of paramount importance. It was recently found evidence that the rate of protein aggregation is related to the mechanical stability of the fibrillar structure, such that the higher the mechanical stability, the faster the fibril is formed. However, this conclusion was supported by a limited dataset. In this work, we expand the previous study to a larger dataset, including the wild type of Aβ42 peptide and its 20 mutants, the aggregation rate of which was measured experimentally. By using all-atom steered molecular dynamics (SMD) simulations we can access the mechanical stability of the fibril structure, which is characterized by the rupture force, pulling work and unbinding free energy barrier. Our result confirms that mechanical stability is indeed related to the aggregation rate. Since estimation of the aggregation rate using all-atom simulations is almost forbidden by the current computational capabilities, our result is useful for predicting it based on information obtained from fast SMD simulations for fibrils.
Collapse
Affiliation(s)
| | - Mai Suan Li
- Theoretical Physics, Institute of Physics, Polish Academy of Sciences, Poland
| |
Collapse
|
31
|
Naeini VF, Baniassadi M, Foroutan M, Rémond Y, George D. Decisive structural elements in water and ion permeation through mechanosensitive channels of large conductance: insights from molecular dynamics simulation. RSC Adv 2022; 12:17803-17816. [PMID: 35765322 PMCID: PMC9201702 DOI: 10.1039/d2ra02284b] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2022] [Accepted: 06/09/2022] [Indexed: 11/21/2022] Open
Abstract
In this paper, a series of equilibrium molecular dynamics simulations (EMD), steered molecular dynamics (SMD), and computational electrophysiology methods are carried out to explore water and ion permeation through mechanosensitive channels of large conductance (MscL). This research aims to identify the pore-lining side chains of the channel in different conformations of MscL homologs by analyzing the pore size. The distribution of permeating water dipole angles through the pore domains enclosed by VAL21 and GLU104 demonstrated that water molecules are oriented toward the charged oxygen headgroups of GLU104 from their hydrogen atoms to retain this interaction in a stabilized fashion. Although, this behavior was not perceived for VAL21. Numerical assessments of the secondary structure clarified that, during the ion permeation, in addition to the secondary structure alterations, the structure of Tb-MscL would also undergo significant conformational changes. It was elucidated that VAL21, GLU104, and water molecules accomplish a fundamental task in ion permeation. The mentioned residues hinder ion permeation so that the pulling SMD force is increased remarkably when the ions permeate through the domains enclosed by VAL21 and GLU102. The hydration level and potassium diffusivity in the hydrophobic gate of the transmembrane domain were promoted by applying the external electric field. Furthermore, the implementation of an external electric field altered the distribution pattern for potassium ions in the system while intensifying the accumulation of Cl− in the vicinity of ARG11 and ARG98. Graphical representation of the most determinant pore-lining side chains of Tb-MscL along with the solid surfaces depicting the spatial shape of the interior pore.![]()
Collapse
Affiliation(s)
- Vahid Fadaei Naeini
- Division of Machine Elements, Luleå University of Technology 97187 Luleå Sweden
| | - Majid Baniassadi
- School of Mechanical Engineering, College of Engineering, University of Tehran Tehran Iran.,University of Strasbourg, ICube Laboratory/CNRS 2 Rue Boussingault 67000 Strasbourg France
| | - Masumeh Foroutan
- Department of Physical Chemistry, School of Chemistry, College of Science, University of Tehran Tehran Iran
| | - Yves Rémond
- University of Strasbourg, ICube Laboratory/CNRS 2 Rue Boussingault 67000 Strasbourg France
| | - Daniel George
- University of Strasbourg, ICube Laboratory/CNRS 2 Rue Boussingault 67000 Strasbourg France
| |
Collapse
|
32
|
Verification of Information Thermodynamics in a Trapped Ion System. ENTROPY 2022; 24:e24060813. [PMID: 35741534 PMCID: PMC9222944 DOI: 10.3390/e24060813] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/22/2022] [Revised: 06/07/2022] [Accepted: 06/08/2022] [Indexed: 02/05/2023]
Abstract
Information thermodynamics has developed rapidly over past years, and the trapped ions, as a controllable quantum system, have demonstrated feasibility to experimentally verify the theoretical predictions in the information thermodynamics. Here, we address some representative theories of information thermodynamics, such as the quantum Landauer principle, information equality based on the two-point measurement, information-theoretical bound of irreversibility, and speed limit restrained by the entropy production of system, and review their experimental demonstration in the trapped ion system. In these schemes, the typical physical processes, such as the entropy flow, energy transfer, and information flow, build the connection between thermodynamic processes and information variation. We then elucidate the concrete quantum control strategies to simulate these processes by using quantum operators and the decay paths in the trapped-ion system. Based on them, some significantly dynamical processes in the trapped ion system to realize the newly proposed information-thermodynamic models is reviewed. Although only some latest experimental results of information thermodynamics with a single trapped-ion quantum system are reviewed here, we expect to find more exploration in the future with more ions involved in the experimental systems.
Collapse
|
33
|
Nguyen H, Li MS. Antibody-nanobody combination increases their neutralizing activity against SARS-CoV-2 and nanobody H11-H4 is effective against Alpha, Kappa and Delta variants. Sci Rep 2022; 12:9701. [PMID: 35690632 PMCID: PMC9188278 DOI: 10.1038/s41598-022-14263-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 06/03/2022] [Indexed: 11/14/2022] Open
Abstract
The global spread of COVID-19 is devastating health systems and economies worldwide. While the use of vaccines has yielded encouraging results, the emergence of new variants of SARS-CoV-2 shows that combating COVID-19 remains a big challenge. One of the most promising treatments is the use of not only antibodies, but also nanobodies. Recent experimental studies revealed that the combination of antibody and nanobody can significantly improve their neutralizing ability through binding to the SARS-CoV-2 spike protein, but the molecular mechanisms underlying this observation remain largely unknown. In this work, we investigated the binding affinity of the CR3022 antibody and H11-H4 nanobody to the SARS-CoV-2 receptor binding domain (RBD) using molecular modeling. Both all-atom steered molecular dynamics simulations and coarse-grained umbrella sampling showed that, consistent with the experiment, CR3022 associates with RBD more strongly than H11-H4. We predict that the combination of CR3022 and H11-H4 considerably increases their binding affinity to the spike protein. The electrostatic interaction was found to control the association strength of CR3022, but the van der Waals interaction dominates in the case of H11-H4. However, our study for a larger set of nanobodies and antibodies showed that the relative role of these interactions depends on the specific complex. Importantly, we showed Beta, Gamma, Lambda, and Mu variants reduce the H11-H4 activity while Alpha, Kappa and Delta variants increase its neutralizing ability, which is in line with experiment reporting that the nanobody elicited from the llama is very promising for fighting against the Delta variant.
Collapse
Affiliation(s)
- Hung Nguyen
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668, Warsaw, Poland
| | - Mai Suan Li
- Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668, Warsaw, Poland. .,Life Science Lab, Institute for Computational Science and Technology, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, Ho Chi Minh City, Vietnam.
| |
Collapse
|
34
|
Sumikama T, Federici Canova F, Gao DZ, Penedo M, Miyazawa K, Foster AS, Fukuma T. Computed Three-Dimensional Atomic Force Microscopy Images of Biopolymers Using the Jarzynski Equality. J Phys Chem Lett 2022; 13:5365-5371. [PMID: 35678499 PMCID: PMC9208010 DOI: 10.1021/acs.jpclett.2c01093] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2022] [Accepted: 05/18/2022] [Indexed: 06/15/2023]
Abstract
Three-dimensional atomic force microscopy (3D-AFM) has resolved three-dimensional distributions of solvent molecules at solid-liquid interfaces at the subnanometer scale. This method is now being extended to the imaging of biopolymer assemblies such as chromosomes or proteins in cells, with the expectation of being able to resolve their three-dimensional structures. Here, we have developed a computational method to simulate 3D-AFM images of biopolymers by using the Jarzynski equality. It is found that some parts of the fiber structure of biopolymers are indeed resolved in the 3D-AFM image. The dependency of 3D-AFM images on the vertical scanning velocity is investigated, and optimum scanning velocities are found. It is also clarified that forces in nonequilibrium processes are measured in 3D-AFM measurements when the dynamics of polymers are slower than the scanning of the probe.
Collapse
Affiliation(s)
- Takashi Sumikama
- PRESTO,
JST, Kawaguchi, Saitama 332-0012, Japan
- Nano
Life Science Institute (WPI-NanoLSI), Kanazawa
University, Kanazawa 920-1192, Japan
| | - Filippo Federici Canova
- Nanolayers
Research Computing Ltd., 1 Granville Court, Granville Road, London N12 0HL, United Kingdom
- Department
of Applied Physics, Aalto University, Aalto 00076, Finland
| | - David Z. Gao
- Nanolayers
Research Computing Ltd., 1 Granville Court, Granville Road, London N12 0HL, United Kingdom
- Department
of Physics, Norwegian University of Science
and Technology (NTNU), 7491 Trondheim, Norway
| | - Marcos Penedo
- Nano
Life Science Institute (WPI-NanoLSI), Kanazawa
University, Kanazawa 920-1192, Japan
- Laboratory
for Bio and Nanoinstrumentation, Institute for Bioengineering, École Polytechnique Fédérale
de Lausanne, Lausanne CH-1015, Switzerland
| | - Keisuke Miyazawa
- Nano
Life Science Institute (WPI-NanoLSI), Kanazawa
University, Kanazawa 920-1192, Japan
- Division
of Electrical Engineering and Computer Science, Kanazawa University, Kanazawa 920-1192, Japan
- Faculty of
Frontier Engineering, Kanazawa University, Kanazawa 920-1192, Japan
| | - Adam S. Foster
- Nano
Life Science Institute (WPI-NanoLSI), Kanazawa
University, Kanazawa 920-1192, Japan
- Department
of Applied Physics, Aalto University, Aalto 00076, Finland
| | - Takeshi Fukuma
- Nano
Life Science Institute (WPI-NanoLSI), Kanazawa
University, Kanazawa 920-1192, Japan
- Division
of Electrical Engineering and Computer Science, Kanazawa University, Kanazawa 920-1192, Japan
| |
Collapse
|
35
|
Abstract
Fluctuation theorems specify the non-zero probability to observe negative entropy production, contrary to a naive expectation from the second law of thermodynamics. For closed particle trajectories in a fluid, Stokes theorem can be used to give a geometric characterization of the entropy production. Building on this picture, we formulate a topological fluctuation theorem that depends only by the winding number around each vortex core and is insensitive to other aspects of the force. The probability is robust to local deformations of the particle trajectory, reminiscent of topologically protected modes in various classical and quantum systems. We demonstrate that entropy production is quantized in these strongly fluctuating systems, and it is controlled by a topological invariant. We demonstrate that the theorem holds even when the probability distributions are non-Gaussian functions of the generated heat. While topology is crucial in complex systems, stochastic thermodynamics uncovers universal constraints for non-equilibrium fluctuations. The authors combine these two areas and formulate a fluctuation theorem for the heat dissipated along closed loops in vortex force fields, which is found to be topologically protected.
Collapse
|
36
|
Louwerse MD, Sivak D. Multidimensional minimum-work control of a 2D Ising model. J Chem Phys 2022; 156:194108. [DOI: 10.1063/5.0086079] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A system's configurational state can be manipulated using dynamic variation of control parameters, such as temperature, pressure, or magnetic field; for finite-duration driving, excess work is required above the equilibrium free-energy change. Minimum-work protocols in multidimensional control-parameter space have potential to significantly reduce work relative to one-dimensional control. By numerically minimizing a linear-response approximation to the excess work, we design protocols in control-parameter spaces of a 2D Ising model that efficiently drive the system from the all-down to all-up configuration. We find that such designed multidimensional protocols take advantage of more flexible control to avoid control-parameter regions of high system resistance, heterogeneously input and extract work to make use of system relaxation, and flatten the energy landscape, making accessible many configurations that would otherwise have prohibitively high energy and thus decreasing spin correlations. Relative to one-dimensional protocols, this speeds up the rate-limiting spin-inversion reaction, thereby keeping the system significantly closer to equilibrium for a wide range of protocol durations, and significantly reducing resistance and hence work.
Collapse
|
37
|
Nguyen H, Lan PD, Nissley DA, O’Brien EP, Li MS. Cocktail of REGN Antibodies Binds More Strongly to SARS-CoV-2 Than Its Components, but the Omicron Variant Reduces Its Neutralizing Ability. J Phys Chem B 2022; 126:2812-2823. [PMID: 35403431 PMCID: PMC9016775 DOI: 10.1021/acs.jpcb.2c00708] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
A promising approach to combat Covid-19 infections is the development of effective antiviral antibodies that target the SARS-CoV-2 spike protein. Understanding the structures and molecular mechanisms underlying the binding of antibodies to SARS-CoV-2 can contribute to quickly achieving this goal. Recently, a cocktail of REGN10987 and REGN10933 antibodies was shown to be an excellent candidate for the treatment of Covid-19. Here, using all-atom steered molecular dynamics and coarse-grained umbrella sampling, we examine the interactions of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein with REGN10987 and REGN10933 separately as well as together. Both computational methods show that REGN10933 binds to RBD more strongly than REGN10987. Importantly, the cocktail binds to RBD (simultaneous binding) more strongly than its components. The dissociation constants of REGN10987-RBD and REGN10933-RBD complexes calculated from the coarse-grained simulations are in good agreement with the experimental data. Thus, REGN10933 is probably a better candidate for treating Covid-19 than REGN10987, although the cocktail appears to neutralize the virus more efficiently than REGN10933 or REGN10987 alone. The association of REGN10987 with RBD is driven by van der Waals interactions, while electrostatic interactions dominate in the case of REGN10933 and the cocktail. We also studied the effectiveness of these antibodies on the two most dangerous variants Delta and Omicron. Consistent with recent experimental reports, our results confirmed that the Omicron variant reduces the neutralizing activity of REGN10933, REGN10987, and REGN10933+REGN10987 with the K417N, N440K, L484A, and Q498R mutations playing a decisive role, while the Delta variant slightly changes their activity.
Collapse
Affiliation(s)
- Hung Nguyen
- Institute
of Physics, Polish Academy of Sciences, al. Lotnikow 32/46, 02-668 Warsaw, Poland
| | - Pham Dang Lan
- Life
Science Lab, Institute for Computational
Science and Technology, Quang Trung Software City, Tan Chanh Hiep Ward, District 12, 729110 Ho Chi Minh City, Vietnam,Faculty
of Physics and Engineering Physics, VNUHCM-University
of Science, 227, Nguyen
Van Cu Street, District 5, 749000 Ho Chi Minh City, Vietnam
| | - Daniel A. Nissley
- Department
of Statistics, University of Oxford, Oxford
Protein Bioinformatics Group, OX1 2JD Oxford, United Kingdom
| | - Edward P. O’Brien
- Department
of Chemistry, Penn State University, University Park, Pennsylvania 16802, United States,Bioinformatics
and Genomics Graduate Program, The Huck Institutes of the Life Sciences, Penn State University, University Park, Pennsylvania 16802, United States,Institute
for Computational and Data Sciences, Penn
State University, University Park, Pennsylvania 16802, United States
| | - Mai Suan Li
- Institute
of Physics, Polish Academy of Sciences, al. Lotnikow 32/46, 02-668 Warsaw, Poland,
| |
Collapse
|
38
|
Sun Z, He Q. Seeding the multi-dimensional nonequilibrium pulling for Hamiltonian variation: indirect nonequilibrium free energy simulations at QM levels. Phys Chem Chem Phys 2022; 24:8800-8819. [PMID: 35352744 DOI: 10.1039/d2cp00355d] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
The combination of free energy simulations in the alchemical and configurational spaces provides a feasible route to access the thermodynamic profiles under a computationally demanding target Hamiltonian. Normally, due to the significant differences between the computational cost of ab initio quantum mechanics (QM) calculations and those of semi-empirical quantum mechanics (SQM) and molecular mechanics (MM), this indirect method could be used to obtain the QM thermodynamics by combining the SQM or MM results and the SQM-to-QM or MM-to-QM corrections. In our previous work, a multi-dimensional nonequilibrium pulling framework for Hamiltonian variations was introduced based on bidirectional pulling and bidirectional reweighting. The method performs nonequilibrium free energy simulations in the configurational space to obtain the thermodynamic profile along the conformational change pathway under a selected computationally efficient Hamiltonian, and uses the nonequilibrium alchemical method to correct or perturb the thermodynamic profile to that under the target Hamiltonian. The BAR-based method is designed to achieve the best generality and transferability and thus leads to modest (∼20 fold) speedup. In this work, we explore the possibility of further accelerating the nonequilibrium free energy simulation by employing unidirectional pulling and using the selection criterion to obtain the initial configurations used to initiate nonequilibrium trajectories following the idea of adaptive steered molecular dynamics (ASMD). A single initial condition is used to seed the whole multi-dimensional nonequilibrium free energy simulation and the sampling is performed fully in the nonequilibrium ensemble. Introducing very short ps-length equilibrium sampling to grab more initial seeds could also be helpful. The ASMD scheme estimates the free energy difference with the unidirectional exponential average (EXP), but it does not follow exactly the requirements of the EXP estimator. Another deficiency of the seeding simulation is the inherently sequential or serial pulling due to the inter-segment dependency, which triggers some problems in the parallelizability of the simulation. Numerical tests are performed to grasp some insights and guidelines for using this selection-criterion-based ASMD scheme. The presented selection-criterion-based multi-dimensional ASMD scheme follows the same perturbation network of the BAR-based method, and thus could be used in various Hamiltonian-variation cases.
Collapse
Affiliation(s)
- Zhaoxi Sun
- Beijing National Laboratory for Molecular Sciences, College of Chemistry and Molecular Engineering, Institute of Theoretical and Computational Chemistry, Peking University, Beijing 100871, China.
| | - Qiaole He
- AI Department of Enzymaster (Ningbo) Bio-Engineering Co., Ltd, North Century Avenue 333, 315100 Ningbo, China
| |
Collapse
|
39
|
Mana T, Kundu J, Lahiri H, Bera S, Kolay J, Sinha S, Mukhopadhyay R. Molecularly resolved, label-free nucleic acid sensing at solid-liquid interface using non-ionic DNA analogues. RSC Adv 2022; 12:9263-9274. [PMID: 35424880 PMCID: PMC8985177 DOI: 10.1039/d2ra00386d] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2022] [Accepted: 03/07/2022] [Indexed: 01/10/2023] Open
Abstract
Nucleic acid-based biosensors, where the capture probe is a nucleic acid, e.g., DNA or its synthetic analogue xeno nucleic acid (XNA), offer interesting ways of eliciting clinically relevant information from hybridization/dehybridization signals. In this respect, the application of XNA probes is attractive since the drawbacks of DNA probes might be overcome. Within the XNA probe repertoire, peptide nucleic acid (PNA) and morpholino (MO) are promising since their backbones are non-ionic. Therefore, in the absence of electrostatic charge repulsion between the capture probe and the target nucleic acid, a stable duplex can be formed. In addition, these are nuclease-resistant probes. Herein, we have tested the molecularly resolved nucleic acid sensing capacity of PNA and MO capture probes using a fluorescent label-free single molecule force spectroscopy approach. As far as single nucleobase mismatch discrimination is concerned, both PNA and MO performed better than DNA, while the performance of the MO probe was the best. We propose that the conformationally more rigid backbone of MO, compared to the conformationally flexible PNA, is an advantage for MO, since the probe orientation can be made more upright on the surface and therefore MO can be more effectively accessed by the target sequences. The performance of the XNA probes has been compared to that of the DNA probe, using fixed nucleobase sequences, so that the effect of backbone variation could be investigated. To our knowledge, this is the first report on molecularly resolved nucleic acid sensing by non-ionic capture probes, here, MO and PNA. Improved nucleic acid sensing in terms of single nucleobase mismatch discrimination, as achieved by the surface-confined non-ionic PNA and MO capture probes, is exemplified by single molecule force spectroscopy.![]()
Collapse
Affiliation(s)
- Tanushree Mana
- School of Biological Sciences, Indian Association for the Cultivation of Science Jadavpur Kolkata 700 032 India +91 33 2473 2805 +91 33 2473 4971 extn 1506
| | - Jayanta Kundu
- School of Applied and Interdisciplinary Sciences, Indian Association for the Cultivation of Science Jadavpur Kolkata 700 032 India
| | - Hiya Lahiri
- School of Biological Sciences, Indian Association for the Cultivation of Science Jadavpur Kolkata 700 032 India +91 33 2473 2805 +91 33 2473 4971 extn 1506
| | - Sudipta Bera
- School of Biological Sciences, Indian Association for the Cultivation of Science Jadavpur Kolkata 700 032 India +91 33 2473 2805 +91 33 2473 4971 extn 1506
| | - Jayeeta Kolay
- School of Biological Sciences, Indian Association for the Cultivation of Science Jadavpur Kolkata 700 032 India +91 33 2473 2805 +91 33 2473 4971 extn 1506
| | - Surajit Sinha
- School of Applied and Interdisciplinary Sciences, Indian Association for the Cultivation of Science Jadavpur Kolkata 700 032 India
| | - Rupa Mukhopadhyay
- School of Biological Sciences, Indian Association for the Cultivation of Science Jadavpur Kolkata 700 032 India +91 33 2473 2805 +91 33 2473 4971 extn 1506
| |
Collapse
|
40
|
Nandi T, Koti Ainavarapu SR. Reconstruction of the Free Energy Profile for SUMO1 from Nonequilibrium Single-Molecule Pulling Experiments. J Phys Chem B 2022; 126:2168-2172. [PMID: 35271281 DOI: 10.1021/acs.jpcb.1c08596] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Free energy profiles form the cornerstone in the study of protein folding and function. In this study, the free energy profile of SUMO1 protein is directly reconstructed using an extension of the Jarzynski equality from atomic force microscope (AFM) based single-molecule force spectroscopy (SMFS) experiments. SUMO1 is a ubiquitin-like posttranslational modifier protein having a β clamp motif in its structure, imparting it with mechanical stability. We use the Jarzynski equality to obtain the equilibrium free energy profile from repeated nonequilibrium single-molecule pulling experiments. Indeed, the free energy values determined by the Jarzynski equality are lesser than the normal work average at all extensions. The free energy profiles constructed for the two velocities (100 and 400 nm/s) overlap with each other. The unfolding free energy barrier is estimated to be ∼7.5 kcal/mol. We anticipate that the Jarzynski equality can be applied in a similar manner to other ubiquitin-like proteins to extract their differences in the free energy profile, and hence, the effect of sequence diversity of structurally homologous proteins on the free energy landscape can be studied.
Collapse
Affiliation(s)
- Tathagata Nandi
- Department of Chemical Sciences, Tata Institute of Fundamental Research Dr. Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Sri Rama Koti Ainavarapu
- Department of Chemical Sciences, Tata Institute of Fundamental Research Dr. Homi Bhabha Road, Colaba, Mumbai 400005, India
| |
Collapse
|
41
|
Chen H, Ogden D, Pant S, Cai W, Tajkhorshid E, Moradi M, Roux B, Chipot C. A Companion Guide to the String Method with Swarms of Trajectories: Characterization, Performance, and Pitfalls. J Chem Theory Comput 2022; 18:1406-1422. [PMID: 35138832 PMCID: PMC8904302 DOI: 10.1021/acs.jctc.1c01049] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The string method with swarms of trajectories (SMwST) is an algorithm that identifies a physically meaningful transition pathway─a one-dimensional curve, embedded within a high-dimensional space of selected collective variables. The SMwST algorithm leans on a series of short, unbiased molecular dynamics simulations spawned at different locations of the discretized path, from whence an average dynamic drift is determined to evolve the string toward an optimal pathway. However conceptually simple in both its theoretical formulation and practical implementation, the SMwST algorithm is computationally intensive and requires a careful choice of parameters for optimal cost-effectiveness in applications to challenging problems in chemistry and biology. In this contribution, the SMwST algorithm is presented in a self-contained manner, discussing with a critical eye its theoretical underpinnings, applicability, inherent limitations, and use in the context of path-following free-energy calculations and their possible extension to kinetics modeling. Through multiple simulations of a prototypical polypeptide, combining the search of the transition pathway and the computation of the potential of mean force along it, several practical aspects of the methodology are examined with the objective of optimizing the computational effort, yet without sacrificing accuracy. In light of the results reported here, we propose some general guidelines aimed at improving the efficiency and reliability of the computed pathways and free-energy profiles underlying the conformational transitions at hand.
Collapse
Affiliation(s)
- Haochuan Chen
- Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Nankai University, Tianjin 300071, China
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche no 7019, Université de Lorraine, B.P. 70239, 54506 Vandœuvre-lès-Nancy Cedex, France
| | - Dylan Ogden
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Shashank Pant
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Wensheng Cai
- Research Center for Analytical Sciences, College of Chemistry, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Nankai University, Tianjin 300071, China
| | - Emad Tajkhorshid
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Department of Biochemistry and Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - Mahmoud Moradi
- Department of Chemistry and Biochemistry, University of Arkansas, Fayetteville, Arkansas 72701, United States
| | - Benoît Roux
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, United States
| | - Christophe Chipot
- Theoretical and Computational Biophysics Group, NIH Center for Macromolecular Modeling and Bioinformatics, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
- Laboratoire International Associé Centre National de la Recherche Scientifique et University of Illinois at Urbana-Champaign, Unité Mixte de Recherche no 7019, Université de Lorraine, B.P. 70239, 54506 Vandœuvre-lès-Nancy Cedex, France
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
42
|
Chekmarev SF. Extraction of kinetics from equilibrium distributions of states using the Metropolis Monte Carlo method. Phys Rev E 2022; 105:034407. [PMID: 35428044 DOI: 10.1103/physreve.105.034407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Accepted: 02/25/2022] [Indexed: 06/14/2023]
Abstract
The Metropolis Monte Carlo (MMC) method is used to extract reaction kinetics from a given equilibrium distribution of states of a complex system. The approach is illustrated by the folding/unfolding reaction for two proteins: a model β-hairpin and α-helical protein α_{3}D. For the β-hairpin, the free energy surfaces (FESs) and free energy profiles (FEPs) are employed as the equilibrium distributions of states, playing a role of the potentials of mean force to determine the acceptance probabilities of new states in the MMC simulations. Based on the FESs and PESs for a set of temperatures that were simulated with the molecular dynamics (MD) method, the MMC simulations are performed to extract folding/unfolding rates. It has been found that the rate constants and first-passage time (FPT) distributions obtained in the MMC simulations change with temperature in good agreement with those from the MD simulations. For α_{3}D protein, whose equilibrium folding/unfolding was studied with the single-molecule FRET method [Chung et al., J. Phys. Chem. A 115, 3642 (2011)1089-563910.1021/jp1009669], the FRET-efficiency histograms at different denaturant concentrations were used as the equilibrium distributions of protein states. It has been found that the rate constants for folding and unfolding obtained in the MMC simulations change with denaturant concentration in reasonable agreement with the constants that were extracted from the photon trajectories on the basis of theoretical models. The simulated FPT distributions are single-exponential, which is consistent with the assumption of two-state kinetics that was made in the theoretical models. The promising feature of the present approach is that it is based solely on the equilibrium distributions of states, without introducing any additional parameters to perform simulations, which suggests its applicability to other complex systems.
Collapse
Affiliation(s)
- Sergei F Chekmarev
- Institute of Thermophysics, SB RAS, 630090 Novosibirsk, Russia and Department of Physics, Novosibirsk State University, 630090 Novosibirsk, Russia
| |
Collapse
|
43
|
Rizzuti B. Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena. BIOCHIMICA ET BIOPHYSICA ACTA. PROTEINS AND PROTEOMICS 2022; 1870:140757. [PMID: 35051666 DOI: 10.1016/j.bbapap.2022.140757] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2021] [Revised: 01/08/2022] [Accepted: 01/10/2022] [Indexed: 12/22/2022]
Abstract
Molecular dynamics simulation is the most popular computational technique for investigating the structural and dynamical behaviour of proteins, in search of the molecular basis of their function. Far from being a completely settled field of research, simulations are still evolving to best capture the essential features of the atomic interactions that govern a protein's inner motions. Modern force fields are becoming increasingly accurate in providing a physical description adequate to this purpose, and allow us to model complex biological systems under fairly realistic conditions. Furthermore, the use of accelerated sampling techniques is improving our access to the observation of progressively larger molecular structures, longer time scales, and more hidden functional events. In this review, the basic principles of molecular dynamics simulations and a number of key applications in the area of protein science are summarized, and some of the most important results are discussed. Examples include the study of the structure, dynamics and binding properties of 'difficult' targets, such as intrinsically disordered proteins and membrane receptors, and the investigation of challenging phenomena like hydration-driven processes and protein aggregation. The findings described provide an overall picture of the current state of this research field, and indicate new perspectives on the road ahead to the upcoming future of molecular simulations.
Collapse
Affiliation(s)
- Bruno Rizzuti
- CNR-NANOTEC, SS Rende (CS), Department of Physics, University of Calabria, 87036 Rende, Italy; Institute for Biocomputation and Physics of Complex Systems (BIFI), Joint Unit GBsC-CSIC-BIFI, University of Zaragoza, 50018 Zaragoza, Spain.
| |
Collapse
|
44
|
Gimenez-Dejoz J, Numata K. Molecular dynamics study of the internalization of cell-penetrating peptides containing unnatural amino acids across membranes. NANOSCALE ADVANCES 2022; 4:397-407. [PMID: 36132688 PMCID: PMC9419563 DOI: 10.1039/d1na00674f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Accepted: 11/10/2021] [Indexed: 06/16/2023]
Abstract
Peptide-based delivery systems that deliver target molecules into cells have been gaining traction. These systems need cell-penetrating peptides (CPPs), which have the remarkable ability to penetrate into biological membranes and help internalize different cargoes into cells through the cell membranes. The molecular internalization mechanism and structure-function relationships of CPPs are not clear, although the incorporation of nonproteinogenic amino acids such as α-aminoisobutyric acid (Aib) has been reported to increase their helicity, biostability and penetration efficiencies. Here, we used molecular dynamics to study two Aib-containing CPPs, poly(LysAibAla)3 (KAibA) and poly(LysAibGly)3 (KAibG), that previously showed high cell internalization efficiency. KAibA and KAibG displayed the lowest internalization energies among the studied CPPs, showing distinct internalization mechanisms depending on the lipid composition of the model membranes. The presence of Aib residues allows these CPPs to adopt amphipathic folding to efficiently penetrate through the membranes. Elucidating how Aib incorporation affects CPP-membrane binding and interactions is beneficial for the design of CPPs for efficient intracellular delivery.
Collapse
Affiliation(s)
- Joan Gimenez-Dejoz
- Biomacromolecules Research Team, RIKEN Center for Sustainable Resource Science Saitama Japan
| | - Keiji Numata
- Biomacromolecules Research Team, RIKEN Center for Sustainable Resource Science Saitama Japan
- Department of Material Chemistry, Graduate School of Engineering, Kyoto University Kyoto Japan
| |
Collapse
|
45
|
Weerakoon D, Carbajo RJ, De Maria L, Tyrchan C, Zhao H. Impact of PROTAC Linker Plasticity on the Solution Conformations and Dissociation of the Ternary Complex. J Chem Inf Model 2022; 62:340-349. [PMID: 35018781 DOI: 10.1021/acs.jcim.1c01036] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The conformational behavior of a small molecule free in solution is important to understand the free energy of binding to its target. This could be of special interest for proteolysis-targeting chimeras (PROTACs) due to their often flexible and lengthy linkers and the need to induce a ternary complex. Here, we report on the molecular dynamics (MD) simulations of two PROTACs, MZ1 and dBET6, revealing different linker conformational behaviors. The simulation of MZ1 in dimethyl sulfoxide (DMSO) agrees well with the nuclear magnetic resonance study, providing strong support for the relevance of our simulations. To further understand the role of linker plasticity in the formation of a ternary complex, the dissociation of the complex von Hippel-Lindau-MZ1-BRD4 is investigated in detail by steered simulations and is shown to follow a two-step pathway. Interestingly, both MZ1 and dBET6 display in water, a tendency toward an intramolecular lipophilic interaction between the two warheads. The hydrophobic contact of the two warheads would prevent them from binding to their respective proteins and might have an effect on the efficacy of induced cellular protein degradation. However, conformations featuring this hydrophobic contact of the two warheads are calculated to be marginally more favorable.
Collapse
Affiliation(s)
- Dhanushka Weerakoon
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg 43183, Sweden
| | - Rodrigo J Carbajo
- Chemistry, Oncology R&D, AstraZeneca, Cambridge CB4 0QA, United Kingdom
| | - Leonardo De Maria
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg 43183, Sweden
| | - Christian Tyrchan
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg 43183, Sweden
| | - Hongtao Zhao
- Medicinal Chemistry, Research and Early Development, Respiratory and Immunology (R&I), BioPharmaceuticals R&D, AstraZeneca, Gothenburg 43183, Sweden
| |
Collapse
|
46
|
Saito T. Preaveraging description of polymer nonequilibrium stretching. Phys Rev E 2022; 105:014502. [PMID: 35193303 DOI: 10.1103/physreve.105.014502] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2021] [Accepted: 11/26/2021] [Indexed: 11/07/2022]
Abstract
This article focuses on a preaveraging description of polymer nonequilibrium stretching, where a single polymer undergoes a transient process from equilibrium to nonequilibrium steady state by pulling one chain end. The preaveraging method combined with mode analysis reduces the original Langevin equation to a simplified form for both a stretched steady state and an equilibrium state, even in the presence of self-avoiding repulsive interactions spanning a long range. However, the transient stretching process exhibits evolution of a hierarchal regime structure, which means a qualitative temporal change in probabilistic distributions assumed in preaveraging. We investigate the preaveraging method for evolution of the regime structure with consideration of the nonequilibrium work relations and deviations from the fluctuation-dissipation relation.
Collapse
Affiliation(s)
- Takuya Saito
- Department of Physical Sciences, Aoyama Gakuin University, Chuo-ku, Sagamihara 252-5258, Japan
| |
Collapse
|
47
|
Kobayashi H, Rohrbach PB, Scheichl R, Wilding NB, Jack RL. Critical point for demixing of binary hard spheres. Phys Rev E 2021; 104:044603. [PMID: 34781560 DOI: 10.1103/physreve.104.044603] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Accepted: 09/22/2021] [Indexed: 11/07/2022]
Abstract
We use a two-level simulation method to analyze the critical point associated with demixing of binary hard-sphere mixtures. The method exploits an accurate coarse-grained model with two- and three-body effective interactions. Using this model within the two-level methodology allows computation of properties of the full (fine-grained) mixture. The critical point is located by computing the probability distribution for the number of large particles in the grand canonical ensemble and matching to the universal form for the 3D Ising universality class. The results have a strong and unexpected dependence on the size ratio between large and small particles, which is related to three-body effective interactions and the geometry of the underlying hard-sphere packings.
Collapse
Affiliation(s)
- Hideki Kobayashi
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Paul B Rohrbach
- DAMTP, University of Cambridge, Centre for Mathematical Sciences, Wilberforce Road, Cambridge CB3 0WA, United Kingdom
| | - Robert Scheichl
- Institute for Applied Mathematics, Heidelberg University, INF 205, 69120 Heidelberg, Germany.,Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, United Kingdom
| | - Nigel B Wilding
- H. H. Wills Physics Laboratory, University of Bristol, Royal Fort, Bristol BS8 1TL, United Kingdom
| | - Robert L Jack
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom.,DAMTP, University of Cambridge, Centre for Mathematical Sciences, Wilberforce Road, Cambridge CB3 0WA, United Kingdom
| |
Collapse
|
48
|
Barden DR, Vashisth H. Water Dynamics in a Peptide-appended Pillar[5]arene Artificial Channel in Lipid and Biomimetic Membranes. Front Chem 2021; 9:753635. [PMID: 34778209 PMCID: PMC8586425 DOI: 10.3389/fchem.2021.753635] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 10/04/2021] [Indexed: 11/21/2022] Open
Abstract
Peptide-appended Pillar[5]arene (PAP) is an artificial water channel that can be incorporated into lipid and polymeric membranes to achieve high permeability and enhanced selectivity for angstrom-scale separations [Shen et al. Nat. Commun.9:2294 (2018)]. In comparison to commonly studied rigid carbon nanotubes, PAP channels are conformationally flexible, yet these channels allow a high water permeability [Y. Liu and H. Vashisth Phys. Chem. Chem. Phys.21:22711 (2019)]. Using molecular dynamics (MD) simulations, we study water dynamics in PAP channels embedded in biological (lipid) and biomimetic (block-copolymer) membranes to probe the effect of the membrane environment on water transport characteristics of PAP channels. We have resolved the free energy surface and local minima for water diffusion within the channel in each type of membrane. We find that water follows single file transport with low free-energy barriers in regions surroundings the central ring of the PAP channel and the single file diffusivity of water correlates with the number of hydrogen bonding sites within the channel, as is known for other sub-nm pore-size synthetic and biological water channels [Horner et al. Sci. Adv.1:e1400083 (2015)].
Collapse
Affiliation(s)
- Daniel Ryan Barden
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, United States
| | - Harish Vashisth
- Department of Chemical Engineering, University of New Hampshire, Durham, NH, United States
| |
Collapse
|
49
|
Werner P, Hartmann AK. Similarity of extremely rare nonequilibrium processes to equilibrium processes. Phys Rev E 2021; 104:034407. [PMID: 34654125 DOI: 10.1103/physreve.104.034407] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 08/31/2021] [Indexed: 11/07/2022]
Abstract
For system coupled to heat baths, typical nonequilibrated processes, e.g., induced by varying an external parameter without waiting for equilibration in between, are very different from the corresponding equilibrium infinitely slow processes. Nevertheless, there are connections between equilibrium and nonequilibrated behaviors, e.g., the theorems of Jarzynski and Crooks, which relate the distribution P(W) of nonequilibrium work to the free energy differences ΔF. Here we study the naturally arising question, whether those relevant but rare trajectories, which exhibit these work values, show a higher degree of similarity to equilibrium. For convenience, we have chosen a simple model of RNA secondary structures (or single-stranded DNA), here modeling a medium-size hairpin structure, under influence of a varying external force. This allows us to measure the work W during the resulting fast unfolding and refolding processes within Monte Carlo simulations, i.e., in nonequilibrium. Also we sample numerically efficiently directly in exact equilibrium, for comparison. Using a sophisticated large-deviation algorithm, we are able to measure work distributions with high precision down to probabilities as small as 10^{-46}, enabling us to verify the Crooks and Jarzynski theorems. Furthermore, we analyze force-extension curves and the configurations of the secondary structures during unfolding and refolding for typical equilibrium processes and nonequilibrated processes. We find that the nonequilibrated processes where the work values are close to those which are most relevant for applying Crooks and Jarzynski theorems, respectively, but which occur with exponential small probabilities, are most and quite similar to the equilibrium processes.
Collapse
Affiliation(s)
- Peter Werner
- Institut für Physik, Universität Oldenburg, 26111 Oldenburg, Germany
| | | |
Collapse
|
50
|
Lbadaoui-Darvas M, Garberoglio G, Karadima KS, Cordeiro MNDS, Nenes A, Takahama S. Molecular simulations of interfacial systems: challenges, applications and future perspectives. MOLECULAR SIMULATION 2021. [DOI: 10.1080/08927022.2021.1980215] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Affiliation(s)
- Mária Lbadaoui-Darvas
- ENAC/IIE; Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne, Switzerland
| | - Giovanni Garberoglio
- European Centre for Theoretical Studies in Nuclear Physics and Related Areas (FBK-ECT*), Trento, Italy
- Trento Institute for Fundamental Physics and Applications (TIFPA-INFN), Trento, Italy
| | - Katerina S. Karadima
- Department of Chemical Engineering, University of Patras, Patras, Greece
- Institute of Chemical Engineering Sciences, Foundation for Research and Technology Hellas(FORTH-ICE/HT), Patras, Greece
| | | | - Athanasios Nenes
- ENAC/IIE; Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne, Switzerland
- Institute of Chemical Engineering Sciences, Foundation for Research and Technology Hellas(FORTH-ICE/HT), Patras, Greece
| | - Satoshi Takahama
- ENAC/IIE; Swiss Federal Institute of Technology Lausanne (EPFL), Lausanne, Switzerland
| |
Collapse
|