Energy Landscapes and Structural Ensembles of Glucagon-like Peptide-1 Monomers.

While GLP-1 and its analogues are important pharmaceutical agents in the treatment of type 2 diabetes and obesity, their susceptibility to aggregate into amyloid fibrils poses a significant safety issue. Many factors may contribute to the aggregation propensity, including pH. While it is known that the monomeric structure of GLP-1 has a strong impact on primary nucleation, probing its diverse structural ensemble is challenging. Here, we investigated the monomer structural ensembles at pH 3, 4, and 7.5 using state-of-the-art computational methods in combination with experimental data. We found significant stabilization of ß-strand structures and destabilization of helical structures at lower pH, correlating with observed aggregation lag times, which are lower under these conditions. We further identified helical defects at pH 4, which led to the fastest observed aggregation, in agreement with our far-UV circular dichroism data. The detailed atomistic structures that result from the computational studies help to rationalize the experimental results on the aggregation propensity of GLP-1. This work provides a new insight into the pH-dependence of monomeric structural ensembles of GLP-1 and connects them to experimental observations.

Going for Gold(-Standard): Attaining Coupled Cluster Accuracy in Oxide-Supported Nanoclusters.

The structure of oxide-supported metal nanoclusters plays an essential role in their sharply enhanced catalytic activity over that of bulk metals. Simulations provide the atomic-scale resolution needed to understand these systems. However, the sensitive mix of metal-metal and metal-support interactions, which govern their structure, puts stringent requirements on the method used, requiring calculations beyond standard density functional theory (DFT). The method of choice is coupled cluster theory [specifically CCSD(T)], but its computational cost has so far prevented its application to these systems. In this work, we showcase two approaches to make CCSD(T) accuracy readily achievable in oxide-supported nanoclusters. First, we leverage the SKZCAM protocol to provide the first benchmarks of oxide-supported nanoclusters, revealing that it is specifically metal-metal interactions that are challenging to capture with DFT. Second, we propose a CCSD(T) correction (ΔCC) to the metal-metal interaction errors in DFT, reaching accuracy comparable to that of the SKZCAM protocol at significantly lower cost. This approach forges a path toward studying larger systems at reliable accuracy, which we highlight by identifying a ground-state structure in agreement with experiments for Au20 on MgO, a challenging system where DFT models have yielded conflicting predictions.

Oscillatory Active State of a Pd Nanocatalyst Identified by In Situ Capture of the Instantaneous Structure-Activity Change at the Atomic Scale.

Identifying the active phase with the highest activity, which is long-believed to be a steady state of the catalyst, is the basis of rational design of heterogeneous catalysis. In this work, we performed detailed in situ investigations, successfully capturing the instantaneous structure-activity change in oscillating Pd nanocatalysts during methane oxidation, which reveals an unprecedented oscillatory active state. Combining in situ quantitative environmental transmission electron microscopy and highly sensitive online mass spectrometry, we identified two distinct phases for the reaction: one where the Pd nanoparticles refill with oxygen, and the other, a period of abrupt pumping of oxygen and boosted methane oxidation within about 1 s. It is the rapid reduction process that shows the highest activity for total oxidation of methane, not a PdO or Pd steady state under the conditions applied here (methane:oxygen = 5:1). This observation challenges the traditional understanding of the active phase and requires a completely different strategy for catalyst optimization.

100 Years of the Lennard-Jones Potential.

It is now 100 years since Lennard-Jones published his first paper introducing the now famous potential that bears his name. It is therefore timely to reflect on the many achievements, as well as the limitations, of this potential in the theory of atomic and molecular interactions, where applications range from descriptions of intermolecular forces to molecules, clusters, and condensed matter.

Multilevel Framework for Analysis of Protein Folding Involving Disulfide Bond Formation.

In this study, a three-layered multicenter ONIOM approach is implemented to characterize the naive folding pathway of bovine pancreatic trypsin inhibitor (BPTI). Each layer represents a distinct level of theory, where the initial layer, encompassing the entire protein, is modeled by a general all-atom force-field GFN-FF. An intermediate electronic structure layer consisting of three multicenter fragments is introduced with the state-of-the-art semiempirical tight-binding method GFN2-xTB. Higher accuracy, specifically addressing the breaking and formation of the three disulfide bonds, is achieved at the innermost layer using the composite DFT method r2SCAN-3c. Our analysis sheds light on the structural stability of BPTI, particularly the significance of interlinking disulfide bonds. The accuracy and efficiency of the multicenter QM/SQM/MM approach are benchmarked using the oxidative formation of cystine. For the folding pathway of BPTI, relative stabilities are investigated through the calculation of free energy contributions for selected intermediates, focusing on the impact of the disulfide bond. Our results highlight the intricate trade-off between accuracy and computational cost, demonstrating that the multicenter ONIOM approach provides a well-balanced and comprehensive solution to describe electronic structure effects in biomolecular systems. We conclude that multiscale energy landscape exploration provides a robust methodology for the study of intriguing biological targets.

Analysis and interpretation of first passage time distributions featuring rare events.

In this contribution we consider theory and associated computational tools to treat the kinetics associated with competing pathways on multifunnel energy landscapes. Multifunnel landscapes are associated with molecular switches and multifunctional materials, and are expected to exhibit multiple relaxation time scales and associated thermodynamic signatures in the heat capacity. Our focus here is on the first passage time distribution, which is encoded in a kinetic transition network containing all the locally stable states and the pathways between them. This network can be renormalised to reduce the dimensionality, while exactly conserving the mean first passage time and approximately conserving the full distribution. The structure of the reduced network can be visualised using disconnectivity graphs. We show how features in the first passage time distribution can be associated with specific kinetic traps, and how the appearance of competing relaxation time scales depends on the starting conditions. The theory is tested for two model landscapes and applied to an atomic cluster and a disordered peptide. Our most important contribution is probably the reconstruction of the full distribution for long time scales, where numerical problems prevent direct calculations. Here we combine accurate treatment of the mean first passage time with the reliable part of the distribution corresponding to faster time scales. Hence we now have a fundamental understanding of both thermodynamic and kinetic signatures of multifunnel landscapes.

Combining experiment and energy landscapes to explore anaerobic heme breakdown in multifunctional hemoproteins.

To survive, many pathogens extract heme from their host organism and break down the porphyrin scaffold to sequester the Fe2+ ion via a heme oxygenase. Recent studies have revealed that certain pathogens can anaerobically degrade heme. Our own research has shown that one such pathway proceeds via NADH-dependent heme degradation, which has been identified in a family of hemoproteins from a range of bacteria. HemS, from Yersinia enterocolitica, is the main focus of this work, along with HmuS (Yersinia pestis), ChuS (Escherichia coli) and ShuS (Shigella dysenteriae). We combine experiments, Energy Landscape Theory, and a bioinformatic investigation to place these homologues within a wider phylogenetic context. A subset of these hemoproteins are known to bind certain DNA promoter regions, suggesting not only that they can catalytically degrade heme, but that they are also involved in transcriptional modulation responding to heme flux. Many of the bacterial species responsible for these hemoproteins (including those that produce HemS, ChuS and ShuS) are known to specifically target oxygen-depleted regions of the gastrointestinal tract. A deeper understanding of anaerobic heme breakdown processes exploited by these pathogens could therefore prove useful in the development of future strategies for disease prevention.

Energy landscapes for proteins described by the UNRES coarse-grained potential.

The self-assembly of proteins is encoded in the underlying potential energy surface (PES), from which we can predict structure, dynamics, and thermodynamic properties. However, the corresponding analysis becomes increasingly challenging with larger protein sizes, due to the computational time required, which grows significantly with the number of atoms. Coarse-grained models offer an attractive approach to reduce the computational cost. In this Feature Article, we describe our implementation of the UNited RESidue (UNRES) coarse-grained potential in the Cambridge energy landscapes software. We have applied this framework to explore the energy landscapes of four proteins that exhibit native states involving different secondary structures. Here we have tested the ability of the UNRES potential to represent the global energy landscape of proteins containing up to 100 amino acid residues. The resulting potential energy landscapes exhibit good agreement with experiment, with low-lying minima close to the PDB geometries and to results obtained using the all-atom AMBER force field. The new program interfaces will allow us to investigate larger biomolecules in future work, using the UNRES potential in combination with all the methodology available in the computational energy landscapes framework.

Global analysis of energy landscapes for materials modeling: A test case for C60.

In this contribution, we employ computational tools from the energy landscape approach to test Gaussian Approximation Potentials (GAPs) for C60. In particular, we apply basin-hopping global optimization and explore the landscape starting from the low-lying minima using discrete path sampling. We exploit existing databases of minima and transition states harvested from previous work using tight-binding potentials. We explore the energy landscape for the full range of structures and pathways spanning from the buckminsterfullerene global minimum up to buckybowls. In the initial GAP model, the fullerene part of the landscape is reproduced quite well. However, there are extensive families of C1@C59 and C2@C58 structures that lie lower in energy. We succeeded in refining the potential to remove these artifacts by simply including two minima from the C2@C58 families found by global landscape exploration. We suggest that the energy landscape approach could be used systematically to test and improve machine learning interatomic potentials.

Energy landscapes and heat capacity signatures for peptides correlate with phase separation propensity.

Phase separation plays an important role in the formation of membraneless compartments within the cell and intrinsically disordered proteins with low-complexity sequences can drive this compartmentalisation. Various intermolecular forces, such as aromatic-aromatic and cation-aromatic interactions, promote phase separation. However, little is known about how the ability of proteins to phase separate under physiological conditions is encoded in their energy landscapes and this is the focus of the present investigation. Our results provide a first glimpse into how the energy landscapes of minimal peptides that contain - and cation- interactions differ from the peptides that lack amino acids with such interactions. The peaks in the heat capacity () as a function of temperature report on alternative low-lying conformations that differ significantly in terms of their enthalpic and entropic contributions. The analysis and subsequent quantification of frustration of the energy landscape suggest that the interactions that promote phase separation lead to features (peaks or inflection points) at low temperatures in . More features may occur for peptides containing residues with better phase separation propensity and the energy landscape is more frustrated for such peptides. Overall, this work links the features in the underlying single-molecule potential energy landscapes to their collective phase separation behaviour and identifies quantities ( and frustration metric) that can be utilised in soft material design.

Exploring energy landscapes for solid-state systems with variable cells at the extended tight-binding level.

The design of novel materials requires a theoretical understanding of dynamical processes in the solid state, including polymorphic transitions and associated pathways. The organization of the potential energy landscape plays a crucial role in such processes, which may involve changes in the periodic boundaries. This study reports the implementation of a general framework for periodic condensed matter systems in our energy landscape analysis software, allowing for variation in both the unit cell and atomic positions. This implementation provides access to basin-hopping global optimization, the doubly nudged elastic band procedure for identifying transition state candidates, the missing connection approach for multi-step pathways, and general tools for the construction and analysis of kinetic transition networks. The computational efficacy of the procedures is explored using the state-of-the-art semiempirical method GFN1-xTB for the first time in this solid-state context. We investigate the effectiveness of this level of theory by characterizing the potential energy and enthalpy landscapes of several systems, including silicon, CdSe, ZnS, and NaCl, and discuss further technical challenges, such as translational permutation of the cell. Despite the expected limitations of the semiempirical method, we find that the resulting energy landscapes provide useful insight into solid-state simulations, which will facilitate detailed analysis of processes such as defect and ion migration, including refinement at higher levels of theory.

Triplet-Encoded Prebiotic RNA Aminoacylation.

The encoding step of translation involves attachment of amino acids to cognate tRNAs by aminoacyl-tRNA synthetases, themselves the product of coded peptide synthesis. So, the question arises─before these enzymes evolved, how were primordial tRNAs selectively aminoacylated? Here, we demonstrate enzyme-free, sequence-dependent, chemoselective aminoacylation of RNA. We investigated two potentially prebiotic routes to aminoacyl-tRNA acceptor stem-overhang mimics and analyzed those oligonucleotides undergoing the most efficient aminoacylation. Overhang sequences do not significantly influence the chemoselectivity of aminoacylation by either route. For aminoacyl-transfer from a mixed anhydride donor strand, the chemoselectivity and stereoselectivity of aminoacylation depend on the terminal three base pairs of the stem. The results support early suggestions of a second genetic code in the acceptor stem.

Energy Landscapes and Heat Capacity Signatures for Monomers and Dimers of Amyloid-Forming Hexapeptides.

Amyloid formation is a hallmark of various neurodegenerative disorders. In this contribution, energy landscapes are explored for various hexapeptides that are known to form amyloids. Heat capacity (CV) analysis at low temperature for these hexapeptides reveals that the low energy structures contributing to the first heat capacity feature above a threshold temperature exhibit a variety of backbone conformations for amyloid-forming monomers. The corresponding control sequences do not exhibit such structural polymorphism, as diagnosed via end-to-end distance and a dihedral angle defined for the monomer. A similar heat capacity analysis for dimer conformations obtained using basin-hopping global optimisation shows clear features in end-to-end distance versus dihedral correlation plots, where amyloid-forming sequences exhibit a preference for larger end-to-end distances and larger positive dihedrals. These results hold true for sequences taken from tau, amylin, insulin A chain, a de novo designed peptide, and various control sequences. While there is a little overall correlation between the aggregation propensity and the temperature at which the low-temperature CV feature occurs, further analysis suggests that the amyloid-forming sequences exhibit the key CV feature at a lower temperature compared to control sequences derived from the same protein.

On the Mechanism of Polaritonic Rate Suppression from Quantum Transition Paths.

Polariton chemistry holds promise for facilitating mode-selective chemical reactions, but the underlying mechanism behind the rate modifications observed under strong vibrational coupling is not well understood. Using the recently developed quantum transition path theory, we have uncovered a mechanism of resonant suppression of a thermal reaction rate in a simple model polaritonic system consisting of a reactive mode in a bath confined to a lossless microcavity with a single photon mode. We observed the formation of a polariton during rate-limiting transitions on reactive pathways and identified the concomitant rate suppression as being due to hybridization between the reactive mode and the cavity mode, which inhibits bath-mediated tunneling. The transition probabilities that define the quantum master equation can be directly translated into a visualization of the corresponding polariton energy landscape. This landscape exhibits a double funnel structure with a large barrier between the initial and final states.

Stable and scalable computation of state visitation probabilities in finite Markov chains.

We report an algorithm based on renormalization to compute the probability that a particular state, or set thereof, is visited along the first passage or transition paths between two endpoint states of a finite Markov chain. The procedure is numerically stable and does not require dense storage of the transition matrix.

Analysing ill-conditioned Markov chains.

Discrete state Markov chains in discrete or continuous time are widely used to model phenomena in the social, physical and life sciences. In many cases, the model can feature a large state space, with extreme differences between the fastest and slowest transition timescales. Analysis of such ill-conditioned models is often intractable with finite precision linear algebra techniques. In this contribution, we propose a solution to this problem, namely partial graph transformation, to iteratively eliminate and renormalize states, producing a low-rank Markov chain from an ill-conditioned initial model. We show that the error induced by this procedure can be minimized by retaining both the renormalized nodes that represent metastable superbasins, and those through which reactive pathways concentrate, i.e. the dividing surface in the discrete state space. This procedure typically returns a much lower rank model, where trajectories can be efficiently generated with kinetic path sampling. We apply this approach to an ill-conditioned Markov chain for a model multi-community system, measuring the accuracy by direct comparison with trajectories and transition statistics. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Molecular Energy Landscapes of Hardware-Efficient Ansätze in Quantum Computing.

Rapid advances in quantum computing have opened up new opportunities for solving the central electronic structure problem in computational chemistry. In the noisy intermediate-scale quantum (NISQ) era, where qubit coherence times are limited, it is essential to exploit quantum algorithms with sufficiently short quantum circuits to maximize qubit efficiency. The procedural construction of hardware-efficient ansätze provides one approach to design such circuits. However, refining the accuracy of the global minimum by increasing circuit depth may lead to a proliferation of local minima that hinders global optimization. To investigate this phenomenon, we explore the energy landscapes of hardware-efficient circuits to identify ground-state energies of the hydrogen, lithium hydride, and beryllium hydride molecules. We also propose a simple dimensionality reduction procedure that reduces quantum gate depth while retaining high accuracy for the global minimum, simplifying the energy landscape, and hence speeding up optimization from both software and hardware perspectives.

Dynamic Diastereomerism on Chiral Surfaces.

Adsorption of chiral molecules on chiral surfaces implies diastereomerism, evident in the adoption of distinct adsorption geometries. We show here that this diastereomerism produces a signature in the motion of chiral molecules desorbing from a chiral surface. The rotations of S- and R-alanine molecules are analyzed upon desorption from R-Cu{531} using first-principles molecular dynamics simulations. S-Ala molecules exhibit a larger angular momentum, with a clear preference for one rotational sense, whereas no such preference is observed for R-Ala molecules upon desorption from this surface. These trends would be reversed for desorption from the S-Cu{531} surface. Possible applications include chiral separation techniques and enantiospecific sensors.

Stochastic paths controlling speed and dissipation.

Natural processes occur in a finite amount of time and dissipate energy, entropy, and matter. Near equilibrium, thermodynamic intuition suggests that fast irreversible processes will dissipate more energy and entropy than slow quasistatic processes connecting the same initial and final states. For small systems, recently discovered thermodynamic speed limits suggest that faster processes will dissipate more than slower processes. Here, we test the hypothesis that this relationship between speed and dissipation holds for stochastic paths far from equilibrium. To analyze stochastic paths on finite timescales, we derive an exact expression for the path probabilities of continuous-time Markov chains from the path summation solution to the master equation. We present a minimal model for a driven system in which relative energies of the initial and target states control the speed, and the nonequilibrium currents of a cycle control the dissipation. Although the hypothesis holds near equilibrium, we find that faster processes can dissipate less under far-from-equilibrium conditions because of strong currents. This model serves as a minimal prototype for designing kinetics to sculpt the nonequilibrium path space so that faster paths produce less dissipation.

Design of self-assembling mesoscopic Goldberg polyhedra.

Palladium ions complexed with nonlinear bidentate ligands have been shown to form hollow, spherical shells with high symmetries. We show that such structures can be reproduced using model anisotropic mesoscale building blocks featuring excluded volume and long-range ionic interactions. A linear building block with a central charged particle, in combination with a bent 'ligand' particle with opposite charges at the ends is sufficient to drive the system towards planar coordination, and the charge ratio determines the coordination number. Similar to the molecular systems, the bend in the 'ligand' particle determines the curvature of the shells that these building blocks prefer. Besides reproducing exotic structures such as M30L60 and M48L96 tetravalent Goldberg polyhedra, we identify highly cooperative single transition state rearrangements between low-energy competing structures as well, corresponding to rotatory motions of a planar subunit within the spherical shell.