Water Freezing and Ice Melting Małolepsza, Edyta; Keyes, Tom
Journal of chemical theory and computation,
2015-Dec-08, Volume:
11, Issue:
12
Journal Article
Peer reviewed
Open access
The generalized replica exchange method (gREM) is designed to sample states with coexisting phases and thereby to describe strong first order phase transitions. The isobaric MD version of the gREM is ...presented and applied to the freezing of liquid water and the melting of hexagonal and cubic ice. It is confirmed that coexisting states are well-sampled. The statistical temperature as a function of enthalpy, T S (H), is obtained. Hysteresis between freezing and melting is observed and discussed. The entropic analysis of phase transitions is applied and equilibrium transition temperatures, latent heats, and surface tensions are obtained for hexagonal ice ↔ liquid and cubic ice ↔ liquid with excellent agreement with published values. A new method is given to assign water molecules among various symmetry types. Pathways for water freezing, ultimately leading to hexagonal ice, are found to contain intermediate layered structures built from hexagonal and cubic ice.
Under ambient conditions, water freezes to either hexagonal ice or a hexagonal/cubic composite ice. The presence of hydrophobic guest molecules introduces a competing pathway: gas hydrate formation, ...with the guests in clathrate cages. Here, the pathways of the phase transitions are sought as sequences of states with coexisting phases, using a generalized replica exchange algorithm designed to sample them in equilibrium, avoiding nonequilibrium processes. For a dilute solution of methane in water under 200 atm, initializing the simulation with the full set of replicas leads to methane trapped in hexagonal/cubic ice, while gradually adding replicas with decreasing enthalpy produces the initial steps of hydrate growth. Once a small amount of hydrate is formed, water rearranges to form empty cages, eventually transforming the remainder of the system to metastable β ice, a scaffolding for hydrates. It is suggested that configurations with empty cages are reaction intermediates in hydrate formation when more guest molecules are available. Free energy profiles show that methane acts as a catalyst reducing the barrier for β ice versus hexagonal/cubic ice formation.
New sets of parameters (maps) for calculating amide I vibrational spectra for proteins through a vibrational exciton model are proposed. The maps are calculated as a function of electric field and ...van der Waals forces on the atoms of peptide bonds, taking into account the full interaction between peptide bonds and the surrounding environment. The maps are designed to be employed using data obtained from standard all-atom molecular simulations without any additional constraints on the system. Six proteins representing a wide range of sizes and secondary structure complexity were chosen as a test set. Spectra calculated for these proteins reproduce experimental data both qualitatively and quantitatively. The proposed maps lead to spectra that capture the weak second peak observed in proteins containing β-sheets, allowing for clear distinction between α-helical and β-sheet proteins. While the parametrization is specific to the CHARMM force field, the methodology presented can be readily applied to any empirical force field.
Combining genetic and cell-type-specific proteomic datasets can generate biological insights and therapeutic hypotheses, but a technical and statistical framework for such analyses is lacking. Here, ...we present an open-source computational tool called Genoppi (lagelab.org/genoppi) that enables robust, standardized, and intuitive integration of quantitative proteomic results with genetic data. We use Genoppi to analyze 16 cell-type-specific protein interaction datasets of four proteins (BCL2, TDP-43, MDM2, PTEN) involved in cancer and neurological disease. Through systematic quality control of the data and integration with published protein interactions, we show a general pattern of both cell-type-independent and cell-type-specific interactions across three cancer cell types and one human iPSC-derived neuronal cell type. Furthermore, through the integration of proteomic and genetic datasets in Genoppi, our results suggest that the neuron-specific interactions of these proteins are mediating their genetic involvement in neurodegenerative diseases. Importantly, our analyses suggest that human iPSC-derived neurons are a relevant model system for studying the involvement of BCL2 and TDP-43 in amyotrophic lateral sclerosis.
A prescription for sampling isobaric generalized ensembles with molecular dynamics is presented and applied to the generalized replica exchange method (gREM), which was designed to simulate ...first-order phase transitions. The properties of the isobaric gREM ensemble are discussed, and a study is presented for the liquid–vapor equilibrium of the guest molecules given for gas hydrate formation with the mW water model. Phase diagrams, critical parameters, and a law of corresponding states are obtained.
Reverse micelles are attractive nanoscale systems used for the confinement of molecules in studies of structure and chemical reactions, including protein folding, and aggregation. The simulation of ...reverse micelles, in which a water “pool” is separated from a nonpolar bulk phase by a surfactant layer, poses significant challenges to empirical force fields due to the diversity of interactions between nonpolar, polar, and charged groups. We have explored the dependence of system density, reverse micelle structure, and water configurational relaxation times as a function of reverse micelle composition, including water:surfactant ratio, absolute number of water molecules, and force field using molecular dynamics simulations. The resulting structures and dynamics are found to depend more on the force field used than on varying interpretations of the water:surfactant ratio in terms of absolute size of the reverse micelle. Substantial deviations from spherical reverse micelle geometries are observed in all unrestrained simulations. Rotational anisotropy decay times and water residence times show a strong dependence on force field and water model used, but power-law relaxation in time is observed independent of the force field. Our results suggest the need for further experimental study of reverse micelles that can provide insight into the distribution and dynamics of shape fluctuations in these complex systems.
Metastable β ice holds small guest molecules in stable gas hydrates, so its solid-liquid equilibrium is of interest. However, aqueous crystal-liquid transitions are very difficult to simulate. A new ...molecular dynamics algorithm generates trajectories in a generalized NPT ensemble and equilibrates states of coexisting phases with a selectable enthalpy. With replicas spanning the range between β ice and liquid water, we find the statistical temperature from the enthalpy histograms and characterize the transition by the entropy, introducing a general computational procedure for first-order transitions.
Knowledge of how intermolecular interactions of amyloid-forming proteins cause protein aggregation and how those interactions are affected by sequence and solution conditions is essential to our ...understanding of the onset of many degenerative diseases. Of particular interest is the aggregation of the amyloid-β (Aβ) peptide, linked to Alzheimer's disease, and the aggregation of the Sup35 yeast prion peptide, which resembles the mammalian prion protein linked to spongiform encephalopathies. To facilitate the study of these important peptides, experimentalists have identified small peptide congeners of the full-length proteins that exhibit amyloidogenic behavior, including the KLVFFAE sub-sequence, Aβ16-22, and the GNNQQNY subsequence, Sup357-13. In this study, molecular dynamics simulations were used to examine these peptide fragments encapsulated in reverse micelles (RMs) in order to identify the fundamental principles that govern how sequence and solution environment influence peptide aggregation. Aβ16-22 and Sup357-13 are observed to organize into anti-parallel and parallel β-sheet arrangements. Confinement in the sodium bis(2-ethylhexyl) sulfosuccinate (AOT) reverse micelles is shown to stabilize extended peptide conformations and enhance peptide aggregation. Substantial fluctuations in the reverse micelle shape are observed, in agreement with earlier studies. Shape fluctuations are found to facilitate peptide solvation through interactions between the peptide and AOT surfactant, including direct interaction between non-polar peptide residues and the aliphatic surfactant tails. Computed amide I IR spectra are compared with experimental spectra and found to reflect changes in the peptide structures induced by confinement in the RM environment. Furthermore, examination of the rotational anisotropy decay of water in the RM demonstrates that the water dynamics are sensitive to the presence of peptide as well as the peptide sequence. Overall, our results demonstrate that the RM is a complex confining environment where substantial direct interaction between the surfactant and peptides plays an important role in determining the resulting ensemble of peptide conformations. By extension the results suggest that similarly complex sequence-dependent interactions may determine conformational ensembles of amyloid-forming peptides in a cellular environment.
The cohesin complex plays an essential role in chromosome maintenance and transcriptional regulation. Recurrent somatic mutations in the cohesin complex are frequent genetic drivers in cancer, ...including myelodysplastic syndromes (MDS) and acute myeloid leukemia (AML). Here, using genetic dependency screens of stromal antigen 2-mutant (STAG2-mutant) AML, we identified DNA damage repair and replication as genetic dependencies in cohesin-mutant cells. We demonstrated increased levels of DNA damage and sensitivity of cohesin-mutant cells to poly(ADP-ribose) polymerase (PARP) inhibition. We developed a mouse model of MDS in which Stag2 mutations arose as clonal secondary lesions in the background of clonal hematopoiesis driven by tet methylcytosine dioxygenase 2 (Tet2) mutations and demonstrated selective depletion of cohesin-mutant cells with PARP inhibition in vivo. Finally, we demonstrated a shift from STAG2- to STAG1-containing cohesin complexes in cohesin-mutant cells, which was associated with longer DNA loop extrusion, more intermixing of chromatin compartments, and increased interaction with PARP and replication protein A complex. Our findings inform the biology and therapeutic opportunities for cohesin-mutant malignancies.