skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Learned mappings for targeted free energy perturbation between peptide conformations
Targeted free energy perturbation uses an invertible mapping to promote configuration space overlap and the convergence of free energy estimates. However, developing suitable mappings can be challenging. Wirnsberger et al. [J. Chem. Phys. 153, 144112 (2020)] demonstrated the use of machine learning to train deep neural networks that map between Boltzmann distributions for different thermodynamic states. Here, we adapt their approach to the free energy differences of a flexible bonded molecule, deca-alanine, with harmonic biases and different spring centers. When the neural network is trained until “early stopping”—when the loss value of the test set increases—we calculate accurate free energy differences between thermodynamic states with spring centers separated by 1 Å and sometimes 2 Å. For more distant thermodynamic states, the mapping does not produce structures representative of the target state, and the method does not reproduce reference calculations.  more » « less
Award ID(s):
2429324
PAR ID:
10523589
Author(s) / Creator(s):
; ;
Publisher / Repository:
AIP Publishing
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
159
Issue:
12
ISSN:
0021-9606
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The mitoribosome translates mitochondrial mRNAs and regulates energy conversion that is a signature of aerobic life forms. We present a 2.2 Å resolution structure of human mitoribosome together with validated mitoribosomal RNA (rRNA) modifications, including aminoacylated CP-tRNAVal. The structure shows how mitoribosomal proteins stabilise binding of mRNA and tRNA helping to align it in the decoding center, whereas the GDP-bound mS29 stabilizes intersubunit communication. Comparison between different states, with respect to tRNA position, allowed us to characterize a non-canonical L1 stalk, and molecular dynamics simulations revealed how it facilitates tRNA transitions in a way that does not require interactions with rRNA. We also report functionally important polyamines that are depleted when cells are subjected to an antibiotic treatment. The structural, biochemical, and computational data illuminate the principal functional components of the translation mechanism in mitochondria and provide a description of the structure and function of the human mitoribosome. 
    more » « less
  2. Abstract Phase diagrams offer substantial predictive power for materials synthesis by identifying the stability regions of target phases. However, thermodynamic phase diagrams do not offer explicit information regarding the kinetic competitiveness of undesired by-product phases. Here we propose a quantitative and computable thermodynamic metric to identify synthesis conditions under which the propensity to form kinetically competing by-products is minimized. We hypothesize that thermodynamic competition is minimized when the difference in free energy between a target phase and the minimal energy of all other competing phases is maximized. We validate this hypothesis for aqueous materials synthesis through two empirical approaches: first, by analysing 331 aqueous synthesis recipes text-mined from the literature; and second, by systematic experimental synthesis of LiIn(IO3)4and LiFePO4across a wide range of aqueous electrochemical conditions. Our results show that even for synthesis conditions that are within the stability region of a thermodynamic Pourbaix diagram, phase-pure synthesis occurs only when thermodynamic competition with undesired phases is minimized. 
    more » « less
  3. Phytochromes (PHYs) are photoreceptor proteins first discovered in plants, where they control a variety of photomorphogenesis events. PHYs as photochromic proteins can reversibly switch between two distinct states: a red light (Pr) and a far-red light (Pfr) absorbing form. The discovery of Bacteriophytochromes (BphPs) in nonphotosynthetic bacteria has opened new frontiers in our understanding of the mechanisms by which these natural photoswitches can control single cell development, although the role of BphPs in vivo remains largely unknown. BphPs are dimeric proteins that consist of a photosensory core module (PCM) and an enzymatic domain, often a histidine kinase. The PCM is composed of three domains (PAS, GAF, and PHY). It holds a covalently bound open-chain tetrapyrrole (biliverdin, BV) chromophore. Upon absorption of light, the double bond between BV rings C and D isomerizes and reversibly switches the protein between Pr and Pfr states. We report crystal structures of the wild-type and mutant (His275Thr) forms of the canonical BphP from the nonphotosynthetic myxobacterium Stigmatella aurantiaca (SaBphP2) in the Pr state. Structures were determined at 1.65 Å and 2.2 Å (respectively), the highest resolution of any PCM construct to date. We also report the room temperature wild-type structure of the same protein determined at 2.1 Å at the SPring-8 Angstrom Compact free electron LAser (SACLA), Japan. Our results not only highlight and confirm important amino acids near the chromophore that play a role in Pr-Pfr photoconversion but also describe the signal transduction into the PHY domain which moves across tens of angstroms after the light stimulus. 
    more » « less
  4. null (Ed.)
    Abstract In the first two decades of the XXI century, corroles have emerged as an important class of porphyrinoids for photonics and biomedical photonics. In comparison with porphyrins, corroles have lower molecular symmetry and higher electron density, which leads to uniquely complementary properties. In macrocycles of free-base corroles, for example, three protons are distributed among four pyrrole nitrogens. It results in distinct tautomers that have different thermodynamic energies. Herein, we focus on the excited-state dynamics of a corrole modified with l -phenylalanine. The tautomerization in the singlet-excited state occurs in the timescales of about 10–100 picoseconds and exhibits substantial kinetic isotope effects. It, however, does not discernably affect nanosecond deactivation of the photoexcited corrole and its basic photophysics. Nevertheless, this excited-state tautomerization dynamics can strongly affect photoinduced processes with comparable or shorter timescales, considering the 100-meV energy differences between the tautomers in the excited state. The effects on the kinetics of charge transfer and energy transfer, initiated prior to reaching the equilibrium thermalization of the excited-state tautomer population, can be indeed substantial. Such considerations are crucially important in the design of systems for artificial photosynthesis and other forms of energy conversion and charge transduction. 
    more » « less
  5. Translocation of proteins is correlated with structural fluctuations that access conformational states higher in free energy than the folded state. We use electric fields at the solid-state nanopore to control the relative free energy and occupancy of different protein conformational states at the single-molecule level. The change in occupancy of different protein conformations as a function of electric field gives rise to shifts in the measured distributions of ionic current blockades and residence times. We probe the statistics of the ionic current blockades and residence times for three mutants of the λ -repressor family in order to determine the number of accessible conformational states of each mutant and evaluate the ruggedness of their free energy landscapes. Translocation becomes faster at higher electric fields when additional flexible conformations are available for threading through the pore. At the same time, folding rates are not correlated with ease of translocation; a slow-folding mutant with a low-lying intermediate state translocates faster than a faster-folding two-state mutant. Such behavior allows us to distinguish among protein mutants by selecting for the degree of current blockade and residence time at the pore. Based on these findings, we present a simple free energy model that explains the complementary relationship between folding equilibrium constants and translocation rates. 
    more » « less