Molecular dynamics simulation on miscibility of trans-1,4,5,8-tetranitro-1,4,5,8 -tetraazadecalin (TNAD) with some propellants

The solubility parameters of TNAD, HMX, RDX, DINA, DNP propellants were predicted by molecular dynamics (MD) simulation in order to evaluate the miscibility of TNAD and the other four propellants. The results show that the order of miscibility is TNAD/DINA > TNAD/DNP > TNAD/RDX > TNAD/HMX from the analysis of miscibility. The densities and binding energies of TNAD/propellants blends were further investigated. The results indicate that the better the miscibility between TNAD and the propellants, the smaller the variation of the density rate. The larger the intermolecular interaction, the better the miscibility between components. The analysis of radial distribution function shows that the main interaction ways between TNAD and other energetic components are short-range interactions.

FT-IR imaging spectroscopy of phase separation in blends of poly(3-hydroxybutyrate) with poly(L-lactic acid) and poly(epsilon-caprolactone)

The detection of phase separation and identification of miscibility in biopolymer blends is an important aspect for the improvement of their physical properties. In this article, the phase separation in blends of poly(3-hydroxybutyrate) (PHB) with poly(L-lactic acid) (PLA) and poly(epsilon-caprolactone) (PCL), respectively, has been studied as a function of the blend composition by FT-IR imaging spectroscopy. For both polymer blend systems, a miscibility gap has been found around the 50:50% (w/w) composition of the two components. Furthermore, the separating phases have been identified as blends of the two polymer components and their compositions could be determined from calibrations based on the spectra of the blends in the compositional range of miscibility. The data derived from FT-IR spectroscopic imaging were corroborated by additional DSC analyses and mechanical stress-strain measurements of polymer blend films, which exhibited a characteristic fracture behavior as a function of PHB composition.     

Molecular dynamics simulations of microstructure and mixing dynamics of cryoprotective solvents in water and in the presence of a lipid membrane

Molecular dynamics (MD) simulation is used to investigate the solubility behavior of cryoprotective (CP) solvents, such as DMSO, ethylene glycol (EG) and glycerol (GL), in pure water and in the presence of a lipid membrane. The MD study is focused on an equilibration timescale required for mixing large CP aggregates with aqueous and aqueous/lipid environments. The MD analysis demonstrates that DMSO mixes rapidly with water, so that all solute molecules are uniformly distributed in the equilibrium aqueous solution. Our investigation of the microstructure of binary EG/water and GL/water systems reveals that, despite the miscibility of both CP solvents with water, they are not ideally mixed in aqueous solutions at the molecular level. The MD simulations show that the mixing dynamics of the large CP cluster and surrounding water is found to be strongly dependent on nature of hydrophilic and hydrophobic interactions acting between cryoprotectant molecules. In particular, a spatial hydrogen-bond network formed between CP molecules plays an important role in the mixing dynamics between CP agents and water. A further analysis on the mixing behavior of the CP solvents with pure water and with aqueous solutions at a lipid membrane interface shows that, due to strong binding of the CP molecules to membrane surface, the equilibration process in the lipid environment becomes very slow, at least of the order of microseconds. The MD results are discussed in the context of the better understanding on the composition of the aqueous mixtures of the EG and GL solvents. Knowledge of the microstructure and the dynamics of these systems helps to develop better cryopreservation protocols and to propose more optimal cooling/warming regimes for cellular cryosolutions.     

Morphology and properties of soy protein and polylactide blends

Blends of soy protein (SP) and a semicrystalline polylactide (PLA) were prepared using a twin-screw extruder. The melt rheology, phase morphology, mechanical properties, water resistance, and thermal and dynamic mechanical properties were investigated on specimens prepared by injection molding of these blends. The melt flowability of soy-based plastics was improved through blending with PLA. Scanning electron microscopy revealed that a co-continuous phase structure existed in the blends with soy protein concentrate (SPC) to PLA ratios ranging from 30:70 to 70:30. SPC/PLA blends showed fine co-continuous phase structures, while soy protein isolate (SPI)/PLA blends presented severe phase coarsening. At the same SP to PLA ratios, SPC/PLA blends demonstrated a higher tensile strength than SPI/PLA blends. The water absorption of soy plastics was greatly reduced by blending with PLA. The compatibility was improved by adding 1-5 phr poly(2-ethyl-2-oxazoline) (PEOX) in the blends, and the resulting blends showed an obvious increase in tensile strength and a reduction in water absorption for SPI/PLA blends. The compatibility between SP and PLA was evaluated by mechanical testing, dynamic mechanical analysis (DMA), water absorption, and scanning electron microscopy (SEM) experiments. Differential scanning calorimetry (DSC) revealed that PLA in the blends was mostly amorphous in the injection molded articles, and SP accelerated the cold crystallization and could increase the final crystallinity of PLA in the blends.     

Conformational dynamics of full-length inducible human Hsp70 derived from microsecond molecular dynamics simulations in explicit solvent

Human 70?kDa heat shock protein (hHsp70) is an ATP-dependent chaperone and is currently an important target for developing new drugs in cancer therapy. Knowledge of the conformations of hHsp70 is central to understand the interactions between its nucleotide-binding domain (NBD) and substrate-binding domain (SBD) and is a prerequisite to design inhibitors. The conformations of ADP-bound (or nucleotide-free) hHsp70 and ATP-bound hHsp70 was investigated by using unbiased all-atom molecular dynamics (MD) simulations of homology models of hHsp70 in explicit solvent on a timescale of .5 and 2.7 μs, respectively. The conformational heterogeneity of hHsp70 was analyzed by computing effective free-energy landscapes (FELs) and distance distribution between selected pair of residues. These theoretical data were compared with those extracted from single-molecule Förster resonance energy transfer (FRET) experiments and to small-angle X-ray scattering (SAXS) data of Hsp70 homologs. The distance between a pair of residues in FRET is systematically larger than the distance computed in MD which is interpreted as an effect of the size and of the dynamics of the fluorescent probes. The origin of the conformational heterogeneity of hHsp70 in the ATP-bound state is due to different binding modes of the helix B of the SBD onto the NBD. In the ADP-bound (or nucleotide-free) state, it arises from the different closed conformations of the SBD and from the different positions of the SBD relative to the NBD. In each nucleotide-binding state, Hsp70 is better represented by an ensemble of conformations on a μs timescale corresponding to different local minima of the FEL.

An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:30     

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.     

Molecular dynamics simulation of the stability of a 22-residue α-helix in water and 30% trifluoroethanol

A molecular dynamics (MD) simulation was performed on the α-helix H8-HC5, the C-terminal part of myoglobin (residue 132–153), under periodic boundary conditions in two different solutions, water and water with 30% (v/v) 2,2,2-trifluoroethanol (TFE), at 300 K to investigate the stability of the helix. In both simulations, the initial configuration was a canonical right-handed α-helix. In the course of the MD trajectory in water (200 ps), the helix clearly destabilized and began to unfold after 100 ps. In the TFE solution, two stable parts of helical regions were observed after 70 ps of a 200-ps MD simulation, supporting the notion that TFE acts as a structure-forming solvent. © 1993 John Wiley & Sons, Inc.     

A study of conformational stability of polyglycine and poly(L-alanine), and polyglycine/poly(L-alanine) blends in the solid state by (13)C cross-polarization/magic angle spinning NMR

13C Cross-Polarization/Magic Angle Spinning nmr and T(1rhoH) experiments of polyglycine (PG), poly(L-alanine) (PLA), and PG/PLA blends prepared from dichloroacetic acid solution have been carried out, in order to elucidate the conformational stability of these polypeptides in the solid state. From these experimental results, it was clarified that the conformations of PG and PLA in their blends are strongly influenced by intermolecular hydrogen-bonding interactions that cause their miscibility at the molecular level.     

Compatibilization effect of poly(epsilon-caprolactone)-b-poly(ethylene glycol) block copolymers and phase morphology analysis in immiscible poly(lactide)/poly(epsilon-caprolactone) blends

The miscibility and phase behavior of two stereoisomer forms of poly(lactide) (PLA: poly (L-lactide) (PLLA) and poly(DL-lactide) (PDLLA)) blends with poly(epsilon-caprolactone)-b-poly(ethylene glycol) (PCL-b-PEG) and PCL-b-monomethoxy-PEG (PCL-b-MPEG) block copolymers have been investigated by differential scanning calorimetry (DSC). The DSC thermal behavior of both the blend systems revealed that PLA is miscible with the PEG segment phase of PCL-b-(M)PEG but is still immiscible with its PCL segment phase although PCL was block-copolymerized with PEG. On the basis of these results, PCL-b-PEG was added as a compatibilizer to PLA/PCL binary blends. The improvement in mechanical properties of PLA/PCL blends was achieved as anticipated upon the addition of PCL-b-PEG. In addition, atomic force microscopy (AFM) measurements have been performed in order to study the compositional synergism to be observed in mechanical tests. AFM observations of the morphological dependency on blend composition indicate that PLA/PCL blends are immiscible but compatible to some extent and that synergism of compatibilizing may be maximized in the compositional blend ratio before apparent phase separation and coarsening.     

Miscibility studies of HPMC/PVA blends in water by viscosity, density, refractive index and ultrasonic velocity method

Hydroxy propyl methyl cellulose (HPMC)/polyvinyl alcohol (PVA) blends are edible polymer films used for food packing and directly in foodstuffs. However they are water-soluble in ordinary temperature and have good mechanical properties. The miscibility of HPMC/PVA blend in water was studied by viscosity, ultrasonic velocity, density and refractive index techniques at 30 and 50 °C. Using viscosity data, the interaction parameters μ and α were calculated. These values revealed that HPMC/PVA blend is miscible when the HPMC content is more than 60% in the blend at 30 and 50 °C. And also the result revealed that the change in temperature has no significant effect on the miscibility of HPMC/PVA polymer blend.     

A study of conformational stability of poly(L-alanine), poly(L-valine), and poly(L-alanine)/poly(L-valine) blends in the solid state by (13)C cross-polarization/magic angle spinning NMR

13C cross-polarization/magic angle spinning (CP/MAS) NMR and (1)H T(1rho) experiments of poly(L-alanine) (PLA), poly(L-valine) (PLV), and PLA/PLV blends have been carried out in order to elucidate the conformational stability of the polypeptides in the solid state. These were prepared by adding a trifluoroacetic acid (TFA) solution of the polymer with a 2.0 wt/wt % of sulfuric acid (H(2)SO(4)) to alkaline water. From these experimental results, it is clarified that the conformations of PLA and PLV in their blends are strongly influenced by intermolecular hydrogen-bonding interactions that cause their miscibility at the molecular level.     

Role of chemical short-range order in atomic dynamics decoupling

Abstract

Using molecular dynamics simulation, α-relaxation times τ_α and self-diffusion coefficients D for Al₉₀Fe₁₀, Al₈₀Fe₂₀, Al₇₀Fe₃₀, Al₆₀Fe₄₀ and Al₈₀Ni₂₀ (as a contrast system) melts have been systematically computed over a wide temperature range (1000–2000 K). The computed results reveal that τ_Fe/τ_Al (or D_Al/D_Fe) for the Al₉₀Fe₁₀ and Al₈₀Fe₂₀ melts exhibit an accelerating increase with cooling at temperatures lower than 1400 K, implying a clear decoupling of dynamics of Al and Fe (here referred to as component decoupling). This component decoupling diminishes in Al₇₀Fe₃₀ melt and disappears in Al₆₀Fe₄₀ melt. We simultaneously checked the relaxation decoupling (i.e. the decoupling between α-relaxation and diffusion). The relaxation decoupling is clear in Al₆₀Fe₄₀ melt, less clear in Al₇₀Fe₃₀ melt and not shown in Al₈₀Fe₂₀ and Al₉₀Fe₁₀ melt. It exhibits a tendency counter to that of component decoupling with changing composition, arguing that relaxation decoupling does not necessarily lead to component decoupling. This finding is contradicted against the conventional view that component decoupling is believed as a result of relaxation decoupling. We further attributed such a contradiction to the difference in the degree of chemical short-range order (CSRO) in melts. The existence of CSRO can increase the cooperativity in dynamics of different components. So it is better to consider component decoupling as a combined effect of relaxation decoupling and CSRO. This work would be helpful in improving our understanding of the relationship between the two kinds of decoupling.     

In situ formation of blends by photopolymerization of poly(ethylene glycol) dimethacrylate and polylactide

Blends of cross-linked poly(ethylene glycol) dimethacrylate (PEGDMA) and poly(d,l-lactide) (PLA) were prepared by mixing photoactive PEGDMA (molecular mass: 875 g/mol) and PLA, and subsequently photopolymerizing the mixture with visible light. The effects of PLA molecular mass and mass fraction on the rheological properties of the PEGDMA/PLA mixtures, and on the degree of methacrylate vinyl conversion (DC), as well as blend miscibility, microstructure, mechanical properties, in vitro swelling behavior, and cell responses were studied. PLA-2K (molecular mass: 2096 g/mol) and PLA-63K (molecular mass: 63 000 g/mol) formed miscible and partially miscible blends with cross-linked PEGDMA, respectively. The addition of the PLA-2K did not affect the immediate or post-cure (>24 h) DC of the PEGDMA upon photopolymerization. However, the addition of PLA-63K decreased the immediate DC of the PEGDMA, which can be increased through extending the curing time or post-curing period. Compared to the cross-linked neat PEGDMA and PLA-2K/PEGDMA blends, PLA-63K/PEGDMA blends were significantly stronger, stiffer, and tougher. Both types of blends and the cross-linked PEGDMA swelled when soaked in a phosphate buffered saline (PBS) solution. The attachment and spreading of MCT3-E1 cells increased with increasing PLA-63K content in the blends. The facile and rapid formation of PEGDMA/PLA blends by photopolymerization represents a simple and efficient approach to a class of biomaterials with a broad spectrum of properties.     

Molecular dynamics simulation of a phospholipid membrane

We present the results of molecular dynamics (MD) simulations of a phospholipid membrane in water, including full atomic detail. The goal of the simulations was twofold: first we wanted to set up a simulation system which is able to reproduce experimental results and can serve as a model membrane in future simulations. This goal being reached it is then further possible to gain insight in to those properties that are experimentally more difficult to access. The system studied is dipalmitoylphosphatidylcholine/water, consisting of 5408 atoms. Using original force field parameters the membrane turned out to approach a gel-like state. With slight changes of the parameters, the system adopted a liquid-crystalline state. Separate 80 ps runs were performed on both the gel and liquid-crystalline systems. Comparison of MD results with reliable experimental data (bilayer repeat distance, surface area per lipid, tail order parameters, atom distributions) showed that our simulations, especially the one in the liquid-crystalline phase, can serve as a realistic model for a phospholipid membrane. Further analysis of the trajectories revealed valuable information on various properties. In the liquid-crystalline phase, the interface turns out to be quite diffuse, with water molecules penetrating into the bilayer to the position of the carbonyl groups. The 10–90% width of the interface turns out to be 1.3 nm and the width of the hydrocarbon interior 3.0 nm. The headgroup dipoles are oriented at a small angle with respect to the bilayer plane. The resulting charge distribution is almost completely cancelled by the water molecules. The electron density distribution shows a large dip in the middle of the membrane. In this part the tails are more flexible. The mean life time between dihedral transitions is 20 ps. The average number of gauche angles per tail is 3.5. The occurrence of kinks is not a significant feature.Abbreviations MD molecular dynamics - DPPC dipalmitoylphosphatidylcholine - SPC simple point charges - DPPE dipalmitoylphosphatidylethanolamine Correspondence to: H. J. C. Berendsen     

The new program OPAL for molecular dynamics simulations and energy refinements of biological macromolecules

Summary A new program for molecular dynamics (MD) simulation and energy refinement of biological macromolecules, OPAL, is introduced. Combined with the supporting program TRAJEC for the analysis of MD trajectories, OPAL affords high efficiency and flexibility for work with diferent force fields, and offers a user-friendly interface and extensive trajectory analysis capabilities. Salient features are computational speeds of up to 1.5 GFlops on vector supercomputers such as the NEC SX-3, ellipsoidal boundaries to reduce the system size for studies in explicit solvents, and natural treatment of the hydrostatic pressure. Practical applications of OPAL are illustrated with MD simulations of pure water, energy minimization of the NMR structure of the mixed disulfide of a mutant E. coli glutaredoxin with glutathione in different solvent models, and MD simulations of a small protein, pheromone Er-2, using either instantaneous or time-averaged NMR restraints, or no restraints.Abbreviations D diffusion constant in cm²/s - Er-2 pheromone 2 from Euplotes raikovi - GFlop one billion floating point operations per second - Grx(C14S)-SG mixed disulfide between a mutant E. coli glutaredoxin, with Cys¹⁴ replaced by Ser, and glutathione - MD molecular dynamics - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - density in g/cm³     

Compatible ternary blends of chitosan/poly(vinyl alcohol)/poly(lactic acid) produced by oil-in-water emulsion processing

Ternary compatible blends of chitosan, poly(vinyl alcohol), and poly(lactic acid) were prepared by an oil-in-water (O/W) emulsion process. Solutions of chitosan in aqueous acetic acid, poly(vinyl alcohol) (PVA) in water, and poly(lactic acid) (PLA) in chloroform were blended with a high-shear mixer. PVA was used as an emulsifier to stabilize the emulsion and to reduce the interfacial tension between the solid polymers in the blends produced. It proved to work very well because the emulsions were stable for periods of days or weeks and compatible blends were obtained when PVA was added. This effect was attributed to a synergistic effect of PVA and chitosan because the binary blends PVA/PLA and chitosan/PLA were completely incompatible. The blends were characterized by scanning electron microscopy (SEM), differential scanning calorimetry (DSC), thermal mechanical analysis (TMA), stress-strain tests, and Fourier transform infrared spectroscopy (FTIR). The results indicated that despite the fact that the system contained distinct phases some degree of molecular miscibility occurred when the three components were present in the blend.     

Molecular dynamics simulations of the lipid bilayer edge

Phospholipid bilayers have been intensively studied by molecular dynamics (MD) simulation in recent years. The properties of bilayer edges are important in determining the structure and stability of pores formed in vesicles and biomembranes. In this work, we use molecular dynamics simulation to investigate the structure, dynamics, and line tension of the edges of bilayer ribbons composed of pure dimyristoylphosphatidylcholine (DMPC) or palmitoyl-oleoylphosphatidylethanolamine (POPE). As expected, we observe a significant reorganization of lipids at and near the edges. The treatment of electrostatic effects is shown to have a qualitative impact on the structure and stability of the edge, and significant differences are observed in the dynamics and structure of edges formed by DMPC and palmitoyl-oleoylphosphatidylethanolamine. From the pressure anisotropy in the simulation box, we calculate a line tension of approximately 10-30 pN for the DMPC edge, in qualitative agreement with experimental estimates for similar lipids.     

New heat shock protein (Hsp90) inhibitors,designed by pharmacophore modeling and virtual screening: synthesis,biological evaluation and molecular dynamics studies

Abstract

Inhibition of heat shock protein 90 (Hsp90) is known to be a significantly effective strategy in cancer therapy. Here, pyrazolopyranopyrimidine derivatives were characterized as new Hsp90 inhibitors. The molecules’ key structure (ZINC02819805) was determined by utilizing a pharmacophore model virtual screening workflow. Structural optimization was then carried out on compound ZINC02819805, pyrazolopyranopyrimidine derivatives were designed and six chosen derivatives were synthesized. The inhibition of Hsp90 ATPase activity of synthesized compounds revealed that para methylphenyl derivative of pyrazolopyranopyrimidine (HM3) was the most potent inhibitor (IC₅₀ = 5.5?µM). The anti-proliferative activity of this compound was evaluated against a panel of cell lines including MCF-7, HeLa and HUVEC (IC₅₀ = 1.28?µM, IC₅₀ = 1.74?µM and IC₅₀ = 61.48?µM respectively) by MTT method. The western blot analysis of treated MCF-7 cells with compound HM3 showed that the expression level of Hsp70 and Her2 proteins changed. The high level of Hsp70 expression and low level of Her2 expression suggest that compound HM3 exhibits inhibitory effect on Hsp90. Finally, the key interactions between HM3 and Hsp90 protein were studied by molecular dynamics simulation and showed that compound HM3 was stable in Hsp90 active cite during 200?ns simulation. Abbreviations Hsp90 Heat shock protein 90

MTT 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide

ATP adenosine triphosphate

MD molecular dynamics simulation

RMSD root-mean-square deviation

RMSF root-mean-square fluctuation

Rg gyration radius

m-SABNPs boehmite nanoparticles-supported sulfamic acid

Communicated by Ramaswamy H. Sarma