The self-diffusion and hydrogen bond interaction in neat liquid alkanols: a molecular dynamic simulation study

Self-diffusion of methanol, ethanol, 1-propanol and 2-propanol has been studied by molecular dynamics simulation in the temperature range between the melting pressure curve and 478 K at pressures up to 300 MPa. The simulation results on self-diffusion of methanol, ethanol and 2-propanol (for 2-propanol, at high temperatures) agree well with experiment, which suggests that the simulation method is a powerful tool to obtain self-diffusion coefficients over wide range of temperature and pressure, under which it is rather difficult for experiments. The local structures of methanol, ethanol and 2-propanol are investigated by calculating the radial distribution functions, H-bond numbers, coordination numbers and the ratios of H-bond number divided by coordination number. The correlation between self-diffusion and structural properties, and the influence of temperature and pressure on them are discussed. The degree of forming H-bond space network in methanol, ethanol and water is higher than that in 2-propanol, and they are all higher than those in ammonia and methylamine. The simulation results demonstrate that the effect of hydrogen bonding on the translational dynamics in methanol and ethanol is more pronounced than that in 2-propanol.     

Calculation of the diffusion coefficient of thiophene oligomers using molecular dynamics

Electropolymerisation is a very useful methodology for conducting polymers synthesis. A total comprehension of this process will help on the designing of new materials with improved optical and electrical properties. In this sense, computational simulations can deliver important information at atomic scale. Within a kinetic Monte Carlo scheme, diffusion rates are crucial to obtain accurate predictions; however, experimental values of this dynamic property for different oligomers are very scarce among literature. In this study, the diffusion coefficient (D) of thiophene oligomers (1Th–6Th) has been calculated using molecular dynamics simulations coupled with the Einstein expression. Results are in the order of experimental values, demonstrating that this methodology is a fast and reliable alternative to calculate diffusion coefficients with low computational costs.     

A niosomal bilayer of sorbitan monostearate in complex with flavones: a molecular dynamics simulation study

Bilayers prepared from sorbitan fatty acid esters (Span) have been frequently used for delivery of drugs including flavonoids. We applied molecular dynamics simulation to characterize the structure of a sorbitan monostearate (Span 60) bilayer in complex with three representative flavones, a subclass of flavonoids. At a low concentration, unsubstituted flavone, the most hydrophobic member, was able to flip over and cross the bilayer with a large diffusion coefficient. At a high concentration, it was accumulated at the bilayer center resulting in a phase separation. The leaflets of the bilayer were pushed in the opposite directions increasing the membrane thickness. Order parameter of the stearate chain of Span 60 was not affected significantly by unsubstituted flavone. In contrast, chrysin with hydroxylated ring A was lined up with the acyl chains of Span 60 with its hydroxyl group facing the membrane surface. Neither flipping nor transbilayer movement were allowed. Diffusion coefficient was only 15–25% of that of unsubstituted flavone and order parameter decreased with the concentration of chrysin. Luteolin, the most hydroxylated member, interacted mainly with the headgroup of Span 60 and assumed many different orientations without crossing the bilayer. Unlike chrysin and unsubstituted flavone the bilayer integrity was disrupted at 50?mol% luteolin. These behaviors and structures of flavones in a Span 60 bilayer can be accounted for by their hydrophobicity and sites of hydroxylation.     

Force-induced unfolding of human telomeric G-quadruplex: A steered molecular dynamics simulation study

We study the unfolding of a parallel G-quadruplex from human telomeric DNA by mechanical stretching using steered molecular dynamics (MD) simulation. We find that the force curves and unfolding processes strongly depend on the pulling sites. With pulling sites located on the sugar-phosphate backbone, the force-extension curve shows a single peak and the unfolding proceeds sequentially. Pulling sites located on the terminal nucleobases lead to a force-extension curve with two peaks and the unfolding is more cooperative. Simulations of the refolding of partially unfolded quadruplexes show very different behavior for the two different pulling modalities. In particular, starting from an unfolded state prepared by nucleobase pulling leads to a long-lived intermediate state whose existence is also corroborated by the free energy profile computed with the Jarzynski equation. Based on this observation, we propose a novel folding pathway for parallel G-quadruplexes with the human telomere sequence.     

Effect of gold nanoparticle on stability of the DNA molecule: A study of molecular dynamics simulation

An understanding of the mechanism of DNA interactions with gold nanoparticles is useful in today medicine applications. We have performed a molecular dynamics simulation on a B-DNA duplex (CCTCAGGCCTCC) in the vicinity of a gold nanoparticle with a truncated octahedron structure composed of 201 gold atoms (diameter ～1.8 nm) to investigate gold nanoparticle (GNP) effects on the stability of DNA. During simulation, the nanoparticle is closed to DNA and phosphate groups direct the particles into the major grooves of the DNA molecule. Because of peeling and untwisting states that are occur at end of DNA, the nucleotide base lies flat on the surface of GNP. The configuration entropy is estimated using the covariance matrix of atom-positional fluctuations for different bases. The results show that when a gold nanoparticle has interaction with DNA, entropy increases. The results of conformational energy and the hydrogen bond numbers for DNA indicated that DNA becomes unstable in the vicinity of a gold nanoparticle. The radial distribution function was calculated for water hydrogen–phosphate oxygen pairs. Almost for all nucleotide, the presence of a nanoparticle around DNA caused water molecules to be released from the DNA duplex and cations were close to the DNA.     

Interfacial water on crystalline silica: a comparative molecular dynamics simulation study

Understanding the properties of interfacial water at solid–liquid interfaces is important in a wide range of applications. Molecular dynamics is becoming a widespread tool for this purpose. Unfortunately, however, the results of such studies are known to strongly depend on the selection of force fields. It is, therefore, of interest to assess the extent by which the implemented force fields can affect the predicted properties of interfacial water. Two silica surfaces, with low and high surface hydroxyl density, respectively, were simulated implementing four force fields. These force fields yield different orientation and flexibility of surface hydrogen atoms, and also different interaction potentials with water molecules. The properties for interfacial water were quantified by calculating contact angles, atomic density profiles, surface density distributions, hydrogen bond density profiles and residence times for water near the solid substrates. We found that at low surface density of hydroxyl groups, the force field strongly affects the predicted contact angle, while at high density of hydroxyl groups, water wets all surfaces considered. From a molecular-level point of view, our results show that the position and intensity of peaks observed from oxygen and hydrogen atomic density profiles are quite different when different force fields are implemented, even when the simulated contact angles are similar. Particularly, the surfaces simulated by the CLAYFF force field appear to attract water more strongly than those simulated by the Bródka and Zerda force field. It was found that the surface density distributions for water strongly depend on the orientation of surface hydrogen atoms. In all cases, we found an elevated number of hydrogen bonds formed between interfacial water molecules. The hydrogen bond density profile does not depend strongly on the force field implemented to simulate the substrate, suggesting that interfacial water assumes the necessary orientation to maximise the number of water–water hydrogen bonds irrespectively of surface properties. Conversely, the residence time for water molecules near the interface strongly depends on the force field and on the flexibility of surface hydroxyl groups. Specifically, water molecules reside for longer times at contact with rigid substrates with high density of hydroxyl groups. These results should be considered when comparisons between simulated and experimental data are attempted.     

In silico approaches towards understanding CALB using molecular dynamics simulation and docking

Candida antarctica lipase B (CALB), a serine protease, is involved in the hydrolysis of substrates at the aqueous lipid interface. There is a significant role played by the helices in serine proteases including acting as a flap covering the active site region. The α5 and α10 helices in the path to the active site of CALB appear to play an important role in the region. This study investigates these helices by mutational studies, docking and molecular dynamics simulations. The mutations were selected based on their proximity to the active site and their presence at the α10-helix in the path of the active site. Molecular dynamics studies reveal the flexibility, stability and hydrogen bonding ability of the α5 helix. The radius of gyration (R _g) clearly showed the compactness of the structure. Docking studies show the changes occurring at the protein's binding site before and after 15 ns of simulation. Results from the study demonstrate the importance of the two helices α5 and α10 in the stability of CALB.     

Complex disruption effect of natural polyphenols on Bcl-2-Bax: molecular dynamics simulation and essential dynamics study

Apoptosis (programmed cell death) is a process by which cells died after completing physiological function or after a severe genetic damage. Apoptosis is mainly regulated by the Bcl-2 family of proteins. Anti apoptotic protein Bcl-2 prevents the Bax activation/oligomerization to form heterodimer which is responsible for release of the cytochrome c from mitochondria to the cytosol in response to death signal. Quercetin and taxifolin (natural polyphenols) efficiently bound to hydrophobic groove of Bcl-2 and altered the structure by inducing conformational changes. Taxifolin was found more efficient when compared to quercetin in terms of interaction energy and collapse of hydrophobic groove. Taxifolin and quercetin were found to dissociate the Bcl-2-Bax complex during 12?ns MD simulation. The effect of taxifolin and quercetin was, further validated by the MD simulation of ligand-unbound Bcl-2-Bax which showed stability during the simulation. Obatoclax (an inhibitor of Bcl-2) had no significant dissociation effect on Bcl-2-Bax during simulation which favored the previous experimental results and disruption effect of taxifolin and quercetin.     

Bovine pancreatic trypsin inhibitor crystals with different morphologies: a molecular dynamics simulation study

A molecular dynamics simulation study is reported for three polymorphic protein crystals (4PTI, 5PTI and 6PTI) of bovine pancreatic trypsin inhibitor (BPTI). The simulated lattice constants are in good agreement with experimental data, indicating the reliability of force field used. The fluctuation patterns of peptide chains in the three crystals are similar, and the protein structures are fairly well maintained during simulation. We observe that water forms a pronounced hydration layer near the protein surface. The diffusion coefficients of water in the three crystals are smaller than in bulk phase, and thus, the activation energies are higher. The porosity, fluctuation of peptide chains and solvent-accessible surface area as well as the diffusion coefficients of water and counterion in 5PTI are the largest among the three crystals. The diffusion of water and counterion is anisotropic, and the degree of anisotropy increases in the order of 4PTI < 5PTI < 6PTI. Despite a slight difference, the structural and diffusion properties in the three BPTI crystals are generally close. This simulation study reveals that crystal polymorphism does not significantly affect microscopic properties in the BPTI crystals with different morphologies.     

New insights into the molecular mechanism of methanol-induced inactivation of Thermomyces lanuginosus lipase: a molecular dynamics simulation study

Methanol intolerance of lipase is a major limitation in lipase-catalysed methanolysis reactions. In this study, to understand the molecular mechanism of methanol-induced inactivation of lipases, we performed molecular dynamics (MD) simulations of Thermomyces lanuginosus lipase (TLL) in water and methanol and compared the observed structural and dynamic properties. The solvent accessibility analysis showed that in methanol, polar residues tended to be buried away from the solvent while non-polar residues tended to be more solvent-exposed in comparison to those in water. Moreover, we observed that in methanol, the van der Waals packing of the core residues in two hydrophobic regions of TLL became weak. Additionally, the catalytically relevant hydrogen bond between Asp²⁰¹ OD2 and His²⁵⁸ ND1 in the active site was broken when enzyme was solvated in methanol. This may affect the stability of the tetrahedral intermediates in the catalytic cycle of TLL. Furthermore, compared to in water, some enzyme surface residues displayed enhanced movement in methanol with higher Cα root-mean-square atomic positional fluctuation values. One of such methanol-affecting surface residues (Ile²⁴¹) was chosen for mutation, and MD simulation of the I241E mutant in methanol was conducted. The structural analysis of the mutant showed that replacing a non-polar surface residue with an acidic one at position 241 contributed to the stabilisation of enzyme structure in methanol. Ultimately, these results, while providing molecular-level insights into the destabilising effect of methanol on TLL, highlight the importance of surface residue redesign to improve the stability of lipases in methanol environments.     

New potential inhibitors of mTOR: a computational investigation integrating molecular docking,virtual screening and molecular dynamics simulation

The mTOR (mammalian or mechanistic Target Of Rapamycin), a complex metabolic pathway that involves multiple steps and regulators, is a major human metabolic pathway responsible for cell growth control in response to multiple factors and that is dysregulated in various types of cancer. The classical inhibition of the mTOR pathway is performed by rapamycin and its analogs (rapalogs). Considering that rapamycin binds to an allosteric site and performs a crucial role in the inhibition of the mTOR complex without causing the deleterious side effects common to ATP-competitive inhibitors, we employ ligand-based drug design strategies, such as virtual screening methodology, computational determination of ADME/Tox properties of selected molecules, and molecular dynamics in order to select molecules with the potential to become non-ATP-competitive inhibitors of the mTOR enzymatic complex. Our findings suggest five novel potential mTOR inhibitors, with similar or better properties than the classic inhibitor complex, rapamycin.     

Effects of cholesterol concentration on the interaction of cytarabine with lipid membranes: a molecular dynamics simulation study

Liposomal cytarabine, DepoCyt, is a chemotherapy agent which is used in cancer treatment. This form of cytarabine has more efficacy and fewer side effects relative to the other forms. Since DepoCyt contains the cytarabine encapsulated within phosphatidylcholine and the sterol molecules, we modeled dioleoylphosphatidylcholine (DOPC)/cholesterol bilayer membrane as a carrier for cytarabine to study drug–bilayer interactions. For this purpose, we performed a series of united-atom molecular dynamics (MD) simulations for 25?ns to investigate the interactions between cytarabine and cholesterol-containing DOPC lipid bilayers. Only the uncharged form of cytarabine molecule was investigated. In this study, different levels of the cholesterol content (0, 20, and 40%) were used. MD simulations allowed us to determine dynamical and structural properties of the bilayer membrane and to estimate the preferred location and orientation of the cytarabine molecule inside the bilayer membrane. Properties such as membrane thickness, area per lipid, diffusion coefficient, mass density, bilayer packing, order parameters, and intermolecular interactions were examined. The results show that by increasing the cholesterol concentration in the lipid bilayers, the bilayer thickness increases and area per lipid decreases. Moreover, in accordance with the experiments, our calculations show that cholesterol molecules have ordering effect on the hydrocarbon acyl chains. Furthermore, the cytarabine molecule preferentially occupies the polar region of the lipid head groups to form specific interactions (hydrogen bonds). Our results fully support the experimental data. Our finding about drug–bilayer interaction is crucial for the liposomal drug design.     

Atomic insight into designed carbamate-based derivatives as acetylcholine esterase (AChE) inhibitors: a computational study by multiple molecular docking and molecular dynamics simulation

Over 100 variants have been designed and studied, using multiple docking methods such as Autodock Vina, ArgusLab, Molegro Virtual Docker, and Hex-Cuda, to study the effect of alteration in the structure of carbamate-based acetylcholyne esterase (AChE) inhibitors. Sixteen selected systems were then subjected to 14 ns molecular dynamics (MD) simulations. Results from all the docking methods are in agreement. Variants that involved biphenyl substituents possess the most negative binding energies in the ?37.64 to ?39.31 kJ mol^?1 range due to their π–π interactions with AChE aromatic residues. The root mean square deviation values showed that all of these components achieved equilibration after 6 ns. Gyration radius (R_g) and solvent accessibility surface area were calculated to further investigate the AChE conformational changes in the presence of these components. MD simulation results suggested that these components might interact with AChE, possibly with no major changes in AChE secondary and tertiary structures.     

Understanding interaction and dynamics of water molecules in the epoxy via molecular dynamics simulation

The robust structural integrity of the epoxy plays an important role in ensuring the long-term service life of its applications, which is affected by the absorbed moisture. In order to understand the mechanism of the moisture effect, the knowledge of the interaction and dynamics of the water molecules inside the epoxy is of great interest. Molecular dynamics simulation is used in this work to investigate the structure and bonding behaviour of the water molecules in the highly cross-linked epoxy network. When the moisture concentration is low, the water molecules are well dispersed in the cross-linked structure and located in the vicinity of the epoxy functional groups, which predominantly form the hydrogen bond (H-bond) with the epoxy network, resulting in the low water mobility in the epoxy. At the high concentration, the water favourably forms the large cluster due to the predominant water–water H-bond interaction, and the water molecules diffuse primarily inside the cluster, which leads to the high water mobility and the accelerated H-bond dynamics. The variation of the bonding behaviour and dynamics of the water molecules reported here could be exploited to understand the material change and predict the long-term performance of the epoxy-based products during the intended service life.     

The mechanism for the complexation and dissociation between siRNA and PMAL: a molecular dynamics simulation study based on a coarse-grained model

Abstract

The capacity of silencing genes makes small interfering RNA (siRNA) becomes potential candidates for curing many fatal diseases. Due to the low stability and delivery efficiency of siRNA, the design of amphiphilic carrier for siRNA delivery is vital for the practical gene therapy. In the present work, we explored how the complexation and dissociation of siRNA with poly (maleic anhydride-alt-1-decene) substituted with 3-(dimethylamino) propylamine (PMAL), which is a recent synthesised amphiphilic polymer and can be used in delivery of siRNAs and proteins, using traditional molecular dynamics simulations, together with steered molecular dynamics simulations. It was shown that the complexation of siRNA with PMALs can spontaneously occur, no matter what unit number of PMAL is. PMALs of different unit numbers form micelle-like structures and interact with siRNA surface. With the increase of unit number, PMAL becomes more flexible and interacts with siRNA from attachment to entanglement. The dissociation of PMAL from siRNA is an energy-consuming process. The free energy difference increases with the unit number of PMAL. The free energy for dissociation involves both the stretch of PMAL and the separation of PMAL from siRNA. Therefore, an optimal unit number of PMAL is critical for the delivery efficiency of siRNA when PMAL is used as carrier. In present work, when the radius of gyration of PMAL approaches to that of siRNA, PMAL gives a favoured both complexation and dissociation between siRNA and PMAL. Finally, we propose the mechanism of complexation and dissociation of PMAL with siRNA. The above simulation established a molecular insight of the interaction between siRNA and PMAL and was helpful for the design and applications of new PMAL-based polymers as siRNA delivery carriers.     

Discovery of promising FtsZ inhibitors by E-pharmacophore, 3D-QSAR,molecular docking study,and molecular dynamics simulation

Filamentous temperature-sensitive protein Z (FtsZ), playing a key role in bacterial cell division, is regarded as a promising target for the design of antimicrobial agent. This study is looking for potential high-efficiency FtsZ inhibitors. Ligand-based pharmacophore and E-pharmacophore, virtual screening and molecular docking were used to detect promising FtsZ inhibitors, and molecular dynamics simulation was used to study the stability of protein-ligand complexes in this paper. Sixty-three inhibitors from published literatures with pIC₅₀ ranging from 2.483 to 5.678 were collected to develop ligand-based pharmacophore model. 4DXD bound with 9PC was selected to develop the E-pharmacophore model. The pharmacophore models validated by test set method and decoy set were employed for virtual screening to exclude inactive compounds against ZINC database. After molecular docking, ADME analysis, IFD docking and MM-GBSA, 8 hits were identified as potent FtsZ inhibitors. A 50?ns molecular dynamics simulation was implemented on the compounds to assess the stability between potent inhibitors and FtsZ. The results indicated that the candidate compounds had a high docking score and were strongly combined with FtsZ by forming hydrogen bonding interactions with key amino acid residues, and van der Waals forces and hydrophobic interactions had significant contribution to the stability of the binding. Molecular dynamics simulation results showed that the protein-ligand compounds performed well in both the stability and flexibility of the simulation process.     

Conformational dynamics of the mitochondrial ADP/ATP carrier: a simulation study

The mitochondrial ADP/ATP carrier is a six helix bundle membrane transport protein, which couples the exit of ATP from the mitochondrial matrix to the entry of ADP. Extended (4×20 ns) molecular dynamics simulations of the carrier, in the presence and absence of bound inhibitor (carboxyatractyloside), have been used to explore the conformational dynamics of the protein in a lipid bilayer environment, in the presence and absence of the carboxyatractyloside inhibitor. The dynamic flexibility (measured as conformational drift and fluctuations) of the protein is reduced in the presence of bound inhibitor. Proline residues in transmembrane helices H1, H3 and H5 appear to form dynamic hinges. Fluctuations in inter-helix salt bridges are also observed over the time course of the simulations. Inhibitor-protein and lipid-protein interactions have been characterised in some detail. Overall, the simulations support a transport mechanism in which flexibility about the proline hinges enables a transition between a ‘closed’ and an ‘open’ pore-like state of the carrier protein.     

Structural stability of myoglobin and glycomyoglobin: a comparative molecular dynamics simulation study

Glycoproteins are formed as the result of enzymatic glycosylation or chemical glycation in the body, and produced in vitro in industrial processes. The covalently attached carbohydrate molecule(s) confer new properties to the protein, including modified stability. In the present study, the structural stability of a glycoprotein form of myoglobin, bearing a glucose unit in the N-terminus, has been compared with its native form by the use of molecular dynamics simulation. Both structures were subjected to temperatures of 300 and 500 K in an aqueous environment for 10 ns. Changes in secondary structures and RMSD were then assessed. An overall higher stability was detected for glycomyoglobin, for which the most stable segments/residues were highlighted and compared with the native form. The simple addition of a covalently bound glucose is suggested to exert its stabilizing effect via increased contacts with surrounding water molecules, as well as a different pattern of interactions with neighbor residues.

Electronic supplementary material

The online version of this article (doi:10.1007/s10867-015-9383-2) contains supplementary material, which is available to authorized users.     

Impact of high-pressure processing on diffusion in polyethylene based on molecular dynamics simulation

Abstract

The impact of high-pressure processing (HPP) on dissusion of antioxidant butylated hydroxytoluene (BHT) in polyethylene (PE) was discussed via the molecular dynamics method. Furthermore, the glass transition temperatures (T_g), the accessible free volumes of PE and the diffusion coefficients of BHT in PE at different HPP treatments were calculated, and the diffusion trajectories of the BHT molecules in PE were also presented. Finally, the diffusion mechanism of BHT in PE under HPP was analyzed based on the aforementioned simulation results. The results show that the T_g of PE increases under high pressure while the fractional free volume (FFV) reduces, and the diffusion coefficient decreases with the pressure on the rise. The diffusion trajectories of BHT in PE under HPP indicate that the BHT molecules are trapped and slowly wriggle in a narrow path among PE molecular chains due to the extreme high pressure. However, the high temperature has an opposite effect on the diffusion behavior of BHT in PE compared with high pressure. As the temperature rises, the FFV of PE and the diffusion coefficient of BHT in PE are elevated. This study is helpful to the research of high-pressure food safety and packaging migration.     

Effect of chirality,length and diameter of carbon nanotubes on the adsorption of 20 amino acids: a molecular dynamics simulation study

We carried out molecular dynamics simulations to study the adsorption of all the 20 amino acids (AAs; aromatic, polar and non-polar) on the surface of chiral, zigzag and armchair single-walled carbon nanotubes. The adsorption was occurring in all systems. In the aromatic AAs, the π–π stacking and the semi-hydrogen bond formation cause a strong interaction with the carbon nanotubes (CNTs). We also investigated the chirality, length and diameter dependencies on adsorption energies. We found that all AAs have more tendency to adsorption on the chiral and zigzag CNTs over the armchair. The results show that increasing both the diameter and the length causes the enhancement of the adsorption energy. But, the effect of the length is more than of the diameter. For example, the adsorption energy of Trp on the surface of CNT (4,1), with 2 nm length, is 20.4 kcal/mol. When the length of CNT becomes twice, the adsorption energy increases by 24 ± 0.3%. But by doubling the diameter, the adsorption energy increased only by 9.8 ± 0.25%.