A Molecular Dynamics Investigation of Compressed Aqueous Alkanoate Monolayers

Molecular dynamics are conducted on a dodecanoic acid monolayer/aqueous surface. Surface pressure is controlled by imposing constant-volume conditions for series of lengths of the square slab constituting the MD cell. The response of the alkanoate chains to the pressure is followed by examining various computed quantities that monitor their conformational order. These include atom-pair radial distribution functions, chain torsional angles, energies, atomic densities perpendicular to the interface, diffusivities and atomic plots. These quantities lead to chain separations which in the range 4-5 Å implying order when the alkanoate chains have a mean area of 0.18 nm 2.     

An Extended Gibbs Ensemble

Abstract

A general extended Gibbs ensemble, obtained by augmenting the standard Gibbs ensemble by intermediate states in the spirit of the scaled particle method of Nezbeda and Kolafa [Molec. Simul., 5, 391 (1991)], is introduced. The intermediate states span the states with different number of particles in the simulation boxes and facilitate the transfer of particles even in such complex systems as e.g., mixtures of very different components, systems of flexible polymeric molecules, or systems at very high densities. A general formulation of the ensemble is given and two implementations are considered in detail. The method is very general and is exemplified by studying the fluid-fluid coexistence in a dense binary mixture of the hard-sphere and square-well fluids. It is found that its efficiency is about by factor three greater in comparison with the standard Gibbs ensemble simulations.     

Calculation of Vapour-Liquid Equilibria of Lennard-Jones Binary Systems by the Gibbs Ensemble Monte Carlo Simulation

Abstract

The Gibbs ensemble Monte Carlo simulation has been used to calculate vapour-liquid equilibria of a Lennard-Jones (LJ) binary mixture. The mixture studied is the LB-2-1 model which has been used in our previous calculations on PVT relation and density-dependent local composition. The P-x-y relation has been established at two different temperatures and used to determine vapour-liquid coexistence region in the PVTx space.     

Did the prion protein become vulnerable to misfolding after an evolutionary divide and conquer event?

Despite high sequence identity among mammalian prion proteins (PrPs), mammals have varying rates of susceptibility to prion disease resulting in a so-called species barrier. The species barrier follows no clear pattern, with closely related species or similar sequences being no more likely to infect each other, and remains an unresolved enigma. Variation of the conformationally flexible regions may alter the thermodynamics of the conformational change, commonly referred to as the conformational conversion, which occurs in the pathogenic process of the mammalian prion protein. A conformational ensemble scenario is supported by the species barrier in prion disease and evidence that there are strains of pathogenic prion with different conformations within species. To study how conformational flexibility has evolved in the prion protein, an investigation was undertaken on the evolutionary dynamics of structurally disordered regions in the mammalian prion protein, non-mammalian prion protein that is not vulnerable to prion disease, and remote homologs Doppel and Shadoo. Structural disorder prediction analyzed in an evolutionary context revealed that the occurrence of increased or altered conformational flexibility in mammalian PrPs coincides with key events among PrP, Doppel, and Shadoo. Comparatively rapid evolutionary dynamics of conformational flexibility in the prion protein suggest that the species barrier is not a static phenomenon. A small number of amino acid substitutions can repopulate the conformational ensemble and have a disproportionately large effect on pathogenesis.     

Conformational restriction approach to BACE1 inhibitors II: SAR study of the isocytosine derivatives fixed with a cis-cyclopropane ring

To improve the efficacy of the conformationally restricted BACE1 inhibitors, structural modifications were investigated using two strategies: (a) modification of the terminal aromatic ring and (b) insertion of a spacer between the aromatic rings. In the latter approach, another type of inhibitor 17 bearing an ethylene spacer between two aromatic rings was found to exhibit good BACE1 inhibitory activity, while the corresponding conformationally unrestricted compound 25 showed no activity. This result revealed an interesting effect of a conformational restriction with a cyclopropane ring.     

Speeding up Monte Carlo molecular simulation by a non-conservative early rejection scheme

Monte Carlo (MC) molecular simulation describes fluid systems with rich information, and it is capable of predicting many fluid properties of engineering interest. In general, it is more accurate and representative than equations of state. On the other hand, it requires much more computational effort and simulation time. For that purpose, several techniques have been developed in order to speed up MC molecular simulations while preserving their precision. In particular, early rejection schemes are capable of reducing computational cost by reaching the rejection decision for the undesired MC trials at an earlier stage in comparison to the conventional scheme. In a recent work, we have introduced a ‘conservative’ early rejection scheme as a method to accelerate MC simulations while producing exactly the same results as the conventional algorithm. In this paper, we introduce a ‘non-conservative’ early rejection scheme, which is much faster than the conservative scheme, yet it preserves the precision of the method. The proposed scheme is tested for systems of structureless Lennard-Jones particles in both canonical and NVT-Gibbs ensembles. Numerical experiments were conducted at several thermodynamic conditions for different number of particles. Results show that at certain thermodynamic conditions, the non-conservative method is capable of doubling the speed of the MC molecular simulations in both canonical and NVT-Gibbs ensembles.