Conformational and dynamics simulation study of antimicrobial peptide hedistin—heterogeneity of its helix–turn–helix motif

Hedistin is an antimicrobial peptide isolated from the coelomocytes of Nereis diversicolor, possessing activity against a large spectrum of bacteria including the methicillin resistant Staphylococcus aureus and Vibrio alginolyticus. The three-dimensional structure of hedistin in both aqueous solution and deuterated dodecylphosphocholine (DPC) micelles was examined using circular dichroism (CD) and nuclear magnetic resonance (NMR) techniques. And, the early events of the antibacterial process of hedistin were simulated using palmitoyl-oleoyl-phophatidylcholine (POPC) lipid bilayers and molecular dynamics (MD) simulation methods. Hedistin lacks secondary structure in aqueous solution, however, in DPC micelles, it features with a heterogeneous helix–turn–helix moiety and exhibits obvious amphipathic nature. The turn region (residues Val9–Thr12) in the moiety is a four-residue hinge, lying in between the first N-terminal α-helix (residues Leu5–Lys8) and the second α-helix (residues Val13–Ala17) regions and causing an ～ 120° angle between the axes of the two helices. The segmental and nonlinear nature of hedistin structure is referred to as the heterogeneity of its helix–turn–helix motif which was found to be corresponding to a kind of discrete dynamics behavior, herein coined as its dynamical heterogeneity, at the early stage (0–50 ns) of the MD simulations. That is, the first helix segment, prior to (at 310 K) or following (at 363 K) the second helix, binds to the lipid head-group region and subsequently permeates into the hydrophobic lipid tail region, and the hinge is the last portion entering the lipid environment. This result implies that hedistin may adopt a “carpet” model action when disrupting bacterial membrane.     

Molecular dynamics simulation of the melting processes of core–shell and pure nanoparticles

Core–shell nanoparticles are nanosized particles that consist of a core and a shell, constructed from different metallic elements. Core–shell nanoparticles have received extensive attention, owing to their various potential applications such as paints, optical films and catalysts. Herein, we investigate the melting behaviours of different core–shell nanoparticles under continuous heating using molecular dynamics simulation. Different metallic elements were examined as core–shell and pure nanoparticles. Five different processes were observed during the melting of core–shell nanoparticles. In contrast, only one process was identified during the melting of pure nanoparticles. These processes were influenced by the nanoparticle size, shell thickness and differences between the lattice constants and melting point temperatures of the metallic elements. Our simulation provides microscopic insights into the melting behaviours of existing and proposed core–shell nanoparticles that would be highly beneficial towards the fabrication of materials with different chemical coatings.     

Molecular dynamics simulation of the gGAPDH–NAD+ complex from Trypanosoma cruzi

A 50-ns molecular dynamics simulation has been used to study the homotetramer of the enzyme glycosomal glyceraldehyde 3-phosphate dehydrogenase (gGAPDH) complexes, from Trypanosoma cruzi, with nicotinamide adenine dinucleotide (NAD⁺) cofactors in aqueous solution. The root mean square deviation indicates that the overall structure of the homotetramer does not undergo significant change. The largest structural change observed was in the NAD⁺ binding domain of subunit (chain) D; as a consequence, the NAD⁺ cofactor was dislocated from its initial position. However, the other subunits were not affected, suggesting that the gGAPDH enzyme exhibits non-cooperative behaviour. Our simulation estimates that the NAD⁺ binding domain rotates about 4.8° relative to the catalytic domain in the apo–holo form transition. The hydrogen bond analysis reveals that the residues R12, I13, D38 and M39 are essential for gGAPDH–NAD⁺ interaction. Furthermore, two promising cavities to be explored in drug design were found: one formed by residues I13, R12, T197, T199, E336 and Y339, and the other by residues C166, H194, R249, I13, R12, T197, T199, E336 and Y339. The results presented in this paper offer new insight into the search for inhibitors of the gGAPDH enzyme of T. cruzi protozoan.     

Molecular dynamics simulation of the graphene–water interface: comparing water models

In this work, different water models (i.e. SPC/E, TIP3P, TIP4P/2005, TIP5P, SPC/Fw, TIP4P/2005f and SWM4_DP) are implemented to simulate water on neutral, negatively charged and positively charged graphene. In all cases ambient conditions are considered. Structural and dynamical properties for water are calculated to quantify the differences among various water models. The results show that SPC/E, TIP4P/2005, SPC/Fw, TIP4P/2005f and SWM4_DP water models yield a similar structure for interfacial water on graphene, whether it is neutral, negatively charged or positively charged. TIP5P is the model whose predictions for the structure of the interface deviate the most from those of the other models. Although qualitatively the results are for the most part similar, a large quantitative variation is observed among the dynamical properties predicted when various water models are implemented. Although experimental data are not available to discriminate the most/least accurate of the model predictions, our results could be useful for comparing results for interfacial water obtained implementing different models. Such critical comparison will benefit practical applications such as the development of energy-storage and water-desalination devices (e.g. electric double-layer capacitors), among others.     

The transmembrane domain of caveolin-1 exhibits a helix–break–helix structure

Caveolin is an integral membrane protein that is found in high abundance in caveolae. Both the N- and C- termini lie on the same side of the membrane, and the transmembrane domain has been postulated to form an unusual intra-membrane horseshoe configuration. To probe the structure of the transmembrane domain, we have prepared a construct of caveolin-1 that encompasses residues 96–136 (the entire intact transmembrane domain). Caveolin-1(96–136) was over-expressed and isotopically labeled in E. coli, purified to homogeneity, and incorporated into lyso-myristoylphosphatidylglycerol micelles. Circular dichroism and NMR spectroscopy reveal that the transmembrane domain of caveolin-1 is primarily α-helical (57–65%). Furthermore, chemical shift indexing reveals that the transmembrane domain has a helix–break–helix structure which could be critical for the formation of the intra-membrane horseshoe conformation predicted for caveolin-1. The break in the helix spans residues 108 to 110, and alanine scanning mutagenesis was carried out to probe the structural significance of these residues. Our results indicate that mutation of glycine 108 to alanine does not disrupt the structure, but mutation of isoleucine 109 and proline 110 to alanine dramatically alters the helix–break–helix structure. To explore the structural determinants further, additional mutagenesis was performed. Glycine 108 can be substituted with other small side chain amino acids (i.e. alanine), leucine 109 can be substituted with other β-branched amino acids (i.e. valine), and proline 110 cannot be substituted without disrupting the helix–break–helix structure.     

Molecular dynamics simulation of diffusion of hydrogen in binary hydrogen–tetrahydrofuran hydrate

The binary structure II hydrogen–tetrahydrofuran (THF) hydrate was studied with molecular dynamics simulation. The simulations were carried out at 300, 310 K and 10.1 MPa, and with various contents of hydrogen and THF. The migrations of hydrogen molecules from cage to cage were observed. The migration process of hydrogen was also analysed, and the diffusion coefficients of hydrogen in the hydrate were calculated. The calculated diffusion coefficients qualitatively agreed with the experimental data. Double and quintet occupancies of hydrogen molecules were observed in the small and large cages, respectively, without changing the hydrate structure.     

Fluid–structure interaction simulation of the brain–skull interface for acute subdural haematoma prediction

Traumatic brain injury is a leading cause of disability and mortality. Finite element-based head models are promising tools for enhanced head injury prediction, mitigation and prevention. The reliability of such models depends heavily on adequate representation of the brain–skull interaction. Nevertheless, the brain–skull interface has been largely simplified in previous three-dimensional head models without accounting for the fluid behaviour of the cerebrospinal fluid (CSF) and its mechanical interaction with the brain and skull. In this study, the brain–skull interface in a previously developed head model is modified as a fluid–structure interaction (FSI) approach, in which the CSF is treated on a moving mesh using an arbitrary Lagrangian–Eulerian multi-material formulation and the brain on a deformable mesh using a Lagrangian formulation. The modified model is validated against brain–skull relative displacement and intracranial pressure responses and subsequently imposed to an experimentally determined loading known to cause acute subdural haematoma (ASDH). Compared to the original model, the modified model achieves an improved validation performance in terms of brain–skull relative motion and is able to predict the occurrence of ASDH more accurately, indicating the superiority of the FSI approach for brain–skull interface modelling. The introduction of the FSI approach to represent the fluid behaviour of the CSF and its interaction with the brain and skull is crucial for more accurate head injury predictions.

     

Molecular dynamics simulation of Axillaridine–A: A potent natural cholinesterase inhibitor

Molecular Dynamics (MD) simulations were carried out for human acetylcholinesterase (hAChE) and its complex with Axillaridine–A, in order to dynamically explore the active site of the protein and the behaviour of the ligand at the peripheral binding site. Simulation of the enzyme alone showed that the active site of AChE is located at the bottom of a deep and narrow cavity whose surface is lined with rings of aromatic residues while Tyr72 is almost perpendicular to the Trp286, which is responsible for stable π -π interactions. The complexation of AChE with Axillaridine-A, results in the reduction of gorge size due to interaction between the ligand and the active site residues. The gorge size was determined by the distance between the center of mass of Glu81 and Trp286. As far as the geometry of the active site is concerned, the presence of ligand in the active site alters its specific conformation, as revealed by stable hydrogen bondings established between amino acids. With the increasing interaction between ligand and the active amino acids, size of the active site of the complex decreases with respect to time. Axillaridine-A, forms stable π -π interactions with the aromatic ring of Tyr124 that results in inhibition of catalytic activity of the enzyme. This π -π interaction keeps the substrate stable at the edge of the catalytic gorge by inhibiting its catalytic activity. The MD results clearly provide an explanation for the binding pattern of bulky steroidal alkaloids at the active site of AChE.     

Molecular dynamics simulation of helium–argon gas mixture under various wall conditions

As computational capabilities increase, molecular dynamics (MD) simulations become important tools of simulating reality. These simulations are especially useful for compressible gas mixture problems. In this study, binary diffusion of helium and argon was examined using a hard-sphere MD simulation method. For the sake of computational speed, low spacing ratios were chosen. Binary mass diffusion of gases in two equally sized halves of a box was simulated for identical initial kinetic energies and number densities. It has been noted that a purely mass diffusion mechanism of different gases is not physically possible. The resultant gas mixtures of several diffusion simulations were used as initial conditions for combined heat transfer – Couette flow, and heating and cooling experiments. The results showed the interesting behaviour of the mixture, which was subjected to various wall conditions. Energy of heavier molecules is found to be more sensitive to the wall velocities and less sensitive to the wall temperatures than lighter molecules. Diffusion, heat transfer, viscosity and heat capacity coefficients are deduced as well.     

Molecular dynamics simulation of diffusion behaviour of gas molecules within oil–paper insulation system

The diffusion behaviour of hydrogen, carbon monoxide, carbon dioxide, methane, acetylene, ethylene and ethane in oil and paper medium was examined using molecular dynamics to reveal the diffusion mechanism of gas molecules in transformer oil–paper insulation system at the microscopic level. These compounds are commonly used in the dissolved gas analysis of power transformers and produced during the ageing process of oil–paper composite insulating material. Two groups of models were constructed using molecular dynamics simulation software to simulate the diffusion behaviour of the aforementioned seven types of small gas molecules in oil and paper. The diffusion coefficients, displacement features, free volume characteristics and interaction energies of the gas molecules were investigated. In particular, the diffusion micro-mechanism of the gas molecules was observed. The differences in diffusion features among the gas molecules were discussed, and the factors influencing the diffusion of the gas molecules were compared. Simulation results indicate that the diffusion coefficients of gas molecules in cellulose is an order of magnitude lower than that in oil, and the diffusion coefficients of these gas molecules in the two types of insulation media have different orders. Free volume of gas molecules is the main factor that influences the diffusion behaviour in oil, whereas intermolecular interaction is the main influencing factor of diffusion behaviour in cellulose.     

Molecular dynamics simulation of the reciprocal fused LiF–KBr mixture: local structure and self-diffusion coefficients

This paper describes the molecular dynamics simulation of the reciprocal fused LiF–KBr mixture, which is located above the critical mixing point, in the temperature range 1280–1450 K. The first coordination sphere is found to form as follows: a smaller ion is formed around a smaller counter-ion, and a larger ion is formed around a larger counter-ion. The calculated concentration dependence of the self-diffusion coefficients and the radial distribution functions of all ion pairs indicate that the degree of association of the Li–F pair increases as the lithium fluoride fraction in the mixture decreases.     

De novo design of transmembrane helix–helix interactions and measurement of stability in a biological membrane

Membrane proteins regulate a large number of cellular functions, and have great potential as tools for manipulation of biological systems. Developing these tools requires a robust and quantitative understanding of membrane protein folding and interactions within the bilayer. With this in mind, we have designed a series of proteins to probe the net thermodynamic contribution of well-known sequence motifs to transmembrane helix-helix association in a biological membrane. The proteins were designed from first principles (de novo) using current knowledge about membrane insertion and stabilizing interaction motifs. A simple poly-Leu “scaffold” was decorated with individual helix interaction motifs (G-XXX-G, polar residues, heptad repeat) to create transmembrane helix–helix interactions of increasing strength. The GALLEX assay, an in vivo assay for measurement of transmembrane helix self-association, was combined with computational methods to characterize the relative strength and mode of interaction for each sequence. In addition, the apparent free energy contribution (ΔΔG^app) of each motif to transmembrane helix self-association was measured in a biological membrane, results that are the first of their kind for these de novo designed sequences, and suggest that the free energy barrier to overcoming weak association is quite small (< 1.4 kcal mol^− 1) in a natural membrane. By quantifying and rationalizing the contribution of key motifs to transmembrane helix association, our work offers a route to direct the design of novel sequences for use in biotechnology or synthetic biology (e.g. molecular switches) and to predict the effects of sequence modification in known transmembrane domains (for control of cellular processes).     

Molecular dynamics simulation on the effect of transition metal binding to the N-terminal fragment of amyloid-β

Abstract

We report molecular dynamics simulations of three possible adducts of Fe(II) to the N-terminal 1–16 fragments of the amyloid-β peptide, along with analogous simulations of Cu(II) and Zn(II) adducts. We find that multiple simulations from different starting points reach pseudo-equilibration within 100–300?ns, leading to over 900?ns of equilibrated trajectory data for each system. The specifics of the coordination modes for Fe(II) have only a weak effect on peptide secondary and tertiary structures, and we therefore compare one of these with analogous models of Cu(II) and Zn(II) complexes. All share broadly similar structural features, with mixture of coil, turn and bend in the N-terminal region and helical structure for residues 11–16. Within this overall pattern, subtle effects due to changes in metal are evident: Fe(II) complexes are more compact and are more likely to occupy bridge and ribbon regions of Ramachandran maps, while Cu(II) coordination leads to greater occupancy of the poly-proline region. Analysis of representative clusters in terms of molecular mechanics energy and atoms-in-molecules properties indicates similarity of four-coordinate Cu and Zn complexes, compared to five-coordinate Fe complex that exhibits lower stability and weaker metal–ligand bonding.

Communicated by Ramaswamy H. Sarma     

Molecular dynamics simulation integrating the inhibition kinetics of hydroxysafflor yellow A on α-glucosidase

Inhibition of α-glucosidase has attracted the attention of researchers due to its connection to type-2 diabetes. Hydroxysafflor yellow A (HSYA) extracted from Carthamus tinctorius L. is a natural antioxidant used in traditional Chinese medicine. In this study, the effect of HSYA on α-glucosidase was evaluated using inhibitory kinetics based on the antioxidant properties of HSYA and by performing computational simulation integration methods. HSYA reversibly inhibited α-glucosidase in a competitive inhibition manner and the evaluated kinetic parameters were IC₅₀ = 1.1 ± 0.22 mM and K_i = 1.04 ± 0.23 mM, respectively. The results of spectrofluorimetry showed that the inner hydrophobic regions of α-glucosidase, which are mostly in the active site, were exposed to the surface with increasing HSYA concentrations, indicating that the inactivation of α-glucosidase by HSYA was accompanied by regional unfolding. The molecular dynamics simulations indicated that the four rings of HSYA interact with four residues such as G217, A278, H279, and G280 at the entrance of the active site. Our study provides insight into the inhibition of α-glucosidase and the accompanying structural changes by HSYA. Based on its α-glucosidase-inhibiting effect and its potential as a natural antioxidant, HSYA is a potential agent for treating α-glucosidase-associated type-2 diabetes.     

Molecular simulation of polymers with a SAFT-γ Mie approach

ABSTRACT

We review the group contribution Statistical Associating Fluid Theory with Mie interaction potentials (SAFT-γ Mie) approach for building coarse-grained models for molecular simulation of polymeric systems. In this top-down method, force field parameters for coarse-grained polymer models can be derived from thermodynamic information on constituent monomer units using the SAFT-γ Mie equation of state (EoS). This strategy can facilitate high-throughput computational screening of polymeric materials, with a corresponding states correlation expediting the force field fitting. Accurate and transferable non-bonded parameters linked to macroscopic thermodynamic data allow for calculation of properties beyond those obtainable from the EoS alone. To overcome limitations of SAFT-γ Mie regarding polymer chain stiffness and branching, hybrid top-down/bottom-up approaches have combined non-bonded parameters from SAFT-γ Mie with bond-stretching and angle-bending potentials from higher-resolution force fields. Our review critically evaluates the performance of recent SAFT-γ Mie polymer models, highlighting the strengths and weaknesses in the context of other equation of state and coarse-graining methods.     

Molecular dynamics approach to investigate the coupling of the hydrophilic–lipophilic balance with the configuration distribution function in biosurfactant-based emulsions

Emulsion stability has been characterized by macroscopic variables such as the hydrophilic–lipophilic balance, with the aim being to predict the surfactant properties of molecules. Nevertheless, this parameter does not take the topology of the molecule into account, as it only considers its lipophilic degree. On the other hand, the classical Derjaguin–Landau–Verwey–Overbeek approach (based on the continuum model), which has been widely utilized to evaluate the stabilities of colloids, polymers, and surfactants, takes some bulk macroscopic parameters such as the shear viscosity coefficient and the dielectric permittivity into account. In the work reported here, molecular dynamics simulations were used to elucidate the mechanism of layer formation and micellar structure for different combinations of valine–aspartic acid peptides in dodecane–water emulsions, as well as their associations with the hydrophilic–lipophilic balance. The peptide–dodecane radial distribution function showed that the first peak intensity was inversely correlated with the hydrophilic–lipophilic balance; moreover, the oscillatory structural forces became increasingly prominent when the hydrophilic–lipophilic balance was decreased. Our results seem to indicate that the radial distribution function could be utilized to evaluate the stabilities of emulsions of peptides via molecular simulations.     

IB–LBM simulation of the haemocyte dynamics in a stenotic capillary

To study the behaviour of a haemocyte when crossing a stenotic capillary, the immersed boundary–lattice Boltzmann method was used to establish a quantitative analysis model. The haemocyte was assumed to be spherical and to have an elastic cell membrane, which can be driven by blood flow to adopt a highly deformable character. In the stenotic capillary, the spherical blood cell was stressed both by the flow and the wall dimension, and the cell shape was forced to be stretched to cross the stenosis. Our simulation investigated the haemocyte crossing process in detail. The velocity and pressure were anatomised to obtain information on how blood flows through a capillary and to estimate the degree of cell damage caused by excessive pressure. Quantitative velocity analysis results demonstrated that a large haemocyte crossing a small stenosis would have a noticeable effect on blood flow, while quantitative pressure distribution analysis results indicated that the crossing process would produce a special pressure distribution in the cell interior and to some extent a sudden change between the cell interior and the surrounding plasma.     

Molecular simulation of solid–solid phase transitions

Applications of dl_poly to solid–solid phase transitions are reviewed, with particular attention to how details of the mechanisms of the transitions may be extracted from molecular dynamics simulations. Two examples in molecular crystals are discussed: the order–disorder transition of p-terphenyl initiated at around 200 K by the unlocking of ring flipping; and the rotator phases of n-alkanes with around 20 carbon atoms per chain, showing distinct molecular mechanisms in the dynamics just below the melting points of odd and even chains. Covalent-ionic materials are represented by an application to an aluminophophate molecular sieve, AlPO₄-5.     

Structural dynamics of a single-stranded RNA–helix junction using NMR

Many regulatory RNAs contain long single strands (ssRNA) that adjoin secondary structural elements. Here, we use NMR spectroscopy to study the dynamic properties of a 12-nucleotide (nt) ssRNA tail derived from the prequeuosine riboswitch linked to the 3′ end of a 48-nt hairpin. Analysis of chemical shifts, NOE connectivity, ¹³C spin relaxation, and residual dipolar coupling data suggests that the first two residues (A25 and U26) in the ssRNA tail stack onto the adjacent helix and assume an ordered conformation. The following U26-A27 step marks the beginning of an A₆-tract and forms an acute pivot point for substantial motions within the tail, which increase toward the terminal end. Despite substantial internal motions, the ssRNA tail adopts, on average, an A-form helical conformation that is coaxial with the helix. Our results reveal a surprising degree of structural and dynamic complexity at the ssRNA–helix junction, which involves a fine balance between order and disorder that may facilitate efficient pseudoknot formation on ligand recognition.     

Molecular dynamics simulation studies of the structural response of an isolated Aβ1–42 monomer localized in the vicinity of the hydrophilic TiO2 surface

We have probed the effect of a model hydrophilic surface, rutile TiO₂, on the full-length amyloid beta (Aβ_1–42) monomer using molecular dynamics simulations. The rutile surface brings about sharp changes in the peptide’s intrinsic behavior in a distance-dependent manner. The intrinsic collapse of the peptide is disrupted, while the β-sheet propensity is sharply enhanced with increased proximity to the surface. The results may have implications for Aβ self-assembly and fibrillogenesis on hydrophilic surfaces and should be taken into consideration in the design of novel nanomaterials for perturbing amyloidogenic behavior.