All-atom Molecular Dynamics Calculations Research Articles

Single-Molecule Unidirectional Processive Movement along a Helical Polymer Chain in Non-aqueous Medium.

In this work, a molecule “walking” along a single chain of a synthetic helical polymer, which is used as a rail on a substrate in an organic solvent at room temperature, is observed. The walking comprises the unidirectional processive movement of a short-chain molecule along a chiral helical chain in 3 nm steps, driven by Brownian motion and a tapping effect of the atomic force microscopy tip based on a flash ratchet mechanism. Furthermore, the rail consists of a long-chain substituted phenylacetylene polymer with pendant cholesteryl groups, along which the short-chain molecule can walk as a result of van der Waals interactions. The macromolecular motion is videoed using a fast-scanning atomic force microscope, and additionally, this phenomenon is also simulated by all-atom molecular dynamics calculations. On the basis of these results, we propose the principle of a polymer molecular motor. This is the first report of a synthetic walking machine of a chiral helical polymer driven by thermal fluctuation as an artificial life function.

Molecular-Scale Description of SPAN80 Desorption from a Squalane-Water Interface.

Extensive all-atom molecular dynamics calculations on the water-squalane interface for nine different loadings with sorbitan monooleate (SPAN80), at T = 300 K, are analyzed for the surface tension equation of state, desorption free-energy profiles as they depend on loading, and to evaluate escape times for adsorbed SPAN80 into the bulk phases. These results suggest that loading only weakly affects accommodation of a SPAN80 molecule by this squalane-water interface. Specifically, the surface tension equation of state is simple through the range of high tension to high loading studied, and the desorption free-energy profiles are weakly dependent on loading here. The perpendicular motion of the centroid of the SPAN80 headgroup ring is well-described by a diffusional model near the minimum of the desorption free-energy profile. Lateral diffusional motion is weakly dependent on loading. Escape times evaluated on the basis of a diffusional model and the desorption free energies are 7 × 10-2 s (into the squalane) and 3 × 102 h (into the water). The latter value is consistent with desorption times of related lab-scale experimental work.

Direct prediction of residual dipolar couplings of small molecules in a stretched gel by stochastic molecular dynamics simulations.

Residual dipolar couplings are highly useful NMR parameters for calculating and refining molecular structures, dynamics, and interactions. For some applications, however, it is inevitable that the preferred orientation of a molecule in an alignment medium is calculated a priori. Several methods have been developed to predict molecular orientations and residual dipolar couplings. Being beneficial for macromolecules and selected small-molecule applications, such approaches lack sufficient accuracy for a large number of organic compounds for which the fine structure and eventually the flexibility of all involved molecules have to be considered or are limited to specific, well-studied liquid crystals. We introduce a simplified model for detailed all-atom molecular dynamics calculations with a polymer strand lined up along the principal axis as a new approach to simulate the preferred orientation of small to medium-sized solutes in polymer-based, gel-type alignment media. As is shown by a first example of strychnine in a polystyrene/CDCl3 gel, the simulations potentially enable the accurate prediction of residual dipolar couplings taking into account structural details and dynamic averaging effects of both the polymer and the solute.

Understanding the structural differences between spherical and rod-shaped human insulin nanoparticles produced by supercritical fluids precipitation.

In this study, the thermal denaturation mechanism and secondary structures of two types of human insulin nanoparticles produced by a process of solution-enhanced dispersion by supercritical fluids using dimethyl sulfoxide (DMSO) and ethanol (EtOH) solutions of insulin are investigated using spectroscopic approaches and molecular dynamics calculations. First, the temperature-dependent IR spectra of spherical and rod-shaped insulin nanoparticles prepared from DMSO and EtOH solution, respectively, are analyzed using principal component analysis (PCA) and 2D correlation spectroscopy to obtain a deeper understanding of the molecular structures and thermal behavior of the two insulin particle shapes. All-atom molecular dynamics (AAMD) calculations are performed to investigate the influence of the solvent molecules on the production of the insulin nanoparticles and to elucidate the geometric differences between the two types of nanoparticles. The results of the PCA, the 2D correlation spectroscopic analysis, and the AAMD calculations clearly reveal that the thermal denaturation mechanisms and the degrees of hydrogen bonding in the spherical and rod-shaped insulin nanoparticles are different. The polarity of the solvent might not alter the structure or function of the insulin produced, but the solvent polarity does influence the synthesis of different shapes of insulin nanoparticles.

All-atom molecular dynamics calculation study of entire poliovirus empty capsids in solution.

Small viruses that belong, for example, to the Picornaviridae, such as poliovirus and foot-and-mouth disease virus, consist simply of capsid proteins and a single-stranded RNA (ssRNA) genome. The capsids are quite stable in solution to protect the genome from the environment. Here, based on long-time and large-scale 6.5 × 10(6) all-atom molecular dynamics calculations for the Mahoney strain of poliovirus, we show microscopic properties of the viral capsids at a molecular level. First, we found equilibrium rapid exchange of water molecules across the capsid. The exchange rate is so high that all water molecules inside the capsid (about 200,000) can leave the capsid and be replaced by water molecules from the outside in about 25 μs. This explains the capsid's tolerance to high pressures and deactivation by exsiccation. In contrast, the capsid did not exchange ions, at least within the present simulation time of 200 ns. This implies that the capsid can function, in principle, as a semipermeable membrane. We also found that, similar to the xylem of trees, the pressure of the solution inside the capsid without the genome was negative. This is caused by coulombic interaction of the solution inside the capsid with the capsid excess charges. The negative pressure may be compensated by positive osmotic pressure by the solution-soluble ssRNA and the counter ions introduced into it.

A Tale of Two Ions and their Membrane Interactions: Clearly the Same or Clearly Different?

A Tale of Two Ions and their Membrane Interactions: Clearly the Same or Clearly Different?

Towards the prediction of order parameters from molecular dynamics simulations in proteins.

A molecular understanding of how protein function is related to protein structure requires an ability to understand large conformational changes between multiple states. Unfortunately these states are often separated by high free energy barriers and within a complex energy landscape. This makes it very difficult to reliably connect, for example by all-atom molecular dynamics calculations, the states, their energies, and the pathways between them. A major issue needed to improve sampling on the intermediate states is an order parameter--a reduced descriptor for the major subset of degrees of freedom--that can be used to aid sampling for the large conformational change. We present a method to combine information from molecular dynamics using non-linear time series and dimensionality reduction, in order to quantitatively determine an order parameter connecting two large-scale conformationally distinct protein states. This new method suggests an implementation for molecular dynamics calculations that may be used to enhance sampling of intermediate states.

Characterizing the unfolded states of proteins using single-molecule FRET spectroscopy and molecular simulations

To obtain quantitative information on the size and dynamics of unfolded proteins we combined single-molecule lifetime and intensity FRET measurements with molecular simulations. We compared the unfolded states of the 64-residue, alpha/beta protein L and the 66-residue, all-beta cold-shock protein CspTm. The average radius of gyration (Rg) calculated from FRET data on freely diffusing molecules was identical for the two unfolded proteins at guanidinium chloride concentrations >3 M, and the FRET-derived Rg of protein L agreed well with the Rg previously measured by equilibrium small-angle x-ray scattering. As the denaturant concentration was lowered, the mean FRET efficiency of the unfolded subpopulation increased, signaling collapse of the polypeptide chain, with protein L being slightly more compact than CspTm. A decrease in Rg with decreasing denaturant was also observed in all-atom molecular dynamics calculations in explicit water/urea solvent, and Langevin simulations of a simplified representation of the polypeptide suggest that collapse can result from either increased interresidue attraction or decreased excluded volume. In contrast to both the FRET and simulation results, previous time-resolved small-angle x-ray scattering experiments showed no collapse for protein L. Analysis of the donor fluorescence decay of the unfolded subpopulation of both proteins gives information about the end-to-end chain distribution and suggests that chain dynamics is slow compared with the donor life-time of approximately 2 ns, whereas the bin-size independence of the small excess width above the shot noise for the FRET efficiency distributions may result from incomplete conformational averaging on even the 1-ms time scale.