Search

Identification of molecules and elucidation of their chemical structure are ubiquitous problems in chemistry. Mass spectrometry (MS) can be used due to its sensitivity and versatility. For detection to occur, analytes must be ionized and transferred to the gas phase. Soft ionization processes such as electrospray ionization are popular; however, resulting microsolvated phases can alter the chemistry of analytes and therefore detection and identification. To understand these processes, we use computational methods to probe the ionization propensity of serine in the gas phase, aqueous microsolvated clusters, and aqueous solution. We show that the tautomeric form of serine is altered by the presence of water, as five water molecules can stabilize the zwitterionic tautomer. Inclusion of cosolutes such as ions can stabilize the zwitterion with as few as one or two water molecules present. We demonstrate that ionization propensity, as measured by gas phase bacisity, can increase by over 100 kJ/mol when placed in a small water-serine cluster, showing the sensitivity of the chemistry of microsolvated analytes. Finally, detailed analysis reveals that small droplets (less than seven water molecules) are extremely sensitive to addition of further water molecules. Beyond this limit, structural and electronic properties change little with droplet size.

The binding of group II metal cations such as Ca2+ and Mg2+ has been largely categorized as electrostatic or ionic using carboxylate symmetric and asymmetric stretching frequency assignments that have been historically used with little regard for the solvation environment of aqueous solutions. However, given the importance of these cations and their binding mechanisms related to biological function and in revealing surface enrichment factors for ocean to marine aerosol transfer, it is imperative that a deeper understanding be sought to include hydration effects. Here, infrared reflection-absorption and Raman spectra for surface and solution phase carboxylate binding information, respectively, are compared against bare (unbound) carboxylate and bidentate Zn2+:carboxylate spectral signatures. Spectral non-coincidence effect analysis, temperature studies, and spectral and potential of mean force calculations result in a concise interpretation of binding motifs that include the role of mediating water molecules, that is, contact and solvent-shared ion pairs. Calcium directly binds to the carboxylate group in contact ion pairs where magnesium rarely does. Moreover, we reveal the dominance of the solvent-shared ion pair of magnesium with carboxylate at the air-water interface and in solution.

The plant acyl-acyl carrier protein (ACP) desaturases are a family of soluble enzymes that convert saturated fatty acyl-ACPs into their cis-monounsaturated equivalents in an oxygen-dependent reaction. These enzymes play a key role in biosynthesis of monounsaturated fatty acids in plants. ACPs are central proteins in fatty acid biosynthesis that deliver acyl chains to desaturases. They have been reported to show a varying degree of local dynamics and structural variability depending on the acyl chain size. It has been suggested that substrate-specific changes in ACP structure and dynamics have a crucial impact on the desaturase enzymatic activity. Using molecular dynamics simulations, we investigated the intrinsic solution structure and dynamics of ACP from spinach with four different acyl chains: capric (C10), myristic (C14), palmitic (C16), and stearic (C18) acids. We found that the fatty acids can adopt two distinct structural binding motifs, which feature different binding free energies and influence the ACP dynamics in a different manner. Docking simulations of ACP to castor Δ 9 -desaturase and ivy Δ 4 -desaturase suggest that ACP desaturase interactions could lead to a preferential selection between the motifs.

Fully synthetic peptoid membranes are known to mimic important features of biological membranes, with several advantages over other biomimetic membranes. A fundamental understanding of how the individual peptoid amphiphiles assemble in solution to form the bilayer membrane is key to unlocking their versatility for application in a broad range of processes. In this study, in situ X-ray scattering and molecular dynamics simulations are used to understand the early stages of assembly of three different peptoids that exhibit distinctly different crystallization kinetics. The in situ measurements reveal that the peptoids aggregate first into a nascent phase that is less crystalline than the assembled peptoid membrane. Anisotropic aromatic interactions are determined to be the dominant driving force in the early stages of membrane formation. These results provide key insights into how the peptoid assembly may be manipulated during the early stages of assembly and nucleation and growth.

Free-standing two-dimensional (2D) organic nanomaterials are highly attractive for biological applications because of their unique structural properties and high biocompatibility. Herein, we design...

Rate theories have found great utility across the chemical sciences by providing a physically transparent way to analyze dynamical processes. Here we demonstrate the benefits of using transition state theory and Marcus theory to study the rate of proton transfer in HCl solutions. By using long

We performed experimental and periodic density functional theory (DFT) IR spectroscopy to investigate the adsorption of acyclic olefins over both acidic and nonacidic zeolites. Two conjugated polye...

Plant plastidial acyl–acyl carrier protein (ACP) desaturases are a soluble class of diiron-containing enzymes that are distinct from the diiron-containing integral membrane desaturases found in plants and other organisms. The archetype of this class is the stearoyl-ACP desaturase which converts stearoyl-ACP into oleoyl (18:1Δ9cis)-ACP. Several variants expressing distinct regioselectivity have been described including a Δ6-16:0-ACP desaturase from black-eyed Susan vine (Thunbergia alata). We solved a crystal structure of the T. alata desaturase at 2.05 Å resolution. Using molecular dynamics (MD) simulations, we identified a low-energy complex between 16:0-ACP and the desaturase that would position C6 and C7 of the acyl chain adjacent to the diiron active site. The model complex was used to identify mutant variants that could convert the T. alata Δ6 desaturase to Δ9 regioselectivity. Additional modeling between ACP and the mutant variants confirmed the predicted regioselectivity. To validate the in-silico predictions, we synthesized two variants of the T. alata desaturase and analyzed their reaction products using gas chromatography-coupled mass spectrometry. Assay results confirmed that mutants designed to convert T. alata Δ6 to Δ9 selectivity exhibited the predicted changes. In complementary experiments, variants of the castor desaturase designed to convert Δ9 to Δ6 selectivity lost some of their Δ9 desaturation ability and gained the ability to desaturate at the Δ6 position. The computational workflow for revealing the mechanistic understanding of regioselectivity presented herein lays a foundation for designing acyl-ACP desaturases with novel selectivities to increase the diversity of monoenes available for bioproduct applications.

Muscovite mica (001) is a widely used model surface for controlling molecular assembly and a common substrate for environmental adsorption processes. The mica (001) surface displays near-trigonal symmetry, but many molecular adsorbates-including water-exhibit unequal probabilities of alignment along its three nominally equivalent lattice directions. Buried hydroxyl groups within the muscovite structure are speculated to be responsible, but direct evidence is lacking. Here, we utilize vibrational sum frequency generation spectroscopy (vSFG) to characterize the orientation and hydrogen-bonding environment of near-surface hydroxyls inside mica. Multiple distinct peaks are detected in the O-H stretch region, which we attribute to Si/Al substitution in the SiO4 tetrahedron and K+ ion adsorption above the hydroxyls based on density functional theory simulations. Our findings demonstrate that vSFG can identify the absolute orientation of -OH groups and, hence, the surface termination at a mica surface, providing a means to investigate how -OH groups influence molecular adsorption and better understand mica stacking-sequences and physical behavior.

Providing a molecular level description of amorphous CaCO\(_3\) nucleation in aqueous solutions necessitates the use of theory and simulation in addition to a direct connection to experiments. Furthermore, the ability to simulate experimental conditions relevant to nucleation requires a systematic reduction in the complexity of a molecular-based model that can be compared to experimental observations and tested for self-consistency. In this review, we further expand on the theoretical details described in (Henzler et al. Sci Adv 4:eaa06283, (2018)), where a molecular-based solution model was constructed based on a solvent-mediated quantum mechanical potential of mean force for Ca\(^{2+}\) and CO\(_3^{2-}\) ion pairing. The statistical mechanics and simulation approaches for solution-phase nucleation are presented and connected to the various CaCO\(_3\) polymorph chemical potentials and interfacial surface energies. The resulting theory presented herein is general and is shown to accurately describe the initial stages of amorphous calcium carbonate nucleation.

We study the prototypical S

We study the prototypical SN2 reaction Cl$^-$ + CH$_3$Cl $\to$ CH$_3$Cl + Cl$^-$ in water using quantum mechanics / molecular mechanics (QM/MM) computer simulations with transition path sampling and inertial likelihood maximization. We have identified a new solvent coordinate to complement the original atom-exchange coordinate used in the classic analysis by Chandrasekhar, Smith, and Jorgensen [Ref1]. The new solvent coordinate quantifies instantaneous solvent induced polarization relative to the equilibrium average charge density at each point along the reaction pathway. On the basis of likelihood scores and committor distributions, the new solvent coordinate improves upon the description of solvent dynamical effects relative to previously proposed solvent coordinates. However, it does not increase the transmission coefficient or the accuracy of a transition state theory rate calculation., 11 pages, 10 figures

When hydrolyzable cations such as aluminum interact with solid-water interfaces, macroscopic interfacial properties (e.g., surface charge and potential) and interfacial phenomena (e.g., particle adhesion) become tightly linked with the microscopic details of ion adsorption and speciation. We use in situ atomic force microscopy to directly image individual aluminum ions at a mica-water interface and show how adsorbate populations change with pH and aluminum activity. Complementary streaming potential measurements then allow us to build a triple layer model (TLM) that links surface potentials to adsorbate populations, via equilibrium binding constants. Our model predicts that hydrolyzed species dominate the mica-water interface, even when unhydrolyzed species dominate the solution. Ab initio molecular dynamics (AIMD) simulations confirm that aluminum hydrolysis is strongly promoted at the interface. The TLM indicates that hydrolyzed adsorbates are responsible for surface-potential inversions, and we find strong correlations between hydrolyzed adsorbates and particle-adhesion forces, suggesting that these species mediate adhesion by chemical bridging.

The ability to reproduce the experimental structure of water around the sodium and potassium ions is a key test of the quality of interaction potentials due to the central importance of these ions in a wide range of important phenomena. Here, we simulate the Na+ and K+ ions in bulk water using three density functional theory functionals: 1) The generalized gradient approximation (GGA) based dispersion corrected revised Perdew, Burke, and Ernzerhof functional (revPBE-D3) 2) The recently developed strongly constrained and appropriately normed (SCAN) functional 3) The random phase approximation (RPA) functional for potassium. We compare with experimental X-ray diffraction (XRD) and X-ray absorption fine structure (EXAFS) measurements to demonstrate that SCAN accurately reproduces key structural details of the hydra- tion structure around the sodium and potassium cations, whereas revPBE-D3 fails to do so. However, we show that SCAN provides a worse description of pure water in comparison with revPBE-D3. RPA also shows an improvement for K+, but slow convergence prevents rigorous comparison. Finally, we analyse cluster energetics to show SCAN and RPA have smaller fluctuations of the mean error of ion-water cluster binding energies compared with revPBE-D3.

Understanding the nature of ionic hydration at a fundamental level has eluded scientists despite intense interest for nearly a century. In particular, the microscopic origins of the asymmetry of ion solvation thermodynamics with respect to the sign of the ionic charge remains a mystery. Here, we determine the response of accurate quantum mechanical water models to strong nanoscale solvation forces arising from excluded volumes and ionic electrostatic fields. This is compared to the predictions of two important limiting classes of classical models of water with fixed point changes, differing in their treatment of "lone pair" electrons. Using the quantum water model as our standard of accuracy, we find that a single fixed classical treatment of lone pair electrons cannot accurately describe solvation of both apolar and cationic solutes, emphasizing the need for a more flexible description of local electronic effects in solvation processes. However, we explicitly show that all water models studied respond to weak long-ranged electrostatic perturbations in a manner that follows macroscopic dielectric continuum models, as would be expected. We emphasize the importance of these findings in the context of realistic ion models, using density functional theory and empirical models, and discuss the implications of our results for quantitatively accurate reduced descriptions of solvation in dielectric media.