Intracellular Aß42 Aggregation Leads to Cellular Thermogenesis.

The aggregation of Aß42 is a hallmark of Alzheimer's disease. It is still not known what the biochemical changes are inside a cell which will eventually lead to Aß42 aggregation. Thermogenesis has been associated with cellular stress, the latter of which may promote aggregation. We perform intracellular thermometry measurements using fluorescent polymeric thermometers to show that Aß42 aggregation in live cells leads to an increase in cell-averaged temperatures. This rise in temperature is mitigated upon treatment with an aggregation inhibitor of Aß42 and is independent of mitochondrial damage that can otherwise lead to thermogenesis. With this, we present a diagnostic assay which could be used to screen small-molecule inhibitors to amyloid proteins in physiologically relevant settings. To interpret our experimental observations and motivate the development of future models, we perform classical molecular dynamics of model Aß peptides to examine the factors that hinder thermal dissipation. We observe that this is controlled by the presence of ions in its surrounding environment, the morphology of the amyloid peptides, and the extent of its hydrogen-bonding interactions with water. We show that aggregation and heat retention by Aß peptides are favored under intracellular-mimicking ionic conditions, which could potentially promote thermogenesis. The latter will, in turn, trigger further nucleation events that accelerate disease progression.

Emergence of Electric Fields at the Water-C12E6 Surfactant Interface.

We study the properties of the interface of water and the surfactant hexaethylene glycol monododecyl ether (C12E6) with a combination of heterodyne-detected vibrational sum frequency generation (HD-VSFG), Kelvin-probe measurements, and molecular dynamics (MD) simulations. We observe that the addition of the hydrogen-bonding surfactant C12E6, close to the critical micelle concentration (CMC), induces a drastic enhancement in the hydrogen bond strength of the water molecules close to the interface, as well as a flip in their net orientation. The mutual orientation of the water and C12E6 molecules leads to the emergence of a broad (∼3 nm) interface with a large electric field of ∼1 V/nm, as evidenced by the Kelvin-probe measurements and MD simulations. Our findings may open the door for the design of novel electric-field-tuned catalytic and light-harvesting systems anchored at the water-surfactant-air interface.

Coupled Local-Mode Approach for the Calculation of Vibrational Spectra: Application to Protonated Water Clusters.

Spectroscopic studies of protonated water clusters (PWCs) have yielded enormous insights into the fundamental nature of the hydrated proton. Here, we introduce a new coupled local-mode (CLM) approach to calculate PWC OH stretch vibrational spectra. The CLM method combines a sampling of representative configurations from density functional theory (DFT)-based ab initio molecular dynamics (AIMD) simulations with DFT calculations of local-mode vibrational frequencies and couplings. Calculations of inhomogeneous OH stretch vibrational spectra for H+(H2O)4 and H+(H2O)21 agree well with experiment and higher-level calculations, and decompositions of the calculated spectra in terms of the coupled modes aids in the interpretation of the spectra. This observation is consistent with the idea that capturing anharmonicity and coupling is as important to accuracy as the underlying level of electronic structure theory. The CLM calculations can easily discern the configuration that dominates the experimental measurement for H+(H2O)5, which can adopt several low-energy conformations.

Understanding Förster Energy Transfer through the Lens of Molecular Dynamics.

A multiscale approach to the dynamics of resonant energy transfer (RET) is presented, combining DFT and TD-DFT results on the energy donor (D) and acceptor (A) moieties with an extensive equilibrium and non-equilibrium molecular dynamics (MD) analysis of a bound D-A pair in solution to build a coarse-grained kinetic model. We demonstrate that a thorough MD study is needed to properly address RET: the enormous configuration space visited by the system cannot be reliably sampled accounting only for a few representative configurations. Moreover, the conformational motion of the RET pair, occurring in a similar time scale as the RET process itself, leads to a sizable increase of the overall process efficiency.

Simulation of Liquids with the Tight-Binding Density-Functional Approach and Improved Atomic Charges.

Theoretical description of liquids, especially liquid water, is an ongoing subject with important implications in various domains such as homogeneous catalysis; solvation of molecular, ionic, and biomolecular species; and reactivity. Various formalisms exist to describe liquids, each one displaying its own balance between accuracy and computational cost that defines its range of applications. The present article revisits the ability of the density-functional-based tight-binding (SCC-DFTB) approach to model liquids by focusing on liquid water and liquid benzene under ambient conditions. To do so, we benchmark a recent correction for the SCC-DFTB atomic charges that allows for a drastic improvement of the pair radial distribution functions of liquid water as compared to both experimental data and density-functional theory results performed in the generalized-gradient approximation. We also report the coupling of the deMonNano and i-PI codes to perform path-integral molecular dynamics. This allows us to rationalize the impact of nuclear quantum effects on the SCC-DFTB description of liquid water. This study evidences the rather good ability of SCC-DFTB to describe liquid water and liquid benzene. As the first example of application, we also present results for a benzene molecule solvated in water with the perspectives of further studies devoted to solvent/water interfaces.

A Data Science Approach to Understanding Water Networks Around Biomolecules: The Case of Tri-Alanine in Liquid Water.

Herein, we use recently developed data science algorithms to illustrate the complexity of the water network surrounding the hydrated peptide tri-alanine extracted from molecular dynamics simulations. We estimate the dimensionality of water variables and show that it is sensitive to the underlying secondary structure of the peptide. We show that water wires threading the peptide encode important information on the secondary structure. Interestingly, the free-energy landscape as revealed by the water wires is very rough for α-configurations and rather smooth for ß-configurations. The structured nature of the free-energy landscape is washed out if one uses more standard collective variables such as the number of hydrogen bonds around the peptide. Our results provide fresh insights into the molecular ingredients behind the coupling of protein and solvent degrees of freedom relevant for many biophysical and chemical processes.

Hierarchical lattice models of hydrogen-bond networks in water.

We develop a graph-based model of the hydrogen-bond network in water, with a view toward quantitatively modeling the molecular-level correlational structure of the network. The networks formed are studied by the constructing the model on two infinite-dimensional lattices. Our models are built bottom up, based on microscopic information coming from atomistic simulations, and we show that the predictions of the model are consistent with known results from ab initio simulations of liquid water. We show that simple entropic models can predict the correlations and clustering of local-coordination defects around tetrahedral waters observed in the atomistic simulations. We also find that orientational correlations between bonds are longer ranged than density correlations, determine the directional correlations within closed loops, and show that the patterns of water wires within these structures are also consistent with previous atomistic simulations. Our models show the existence of density and compressibility anomalies, as seen in the real liquid, and the phase diagram of these models is consistent with the singularity-free scenario previously proposed by Sastry and coworkers [Phys. Rev. E 53, 6144 (1996)1063-651X10.1103/PhysRevE.53.6144].

Quantum mechanical effects in zwitterionic amino acids: The case of proline, hydroxyproline, and alanine in water.

In this work, we use ab initio molecular dynamics simulations to elucidate the electronic properties of three hydrated zwitterionic amino acids, namely proline, hydroxyproline, and alanine, the former two forming an important constituent of collagen. In all three systems, we find a substantial amount of charge transfer between the amino acids and surrounding solvent, which, rather surprisingly, also involves the reorganization of electron density near the hydrophobic non-polar groups. Water around proline appears to be slightly more polarized, as reflected by the enhanced water dipole moment in its hydration shell. This observation is also complemented by an examination of the IR spectra of the three systems where there is a subtle red and blue shift in the O-H stretch and bend regions, respectively, for proline. We show that polarizability of these amino acids as revealed by a dipole moment analysis involves a significant enhancement from the solvent and that this also involves non-polar groups. Our results suggest that quantum mechanical effects are likely to be important in understanding the coupling between biomolecules and water in general and in hydrophobic interactions.

Decomposition of the Experimental Raman and Infrared Spectra of Acidic Water into Proton, Special Pair, and Counterion Contributions.

Textbooks describe excess protons in liquid water as hydronium (H3O+) ions, although their true structure remains lively debated. To address this question, we have combined Raman and infrared (IR) multivariate curve resolution spectroscopy with ab initio molecular dynamics and anharmonic vibrational spectroscopic calculations. Our results are used to resolve, for the first time, the vibrational spectra of hydrated protons and counterions and reveal that there is little ion-pairing below 2 M. Moreover, we find that isolated excess protons are strongly IR active and nearly Raman inactive (with vibrational frequencies of ∼1500 ± 500 cm-1), while flanking water OH vibrations are both IR and Raman active (with higher frequencies of ∼2500 ± 500 cm-1). The emerging picture is consistent with Georg Zundel's seminal work, as well as recent ultrafast dynamics studies, leading to the conclusion that protons in liquid water are primarily hydrated by two flanking water molecules, with a broad range of proton hydrogen bond lengths and asymmetries.

The excess proton at the air-water interface: The role of instantaneous liquid interfaces.

The magnitude of the pH of the surface of water continues to be a contentious topic in the physical chemistry of aqueous interfaces. Recent theoretical studies have shown little or no preference for the proton to be at the surface compared to the bulk. Using ab initio molecular dynamics simulations, we revisit the propensity of the excess proton for the air-water interface with a particular focus on the role of instantaneous liquid interfaces. We find a more pronounced presence for the proton to be at the air-water interface. The enhanced water structuring around the proton results in the presence of proton wires that run parallel to the surface as well as a hydrophobic environment made up of under-coordinated topological defect water molecules, both of which create favorable conditions for proton confinement at the surface. The Grotthuss mechanism within the structured water layer involves a mixture of both concerted and closely spaced stepwise proton hops. The proton makes excursions within the first solvation layer either in proximity to or along the instantaneous interface.

Biomolecules at the amorphous silica/water interface: Binding and fluorescence anisotropy of peptides.

We investigate binding of the tripeptides Lys-Trp-Lys (KWK) and Glu-Trp-Glu (EWE) to the amorphous silica surface using atomistic simulations. These peptides were chosen because they were previously utilized in experiments measuring binding affinity and steady-state fluorescence anisotropy from the indole chromophore of the tryptophan residue. Our simulations were performed using silica with surface change density of -0.8 elementary charges per square nanometer, which is expected at neutral pH. Even though positive charged KWK binds more strongly to the negatively charged silica surface, EWE also forms bound complexes with the surface that are stable for at least 15ns of simulation, in agreement with the experiments which revealed evidence for binding of both KWK and EWE to silica. Binding mechanisms include a wide variety of electrostatic interactions, as well as hydrophobic interactions between the indole group and hydrophobic areas of the heterogeneous silica surface. The long-time limit of the fluorescence anisotropy of tryptophan is calculated from the simulations in order to help interpret the recent experiments. We identify several factors which control the magnitude of the fluorescence anisotropy for each binding configuration.

A computational study on how structure influences the optical properties in model crystal structures of amyloid fibrils.

Amyloid fibrils have been shown to have peculiar optical properties since they can exhibit fluorescence in the absence of aromatic residues. In a recent study, we have shown that proton transfer (PT) events along hydrogen bonds (HBs) are coupled to absorption in the near UV range. Here, we gain more insights into the different types of hydrogen bonding interactions that occur in our model systems and the molecular factors that control the susceptibility of the protons to undergo PT and how this couples to the optical properties. In the case of the strong N-C termini interactions, a nearby methionine residue stabilizes the non-zwitterionic NH2-COOH pair, while zwitterionic NH3+-COO- is stabilized by the proximity of nearby crystallographic water molecules. Proton motion along the hydrogen bonds in the fibril is intimately coupled to the compression of the heavier atoms, similar to what is observed in bulk water. Small changes in the compression of the hydrogen bonds in the protein can lead to significant changes in both the ground and excited state potential energy surfaces associated with PT. Finally, we also reinforce the importance of nuclear quantum fluctuations of protons in the HBs of the amyloid proteins.

Probing Defects and Correlations in the Hydrogen-Bond Network of ab Initio Water.

The hydrogen-bond network of water is characterized by the presence of coordination defects relative to the ideal tetrahedral network of ice, whose fluctuations determine the static and time-dependent properties of the liquid. Because of topological constraints, such defects do not come alone but are highly correlated coming in a plethora of different pairs. Here we discuss in detail such correlations in the case of ab initio water models and show that they have interesting similarities to regular and defective solid phases of water. Although defect correlations involve deviations from idealized tetrahedrality, they can still be regarded as weaker hydrogen bonds that retain a high degree of directionality. We also investigate how the structure and population of coordination defects is affected by approximations to the interatomic potential, finding that, in most cases, the qualitative features of the hydrogen-bond network are remarkably robust.

Ab Initio Quality NMR Parameters in Solid-State Materials Using a High-Dimensional Neural-Network Representation.

Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful experimental tools to probe the local atomic order of a wide range of solid-state compounds. However, due to the complexity of the related spectra, in particular for amorphous materials, their interpretation in terms of structural information is often challenging. These difficulties can be overcome by combining molecular dynamics simulations to generate realistic structural models with an ab initio evaluation of the corresponding chemical shift and quadrupolar coupling tensors. However, due to computational constraints, this approach is limited to relatively small system sizes which, for amorphous materials, prevents an adequate statistical sampling of the distribution of the local environments that is required to quantitatively describe the system. In this work, we present an approach to efficiently and accurately predict the NMR parameters of very large systems. This is achieved by using a high-dimensional neural-network representation of NMR parameters that are calculated using an ab initio formalism. To illustrate the potential of this approach, we applied this neural-network NMR (NN-NMR) method on the (17)O and (29)Si quadrupolar coupling and chemical shift parameters of various crystalline silica polymorphs and silica glasses. This approach is, in principal, general and has the potential to be applied to predict the NMR properties of various materials.

Role of Quantum Vibrations on the Structural, Electronic, and Optical Properties of 9-Methylguanine.

In this work, we report theoretical predictions of the UV-absorption spectra of 9-methylguanine using time dependent density functional theory (TDDFT). Molecular dynamics simulations of the hydrated DNA base are peformed using an empirical force field, Born-Oppenheimer ab initio molecular dynamics (AIMD), and finally path-integral AIMD to understand the role of the underlying electronic potential, solvation, and nuclear quantum vibrations on the absorption spectra. It is shown that the conformational distributions, including hydrogen bonding interactions, are perturbed by the inclusion of nuclear quantum effects, leading to significant changes in the total charge and dipole fluctuations of the DNA base. The calculated absorption spectra using the different sampling protocols shows that the inclusion of nuclear quantum effects causes a significant broadening and red shift of the spectra bringing it into closer agreement with experiments.

Structure and Dynamics of the Instantaneous Water/Vapor Interface Revisited by Path-Integral and Ab Initio Molecular Dynamics Simulations.

The structure and dynamics of the water/vapor interface is revisited by means of path-integral and second-generation Car-Parrinello ab initio molecular dynamics simulations in conjunction with an instantaneous surface definition [Willard, A. P.; Chandler, D. J. Phys. Chem. B 2010, 114, 1954]. In agreement with previous studies, we find that one of the OH bonds of the water molecules in the topmost layer is pointing out of the water into the vapor phase, while the orientation of the underlying layer is reversed. Therebetween, an additional water layer is detected, where the molecules are aligned parallel to the instantaneous water surface.

DNA Binding to the Silica Surface.

We investigate the DNA-silica binding mechanism using molecular dynamics simulations. This system is of technological importance, and also of interest to explore how negatively charged DNA can bind to a silica surface, which is also negatively charged at pH values above its isoelectric point near pH 3. We find that the two major binding mechanisms are attractive interactions between DNA phosphate and surface silanol groups and hydrophobic bonding between DNA base and silica hydrophobic region. Umbrella sampling and the weighted histogram analysis method (WHAM) are used to calculate the free energy surface for detachment of DNA from a binding configuration to a location far from the silica surface. Several factors explain why single-stranded DNA (ssDNA) has been observed to be more strongly attracted to silica than double-stranded (dsDNA): (1) ssDNA is more flexible and therefore able to maximize the number of binding interactions. (2) ssDNA has free unpaired bases to form hydrophobic attachment to silica while dsDNA has to break hydrogen bonds with base partners to get free bases. (3) The linear charge density of dsDNA is twice that of ssDNA. We devise a procedure to approximate the atomic forces between biomolecules and amorphous silica to enable large-scale biomolecule-silica simulations as reported here.

Ab initio molecular dynamics study of the mechanism of proton recombination with a weak base.

Despite its fundamental nature, many of the microscopic features of acid­base recombination remain poorly understood. In this work, we use ab initio molecular dynamics simulations to study the recombination of the proton with a weak base, the carbonate ion CO3(2­). Our simulations elucidate the network structure around CO3(2­) that provides a distribution of pathways over which recombination can occur. We observe that the penultimate neutralization step involves a correlated behavior of the transferred protons that is mediated by the water wires decorating the carbonate. These concerted proton transfers are coupled to collective compressions of these water wires. We show further that these processes are dynamically coupled to the reorganization of the water molecules hydrating the CO3(2­) ion. The insights from these simulations help to bridge the structural and dynamical complexity of the microscopic mechanisms with those of phenomenological models invoked by experiments in this field.

The role of quantum effects on structural and electronic fluctuations in neat and charged water.

In this work, we revisit the role of nuclear quantum effects on the structural and electronic properties of the excess proton in bulk liquid water using advanced molecular dynamics techniques. The hydronium ion is known to be a weak acceptor of a hydrogen bond which gives it some hydrophobic character. Quantum effects reduce the degree of this hydrophobicity which facilitates the fluctuations of the protons along the wires compared to the classical proton. Although the Eigen and Zundel species still appear to be dominant motifs, quantum fluctuations result in rather drastic events where both transient autoionization and delocalization over extended proton wires can simultaneously occur. These wild fluctuations also result in a significant change of the electronic properties of the system such as the broadening of the electronic density of states. An analysis of the Wannier functions indicate that quantum fluctuations of neat water molecules result in transient charging with subtle similarities and differences to that of the excess proton.

Aqueous solutions: state of the art in ab initio molecular dynamics.

The simulation of liquids by ab initio molecular dynamics (AIMD) has been a subject of intense activity over the last two decades. The significant increase in computational resources as well as the development of new and efficient algorithms has elevated this method to the status of a standard quantum mechanical tool that is used by both experimentalists and theoreticians. As AIMD computes the electronic structure from first principles, it is free of ad hoc parametrizations and has thus been applied to a large variety of physical and chemical problems. In particular, AIMD has provided microscopic insight into the structural and dynamical properties of aqueous solutions which are often challenging to probe experimentally. In this review, after a brief theoretical description of the Born-Oppenheimer and Car-Parrinello molecular dynamics formalisms, we show how AIMD has enhanced our understanding of the properties of liquid water and its constituent ions: the proton and the hydroxide ion. Thereafter, a broad overview of the application of AIMD to other aqueous systems, such as solvated organic molecules and inorganic ions, is presented. We also briefly describe the latest theoretical developments made in AIMD, such as methods for enhanced sampling and the inclusion of nuclear quantum effects.