Efficient tuning of potential parameters for liquid–solid interactions

Spherical and cylindrical water droplets on silicon surface are studied to tune the silicon–oxygen interaction. We use molecular dynamics simulations to estimate the contact angle of two different shaped droplets. We found that the cylindrical droplets are independent of the line tension as their three phases curvature is equal zero. Additionally, we compare an analytical model, taking into account or not the Tolman length and we show that for spherical small size droplets, this length is important to be included, in contrast to cylindrical droplets in which the influence of the Tolman length is negligible. We demonstrate that the usual convenient way to exclude linear tension in the general case can give wrong results. Here, we consider cylindrical droplets, since their contact angle does not depend on the droplet size in the range of few to 10ths of nanometres. The droplets are stabilised due to the periodic boundary conditions. This allows us to propose a new parameterisation for nanoscale droplets, which is independent the size of the droplets or its shape, minimising at the same time the calculation procedure. With the proposed methodology, we can extract the epsilon parameter of the interaction potential between a liquid and a solid from the nanoscaled molecular simulation with only as input the macrosized experimental wetting angle for a given temperature.     

Antimicrobial peptides in toroidal and cylindrical pores

Antimicrobial peptides (AMPs) are small, usually cationic peptides, which permeabilize biological membranes. Their mechanism of action is still not well understood. Here we investigate the preference of alamethicin and melittin for pores of different shapes, using molecular dynamics (MD) simulations of the peptides in pre-formed toroidal and cylindrical pores. When an alamethicin hexamer is initially embedded in a cylindrical pore, at the end of the simulation the pore remains cylindrical or closes if glutamines in the N-termini are not located within the pore. On the other hand, when a melittin tetramer is embedded in toroidal pore or in a cylindrical pore, at the end of the simulation the pore is lined both with peptides and lipid headgroups, and, thus, can be classified as a toroidal pore. These observations agree with the prevailing views that alamethicin forms barrel-stave pores whereas melittin forms toroidal pores. Both alamethicin and melittin form amphiphilic helices in the presence of membranes, but their net charge differs; at pH ∼ 7, the net charge of alamethicin is − 1 whereas that of melittin is + 5. This gives rise to stronger electrostatic interactions of melittin with membranes than those of alamethicin. The melittin tetramer interacts more strongly with lipids in the toroidal pore than in the cylindrical one, due to more favorable electrostatic interactions.     

On the Application of Widom's Test Particle Method to Homogeneous and Inhomogeneous Fluids

Abstract

Chemical potentials of a homogeneous and an inhomogeneous Lennard-Jones fluid have been determined by molecular dynamics simulations on the vector computer CYBER 205 by applying essentially the fictitious test particle method of Widom. For the homogeneous fluid we find, contrary to the previous result of Guillot and Guissani, that the simulated chemical potential is independent of the particle number. The crucial point, however, is a sufficiently large cut-off radius in the evaluation of the Boltzmann factor. Comparing with our WCA-type perturbation theory, we get agreement in the chemical potentials within 0.1 kT up to the density n[sgrave]³ = 0.80 and a difference of 0.2 kT at n[sgrave]³ = 0.85. For the inhomogeneous case we consider a fluid in a cylindrical pore and integrate Widom's equation over a certain probe volume as suggested earlier by us. Chemical potentials are then calculated independently in five different probe volumes, which are cylindrical shells. The results agree well from the second to the fourth shell. Inaccuracies in the innermost cylinder can be easily explained by bad statistics. In the shell close to the wall the extremely high local density is responsible for the inaccuracies. Extending the probe volume over all cylindrical shells besides the one closest to the wall is thought to yield rather reliable results for the chemical potential. As a by-product of the simulations we also obtained diffusion coefficients, which are given in an appendix.     

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand–solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand–DNA binding.     

Function-related asymmetry of the specific cardiolipin binding sites on the mitochondrial ADP/ATP carrier

The ADP/ATP carrier (AAC) transports matrix ATP and cytosolic ADP across the inner mitochondrial membrane (IMM). It is well known that cardiolipin (CL) plays an important role in regulating the function of AAC, yet the underlying mechanism still remains elusive. AAC is composed of three homologous domains, and three specific CL binding sites are located at the domain-domain interfaces near the matrix side. Here we report an in-depth investigation on the dynamic properties of the bound CL within the three specific sites through all-atom molecular dynamics simulations of up to 13 μs in total. Our results highlight the importance of the basic and polar residues in CL binding. The basic residues from the linker helix and/or the [Y/W/F][K/R]G motif enable the bound CL to form an intra-domain binding mode, and the canonical inter-domain binding mode only forms when these basic residues are occupied by an additional phospholipid. Of special significance, differences in the basic and polar residues lead to remarkable asymmetry among the three specific CL binding sites. We found that the bound CL at the interface of domains 2 and 3 predominantly adopts inter-domain binding mode, while CLs at the other two sites have much more intra-domain populations. This is consistent with the asymmetric crystal structure of the matrix state (m-state) AAC which implies an asymmetric transport mechanism. The dynamic equilibrium between the inter-domain and intra-domain binding modes observed in our simulations could be highly important for the bound CLs to adapt to the movements during state transitions.     

Non-Ewald methods: theory and applications to molecular systems

Several non-Ewald methods for calculating electrostatic interactions have recently been developed, such as the Wolf method, the reaction field method, the pre-averaging method, and the zero-dipole summation method, for molecular dynamics simulations of various physical systems, including biomolecular systems. We review the theories of these approaches and their potential applications to molecular simulations, and discuss their relationships.     

A comparison of structural and dynamic properties of different simulation methods applied to SH3.

The dynamic and static properties of molecular dynamics simulations using various methods for treating solvent were compared. The SH3 protein domain was chosen as a test case because of its small size and high surface-to-volume ratio. The simulations were analyzed in structural terms by examining crystal packing, distribution of polar residues, and conservation of secondary structure. In addition, the "essential dynamics" method was applied to compare each of the molecular dynamics trajectories with a full solvent simulation. This method proved to be a powerful tool for the comparison of large concerted atomic motions in SH3. It identified methods of simulation that yielded significantly different dynamic properties compared to the full solvent simulation. Simulating SH3 using the stochastic dynamics algorithm with a vacuum (reduced charge) force field produced properties close to those of the full solvent simulation. The application of a recently described solvation term did not improve the dynamic properties. The large concerted atomic motions in the full solvent simulation as revealed by the essential dynamics method were analyzed for possible biological implications. Two loops, which have been shown to be involved in ligand binding, were seen to move in concert to open and close the ligand-binding site.     

Chemiluminescence of off‐line and on‐line gold nanoparticle‐catalyzed luminol system in the presence of flavonoid

It was found that flavonoids could remarkably inhibit the chemiluminescence (CL) intensity of an off‐line gold nanoparticle (AuNP)‐catalyzed luminol–H₂O₂ CL system. By contrast, flavonoids enhanced the CL intensity of an on‐line AuNP‐catalyzed luminol–H₂O₂ CL system. In the off‐line system, the AuNPs were prepared beforehand, whereas in the on‐line system, AuNPs were produced by on‐line mixing of luminol prepared in a buffer solution of NaHCO₃ ? Na₂CO₃ and HAuCl₄ with no need for the preliminary preparation of AuNPs. The on‐line system had prominent advantages over the off‐line system, namely a lowering of the background noise and improvements in the stability of the CL system. The results show that differences in the signal suppression effect of flavonoids on the off‐line AuNP‐catalyzed CL system are influenced by the combined action of a free radical scavenging effect and occupy‐sites function; the latter was proved to be predominant using controlled experiments. Enhancement of the on‐line system was ascribed to the presence of flavonoids promoting the on‐line formation of AuNPs, which better catalyzed the luminol–H₂O₂ CL reaction, and the enhancement activity of the six flavonoids increased with the increase in reducibility. This work broadens the scope of practical applications of an AuNP‐catalyzed CL system.     

The effects of cardiolipin on the structural dynamics of the mitochondrial ADP/ATP carrier in its cytosol-open state

Cardiolipin (CL) has been shown to play a crucial role in regulating the function of proteins in the inner mitochondrial membrane. As the most abundant protein of the inner mitochondrial membrane, the ADP/ATP carrier (AAC) has long been the model of choice to study CL-protein interactions, and specifically bound CLs have been identified in a variety of crystal structures of AAC. However, how CL binding affects the structural dynamics of AAC in atomic detail remains largely elusive. Here we compared all-atom molecular dynamics simulations on bovine AAC1 in lipid bilayers with and without CLs. Our results show that on the current microsecond simulation time scale: 1) CL binding does not significantly affect overall stability of the carrier or structural symmetry at the matrix-gate level; 2) pocket volumes of the carrier and interactions involved in the matrix-gate network become more heterogeneous in parallel simulations with membranes containing CLs; 3) CL binding consistently strengthens backbone hydrogen bonds within helix H2 near the matrix side; and 4) CLs play a consistent stabilizing role on the domain 1-2 interface through binding with the R30:R71:R151 stacking structure and fixing the M2 loop in a defined conformation. CL is necessary for the formation of this stacking structure, and this structure in turn forms a very stable CL binding site. Such a delicate equilibrium suggests the strictly conserved R30:R71:R151stacking structure of AACs could function as a switch under regulation of CLs. Taken together, these results shed new light on the CL-mediated modulation of AAC function.     

Systematic probing of an atomic charge set of sialic acid disaccharides for the rational molecular modeling of avian influenza virus based on molecular dynamics simulations

A systematic searching approach for an atomic charge set through molecular dynamics simulations is introduced to calculate a reasonable sialic acid carbohydrate conformation with respect to the experimentally observed structures. The present molecular dynamics simulation study demonstrated that B3LYP/6-31G is the most suitable basis set for the sialic acid disaccharides, attaining good agreement with experimental data.     

Multiscale Coarse-Grained Approach to Investigate Self-Association of Antibodies

Self-association of therapeutic monoclonal antibodies (mabs) are thought to modulate the undesirably high viscosity observed in their concentrated solutions. Computational prediction of such a self-association behavior is advantageous early during mab drug candidate selection when material availability is limited. Here, we present a coarse-grained (CG) simulation method that enables microsecond molecular dynamics simulations of full-length antibodies at high concentrations. The proposed approach differs from others in two ways: first, charges are assigned to CG beads in an effort to reproduce molecular multipole moments and charge asymmetry of full-length antibodies instead of only localized charges. This leads to great improvements in the agreement between CG and all-atom electrostatic fields. Second, the distinctive hydrophobic character of each antibody is incorporated through empirical adjustments to the short-range van der Waals terms dictated by cosolvent all-atom molecular dynamics simulations of antibody variable regions. CG simulations performed on a set of 15 different mabs reveal that diffusion coefficients in crowded environments are markedly impacted by intermolecular interactions. Diffusion coefficients computed from the simulations are in correlation with experimentally measured observables, including viscosities at a high concentration. Further, we show that the evaluation of electrostatic and hydrophobic characters of the mabs is useful in predicting the nonuniform effect of salt on the viscosity of mab solutions. This CG modeling approach is particularly applicable as a material-free screening tool for selecting antibody candidates with desirable viscosity properties.     

Parametric dependence on shear viscosity of SPC/E water from equilibrium and non-equilibrium molecular dynamics simulations

A parametric dependent study is crucial for the accurate determination of transport coefficients such as shear viscosity. In this study, we calculate the shear viscosity of extended simple point charge water using a transverse current auto-correlation function (TCAF) from equilibrium molecular dynamics (EMD) and the periodic perturbation method from non-equilibrium molecular dynamics (NEMD) simulations for varying coupling time and system sizes. Results show that the shear viscosity calculated using EMD simulations with different thermostats varies significantly with coupling times and system size. The use of Berendsen and velocity-rescale thermostats in NEMD simulations generates a significant drift from the target temperature and results in an inconsistent shear viscosity with coupling time and system size. The use of Nosé–Hoover thermostat in NEMD simulations offers thermodynamic stability which results in a consistent shear viscosity for various coupling times and system sizes.     

The investigation of the effects of counterions in protein dynamics simulations.

Molecular simulations able to exactly represent solvated charged proteins are helpful in understanding protein dynamics, structure and function. In the present study we have used two different starting structures of papain (a typical, stable, globular protein of intermediate net charge) and different modeling procedures to evaluate some effects of counterions in simulations. A number of configurations have been generated and relaxed for each system by various combinations of constrained simulated annealing and molecular dynamics procedures, using the AMBER force field. The analysis of trajectories shows that the simulations of solvated proteins are moderately sensitive to the presence of counterions. However, this sensitivity is highly dependent on the starting model and different procedures of equilibration used. The neutralized systems tend to evince smaller root mean square deviations regardless of the system investigated and the simulation procedure used. The results of parameterized fitting of the simulated structures to the crystallographic data, giving quantitative measure of the total charge influence on the stability of various elements of the secondary structure, revealed a clear scatter of different reactions of various systems' secondary structures to counterions addition: some systems apparently were stabilized when neutralized, while the others were not. Thus, one cannot unequivocally state, despite consideration of specific simulation conditions, whether protein secondary structures are more stable when they have neutralized charges. This suggests that caution should be taken when claiming the stabilizing effect of counterions in simulations other than those involving small, unstable polypeptides or highly charged proteins.     

The impact of line tension on the contact angle of nanodroplets

The line tension for a Lennard–Jones (LJ) fluid on a (9, 3) solid of varying strength was calculated using Monte Carlo simulations. A new perturbation method was used to determine the interfacial tension between liquid–vapour, solid–liquid and solid–vapour phases for this system to determine the Young's equation contact angle. Cylindrical and spherical nanodroplets were simulated for comparison. The contact angles from the cylindrical drops and Young's equation agree very well over the range of surface strengths and cylindrical drop sizes, except on a very weak surface. Tolman length effects were not observable for cylindrical drops. This shows that quite small systems can reproduce macroscopic contact angles. For spherical droplets, a deviation between the contact angle of spherical droplets and Young's equation was evident, but decreased with increasing interaction strengths to be negligible for contact angles less than 90°. Linear fitting of the contact angle data for varying droplet sizes showed no clear effect by line tension on contact angle. All calculated line tension values have a magnitude less than 4 × 10^? 12 J/m with both negative and positive signs. The best estimate of line tension for this system of LJ droplets was 1 × 10^? 13 J/m, which is smaller than the reported estimations in the literature, and is too small to be conclusively positive or negative in value.     

TTAPE-Me dye is not selective to cardiolipin and binds to common anionic phospholipids nonspecifically

Identification, visualization, and quantitation of cardiolipin (CL) in biological membranes is of great interest because of the important structural and physiological roles of this lipid. Selective fluorescent detection of CL using noncovalently bound fluorophore 1,1,2,2-tetrakis[4-(2-trimethylammonioethoxy)-phenylethene (TTAPE-Me) has been recently proposed. However, this dye was only tested on wild-type mitochondria or liposomes containing negligible amounts of other anionic lipids, such as phosphatidylglycerol (PG) and phosphatidylserine (PS). No clear preference of TTAPE-Me for binding to CL compared to PG and PS was found in our experiments on artificial liposomes, Escherichia coli inside-out vesicles, or Saccharomyces cerevisiae mitochondria in vitro or in situ, respectively. The shapes of the emission spectra for these anionic phospholipids were also found to be indistinguishable. Thus, TTAPE-Me is not suitable for detection, visualization, and localization of CL in the presence of other anionic lipids present in substantial physiological amounts. Our experiments and complementary molecular dynamics simulations suggest that fluorescence intensity of TTAPE-Me is regulated by dynamic equilibrium between emitting dye aggregates, stabilized by unspecific but thermodynamically favorable electrostatic interactions with anionic lipids, and nonemitting dye monomers. These results should be taken into consideration when interpreting past and future results of CL detection and localization studies with this probe in vitro and in vivo. Provided methodology emphasizes minimal experimental requirements, which should be considered as a guideline during the development of novel lipid-specific probes.     

A molecular dynamics study of the formation, stability, and oligomerization state of two designed coiled coils: possibilities and limitations

The formation, relative stability, and possible stoichiometries of two (self-)complementary peptide sequences (B and E) designed to form either a parallel homodimeric (B + B) or an antiparallel heterodimeric (B + E) coiled coil have been investigated. Peptide B shows a characteristic coiled coil pattern in circular dichroism spectra at pH 7.4, whereas peptide E is apparently random coiled under these conditions. The peptides are complementary to each other, with peptide E forming a coiled coil when mixed with peptide B. Molecular dynamics simulations show that combinations of B + B and B + E readily form a dimeric coiled coil, whereas E + E does not fall in line with the experimental data. However, the simulations strongly suggest the preferred orientation of the helices in the homodimeric coiled coil is antiparallel, with interactions at the interface quite different to that of the idealized model. In addition, molecular dynamics simulations suggest equilibrium between dimers, trimers, and tetramers of alpha-helices for peptide B.     

Cation-mediated interplay of loops in chaperonin-10

The ubiquitously occurring chaperonins consist of a large tetradecameric Chaperonin-60, forming a cylindrical assembly, and a smaller heptameric Chaperonin-10. For a functional protein folding cycle, Chaperonin-10 caps the cylindrical Chaperonin-60 from one end forming an asymmetric complex. The oligomeric assembly of Chaperonin-10 is known to be highly plastic in nature. In Mycobacterium tuberculosis, the plasticity has been shown to be modulated by reversible binding of divalent cations. Binding of cations confers rigidity to the metal binding loop, and also promotes stability of the oligomeric structure. We have probed the conformational effects of cation binding on the Chaperonin-10 structure through fluorescence studies and molecular dynamics simulations. Fluorescence studies show that cation binding induces reduced exposure and flexibility of the dome loop. The simulations corroborate these results and further indicate a complex landscape of correlated motions between different parts of the molecule. They also show a fascinating interplay between two distantly spaced loops, the metal binding "dome loop" and the GroEL-binding "mobile loop", suggesting an important cation-mediated role in the recognition of Chaperonin-60. In the presence of cations the mobile loop appears poised to dock onto the Chaperonin-60 structure. The divalent metal ions may thus act as key elements in the protein folding cycle, and trigger a conformational switch for molecular recognition.     

Atomic Level Rendering of DNA-Drug Encounter

Computer simulations have been demonstrated to be important for unraveling atomic mechanisms in biological systems. In this study, we show how combining unbiased molecular dynamic simulations with appropriate analysis tools can successfully describe metal-based drug interactions with DNA. To elucidate the noncovalent affinity of cisplatin’s family to DNA, we performed extensive all-atom molecular dynamics simulations (3.7 μs total simulation length). The results show that the parent drug, cisplatin, has less affinity to form noncovalent adducts in the major groove than its aquo complexes. Furthermore, the relative position in which the drugs enter the major groove is dependent on the compound’s net charge. Based on the simulations, we estimated noncovalent binding free energies through the use of Markov state models. In addition, and to overcome the lack of experimental information, we employed two additional methods: Molecular Mechanics Poisson-Boltzmann Surface Area (MMPB-SA) and steered molecular dynamics with the Jarzynski estimator, with an overall good agreement between the three methods. All complexes show interaction energies below 3 kcal/mol with DNA but the charged hydrolysis products have slightly more favorable binding free energies than the parent drug. Moreover, this study sets the precedent for future unbiased DNA-ligand simulations of more complex binders.     

Reaction Field Effects on the Simulated Properties of Liquid Water

Abstract

The effects of including a reaction field contribution on the structure and dynamics of liquid water have been investigated using molecular dynamics simulations. Reaction field effects are determined for two models of water, the simple point charge (SPC) model and the extended simple point charge (SPC/E) model, and at two temperatures (277 K and 300 K). Inclusion of the reaction field leads to a reduced system density, an increase in translational diffusion, which is model dependent, an increase in internal energy, and an increase in rotational diffusion rates, in addition to the large (known) changes in the dielectric properties of liquid water. It is concluded that continued use of the reaction field technique should involve a reparameterization of the water model and not merely a merging with the original model parameters.