Higher Energy Conformers Research Articles

State-specific morphological deformations of the lipid bilayer explain mechanosensitive gating of MSCs ion channels.

The so-called force-from-lipids hypothesis of cellular mechanosensation posits that ion channel gating directly responds to changes in the physical state of the membrane induced for example by lateral tension. Here, we investigate the molecular basis for this transduction mechanism by studying the mechanosensitive ion channel MscS from Escherichia coli and its eukaryotic homolog, MSL1 from Arabidopsis thaliana. First, we use single-particle cryo-EM to determine the structure of a novel open conformation of wild-type MscS, stabilize in thinned lipid nanodisc. Compared with the closed state, the structure shows a reconfiguration of helices TM1, TM2 and TM3a, resulting in widening of the central pore. Based on these structures, we examined how the morphology of the lipid bilayer is altered upon gating, using molecular dynamics simulations. The simulations reveal that, in the absence of extrinsic forces, closed-state MscS causes drastic protrusions in the inner leaflet of the lipid bilayer. These deformations develop to provide adequate solvation to hydrophobic features of the protein surface, and clearly reflect a high energy conformation for the membrane. Strikingly, the protrusions are largely eradicated upon channel opening. An analogous computational study of open and closed states of MSL1 recapitulates these findings. The gating equilibrium of MscS channels thus appears to be dictated by two opposing conformational preferences, namely those of the lipid membrane and of the protein structure. We theorize that membrane tension controls this balance because it increases the energetic cost of the membrane deformations that specifically stabilize the closed channel, shifting the gating equilibrium towards the conductive form.

State-specific morphological deformations of the lipid bilayer explain mechanosensitive gating of MscS ion channels.

The force-from-lipids hypothesis of cellular mechanosensation posits that membrane channels open and close in response to changes in the physical state of the lipid bilayer, induced for example by lateral tension. Here, we investigate the molecular basis for this transduction mechanism by studying the mechanosensitive ion channel MscS from Escherichia coli and its eukaryotic homolog MSL1 from Arabidopsis thaliana. First, we use single-particle cryo-electron microscopy to determine the structure of a novel open conformation of wild-type MscS, stabilized in a thinned lipid nanodisc. Compared with the closed state, the structure shows a reconfiguration of helices TM1, TM2, and TM3a, and widening of the central pore. Based on these structures, we examined how the morphology of the membrane is altered upon gating, using molecular dynamics simulations. The simulations reveal that closed-state MscS causes drastic protrusions in the inner leaflet of the lipid bilayer, both in the absence and presence of lateral tension, and for different lipid compositions. These deformations arise to provide adequate solvation to hydrophobic crevices under the TM1-TM2 hairpin, and clearly reflect a high-energy conformation for the membrane, particularly under tension. Strikingly, these protrusions are largely eradicated upon channel opening. An analogous computational study of open and closed MSL1 recapitulates these findings. The gating equilibrium of MscS channels thus appears to be dictated by opposing conformational preferences, namely those of the lipid membrane and of the protein structure. We propose a membrane deformation model of mechanosensation, which posits that tension shifts the gating equilibrium towards the conductive state not because it alters the mode in which channel and lipids interact, but because it increases the energetic cost of the morphological perturbations in the membrane required by the closed state.

Energy landscapes from cryo-EM snapshots: a benchmarking study

Biomolecules undergo continuous conformational motions, a subset of which are functionally relevant. Understanding, and ultimately controlling biomolecular function are predicated on the ability to map continuous conformational motions, and identify the functionally relevant conformational trajectories. For equilibrium and near-equilibrium processes, function proceeds along minimum-energy pathways on one or more energy landscapes, because higher-energy conformations are only weakly occupied. With the growing interest in identifying functional trajectories, the need for reliable mapping of energy landscapes has become paramount. In response, various data-analytical tools for determining structural variability are emerging. A key question concerns the veracity with which each data-analytical tool can extract functionally relevant conformational trajectories from a collection of single-particle cryo-EM snapshots. Using synthetic data as an independently known ground truth, we benchmark the ability of four leading algorithms to determine biomolecular energy landscapes and identify the functionally relevant conformational paths on these landscapes. Such benchmarking is essential for systematic progress toward atomic-level movies of continuous biomolecular function.

The conformational behavior of N‑ethylformamide as observed by rotational spectroscopy and quantum chemistry

A peptide linkage (CO)NH containing molecule, N-ethylformamide, was investigated by rotational spectroscopy in order to determine the molecular constants of its highest-energy conformer, cis-ac. Its rotational spectrum was observed in two different frequency ranges, in the 4–26 GHz frequency region using a Fourier transform microwave spectrometer and at millimeter wave frequencies between 75 and 116 GHz, employing a broadband high-resolution rotational spectrometer. The measurements at low frequencies allowed to resolve the hyperfine structure components due to nitrogen nuclear quadrupole coupling while the data at higher frequencies provided spectroscopic information about high order centrifugal effects. From a merged fit using all the observational data we have determined a total of thirteen molecular constants that provide a more accurate spectral modelling of the cis-ac conformer and serves a basis for their astronomical search. We have also observed spectra of five singly substituted isotopologues for the cis-ac conformer, three 13C and one for each of 15N and the deuterated species on the N–D position, from which we derived a partial r0 structure, in fair agreement with an ab initio result. In addition, the rotational transitions of the deuterated species of the most stable trans-sc conformer were observed and assigned and three rotational, five centrifugal distortion constants and nuclear quadrupole coupling constants of the nitrogen and deuterium nuclei were determined.

Infrared Spectra and Phototransformations of meta-Fluorophenol Isolated in Argon and Nitrogen Matrices

Monomers of meta-fluorophenol (mFP) were trapped from the gas phase into cryogenic argon and nitrogen matrices. The estimated relative energies of the two conformers are very close, and in the gas phase they have nearly equal populations. Due to the similarity of their structures (they only differ in the orientation of the OH group), the two conformers have also similar predicted vibrational signatures, which makes the vibrational characterization of the individual rotamers challenging. In the present work, it has been established that in an argon matrix only the most stable trans conformer of mFP exists (the OH group pointing away from the fluorine atom). On the other hand, the IR spectrum of mFP in a nitrogen matrix testifies to the simultaneous presence in this matrix of both the trans conformer and of the higher-energy cis conformer (the OH group pointing toward the fluorine atom), which is stabilized by interaction with the matrix gas host. We found that the exposition of the cryogenic N2 matrix to the Globar source of the infrared spectrometer affects the conformational populations. By collecting experimental spectra, either in the full mid-infrared range or only in the range below 2200 cm-1, we were able to reliably distinguish two sets of experimental bands originating from individual conformers. A comparison of the two sets of experimental bands with computed infrared spectra of the conformers allowed, for the first time, the unequivocal vibrational identification of each of them. The joint implementation of computational vibrational spectroscopy and matrix-isolation infrared spectroscopy proved to be a very accurate method of structural analysis. Some mechanistic insights into conformational isomerism (the quantum tunneling of hydrogen atom and vibrationally-induced conformational transformations) have been addressed. Finally, we also subjected matrix-isolated mFP to irradiations with UV light, and the phototransformations observed in these experiments are also described.

Breaking and Formation of Intramolecular Hydrogen Bonds in Dihydroxybenzaldehydes through UV-Induced Conformational Changes in a Low-Temperature Matrix.

Hydrogen bonding (HB) has been receiving attention from both experimental and theoretical researchers since its discovery in the 1920s due to its impact on many chemical and biological processes. However, despite the large number of investigations conducted on this topic, the nature of the HBs and, in particular, the estimation of intramolecular HB energies are still very active subjects of research. In this context, here we report a matrix isolation infrared spectroscopy study of 2,3-dihydroxybenzaldehyde (2,3-DHBA) and 2,4-dihydroxybenzaldehyde (2,4-DHBA), which contain two [one resonance-assisted HB (RAHB) and one conventional HB] and one (RAHB) intramolecular hydrogen bonds, respectively, in their most stable conformer. After isolation of the compounds in cryogenic (15 K) krypton matrices, ultraviolet irradiation led to the formation of higher-energy conformers (by a 180° rotation of the hydroxyl and aldehyde groups), which implies breaking of the intramolecular HBs initially existing in the isolated species and, in the case of 2,3-DHBA, to the formation of a new intramolecular HB. In this way, we were able to manipulate the structure of the molecules, allowing to characterize a diversity of intramolecular HBs in which the OH groups participate (from strong intramolecular RAHBs to weaker conventional HBs, and also no intramolecular HBs) through the corresponding vibrational signatures. The spectroscopic studies were complemented by natural bond orbital analysis and the molecular tailoring approach method, in order to estimate the relative intramolecular HB energies and explore the substitution effects on HB strength.

Structure, conformational landscape and UV-induced selective rotamerization of 2-aminoacetanilide

The molecular structure, conformational space and UV-induced rotamerization and decarbonylation of 2-aminoacetanilide (2AA) were investigated by matrix isolation infrared spectroscopy and quantum chemical electronic structure calculations. The compound was isolated in a cryogenic (15 K) argon matrix, and four conformers were spectroscopically identified: two trans-amide (OCNH dihedral ∼ 180°) conformers (trans-I and trans-II), and two higher-energy cis-amide (OCNH dihedral ∼ 0°) conformers (cis-I and cis-II). These four conformers were found to be present in the as-deposited matrix of 2AA in a population ratio matching well that predicted by the calculations for the gas-phase equilibrium at the temperature of the vapor of the compound before deposition (60 °C = 333.15 K). The calculations predicted also the existence of an additional higher-energy trans-amide conformer (trans-III) which was not observed experimentally. In situ broadband UV irradiation (λ > 235 nm) of the matrix-isolated compound was found to induce selective conversion of conformer trans-I into conformer trans-II, in a few minutes (55 % after 2 min. of irradiation; 70 % after 10 min.), while the populations of the cis-amide conformers did not change. Prolonged UV irradiation was found to result in decomposition of the compound, leading to generation of carbon monoxide and 2-amine-N-methylaniline. The infrared spectra of the experimentally relevant conformers of 2AA were interpreted and assigned with help of normal coordinate analysis, and the different behavior of the conformers of 2AA upon UV-irradiation of the argon matrix was explained with help of time-dependent DFT (TD-DFT) calculations.

Studies on the Selectivity Mechanism of Wild-Type E. coli Thioesterase ‘TesA and Its Mutants for Medium- and Long-Chain Acyl Substrates

E. coli thioesterase ‘TesA is an important enzyme in fatty acid production. Medium-chain fatty acids (MCFAs, C6-C10) are of great interest due to their similar physicochemical properties to petroleum-based oleo-chemicals. It has been shown that wild-type ‘TesA had better selectivity for long-chain acyl substrates (≥C16), while the two mutants ‘TesAE142D/Y145G and ‘TesAM141L/E142D/Y145G had better selectivity for medium-chain acyl substrates. However, it is difficult to obtain the selectivity mechanism of substrates for proteins by traditional experimental methods. In this study, in order to obtain more MCFAs, we analyzed the binding mode of proteins (‘TesA, ‘TesAE142D/Y145G and ‘TesAM141L/E142D/Y145G) and substrates (C16/C8-N-acetylcysteamine analogs, C16/C8-SNAC), the key residues and catalytic mechanisms through molecular docking, molecular dynamics simulations and the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA). The results showed that several main residues related to catalysis, including Ser10, Asn73 and His157, had a strong hydrogen bond interaction with the substrates. The mutant region (Met141-Tyr146) and loop107–113 were mainly dominated by Van der Waals contributions to the substrates. For C16-SNAC, except for ‘TesAM141L/E142D/Y145G with large conformational changes, there were strong interactions at both head and tail ends that distorted the substrate into a more favorable high-energy conformation for the catalytic reaction. For C8-SNAC, the head and tail found it difficult to bind to the enzyme at the same time due to insufficient chain length, which made the substrate binding sites more variable, so ‘TesAM141L/E142D/Y145G with better binding sites had the strongest activity, and ‘TesA had the weakest activity, conversely. In short, the matching substrate chain and binding pocket length are the key factors affecting selectivity. This will be helpful for the further improvement of thioesterases.

Lowering of Reaction Rates by Energetically Favorable Hydrogen Bonding in the Transition State. Degradation of Biofuel Ketohydroperoxides by OH.

Ketohydroperoxides (KHPs) are oxygenates with carbonyl and hydroperoxy functional groups, and they are generated under combustion and atmospheric conditions. Their fate is crucial for secondary organic aerosol formation in the troposphere and for the ignition processes of biofuels in advanced combustion engines. We investigated the thermodynamics and kinetics of nine hydrogen abstraction reactions from four ether KHPs by OH. We find that the rate constants are strongly affected by entropic effects whose estimation requires a consideration of higher-energy conformers of the transition state. A density functional was selected for these reactions by comparison to coupled cluster calculations, and it was used for calculations by multistructural canonical transition-state theory with multidimensional tunneling over the temperature range of 200-2000 K. We find that the effect of multistructural torsional anharmonicity is very large and quite different for the various ether KHP reactions. A leading cause of the structural dependence is the dominance of entropic factors due to the lack of hydrogen bonding in some of the higher-energy conformers of the transition states. Four of the reactions involve abstraction from the α-carbon (the carbon vicinal to the hydroperoxide group); they exhibit nonmonotonic temperature dependence with complex fuel-specific dependence. The rate constants for abstraction from a non-α-carbon of a given KHP can be faster than the ones for abstraction from an α-carbon; in two cases, this is due to entropy, and in one case, the non-α-carbon abstraction has a lower energy barrier. Tunneling and recrossing effects are also found to be important.

Cost-Effective Implementation of Multiconformer Transition State Theory for Alkoxy Radical Unimolecular Reactions.

Alkoxy radicals are important intermediates in the gas-phase oxidation of volatile organic compounds (VOCs) determining the nature of the first-generation products. An accurate description of their chemistry under atmospheric conditions is essential for understanding the atmospheric oxidation of VOCs. Unfortunately, experimental measurements of the rate coefficients of unimolecular alkoxy radical reactions are scarce, especially for larger systems. As has previously been done for peroxy radical hydrogen shift reactions, we present a cost-effective approach to the practical implementation of multiconformer transition state theory (MC-TST) for alkoxy radical unimolecular (H-shift and decomposition) reactions. Specifically, we test the optimal approach for the conformational sampling as well as the best value for a cutoff of high-energy conformers. In order to obtain accurate rate coefficients at a reduced computational cost, an energy cutoff is employed to reduce the required number of high-level calculations. The rate coefficients obtained with the developed theoretical approach are compared to available experimental rate coefficients for both 1,5 H-shifts and decomposition reactions. For all but one of the reactions tested, the calculated MC-TST rate coefficients agree with experimental results to within a factor of 7. The discrepancy for the final reaction is about a factor of 15, but part of the discrepancy is caused by pressure effects, which are not included in MC-TST. Thus, for the fastest alkoxy reactions, deviation from the high-pressure limit even at 1 bar should be considered.

Markov state models elucidate the stability of DNA influenced by the chiral 5S-Tg base.

The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol−1 through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol−1 of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6HO6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson–Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 KBT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.

Low-Temperature Infrared Spectra and Ultraviolet-Induced Rotamerization of 5-Chlorosalicylaldehyde.

5-Chlorosalicylaldehyde (abbreviated as 5CSA) is an important chemical used in the synthesis of fragrances, dyes, and pharmaceuticals. In this investigation, 5CSA isolated in solid N2, at 10 K, and in its neat amorphous and crystalline phases, at 50 and 190 K, respectively, were investigated by infrared spectroscopy and DFT(B3LYP)/6-311++G(d,p) calculations. The systematic theoretical analysis of the 5CSA conformational landscape showed that the compound exhibits four different conformers, which were structurally characterized in detail. In the as-deposited low-temperature matrices of 5CSA, only the most stable conformer, the intramolecularly hydrogen-bonded form I, was found. The same was observed in the case of the investigated low-temperature amorphous and crystalline phases of 5CSA. Conformer I was successfully converted into a higher-energy conformer(II), where both aldehyde and hydroxyl groups are rotated by 180° relative to their position in the initial conformer, through narrowband ultraviolet (UV) (λ = 308 nm) in situ irradiation of the as-deposited N2 matrix of 5CSA. The infrared spectra of both matrix-isolated conformers, as well as those of the neat amorphous and crystalline phases of 5CSA, were assigned and interpreted in comparative terms, allowing us to elucidate structurally and vibrationally relevant effects of the main intra- and intermolecular interactions operating in the different studied phases. Very interestingly, the observed UV-induced I → II rotamerization was found to take place in an exclusive basis, with no other photochemical processes being observed to occur upon UV irradiation, under the experimental conditions used in the present investigation.

The Local Topological Free Energy of the SARS-CoV-2 Spike Protein.

The novel coronavirus SARS-CoV-2 infects human cells using a mechanism that involves binding and structural rearrangement of its Spike protein. Understanding protein rearrangement and identifying specific amino acids where mutations affect protein rearrangement has attracted much attention for drug development. In this manuscript, we use a mathematical method to characterize the local topology/geometry of the SARS-CoV-2 Spike protein backbone. Our results show that local conformational changes in the FP, HR1, and CH domains are associated with global conformational changes in the RBD domain. The SARS-CoV-2 variants analyzed in this manuscript (alpha, beta, gamma, delta Mink, G614, N501) show differences in the local conformations of the FP, HR1, and CH domains as well. Finally, most mutations of concern are either in or in the vicinity of high local topological free energy conformations, suggesting that high local topological free energy conformations could be targets for mutations with significant impact of protein function. Namely, the residues 484, 570, 614, 796, and 969, which are present in variants of concern and are targeted as important in protein function, are predicted as such from our model.

Diethyl ether oxidation: Revisiting the kinetics of the intramolecular hydrogen abstraction reactions of its primary and secondary peroxy radicals

Diethyl ether is an attractive alternative to diesel fuel for its high autoignition propensity, high energy density, and potential to reduce pollutant emissions. Low-temperature oxidation kinetics play an important role governing fuel autoignition in engines. Peroxy radical intramolecular hydrogen atom abstraction reactions are key isomerization steps in the chain branching reaction sequence at low temperatures. In this study, the kinetics of intramolecular hydrogen atom abstraction reactions of the primary and secondary peroxy radicals were revisited using the multistructural torsional variational transition state theory with small curvature tunneling corrections. High-pressure limit rate constants within a broad temperature range were obtained, which differ from those previously reported for the investigated reactive systems due to the differences in the estimated barrier heights, the role played by multistructural torsional anharmonicity, and the effects of tunneling. We also observed that some higher energy conformers of the saddle point species, some of which with diastereomers, play a surprisingly large role. Similarly, tunneling is pronounced at low temperature combustion conditions, claiming for robust methods to estimate this effect. Different kinetic models for diethyl ether were used to test our calculated rate constants and NASA polynomials for both diethyl ether peroxy radicals, and the results showed an overall slightly better performance in the prediction of ignition delay times at high pressures.

Matrix Isolation Study of Fumaric and Maleic Acids in Solid Nitrogen.

Fumaric and maleic acids ((E)- and (Z)-HOOC-CH═CH-COOH, FA and MA) were studied experimentally by infrared spectroscopy in nitrogen matrixes and theoretically by quantum chemical calculations. The calculations, carried out at the DFT(B3LYP) and MP2 levels of theory, predicted the existence of at least 5 conformers of maleic acid and 10 conformers of fumaric acid. After the deposition of the matrixes, two conformers of maleic acid (I and II) and three conformers (I-III) of fumaric acid were observed and characterized vibrationally. Selective narrowband near-infrared (NIR) excitation of the first OH stretching overtones of the different conformers of maleic and fumaric acids initially present in the matrixes allowed the generation of higher-energy forms, never before observed experimentally. In the case of maleic acid, conformers I (a cis-trans form, where cis and trans designate the conformation of the carboxylic groups of the molecule) and II (cis-cis) were found to generate the novel conformers VI (trans-trans) and VII (cis-trans), respectively. The conversion of conformer II into the most stable conformer I was also observed. For fumaric acid, the cis-cis conformers I-III were found to give rise to the new cis-trans conformers IV-VII, respectively. The tunneling decay of the new conformers produced upon NIR excitation of the lowest-energy conformers initially trapped in the matrixes was observed, and their lifetimes in solid N2 were determined. The increased stability of all of the observed high-energy conformers of the studied acids in the N2 matrix, compared to the argon matrix, where they could not be observed experimentally, demonstrates the stabilizing effect of the interaction between the OH groups of the acids with the matrix N2 molecules, in line with previous observations for other carboxylic acids. In addition, the photochemistry of matrix-isolated maleic and fumaric acids upon broad-band UV irradiation (λ > 235 nm) was also investigated. UV-induced isomerization of both acids around the C═C double bond was observed, together with their decarboxylation to acrylic acid.

Structural basis for the synthesis of the core 1 structure by C1GalT1

C1GalT1 is an essential inverting glycosyltransferase responsible for synthesizing the core 1 structure, a common precursor for mucin-type O-glycans found in many glycoproteins. To date, the structure of C1GalT1 and the details of substrate recognition and catalysis remain unknown. Through biophysical and cellular studies, including X-ray crystallography of C1GalT1 complexed to a glycopeptide, we report that C1GalT1 is an obligate GT-A fold dimer that follows a SN2 mechanism. The binding of the glycopeptides to the enzyme is mainly driven by the GalNAc moiety while the peptide sequence provides optimal kinetic and binding parameters. Interestingly, to achieve glycosylation, C1GalT1 recognizes a high-energy conformation of the α-GalNAc-Thr linkage, negligibly populated in solution. By imposing this 3D-arrangement on that fragment, characteristic of α-GalNAc-Ser peptides, C1GalT1 ensures broad glycosylation of both acceptor substrates. These findings illustrate a structural and mechanistic blueprint to explain glycosylation of multiple acceptor substrates, extending the repertoire of mechanisms adopted by glycosyltransferases.

Proteins' Evolution upon Point Mutations.

As the reader must be already aware, state-of-the-art protein folding prediction methods have reached a smashing success in their goal of accurately determining the three-dimensional structures of proteins. Yet, a solution to simple problems such as the effects of protein point mutations on their (i) native conformation; (ii) marginal stability; (iii) ensemble of high-energy nativelike conformations; and (iv) metamorphism propensity and, hence, their evolvability, remains as an unsolved problem. As a plausible solution to the latter, some properties of the amide hydrogen-deuterium exchange, a highly sensitive probe of the structure, stability, and folding of proteins, are assessed from a new perspective. The preliminary results indicate that the protein marginal stability change upon point mutations provides the necessary and sufficient information to estimate, through a Boltzmann factor, the evolution of the amide hydrogen exchange protection factors and, consequently, that of the ensemble of folded conformations coexisting with the native state. This work contributes to our general understanding of the effects of point mutations on proteins and may spur significant progress in our efforts to develop methods to determine the appearance of new folds and functions accurately.

Geometry meta-optimization.

Recent work has demonstrated the promise of using machine-learned surrogates, in particular, Gaussian process (GP) surrogates, in reducing the number of electronic structure calculations (ESCs) needed to perform surrogate model based (SMB) geometry optimization. In this paper, we study geometry meta-optimization with GP surrogates where a SMB optimizer additionally learns from its past "experience" performing geometry optimization. To validate this idea, we start with the simplest setting where a geometry meta-optimizer learns from previous optimizations of the same molecule with different initial-guess geometries. We give empirical evidence that geometry meta-optimization with GP surrogates is effective and requires less tuning compared to SMB optimization with GP surrogates on the ANI-1 dataset of off-equilibrium initial structures of small organic molecules. Unlike SMB optimization where a surrogate should be immediately useful for optimizing a given geometry, a surrogate in geometry meta-optimization has more flexibility because it can distribute its ESC savings across a set of geometries. Indeed, we find that GP surrogates that preserve rotational invariance provide increased marginal ESC savings across geometries. As a more stringent test, we also apply geometry meta-optimization to conformational search on a hand-constructed dataset of hydrocarbons and alcohols. We observe that while SMB optimization and geometry meta-optimization do save on ESCs, they also tend to miss higher energy conformers compared to standard geometry optimization. We believe that further research into characterizing the divergence between GP surrogates and potential energy surfaces is critical not only for advancing geometry meta-optimization but also for exploring the potential of machine-learned surrogates in geometry optimization in general.

Rapid Allylic 1,6 H-Atom Transfer in an Unsaturated Criegee Intermediate.

A novel allylic 1,6 hydrogen-atom-transfer mechanism is established through infrared activation of the 2-butenal oxide Criegee intermediate, resulting in very rapid unimolecular decay to hydroxyl (OH) radical products. A new precursor, Z/E-1,3-diiodobut-1-ene, is synthesized and photolyzed in the presence of oxygen to generate a new four-carbon Criegee intermediate with extended conjugation across the vinyl and carbonyl oxide groups that facilitates rapid allylic 1,6 H-atom transfer. A low-energy reaction pathway involving isomerization of 2-butenal oxide from a lower-energy (tZZ) conformer to a higher-energy (cZZ) conformer followed by 1,6 hydrogen transfer via a seven-membered ring transition state is predicted theoretically and shown experimentally to yield OH products. The low-lying (tZZ) conformer of 2-butenal oxide is identified based on computed anharmonic frequencies and intensities of its conformers. Experimental IR action spectra recorded in the fundamental CH stretch region with OH product detection by UV laser-induced fluorescence reveal a distinctive IR transition of the low-lying (tZZ) conformer at 2996 cm-1 that results in rapid unimolecular decay to OH products. Statistical RRKM calculations involving a combination of conformational isomerization and unimolecular decay via 1,6 H-transfer yield an effective decay rate keff(E) on the order of 108 s-1 at ca. 3000 cm-1 in good accord with the experiment. Unimolecular decay proceeds with significant enhancement due to quantum mechanical tunneling. A rapid thermal decay rate of ca. 106 s-1 is predicted by master-equation modeling of 2-butenal oxide at 298 K, 1 bar. This novel unimolecular decay pathway is expected to increase the nonphotolytic production of OH radicals upon alkene ozonolysis in the troposphere.

Bound nucleotide can control the dynamic architecture of monomeric actin

Polymerization of actin into cytoskeletal filaments is coupled to its bound adenine nucleotides. The mechanism by which nucleotide modulates actin functions has not been evident from analyses of ATP- and ADP-bound crystal structures of the actin monomer. We report that NMR chemical shift differences between the two forms are globally distributed. Furthermore, microsecond–millisecond motions are spread throughout the molecule in the ATP form, but largely confined to subdomains 1 and 2, and the nucleotide binding site in the ADP form. Through these motions, the ATP- and ADP-bound forms sample different high-energy conformations. A deafness-causing, fast-nucleating actin mutant populates the high-energy conformer of ATP-actin more than the wild-type protein, suggesting that this conformer may be on the pathway to nucleation. Together, the data suggest a model in which differential sampling of a nucleation-compatible form of the actin monomer may contribute to control of actin filament dynamics by nucleotide.