Dissipative particle dynamics of non-spherical particles using a Gaussian density model

The Gaussian density molecular model has been adapted for dissipative particle dynamics. The model, when combined with a soft potential, is shown to be a very flexible mesoscale model exhibiting a wide range of phase behaviour. The soft potential allows relatively large time steps to be used and hence a more rapid equilibration. In addition, the model can be used to study both uniaxial and biaxial systems. We have undertaken a number of pilot studies and have demonstrated that the Gaussian model is able to identify nematic–isotropic phase transitions in liquid crystals and the formation of ordered discotic phases.     

Irreversibility and Interatomic Potentials in One-dimensional Lattice Models

Abstract

Since the computational work on the relaxation process toward the equilibrium state by Fermi, Pasta, and Ulam (FPU) [1], the stochastic behaviors in one dimensional lattice models have been studied by many physicists numerically and theoretically. The objective of those studies is to understand the origin of thermodynamical properties from the microscopic particle motions which are determined by the deterministic equations of motion. In the FPU work, since the given energy was too small and anharmonicity was too weak, thermodynamic irreversibility phenomena, which lead to the equipartition of energy among modes, were not observed. Many studies after this work confirmed the existence of the threshold for the chaotic motions [2] [3]. However, the relation of the chaotic motion in the system with many degrees of freedom to the thermodynamical properties is not well understood [4]. The steady lattice thermal conduction in one dimensional lattice poses an interesting problem in this; respect. A clear linear internal temperature gradient was not observed in the FPU model [5], while it was observed in the the ding-a-ling model by Casati et al. [6] and the diatomic Toda lattice (DTL) by Mokross and Büttner [7]. Recently, the system size dependence of the coefficient of thermal conductivity has been studied [8][9] and the Fourier's law of heat conduction is begining to be confirmed in the DTL.     

Accelerating all-atom protein folding simulations through reduced dihedral barriers

Protein folding requires extensive changes of backbone and sidechain dihedral angles, whose energy barriers constitute obstacles for the folding kinetics. Folding of small proteins is furthermore thought to be path-independent. Here, we propose that time-consuming all-atom protein folding simulations may be accelerated through a reduction of the dihedral barriers of the force field. In order to investigate this hypothesis, we performed various folding simulations of two small proteins. We report an acceleration towards smaller root-mean-square deviations from the native protein structure using our proposed method.     

A novel model for the molecular dynamics simulation study on mechanical properties of HMX/F2311 polymer-bonded explosive

The ‘insert’ model for β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX)-based polymer-bonded explosive (PBX) was proposed for finding the relation of temperatures with mechanical properties. This model was simulated by using molecular dynamics models. The elastic constants and the effective moduli were calculated with static analysis method. Cauchy pressure was also calculated. It is found that the rigidity is weakened and the ductibility is improved by adding a small amount of F₂₃₁₁ in the crystalline HMX. The rigidity is also weakened with increasing temperature. However, the ductibility of HMX/F₂₃₁₁ PBX changes as a parabola with increasing temperature duo to the enhancement of F₂₃₁₁ molecular chain movement and simultaneously the increment of high energy conformation ratio in this molecular chain, i.e. the increment of the molecular chain rigidity.     

Molecular docking and pharmacophore model studies of Rho kinase inhibitors

Molecular docking and pharmacophore model approaches were used to characterise the binding features of four different series of Rho kinase (ROCK) inhibitors. Docking simulation of 20 inhibitors with ROCK was performed. The binding conformations and binding affinities of these inhibitors were obtained using AutoDock 4.0 software. The predicted binding affinities correlate well with the activities of these inhibitors (R ² = 0.904). 3D pharmacophore models were generated for ROCK based on highly active inhibitors implemented in Catalyst 4.11 program. The best pharmacophore model consists of one hydrogen bond acceptor feature and two hydrophobic features, and they all seemed to be essential for inhibitors in terms of their binding activities. It is anticipated that the findings reported in this paper may provide very useful information for designing new ROCK inhibitors.     

Selenium and selenium-antagonistic elements in nutritional cancer prevention

Selenium is an essential trace element with antioxidative, antimutagenic, antiviral and anticarcinogenic properties. There is increasing evidence that the dietary selenium intakes are sub-optimal in the populations of many countries and that human cancer mortalities would significantly decline if additional selenium was made available either through supplementation or the fortification of certain foods. An important property of selenium is its interaction with other elements that may be present in foods, the water, the workplace and the environment, e.g. As, Cu, Ni, Co, Cr, Mn, Zn, Cd, Sn, Pb, Hg, Bi, Mo, Ag, Au, etc. The sequestration of elements by selenium represents an efficient natural detoxification mechanism for some of these elements but also results in the physiological inactivation of selenium. Animal experiments confirm that the chronic exposure to low levels of these elements abolishes the cancer-protective effect of selenium. Human cancer is likewise significantly determined by the interactions of selenium with other elements, as evidenced by epidemiological, ecological and case-control studies. Cadmium, for example, is a key risk-increasing element for prostate cancer; for breast cancer, Cd, Cr, Zn are mainly contributing; for bronchial cancer (in smelter workers), Cd, As, Cr, Sb, Co, La, all these elements are in a reciprocal relationship with Se. While selenium remains the key cancer-protective trace element, the interpretation of its mode of action necessitates consideration of the effects of selenium antagonistic elements.     

Microsecond Scale Replica Exchange Molecular Dynamic Simulation of Villin Headpiece: An Insight into the Folding Landscape

Abstract

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 μs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins     

Identification of structural motifs in the E2 glycoprotein of Chikungunya involved in virus–host interaction

Chikungunya fever is one of the reemerging vector-borne diseases. It has become a major global health problem especially in the developing countries. There are no vaccines or specific antiviral drugs available to date. This study reports small molecule inhibitors of envelope glycoprotein 2 (E2 glycoprotein) which are predicted based on Chikungunya virus–host interactions. E2 glycoprotein of Chikungunya virus interacts at 216 residue of the host receptor protein which plays a vital role in initiating infection. Understanding the structural aspects of E2 glycoprotein is crucial to develop specific inhibitors to prevent the virus binding from host receptors. In silico method was adopted to predict the sequence motifs of envelope protein, as the method like yeast two hybrid system is laborious, time consuming, and costly. The E2 glycoprotein structure of the Indian isolate was modeled using two templates (2XFC and 3JOC) and then validated. The class III PDZ domain binding motif was found to be identified at 213–216 amino acids. The corresponding peptide structures which recognize the PDZ domain binding motif were identified by the literature search and were used for generating five point pharmacophore model (ADDDR) containing acceptor, donor and aromatic ring features. Databases such as Asinex, TosLab and Maybridge were searched for the matches for the predicted pharmacophore model. Two compounds were identified as lead molecules as their glide score is?>?5?kcal/mol. Since the pharmacophore model is developed based on Chikungunya virus–host interaction, it can be used for designing promising antiviral lead compounds for the treatment of Chikungunya fever.An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:21     

Protein Sequence Modules

Abstract

Conserved protein sequence segments are commonly believed to correspond to functional sites in the protein sequence. A novel approach is proposed to profile the changing degree of conservation along the protein sequence, by evaluating the occurrence frequencies of all short oligopeptides of the given sequence in a large proteome database. Thus, a protein sequence conservation profile can be plotted for every protein. The profile indicates where along the sequences the potential functional (conserved) sites are located. The corresponding oligopeptides belonging to the sites are very frequent across many prokaryotic species. Analysis of a representative set of such profiles reveals a common feature of all examined proteins: they consist of sequence modules represented by the peaks of conservation. Typical size of the modules (peak-to-peak distance) is 25–30 amino acid residues.