Molecular dynamics simulation of brittle fracture in silicon

The fracture process involves converting potential energy from a strained body into surface energy, thermal energy, and the energy needed to create lattice defects. In dynamic fracture, energy is also initially converted into kinetic energy. This paper uses molecular dynamics (MD) to simulate brittle frcture in silicon and determine how energy is converted from potential energy (strain energy) into other forms.

Molecular dynamics simulation of methane in sodium montmorillonite clay hydrates at elevated pressures and temperatures

Computer simulation has been used to study the structure and dynamics of methane in hydrated sodium montmorillonite clays under conditions encountered in sedimentary basins. Systems containing approximately one, two, three and four molecular layers of water have followed gradients of 150 bar km-1 and 30Kkm-1, to a maximum burial depth of 6 km (900 bar and 460 K). Methane is coordinated to approximately 19 oxygen atoms, of which typically 6 are provided by the clay surface. Only in the three- and four-layer hydrates is methane able to leave the clay surface. Diffusion depends strongly on the porosity (water content) and burial depth: self-diffusion coefficients are in the range 0.12 × 10-9m2s-1 for water and 0.04 × 10−9m2s−1 < D < 8.64 × 10−9m2s−1 for methane. Bearing in mind that porosity decreases with burial depth, it is estimated that maximum diffusion occurs at around 3 km. This is in good agreement with the known location of methane reservoirs in sedimentary basins.

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

Molecular Dynamics implementation of BN2D or 'Mercedes Benz' water model

Two-dimensional 'Mercedes Benz' (MB) or BN2D water model (Naim, 1971) is implemented in Molecular Dynamics. It is known that the MB model can capture abnormal properties of real water (high heat capacity, minima of pressure and isothermal compressibility, negative thermal expansion coefficient) (Silverstein et al., 1998). In this work formulas for calculating the thermodynamic, structural and dynamic properties in microcanonical (NVE) and isothermal-isobaric (NPT) ensembles for the model from Molecular Dynamics simulation are derived and verified against known Monte Carlo results. The convergence of the thermodynamic properties and the system's numerical stability are investigated. The results qualitatively reproduce the peculiarities of real water making the model a visually convenient tool that also requires less computational resources, thus allowing simulations of large (hydrodynamic scale) molecular systems. We provide the open source code written in C/C++ for the BN2D water model implementation using Molecular Dynamics.

Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

Multiscale molecular dynamics/hydrodynamics implementation of two dimensional “Mercedes Benz” water model

A multiscale Molecular Dynamics/Hydrodynamics implementation of the 2D Mercedes Benz (MB or BN2D) [1] water model is developed and investigated. The concept and the governing equations of multiscale coupling together with the results of the two-way coupling implementation are reported. The sensitivity of the multiscale model for obtaining macroscopic and microscopic parameters of the system, such as macroscopic density and velocity fluctuations, radial distribution and velocity autocorrelation functions of MB particles, is evaluated. Critical issues for extending the current model to large systems are discussed.

A hybrid molecular dynamics/fluctuating hydrodynamics method for modelling liquids at multiple scales in space and time

A new 3D implementation of a hybrid model based on the analogy with two-phase hydrodynamics has been developed for the simulation of liquids at microscale. The idea of the method is to smoothly combine the atomistic description in the molecular dynamics zone with the Landau-Lifshitz fluctuating hydrodynamics representation in the rest of the system in the framework of macroscopic conservation laws through the use of a single "zoom-in" user-defined function s that has the meaning of a partial concentration in the two-phase analogy model. In comparison with our previous works, the implementation has been extended to full 3D simulations for a range of atomistic models in GROMACS from argon to water in equilibrium conditions with a constant or a spatially variable function s. Preliminary results of simulating the diffusion of a small peptide in water are also reported.

Self-healing of cracks during ductile regime machining of silicon: Insights from molecular dynamics simulation

During nanoindentation and ductile-regime machining of silicon, a phenomenon known as “self-healing” takes place in that the microcracks, microfractures, and small spallings generated during the machining are filled by the plastically flowing ductile phase of silicon. However, this phenomenon has not been observed in simulation studies. In this work, using a long-range potential function, molecular dynamics simulation was used to provide an improved explanation of this mechanism. A unique phenomenon of brittle cracking was discovered, typically inclined at an angle of 45° to 55° to the cut surface, leading to the formation of periodic arrays of nanogrooves being filled by plastically flowing silicon during cutting. This observation is supported by the direct imaging. The simulated X-ray diffraction analysis proves that in contrast to experiments, Si-I to Si-II (beta tin) transformation during ductile-regime cutting is highly unlikely and solid-state amorphisation of silicon caused solely by the machining stress rather than the cutting temperature is the key to its brittle-ductile transition observed during the MD simulations

Analysis of Molecular Dynamics Simulations of Protein Folding

Microsecond long Molecular Dynamics (MD) trajectories of biomolecular processes are now possible due to advances in computer technology. Soon, trajectories long enough to probe dynamics over many milliseconds will become available. Since these timescales match the physiological timescales over which many small proteins fold, all atom MD simulations of protein folding are now becoming popular. To distill features of such large folding trajectories, we must develop methods that can both compress trajectory data to enable visualization, and that can yield themselves to further analysis, such as the finding of collective coordinates and reduction of the dynamics. Conventionally, clustering has been the most popular MD trajectory analysis technique, followed by principal component analysis (PCA). Simple clustering used in MD trajectory analysis suffers from various serious drawbacks, namely, (i) it is not data driven, (ii) it is unstable to noise and change in cutoff parameters, and (iii) since it does not take into account interrelationships amongst data points, the separation of data into clusters can often be artificial. Usually, partitions generated by clustering techniques are validated visually, but such validation is not possible for MD trajectories of protein folding, as the underlying structural transitions are not well understood. Rigorous cluster validation techniques may be adapted, but it is more crucial to reduce the dimensions in which MD trajectories reside, while still preserving their salient features. PCA has often been used for dimension reduction and while it is computationally inexpensive, being a linear method, it does not achieve good data compression. In this thesis, I propose a different method, a nonmetric multidimensional scaling (nMDS) technique, which achieves superior data compression by virtue of being nonlinear, and also provides a clear insight into the structural processes underlying MD trajectories. I illustrate the capabilities of nMDS by analyzing three complete villin headpiece folding and six norleucine mutant (NLE) folding trajectories simulated by Freddolino and Schulten [1]. Using these trajectories, I make comparisons between nMDS, PCA and clustering to demonstrate the superiority of nMDS. The three villin headpiece trajectories showed great structural heterogeneity. Apart from a few trivial features like early formation of secondary structure, no commonalities between trajectories were found. There were no units of residues or atoms found moving in concert across the trajectories. A flipping transition, corresponding to the flipping of helix 1 relative to the plane formed by helices 2 and 3 was observed towards the end of the folding process in all trajectories, when nearly all native contacts had been formed. However, the transition occurred through a different series of steps in all trajectories, indicating that it may not be a common transition in villin folding. The trajectories showed competition between local structure formation/hydrophobic collapse and global structure formation in all trajectories. Our analysis on the NLE trajectories confirms the notion that a tight hydrophobic core inhibits correct 3-D rearrangement. Only one of the six NLE trajectories folded, and it showed no flipping transition. All the other trajectories get trapped in hydrophobically collapsed states. The NLE residues were found to be buried deeply into the core, compared to the corresponding lysines in the villin headpiece, thereby making the core tighter and harder to undo for 3-D rearrangement. Our results suggest that the NLE may not be a fast folder as experiments suggest. The tightness of the hydrophobic core may be a very important factor in the folding of larger proteins. It is likely that chaperones like GroEL act to undo the tight hydrophobic core of proteins, after most secondary structure elements have been formed, so that global rearrangement is easier. I conclude by presenting facts about chaperone-protein complexes and propose further directions for the study of protein folding.

Data of the molecular dynamics simulations of mutations in the human connexin46 docking interface

The structure of hCx26 derived from the X-ray analysis was used to generate a homology model for hCx46. Interacting connexin molecules were used as starting model for the molecular dynamics (MD) simulation using NAMD and allowed us to predict the dynamic behavior of hCx46wt and the cataract related mutant hCx46N188T as well as two artificial mutants hCx46N188Q and hCx46N188D. Within the 50 ns simulation time the docked complex composed of the mutants dissociate while hCx46wt remains stable. The data indicates that one hCx46 molecule forms 5-7 hydrogen bonds (HBs) with the counterpart connexin of the opposing connexon. These HBs appear essential for a stable docking of the connexons as shown by the simulation of an entire gap junction channel and were lost for all the tested mutants. The data described here are related to the research article entitled "The cataract related mutation N188T in human connexin46 (hCx46) revealed a critical role for residue N188 in the docking process of gap junction channels" (Schadzek et al., 2015) [1].

Distribution of neutral prilocaine in a phospholipid bilayer: Insights from molecular dynamics simulations

In this work, we report a 20-ns constant pressure molecular dynamics simulation of prilocaine (PLC), in amine-amide local anesthetic, in a hydrated liquid crystal bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine. The partition of PLC induces the lateral expansion of the bilayer and a concomitant contraction in its thickness. PLC molecules are preferentially found in the hydrophobic acyl chains region, with a maximum probability at similar to 12 angstrom from the center of the bilayer (between the C(4) and C(5) methylene groups). A decrease in the acyl chain segmental order parameter, vertical bar S-CD vertical bar, compared to neat bilayers, is found, in good agreement with experimental H-2-NMR studies. The decrease in vertical bar S-CD vertical bar induced by PLC is attributed to a larger accessible volume per lipid in the acyl chain region. (C) 2008 Wiley Periodicals, Inc.

Preferential location of prilocaine and etidocaine in phospholipid bilayers: A molecular dynamics study

In this work, we report a 20-ns constant pressure molecular dynamics simulation of the uncharged form of two amino-amide local anesthetics (LA). etidocaine and prilocaine, present at 1:3 LA:lipid, molar ratio inside the membrane, in the hydrated liquid crystal bilayer phase of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC). Both LAs induced lateral expansion and a concomitant contraction in the bilayer thickness. A decrease in the acyl chain segment order parameter, -S(CD), compared to neat bilayers, was also observed. Besides, both LA molecules got preferentially located in the hydrophobic acyl chains region, with a maximum probability at similar to 12 and similar to 10 angstrom from the center of the bilayer for prilocaine and etidocaine, respectively. (C) 2009 Elsevier B.V. All rights reserved.

PARAMETRIC STUDY OF REAXFF SIMULATION PARAMETERS FOR MOLECULAR DYNAMICS MODELING OF REACTIVE CARBON GASES

Abstract The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although the Reax Force Field (ReaxFF) can be used to simulate the chemical behavior of carbon-based systems, the simulation settings required for accurate predictions have not been fully explored. Using the ReaxFF, molecular dynamics (MD) simulations are used to simulate the chemical behavior of pure carbon and hydrocarbon reactive gases that are involved in the formation of carbon structures such as graphite, buckyballs, amorphous carbon, and carbon nanotubes. It is determined that the maximum simulation time step that can be used in MD simulations with the ReaxFF is dependent on the simulated temperature and selected parameter set, as are the predicted reaction rates. It is also determined that different carbon-based reactive gases react at different rates, and that the predicted equilibrium structures are generally the same for the different ReaxFF parameter sets, except in the case of the predicted formation of large graphitic structures with the Chenoweth parameter set under specific conditions.

Diphenylhexatriene membrane probes DPH and TMA-DPH: A comparative molecular dynamics simulation study

Fluorescence spectroscopy andmicroscopy have been utilized as tools in membrane biophysics for decades now. Because phospholipids are non-fluorescent, the use of extrinsic membrane probes in this context is commonplace. Among the latter, 1,6-diphenylhexatriene (DPH) and its trimethylammonium derivative (TMA-DPH) have been extensively used. It is widely believed that, owing to its additional charged group, TMA-DPH is anchored at the lipid/water interface and reports on a bilayer region that is distinct from that of the hydrophobic DPH. In this study, we employ atomistic MD simulations to characterize the behavior of DPH and TMA-DPH in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and POPC/cholesterol (4:1) bilayers. We show that although the dynamics of TMA-DPH in thesemembranes is noticeably more hindered than that of DPH, the location of the average fluorophore of TMA-DPH is only ~3–4 Å more shallow than that of DPH. The hindrance observed in the translational and rotational motions of TMA-DPH compared to DPH is mainly not due to significant differences in depth, but to the favorable electrostatic interactions of the former with electronegative lipid atoms instead. By revealing detailed insights on the behavior of these two probes, our results are useful both in the interpretation of past work and in the planning of future experiments using themasmembrane reporters.

Diphenylhexatriene membrane probes DPH and TMA-DPH: A comparative molecular dynamics simulation study

Fluorescence spectroscopy andmicroscopy have been utilized as tools in membrane biophysics for decades now. Because phospholipids are non-fluorescent, the use of extrinsic membrane probes in this context is commonplace. Among the latter, 1,6-diphenylhexatriene (DPH) and its trimethylammonium derivative (TMA-DPH) have been extensively used. It is widely believed that, owing to its additional charged group, TMA-DPH is anchored at the lipid/water interface and reports on a bilayer region that is distinct from that of the hydrophobic DPH. In this study, we employ atomistic MD simulations to characterize the behavior of DPH and TMA-DPH in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and POPC/cholesterol (4:1) bilayers. We show that although the dynamics of TMA-DPH in thesemembranes is noticeably more hindered than that of DPH, the location of the average fluorophore of TMA-DPH is only ~3–4 Å more shallow than that of DPH. The hindrance observed in the translational and rotational motions of TMA-DPH compared to DPH is mainly not due to significant differences in depth, but to the favorable electrostatic interactions of the former with electronegative lipid atoms instead. By revealing detailed insights on the behavior of these two probes, our results are useful both in the interpretation of past work and in the planning of future experiments using themasmembrane reporters.