Search

The principle of microscopic reversibility states that at equilibrium the number of molecules entering a state by a given path must equal those exiting the state via the same path under identical conditions or, in structural terms, that the conformations along the two pathways are the same. There has been some indirect evidence indicating that protein folding is such a process, but there have been few conclusive findings. In this study, we performed molecular dynamics simulations of an ultrafast unfolding and folding protein at its melting temperature to observe, on an atom-by-atom basis, the pathways the protein followed as it unfolded and folded within a continuous trajectory. In a total of 0.67 µs of simulation in water, we found six transient denaturing events near the melting temperature (323 and 330 K) and an additional refolding event following a previously identified unfolding event at a high temperature (373 K). In each case, unfolding and refolding transition state ensembles were identified, and they agreed well with experiment on the basis of a comparison of S and Φ values. On the basis of several structural properties, these 13 transition state ensembles agreed very well with each other and with four previously identified transition states from high-temperature denaturing simulations. Thus, not only were the unfolding and refolding transition states part of the same ensemble, but in five of the seven cases, the pathway the protein took as it unfolded was nearly identical to the subsequent refolding pathway. These events provide compelling evidence that protein folding is a microscopically reversible process. In the other two cases, the folding and unfolding transition states were remarkably similar to each other but the paths deviated. [ABSTRACT FROM AUTHOR]

The correct treatment of van der Waals and electrostatic nonbonded interactions in molecular force fields is essential for performing realistic molecular dynamics (MD) simulations of solvated polypeptides. The most computationally tractable treatment of nonbonded interactions in MD utilizes a spherical distance cutoff (typically, 8-12 A) to reduce the number of pairwise interactions. In this work, we assess three spherical atom-based cutoff approaches for use with all-atom explicit solvent MD: abrupt truncation, a CHARMM-style electrostatic shift truncation, and our own force-shifted truncation. The chosen system for this study is an end-capped 17-residue alanine-based a-helical peptide, selected because of its use in previous computational and experimental studies. We compare the time-averaged helical content calculated from these MD trajectories with experiment. We also examine the effect of varying the cutoff treatment and distance on energy conservation. We find that the abrupt truncation approach is pathological in its inability to conserve energy. The CHARMM-style shift truncation performs quite well but suffers from energetic instability. On the other hand, the force-shifted spherical cutoff method conserves energy, correctly predicts the experimental helical content, and shows convergence in simulation statistics as the cutoff is increased. This work demonstrates that by using proper and rigorous techniques, it is possible to correctly model polypeptide dynamics in solution with a spherical cutoff. The inherent computational advantage of spherical cutoffs over Ewald summation (and related) techniques is essential in accessing longer MD time scales. [ABSTRACT FROM AUTHOR]

Halide perovskites have the potential to disrupt the photovoltaics market based on their high performance and low cost. However, the decomposition of perovskites under moisture, oxygen, and light raises concerns about service lifetime, especially because degradation mechanisms and the corresponding rate laws that fit the observed data have thus far eluded researchers. Here, we report a water-accelerated photooxidation mechanism dominating the degradation kinetics of archetypal perovskite CH 3 NH 3 PbI 3 in air under >1% relative humidity at 25 °C. From this mechanism, we develop a kinetic model that quantitatively predicts the degradation rate as a function of temperature, ambient O 2 and H 2 O levels, and illumination. Because water is a possible product of dry photooxidation, these results highlight the need for encapsulation schemes that rigorously block oxygen ingress, as product water may accumulate beneath the encapsulant and initiate the more rapid water-accelerated photooxidative decomposition.

Engineering microorganisms into biological factories that convert renewable feedstocks into valuable materials is a major goal of synthetic biology; however, for many nonmodel organisms, we do not yet have the genetic tools, such as suites of strong promoters, necessary to effectively engineer them. In this work, we developed a computational framework that can leverage standard RNA-seq data sets to identify sets of constitutive, strongly expressed genes and predict strong promoter signals within their upstream regions. The framework was applied to a diverse collection of RNA-seq data measured for the methanotroph Methylotuvimicrobium buryatense 5GB1 and identified 25 genes that were constitutively, strongly expressed across 12 experimental conditions. For each gene, the framework predicted short (27-30 nucleotide) sequences as candidate promoters and derived -35 and -10 consensus promoter motifs (TTGACA and TATAAT, respectively) for strong expression in M. buryatense . This consensus closely matches the canonical E. coli sigma-70 motif and was found to be enriched in promoter regions of the genome. A subset of promoter predictions was experimentally validated in a XylE reporter assay, including the consensus promoter, which showed high expression. The pmoC , pqqA , and ssrA promoter predictions were additionally screened in an experiment that scrambled the -35 and -10 signal sequences, confirming that transcription initiation was disrupted when these specific regions of the predicted sequence were altered. These results indicate that the computational framework can make biologically meaningful promoter predictions and identify key pieces of regulatory systems that can serve as foundational tools for engineering diverse microorganisms for biomolecule production.

Designing new ionic liquids (ILs) is of crucial importance for various industrial applications. However, this always leads to a daunting challenge, as the number of possible combinations of cation and anion are very high and it is impossible to experimentally propose and screen a wide pool of potential candidates. However, recent applications of machine learning (ML) models have greatly improved the overall chemical discovery pipeline. In this study, we compare different generative methods for producing ionic liquids. In this comparison, we show the following: (1) when training data is scarce, a transfer learning approach can be applied to variational autoencoders (VAEs) to generate molecular structures of the target molecule type; (2) in a VAE-like structure, separate latent spaces for the cationic and anionic moieties can result in meaningful representations for their combinative, macroscopic properties; (3) interpolating between ILs with desired properties can result in a new IL with attributes similar to the two structural end points.