Impact of Vibrational Entropy on the Stability of Unsolvated Peptide Helices with Increasing Length

Ultralight Ultrafast Enzymes

Inorganic materials depleted of heavy stable isotopes are known to deviate strongly in some physicochemical properties from their isotopically natural counterparts. Here we explored for the first time the effect of simultaneous depletion of the heavy carbon, hydrogen, oxygen and nitrogen isotopes on the bacterium E. coli and the enzymes expressed in it. Bacteria showed faster growth, with most proteins exhibiting higher thermal stability, while for recombinant enzymes expressed in depleted media, faster kinetics was discovered. At room temperature, luciferase, thioredoxin and dihydrofolate reductase and Pfu DNA polymerase showed up to a 250 % increase in activity compared to the native counterparts, with an additional ∼50 % increase at 10 °C. Diminished conformational and vibrational entropy is hypothesized to be the cause of the accelerated kinetics. Ultralight enzymes may find an application where extreme reaction rates are required.

Exploring the Conformers of an Organic Molecule on a Metal Cluster with Bayesian Optimization

Finding low-energy conformers of organic molecules is a complex problem due to the flexibilities of the molecules and the high dimensionality of the search space. When such molecules are on nanoclusters, the search complexity is exacerbated by constraints imposed by the presence of the cluster and other surrounding molecules. To address this challenge, we modified our previously developed active learning molecular conformer search method based on Bayesian optimization and density functional theory. Especially, we have developed and tested strategies to avoid steric clashes between a molecule and a cluster. In this work, we chose a cysteine molecule on a well-studied gold-thiolate cluster as a model system to test and demonstrate our method. We found that cysteine conformers in a cluster inherit the hydrogen bond types from isolated conformers. However, the energy rankings and spacings between the conformers are reordered.

Propensities of Some Amino Acid Pairings in α-Helices Vary with Length

The results of secondary structure prediction methods are widely used in applications in biotechnology and bioinformatics. However, the accuracy limit of these methods could be improved up to 92%. One approach to achieve this goal is to harvest information from the primary structure of the peptide. This study aims to contribute to this goal by investigating the variations in propensity of amino acid pairings to α-helices in globular proteins depending on helix length. (n):(n + 4) residue pairings were determined using a comprehensive peptide data set according to backbone hydrogen bond criterion which states that backbone hydrogen bond is the dominant driving force of protein folding. Helix length is limited to 13 to 26 residues. Findings of this study show that propensities of ALA:GLY and GLY:GLU pairings to α-helix in globular protein increase with increasing helix length but of ALA:ALA and ALA:VAL decrease. While the frequencies of ILE:ALA, LEU:ALA, LEU:GLN, LEU:GLU, LEU:LEU, MET:ILE and VAL:LEU pairings remain roughly constant with length, the 25 residue pairings have varying propensities in narrow helix lengths. The remaining pairings have no prominent propensity to α-helices.

Kinetically Controlled Structural Transitions in Layered Halide-Based Perovskites: An Approach to Modulate Spin Splitting

Two-dimensional hybrid organic-inorganic perovskite (HOIP) semiconductors with pronounced spin splitting, mediated by strong spin-orbit coupling and inversion symmetry breaking, offer the potential for spin manipulation in future spintronic applications. However, HOIPs exhibiting significant conduction/valence band splitting are still relatively rare, given the generally observed preference for (near)centrosymmetric inorganic (especially lead-iodide-based) sublattices, and few approaches are available to control this symmetry breaking within a given HOIP. Here, we demonstrate, using (S-2-MeBA)₂PbI₄ (S-2-MeBA = (S)-(-)-2-methylbutylammonium) as an example, that a temperature-induced structural transition (at ∼180 K) serves to change the degree of chirality transfer to and inversion symmetry breaking within the inorganic layer, thereby enabling modulation of HOIP structural and electronic properties. The cooling rate is shown to dictate whether the structural transition occurs─i.e., slow cooling induces the transition while rapid quenching inhibits it. Ultrafast calorimetry indicates a minute-scale structural relaxation time at the transition temperature, while quenching to lower temperatures allows for effectively locking in the metastable room-temperature phase, thus enabling kinetic control over switching between distinct states with different degrees of structural distortions within the inorganic layers at these temperatures. Density functional theory further highlights that the low-temperature phase of (S-2-MeBA)₂PbI₄ shows more significant spin splitting relative to the room-temperature phase. Our work opens a new pathway to use kinetic control of crystal-to-crystal transitions and thermal cycling to modulate spin splitting in HOIPs for future spintronic applications, and further points to using such "sluggish" phase transitions for switching and control over other physical phenomena, particularly those relying on structural distortions and lattice symmetry.

Molecular Conformer Search with Low-Energy Latent Space

Identifying low-energy conformers with quantum mechanical accuracy for molecules with many degrees of freedom is challenging. In this work, we use the molecular dihedral angles as features and explore the possibility of performing molecular conformer search in a latent space with a generative model named variational auto-encoder (VAE). We bias the VAE towards low-energy molecular configurations to generate more informative data. In this way, we can effectively build a reliable energy model for the low-energy potential energy surface. After the energy model has been built, we extract local-minimum conformations and refine them with structure optimization. We have tested and benchmarked our low-energy latent-space (LOLS) structure search method on organic molecules with 5-9 searching dimensions. Our results agree with previous studies.

Efficient Amino Acid Conformer Search with Bayesian Optimization

Finding low-energy molecular conformers is challenging due to the high dimensionality of the search space and the computational cost of accurate quantum chemical methods for determining conformer structures and energies. Here, we combine active-learning Bayesian optimization (BO) algorithms with quantum chemistry methods to address this challenge. Using cysteine as an example, we show that our procedure is both efficient and accurate. After only 1000 single-point calculations and approximately 80 structure relaxations, which is less than 10% computational cost of the current fastest method, we have found the low-energy conformers in good agreement with experimental measurements and reference calculations. To test the transferability of our method, we also repeated the conformer search of serine, tryptophan, and aspartic acid. The results agree well with previous conformer search studies.

On the Role of Entropy in the Stabilization of α-Helices

Protein folding evolves by exploring the conformational space with a subtle balance between enthalpy and entropy changes which eventually leads to a decrease of free energy upon reaching the folded structure. A complete understanding of this process requires, therefore, a deep insight into both contributions to free energy. In this work, we clarify the role of entropy in favoring the stabilization of folded structures in polyalanine peptides with up to 12 residues. We use a novel method referred to as K2V that allows us to obtain the potential-energy landscapes in terms of residue conformations extracted from molecular dynamics simulations at conformational equilibrium and yields folding thermodynamic magnitudes, which are in agreement with the experimental data available. Our results demonstrate that the folded structures of the larger polyalanine chains are stabilized with respect to the folded structures of the shorter chains by both an energetic contribution coming from the formation of the intramolecular hydrogen bonds and an entropic contribution coming from an increase of the entropy of the solvent with approximate weights of 60 and 40%, respectively, thus unveiling a key piece in the puzzle of protein folding. In addition, the ability of the K2V method to provide the enthalpic and entropic contributions for individual residues along the peptide chain makes it clear that the energetic and entropic stabilizations are basically governed by the nearest neighbor residue conformations, with the folding propensity being rationalized in terms of triads of residues.

Entropy of Proteins Using Multiscale Cell Correlation

A new multiscale method is presented to calculate the entropy of proteins from molecular dynamics simulations. Termed Multiscale Cell Correlation (MCC), the method decomposes the protein into sets of rigid-body units based on their covalent-bond connectivity at three levels of hierarchy: molecule, residue, and united atom. It evaluates the vibrational and topographical entropy from forces, torques, and dihedrals at each level, taking into account correlations between sets of constituent units that together make up a larger unit at the coarser length scale. MCC gives entropies in close agreement with normal-mode analysis and smaller than those using quasiharmonic analysis as well as providing much faster convergence. Moreover, MCC provides an insightful decomposition of entropy at each length scale and for each type of amino acid according to their solvent exposure and whether they are terminal residues. While the residue entropy depends weakly on solvent exposure, there is greater variation in entropy components for larger, more polar amino acids, which have increased conformational entropy but reduced vibrational entropy with greater solvent exposure.

Probing the conformational landscape and thermochemistry of DNA dinucleotide anions via helium nanodroplet infrared action spectroscopy

Isolation of biomolecules in vacuum facilitates characterization of the intramolecular interactions that determine three-dimensional structure, but experimental quantification of conformer thermochemistry remains challenging. Infrared spectroscopy of molecules trapped in helium nanodroplets is a promising methodology for the measurement of thermochemical parameters. When molecules are captured in a helium nanodroplet, the rate of cooling to an equilibrium temperature of ca. 0.4 K is generally faster than the rate of isomerization, resulting in "shock-freezing" that kinetically traps molecules in local conformational minima. This unique property enables the study of temperature-dependent conformational equilibria via infrared spectroscopy at 0.4 K, thereby avoiding the deleterious effects of spectral broadening at higher temperatures. Herein, we demonstrate the first application of this approach to ionic species by coupling electrospray ionization mass spectrometry (ESI-MS) with helium nanodroplet infrared action spectroscopy to probe the structure and thermochemistry of deprotonated DNA dinucleotides. Dinucleotide anions were generated by ESI, confined in an ion trap at temperatures between 90 and 350 K, and entrained in traversing helium nanodroplets. The infrared action spectra of the entrained ions show a strong dependence on pre-pickup ion temperature, consistent with the preservation of conformer population upon cooling to 0.4 K. Non-negative matrix factorization was utilized to identify component conformer infrared spectra and determine temperature-dependent conformer populations. Relative enthalpies and entropies of conformers were subsequently obtained from a van't Hoff analysis. IR spectra and conformer thermochemistry are compared to results from ion mobility spectrometry (IMS) and electronic structure methods. The implementation of ESI-MS as a source of dopant molecules expands the diversity of molecules accessible for thermochemical measurements, enabling the study of larger, non-volatile species.

Quantum mechanics of proteins in explicit water: The role of plasmon-like solute-solvent interactions

Quantum-mechanical van der Waals dispersion interactions play an essential role in intraprotein and protein-water interactions-the two main factors affecting the structure and dynamics of proteins in water. Typically, these interactions are only treated phenomenologically, via pairwise potential terms in classical force fields. Here, we use an explicit quantum-mechanical approach of density-functional tight-binding combined with the many-body dispersion formalism and demonstrate the relevance of many-body van der Waals forces both to protein energetics and to protein-water interactions. In contrast to commonly used pairwise approaches, many-body effects substantially decrease the relative stability of native states in the absence of water. Upon solvation, the protein-water dispersion interaction counteracts this effect and stabilizes native conformations and transition states. These observations arise from the highly delocalized and collective character of the interactions, suggesting a remarkable persistence of electron correlation through aqueous environments and providing the basis for long-range interaction mechanisms in biomolecular systems.

Decisive role of nuclear quantum effects on surface mediated water dissociation at finite temperature

Water molecules adsorbed on inorganic substrates play an important role in several technological applications. In the presence of light atoms in adsorbates, nuclear quantum effects (NQEs) influence the structural stability and the dynamical properties of these systems. In this work, we explore the impact of NQEs on the dissociation of water wires on stepped Pt(221) surfaces. By performing ab initio molecular dynamics simulations with van der Waals corrected density functional theory, we note that several competing minima for both intact and dissociated structures are accessible at finite temperatures, making it important to assess whether harmonic estimates of the quantum free energy are sufficient to determine the relative stability of the different states. We thus perform ab initio path integral molecular dynamics (PIMD) in order to calculate these contributions taking into account the conformational entropy and anharmonicities at finite temperatures. We propose that when adsorption is weak and NQEs on the substrate are negligible, PIMD simulations can be performed through a simple partition of the system, resulting in considerable computational savings. We then calculate the full contribution of NQEs to the free energies, including also anharmonic terms. We find that they result in an increase of up to 20% of the quantum contribution to the dissociation free energy compared with the harmonic estimates. We also find that the dissociation process has a negligible contribution from tunneling but is dominated by zero point energies, which can enhance the rate of dissociation by three orders of magnitude. Finally we highlight how both temperature and NQEs indirectly impact dipoles and the redistribution of electron density, causing work function changes of up to 0.4 eV with respect to static estimates. This quantitative determination of the change in the work function provides a possible approach to determine experimentally the most stable configurations of water oligomers on the stepped surfaces.

Testing the mutual information expansion of entropy with multivariate Gaussian distributions

The mutual information expansion (MIE) represents an approximation of the configurational entropy in terms of low-dimensional integrals. It is frequently employed to compute entropies from simulation data of large systems, such as macromolecules, for which brute-force evaluation of the full configurational integral is intractable. Here, we test the validity of MIE for systems consisting of more than m = 100 degrees of freedom (dofs). The dofs are distributed according to multivariate Gaussian distributions which were generated from protein structures using a variant of the anisotropic network model. For the Gaussian distributions, we have semi-analytical access to the configurational entropy as well as to all contributions of MIE. This allows us to accurately assess the validity of MIE for different situations. We find that MIE diverges for systems containing long-range correlations which means that the error of consecutive MIE approximations grows with the truncation order n for all tractable n ≪ m. This fact implies severe limitations on the applicability of MIE, which are discussed in the article. For systems with correlations that decay exponentially with distance, MIE represents an asymptotic expansion of entropy, where the first successive MIE approximations approach the exact entropy, while MIE also diverges for larger orders. In this case, MIE serves as a useful entropy expansion when truncated up to a specific truncation order which depends on the correlation length of the system.

Probing chirality recognition of protonated glutamic acid dimers by gas-phase vibrational spectroscopy and first-principles simulations

We characterize stereospecific aspects of homochiral and heterochiral dimers of glutamic acid by infrared spectroscopy and first-principles molecular dynamics simulations.

Vibrational entropy of a protein: large differences between distinct conformations

In this article, it is investigated whether vibrational entropy (VE) is an important contribution to the free energy of globular proteins at ambient conditions. VE represents the major configurational-entropy contribution of these proteins. By definition, it is an average of the configurational entropies of the protein within single minima of the energy landscape, weighted by their occupation probabilities. Its large part originates from thermal motion of flexible torsion angles giving rise to the finite peak widths observed in torsion angle distributions. While VE may affect the equilibrium properties of proteins, it is usually neglected in numerical calculations as its consideration is difficult. Moreover, it is sometimes believed that all well-packed conformations of a globular protein have similar VE anyway. Here, we measure explicitly the VE for six different conformations from simulation data of a test protein. Estimates are obtained using the quasi-harmonic approximation for three coordinate sets, Cartesian, bond-angle-torsion (BAT), and a new set termed rotamer-degeneracy lifted BAT coordinates by us. The new set gives improved estimates as it overcomes a known shortcoming of the quasi-harmonic approximation caused by multiply populated rotamer states, and it may serve for VE estimation of macromolecules in a very general context. The obtained VE values depend considerably on the type of coordinates used. However, for all coordinate sets we find large entropy differences between the conformations, of the order of the overall stability of the protein. This result may have important implications on the choice of free energy expressions used in software for protein structure prediction, protein design, and NMR refinement.

Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol

Molecular crystals often exist in multiple competing polymorphs, showing significantly different physicochemical properties. Computational crystal structure prediction is key to interpret and guide the search for the most stable or useful form, a real challenge due to the combinatorial search space, and the complex interplay of subtle effects that work together to determine the relative stability of different structures. Here we take a comprehensive approach based on different flavors of thermodynamic integration in order to estimate all contributions to the free energies of these systems with density-functional theory, including the oft-neglected anharmonic contributions and nuclear quantum effects. We take the two main stable forms of paracetamol as a paradigmatic example. We find that anharmonic contributions, different descriptions of van der Waals interactions, and nuclear quantum effects all matter to quantitatively determine the stability of different phases. Our analysis highlights the many challenges inherent in the development of a quantitative and predictive framework to model molecular crystals. However, it also indicates which of the components of the free energy can benefit from a cancellation of errors that can redeem the predictive power of approximate models, and suggests simple steps that could be taken to improve the reliability of ab initio crystal structure prediction.

First-principles data set of 45,892 isolated and cation-coordinated conformers of 20 proteinogenic amino acids

We present a structural data set of the 20 proteinogenic amino acids and their amino-methylated and acetylated (capped) dipeptides. Different protonation states of the backbone (uncharged and zwitterionic) were considered for the amino acids as well as varied side chain protonation states. Furthermore, we studied amino acids and dipeptides in complex with divalent cations (Ca²⁺, Ba²⁺, Sr²⁺, Cd²⁺, Pb²⁺, and Hg²⁺). The database covers the conformational hierarchies of 280 systems in a wide relative energy range of up to 4 eV (390 kJ/mol), summing up to a total of 45,892 stationary points on the respective potential-energy surfaces. All systems were calculated on equal first-principles footing, applying density-functional theory in the generalized gradient approximation corrected for long-range van der Waals interactions. We show good agreement to available experimental data for gas-phase ion affinities. Our curated data can be utilized, for example, for a wide comparison across chemical space of the building blocks of life, for the parametrization of protein force fields, and for the calculation of reference spectra for biophysical applications.

Beyond static structures: Putting forth REMD as a tool to solve problems in computational organic chemistry

Computational studies of organic systems are frequently limited to static pictures that closely align with textbook style presentations of reaction mechanisms and isomerization processes. Of course, in reality chemical systems are dynamic entities where a multitude of molecular conformations exists on incredibly complex potential energy surfaces (PES). Here, we borrow a computational technique originally conceived to be used in the context of biological simulations, together with empirical force fields, and apply it to organic chemical problems. Replica-exchange molecular dynamics (REMD) permits thorough exploration of the PES. We combined REMD with density functional tight binding (DFTB), thereby establishing the level of accuracy necessary to analyze small molecular systems. Through the study of four prototypical problems: isomer identification, reaction mechanisms, temperature-dependent rotational processes, and catalysis, we reveal new insights and chemistry that likely would be missed using static electronic structure computations. The REMD-DFTB methodology at the heart of this study is powered by i-PI, which efficiently handles the interface between the DFTB and REMD codes.

Probing the role of electrostatics of polypeptide main-chain in protein folding by perturbing N-terminal residue stereochemistry: DFT study with oligoalanine models

Energetics of folding (ΔH_E→F, in kcal mol⁻¹) from the extended (E) structure to the folded (F) structure for Ia and Ib critically depend on the geometrical relationship between the backbone peptide units of the polypeptide structure.

Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.

Stability of Complex Biomolecular Structures: van der Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects

Biomolecules are complex systems stabilized by a delicate balance of weak interactions, making it important to assess all energetic contributions in an accurate manner. However, it is a priori unclear which contributions make more of an impact. Here, we examine stacked polyglutamine (polyQ) strands, a peptide repeat often found in amyloid aggregates. We investigate the role of hydrogen bond (HB) cooperativity, van der Waals (vdW) dispersion interactions, and quantum contributions to free energies, including anharmonicities through density functional theory and ab initio path integral simulations. Of these various factors, we find that the largest impact on structural stabilization comes from vdW interactions. HB cooperativity is the second largest contribution as the size of the stacked chain grows. Competing nuclear quantum effects make the net quantum contribution small but very sensitive to anharmonicities, vdW, and the number of HBs. Our results suggest that a reliable treatment of these systems can only be attained by considering all of these components.

Native like helices in a specially designed β peptide in the gas phase

First principles simulations and gas phase spectroscopy suggest equilibrium of helices for an oligomer of open chain β amino acids.

Exploring the conformational preferences of 20-residue peptides in isolation: Ac-Ala19-Lys + H+vs. Ac-Lys-Ala19 + H+ and the current reach of DFT

Using a high-level density functional and an exhaustive search of conformation space, the predicted conformation of a 20-amino acid peptide explains two seemingly contradictory experiments.