Vapour-Liquid Phase Equilibria of n-alkanes by Direct Monte Carlo Simulations

Abstract

We report results of direct Monte Carlo simulations of n-pentane and n-decane at the liquidvapour interface for a number of temperatures. The intermolecular interactions are modeled using the last version of the anisotropic united atom model (AUA4). We have used the local long range correction energy and an algorithm allowing to select randomly with equal probability two different displacements. The liquid and vapour densities are in excellent agreement with experimental data and with those previously calculated using the GEMC method.     

A Leap-frog Algorithm for Stochastic Dynamics

Abstract

A third-order algorithm for stochastic dynamics (SD) simulations is proposed, identical to the powerful molecular dynamics leap-frog algorithm in the limit of infinitely small friction coefficient γ. It belongs to the class of SD algorithms, in which the integration time step Δt is not limited by the condition Δt ≤ γ^?1, but only by the properties of the systematic force. It is shown how constraints, such as bond length or bond angle constraints, can be incorporated in the computational scheme. It is argued that the third-order Verlet-type SD algorithm proposed earlier may be simplified without loosing its third-order accuracy. The leap-frog SD algorithm is proven to be equivalent to the verlet-type SD algorithm. Both these SD algorithms are slightly more economical on computer storage than the Beeman-type SD algorithm.     

The Force-Biased Algorithm for the Irregular Close Packing of Equal Hard Spheres

Abstract

We present a “force-biased” algorithm for generating the irregular close packing of hard spheres. The algorithm is partly based on Jodrey and Tory's ideas [9] and incorporates methods from Molecular Dynamics. Packings generated by means of the two algorithms are consistent up to final packing fraction of 0.65, which seems to be the limit density of Jodrey and Tory's method. Significantly higher densities (up to 0.71) can be achieved for small numbers of spheres by the force-biased algorithm. However the shape of the radial and angle distribution functions implies that a partial short-range ordering occurs in packings of those densities.     

Molecular dynamics simulation of the Al2O3 film structure during atomic layer deposition

Based on an understanding of atomic layer deposition (ALD) from prior experimental and computational results, all-atom molecular dynamics (MD) simulations are used to model the Al₂O₃ film structure and composition during ALD processing. By separating the large time-scale surface reactions from the small time-scale structural relaxation, we have focused on the growth dynamics of amorphous Al₂O₃ films at the atomic scale. The simulations are able to reproduce some important properties and growth mechanisms of Al₂O₃ ALD films, and hence provide a bridge between atomic-level information and experimental measurements. Information about the evolution of the microscopic structures of the Al₂O₃ films is generated, and the influence of operation parameters on the Al₂O₃ ALD process. The simulations predict a strong influence of the initial surface composition and process temperature on the surface roughness, growth rate and growth mode of the deposited films.     

Understanding the specificity of serpin–protease complexes through interface analysis

Serpins such as antithrombin, heparin cofactor II, plasminogen activator inhibitor, antitrypsin, antichymotrypsin, and neuroserpin are involved in important biological processes by inhibiting specific serine proteases. Initially, the protease recognizes the mobile reactive loop of the serpin eliciting conformational changes, where the cleaved loop together with the protease inserts into β-sheet A, translocating the protease to the opposite side of inhibitor leading to its inactivation. Serpin interaction with proteases is governed mainly by the reactive center loop residues (RCL). However, in some inhibitory serpins, exosite residues apart from RCL have been shown to confer protease specificity. Further, this forms the basis of multi-specificity of some serpins, but the residues and their dimension at interface in serpin-protease complexes remain elusive. Here, we present a comprehensive structural analysis of the serpin-protease interfaces using bio COmplexes COntact MAPS (COCOMAPS), PRotein Interface Conservation and Energetics (PRICE), and ProFace programs. We have carried out interface, burial, and evolutionary analysis of different serpin-protease complexes. Among the studied complexes, non-inhibitory serpins exhibit larger interface region with greater number of residue involvement as compared to the inhibitory serpins. On comparing the multi-specific serpins (antithrombin and antitrypsin), a difference in the interface area and residue number was observed, suggestive of a differential mechanism of action of these serpins in regulating their different target proteases. Further, detailed study of these multi-specific serpins listed few essential residues (common in all the complexes) and certain specificity (unique to each complex) determining residues at their interfaces. Structural mapping of interface residues suggested that individual patches with evolutionary conserved residues in specific serpins determine their specificity towards a particular protease.     

Structural analysis of the envelope gp120 V3 loop for some HIV-1 variants circulating in the countries of Eastern Europe

The V3 loop on gp120 from human immunodeficiency virus type 1 (HIV-1) is a focus of many research groups involved in anti-AIDS drug development because this region of the protein is a principal target for neutralizing antibodies and a major determinant for cell tropism and syncytium formation. In this study, the nucleotide sequences of the env gene region coding the V3 loop were determined by DNA sequencing methods for four novel HIV-1 strains that circulate in the countries of Eastern Europe, such as Russia, Belarus, Ukraine, etc. Based on the empirical data obtained, the 3D structures of the V3 loops associated with these viral modifications were generated by computer modeling and then compared to discover similarities in the spatial arrangement of this functionally important site of gp120. Despite the HIV-1 genetic variety, several regions of the V3 loop that contain residues critical for cell tropism were shown to be structurally invariant, which may explain its exceptional role in a co-receptor usage. These data together with those on the biological activity of the V3 individual residues clearly show that these conserved structural motifs of gp120 represent potential HIV-1 weak points most suitable for therapeutic intervention.     

Transmembrane helix 6 observed at the interface of β2AR homodimers in blind docking studies

Peptide– and protein–protein dockings were carried out on β₂-adrenergic receptor (β₂AR) to confirm the presence of transmembrane helix 6 (TM6) at the interface region between two β₂AR monomers, thereby its possible role in dimerization as suggested in numerous experimental and computational studies. Initially, a portion of TM6 was modeled as a peptide consisting of 23 residues and blindly docked to β₂AR monomer using a rigid body approach. Interestingly, all highest score conformations preferred to be near TM5 and TM6 regions of the receptor. Furthermore, longer peptides generated from a whole TM region were blindly docked to β₂AR using the same rigid body approach. This yielded a total of seven docked peptides, each derived from one TM helix. Most interestingly, for each peptide, TM6 was among the most preferred binding site region in the receptor. Besides the peptide dockings, two β₂AR monomers were blindly docked to each other using a full rigid-body search of docking orientations, which yielded a total of 16,000 dimer conformations. Each dimer was then filtered according to a fitness value based on the membrane topology. Among 149 complexes that met the topology requirements, 102 conformers were composed of two monomers oriented in opposite directions, whereas in the remaining 47, the monomers were arranged in parallel. Lastly, all 149 conformers were clustered based on a root mean-squared distance value of 6 Å. In agreement with the peptide results, the clustering yielded the largest population of conformers with the highest Z-score value having TM6 at the interface region.     

Electrotactile stimulation on the tongue: Intensity perception,discrimination, and cross-modality estimation

Due to its high sensitivity and conductivity, electrotactile stimulation (ETS) on the tongue has proven to be a useful and technically convenient tool to substitute and/or augment sensory capabilities. However, most of its applications have only provided spatial attributes and little is known about (a) the ability of the tongue's sensory system to process electrical stimuli of varying magnitudes and (b) how modulation of ETS intensity affects subjects’ ability to decode stimulus intensity. We addressed these questions by quantifying: (1) the magnitude of the dynamic range (DR; maximal comfortable intensity/perception threshold) and its sensitivity to prolonged exposure; (2) subjects’ ability to perceive intensity changes; and (3) subjects’ ability to associate intensity with angular excursions of a protractor's handle. We found that the average DR (17 dB) was generally large in comparison with other tactile loci and of a relatively constant magnitude among subjects, even after prolonged exposure, despite a slight but significant upward drift (p < 0.001). Additionally, our results showed that as stimulus intensity increased, subjects’ ability to discriminate ETS stimuli of different intensities improved (p < 0.05) while estimation accuracy, in general, slightly decreased (increasing underestimation). These results suggest that higher ETS intensity may increase recruitment of rapidly adapting mechanoreceptor fibers, as these are specialized for coding stimulus differences rather than absolute intensities. Furthermore, our study revealed that the tongue's sensory system can effectively convey electrical stimuli despite minimal practice and when information transfer is limited by memory and DR drift.     

Asymmetric structure and domain binding interfaces of human tyrosyl‐tRNA synthetase studied by molecular dynamics simulations

Human tyrosyl‐tRNA synthetase (HsTyrRS) is composed of two structural modules: N‐terminal catalytic core and an EMAP II‐like C‐terminal domain. The structures of these modules are known, but no crystal structure of the full‐length HsTyrRS is currently available. An all‐atom model of the full‐length HsTyrRS was developed in this work. The structure, dynamics, and domain binding interfaces of HsTyrRS were investigated by extensive molecular dynamics (MD) simulations. Our data suggest that HsTyrRS in solution consists of a number of compact asymmetric conformations, which differ significantly by their rigidity, internal mobility, orientation of C‐terminal modules, and the strength of interdomain binding. Interfaces of domain binding obtained in MD simulations are in perfect agreement with our previous coarse‐grained hierarchical rotations technique simulations. Formation of the hydrogen bonds between R93 residue of the ELR cytokine motif and the residues A340 and E479 in the C‐module was observed. This observation supports the idea that the lack of cytokine activity in the full‐length HsTyrRS is explained by interactions between N‐modules and C‐modules, which block the ELR motif. Copyright © 2013 John Wiley & Sons, Ltd.