Search

Water and ice are routinely studied with X-rays to reveal their diverse structures and anomalous properties. We employ a hybrid collisional-radiative/molecular-dynamics method to explore how femtosecond X-ray pulses interact with hexagonal ice. We find that ice makes a phase transition into a crystalline plasma where its initial structure is maintained up to tens of femtoseconds. The ultrafast melting process occurs anisotropically, where different geometric configurations of the structure melt on different time scales. The transient state and anisotropic melting of crystals can be captured by X-ray diffraction, which impacts any study of crystalline structures probed by femtosecond X-ray lasers., Competing Interests: The authors declare no competing financial interest., (© 2024 The Authors. Published by American Chemical Society.)

We demonstrate site-specific X-ray induced fragmentation across the sulfur L-edge of protonated cystine, the dimer of the amino acid cysteine. Ion yield NEXAFS were performed in the gas phase using electrospray ionization (ESI) in combination with an ion trap. The interpretation of the sulfur L-edge NEXAFS spectrum is supported by Restricted Open-Shell Configuration Interaction (ROCIS) calculations. The fragmentation pathway of triply charged cystine ions was modeled by Molecular Dynamics (MD) simulations. We have deduced a possible pathway of fragmentation upon excitation and ionization of S 2p electrons. The disulfide bridge breaks for resonant excitation at lower photon energies but remains intact upon higher energy resonant excitation and upon ionization of S 2p. The larger fragments initially formed subsequently break into smaller fragments.

A nanopore is a tool in single-molecule sensing biotechnology that offers label-free identification with high throughput. Nanopores have been successfully applied to sequence DNA and show potential in the study of proteins. Nevertheless, the task remains challenging due to the large variability in size, charges, and folds of proteins. Miniproteins have a small number of residues, limited secondary structure, and stable tertiary structure, which can offer a systematic way to reduce complexity. In this computational work, we theoretically evaluated sensing two miniproteins found in the human body using a silicon nitride nanopore. We employed molecular dynamics methods to compute occupied-pore ionic current magnitudes and electronic structure calculations to obtain interaction strengths between pore wall and miniprotein. From the interaction strength, we derived dwell times using a mix of combinatorics and numerical solutions. This latter approach circumvents typical computational demands needed to simulate translocation events using molecular dynamics. We focused on two miniproteins potentially difficult to distinguish owing to their isotropic geometry, similar number of residues, and overall comparable structure. We found that the occupied-pore current magnitudes not to vary significantly, but their dwell times differ by 1 order of magnitude. Together, these results suggest a successful identification protocol for similar miniproteins., Competing Interests: The authors declare no competing financial interest., (© 2021 The Authors. Published by American Chemical Society.)

Atmospheric aerosols contain a variety of compounds, among them free amino acids and salt ions. The pH of the aerosol droplets depends on their origin and environment. Consequently, compounds like free amino acids found in the droplets will be at different charge states, since these states to a great extent depend on the surrounding pH condition. In droplets of marine origin, amino acids are believed to drive salt ions to the water surface and a pH-dependent amino acid surface propensity will, therefore, indirectly affect many processes in atmospheric chemistry and physics such as for instance cloud condensation. To understand the surface propensity of glycine, valine, and phenylalanine at acidic, neutral, and basic pH, we used molecular dynamics (MD) simulations to investigate them at three different charge states in water. Their respective surface propensities were obtained by the means of a potential of mean force (PMF) in an umbrella sampling approach. Glycine was found to have no preference for the surface, while both valine and phenylalanine showed high propensities. Among the charge states of the surface-enriched ones, the cation, representing the amino acids at low pH, was found to have the highest affinity. Free energy decomposition revealed that the driving forces depend strongly on the nature of the amino acid and its charge state. In phenylalanine, the main factor was found to be a substantial entropy gain, likely related to the side chain, whereas in valine, hydrogen bonding to the functional groups leads to favorable energies and, in turn, affects the surface propensity. A significant gain in water-water enthalpy was seen for both valine and phenylalanine.

One of the challenges facing single particle imaging with ultrafast X-ray pulses is the structural heterogeneity of the sample to be imaged. For the method to succeed with weakly scattering samples, the diffracted images from a large number of individual proteins need to be averaged. The more the individual proteins differ in structure, the lower the achievable resolution in the final reconstructed image. We use molecular dynamics to simulate two globular proteins in vacuum, fully desolvated as well as with two different solvation layers, at various temperatures. We calculate the diffraction patterns based on the simulations and evaluate the noise in the averaged patterns arising from the structural differences and the surrounding water. Our simulations show that the presence of a minimal water coverage with an average 3 Å thickness will stabilize the protein, reducing the noise associated with structural heterogeneity, whereas additional water will generate more background noise.

Auger spectroscopy has previously been used to study changes in the hydrogen bond network in liquid water, but to the best of our knowledge it has not been used to track such changes as a function of temperature. We show Auger spectroscopy to reflect the weakening of the hydrogen bond network upon heating. This shows that the radiation response of water, i.e., the relative propensity of the different processes occurring after radiation exposure, including femtosecond proton dynamics, depends on the temperature of the system. This proof-of-principle study further demonstrates the suitability of the technique to help elucidate information on the intermolecular structure of liquids such as water, opening the door to further temperature-dependent studies.

Gas-phase electric properties of molecules can be computed routinely using wave function methods or density functional theory (DFT). However, these methods remain computationally expensive for high-throughput screening of the vast chemical space of virtual compounds. Therefore, empirical force fields are a more practical choice in many cases, particularly since force field methods allow one to routinely predict the physicochemical properties in the condensed phases. This work presents Drude polarizable models, to increase the physical realism in empirical force fields, where the core particle is treated as a point charge and the Drude particle is treated either as a 1 s-Gaussian or a ns-Slater ( n = 1, 2, 3) charge density. Systematic parametrization to large high-quality quantum chemistry data obtained from the open access Alexandria Library ( https://doi.org/10.5281/zenodo.1004711 ) ensures the transferability of these parameters. The dipole moments and isotropic polarizabilities of the isolated molecules predicted by the proposed Drude models are in agreement with experiment with accuracy similar to DFT calculations at the B3LYP/aug-cc-pVTZ level of theory. The results show that the inclusion of explicit polarization into the models reduces the root-mean-square deviation with respect to DFT calculations of the predicted dipole moments of 152 dimers and clusters by more than 50%. Finally, we show that the accuracy of the electrostatic interaction energy of the water dimers can be improved systematically by the introduction of polarizable smeared charges as a model for charge penetration.

Single-particle imaging using X-ray free-electron lasers is an emerging technique that could provide high-resolution structures of macromolecules in the gas phase. One of the largest difficulties in realizing this goal is the unknown orientation of the individual sample molecules at the time of exposure. Preorientation of the molecules has been identified as a possible solution to this problem. Using molecular dynamics simulations, we identify a range of electric field strengths where proteins become oriented without losing their structure. For a number of experimentally relevant cases we show that structure determination is possible only when orientation information is included in the orientation-recovery process. We conclude that nondestructive field orientation of intact proteins is feasible and that it enables a range of new structural investigations with single-particle imaging.

The chemical composition of small organic molecules is often very similar to amino acid side chains or the bases in nucleic acids, and hence there is no a priori reason why a molecular mechanics force field could not describe both organic liquids and biomolecules with a single parameter set. Here, we devise a benchmark for force fields in order to test the ability of existing force fields to reproduce some key properties of organic liquids, namely, the density, enthalpy of vaporization, the surface tension, the heat capacity at constant volume and pressure, the isothermal compressibility, the volumetric expansion coefficient, and the static dielectric constant. Well over 1200 experimental measurements were used for comparison to the simulations of 146 organic liquids. Novel polynomial interpolations of the dielectric constant (32 molecules), heat capacity at constant pressure (three molecules), and the isothermal compressibility (53 molecules) as a function of the temperature have been made, based on experimental data, in order to be able to compare simulation results to them. To compute the heat capacities, we applied the two phase thermodynamics method (Lin et al. J. Chem. Phys.2003, 119, 11792), which allows one to compute thermodynamic properties on the basis of the density of states as derived from the velocity autocorrelation function. The method is implemented in a new utility within the GROMACS molecular simulation package, named g_dos, and a detailed exposé of the underlying equations is presented. The purpose of this work is to establish the state of the art of two popular force fields, OPLS/AA (all-atom optimized potential for liquid simulation) and GAFF (generalized Amber force field), to find common bottlenecks, i.e., particularly difficult molecules, and to serve as a reference point for future force field development. To make for a fair playing field, all molecules were evaluated with the same parameter settings, such as thermostats and barostats, treatment of electrostatic interactions, and system size (1000 molecules). The densities and enthalpy of vaporization from an independent data set based on simulations using the CHARMM General Force Field (CGenFF) presented by Vanommeslaeghe et al. (J. Comput. Chem.2010, 31, 671) are included for comparison. We find that, overall, the OPLS/AA force field performs somewhat better than GAFF, but there are significant issues with reproduction of the surface tension and dielectric constants for both force fields.

Structural studies of biological macromolecules are severely limited by radiation damage. Traditional crystallography curbs the effects of damage by spreading damage over many copies of the molecule of interest in the crystal. X-ray lasers offer an additional opportunity for limiting damage by out-running damage processes with ultrashort and very intense X-ray pulses. Such pulses may allow the imaging of single molecules, clusters, or nanoparticles. Coherent flash imaging will also open up new avenues for structural studies on nano- and microcrystalline substances. This paper addresses the theoretical potentials and limitations of nanocrystallography with extremely intense coherent X-ray pulses. We use urea nanocrystals as a model for generic biological substances and simulate the primary and secondary ionization dynamics in the crystalline sample. The results establish conditions for ultrafast single-shot nanocrystallography diffraction experiments as a function of X-ray fluence, pulse duration, and the size of nanocrystals. Nanocrystallography using ultrafast X-ray pulses has the potential to open up a new route in protein crystallography to solve atomic structures of many systems that remain inaccessible using conventional X-ray sources.