RESUMO
We present the results for energy landscapes of hexapeptides obtained using interfaces to the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) program. We have used basin-hopping global optimization and discrete path sampling to explore the landscapes of hexapeptide monomers, dimers, and oligomers containing 10, 100, and 200 monomers modeled using a residue-level coarse-grained potential, Mpipi, implemented in LAMMPS. We find that the dimers of peptides containing amino acid residues that are better at promoting phase separation, such as tyrosine and arginine, have melting peaks at higher temperature in their heat capacity compared to phenylalanine and lysine, respectively. This observation correlates with previous work on the same uncapped hexapeptide monomers modeled using atomistic potential. For oligomers, we compare the variation in monomer conformations with radial distance and observe trends for selected angles calculated for each monomer. The LAMMPS interfaces to the GMIN and OPTIM programs for landscape exploration offer new opportunities to investigate larger systems and provide access to the coarse-grained potentials implemented within LAMMPS.
RESUMO
We construct and analyze disconnectivity graphs to provide the first graphical representation of the vibrational energy landscape of a protein, in this study ß2AR, a G-protein coupled receptor (GPCR), in active and inactive states. The graphs, which indicate the relative free energy of each residue and the minimum free energy barriers for energy transfer between them, reveal important composition, structural and dynamic properties that mediate the flow of energy. Prolines and glycines, which contribute to GPCR plasticity and function, are identified as bottlenecks to energy transport along the backbone from which alternative pathways for energy transport via nearby noncovalent contacts emerge, seen also in the analysis of first passage time (FPT) distributions presented here. Striking differences between the disconnectivity graphs and FPT distributions for the inactive and active states of ß2AR are found where structural and dynamic changes occur upon activation, contributing to allosteric regulation.
Assuntos
Termodinâmica , Vibração , Receptores Adrenérgicos beta 2/química , Receptores Adrenérgicos beta 2/metabolismo , Receptores Acoplados a Proteínas G/metabolismo , Receptores Acoplados a Proteínas G/química , Transferência de Energia , Simulação de Dinâmica MolecularRESUMO
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.
Assuntos
Peptídeo 1 Semelhante ao Glucagon , Peptídeo 1 Semelhante ao Glucagon/química , Peptídeo 1 Semelhante ao Glucagon/metabolismo , Concentração de Íons de Hidrogênio , Termodinâmica , Simulação de Dinâmica Molecular , Agregados ProteicosRESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
Assuntos
Dissulfetos , Dobramento de Proteína , Animais , Bovinos , Aprotinina/química , Cistina/química , Dissulfetos/química , ProteínasRESUMO
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.
RESUMO
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.
Assuntos
Hemeproteínas , Anaerobiose , Filogenia , Hemeproteínas/metabolismo , Heme/metabolismo , Escherichia coli/metabolismoRESUMO
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.
Assuntos
Proteínas , Software , Conformação Proteica , Proteínas/química , Estrutura Secundária de Proteína , Termodinâmica , Simulação de Dinâmica MolecularRESUMO
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.
RESUMO
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.
RESUMO
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.
RESUMO
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.
Assuntos
Temperatura Alta , Doenças Neurodegenerativas , Humanos , Amiloide/química , Proteínas Amiloidogênicas , TemperaturaRESUMO
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.
RESUMO
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.
Assuntos
Aminoacil-tRNA Sintetases , RNA , RNA/metabolismo , Aminoacilação , Sequência de Bases , Código Genético , RNA de Transferência/química , Aminoacil-tRNA Sintetases/metabolismo , Conformação de Ácido NucleicoRESUMO
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.
RESUMO
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'.
RESUMO
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.
RESUMO
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.