Identification and structural characterization of deleterious non-synonymous single nucleotide polymorphisms in the human SKP2 gene

In SCF (Skp, Cullin, F-box) ubiquitin-protein ligase complexes, S-phase kinase 2 (SKP2) is one of the major players of F-box family, that is responsible for the degradation of several important cell regulators and tumor suppressor proteins. Despite of having significant evidence for the role of SKP2 on tumorgenesis, there is a lack of available data regarding the effect of non-synonymous polymorphisms. In this communication, the structural and functional consequences of non-synonymous single nucleotide polymorphisms (nsSNPs) of SKP2 have been reported by employing various computational approaches and molecular dynamics simulation. Initially, several computational tools like SIFT, PolyPhen-2, PredictSNP, I-Mutant 2.0 and ConSurf have been implicated in this study to explore the damaging SNPs. In total of 172 nsSNPs, 5 nsSNPs were identified as deleterious and 3 of them were predicted to be decreased the stability of protein. Guided from ConSurf analysis, P101L (rs761253702) and Y346C (rs755010517) were categorized as the highly conserved and functional disrupting mutations. Therefore, these mutations were subjected to three dimensional model building and molecular dynamics simulation study for the detailed structural consequences upon the mutations. The study revealed that P101L and Y346C mutations increased the flexibility and changed the structural dynamics. As both these mutations are located in the most functional regions of SKP2 protein, these computational insights might be helpful to consider these nsSNPs for wet-lab confirmatory analysis as well as in rationalizing future population based studies and structure based drug design against SKP2.     

Computational prediction of the effects of non-synonymous single nucleotide polymorphisms on the GPI-anchor transamidase subunit GPI8p of Plasmodium falciparum

Drug resistance is increasingly evolving in malaria parasites; hence, it is important to discover and establish alternative drug targets. In this context, GPI-anchor transamidase (GPI-T) is a potential drug target primarily of its crucial role in the development and survival of the parasite in the GPI anchor biosynthesis pathway. The present investigation was undertaken to explore the plausible effects of nsSNP on the structure and functions of GPI-T subunit GPI8p of Plasmodium falciparum. The GPI8p (PF3D7_1128700) was analyzed using various sequence-based and structure-based computational tools such as SIFT, PROVEAN, PredictSNP, SNAP2, I-Mutant, MuPro, ConSurf, NetSurfP, MUSTER, COACH server and STRING server. Of the 34 nsSNPs submitted for functional analysis, 18 nsSNPs (R124 L, N143 K, Y145 F, V157I, T195S, K379E, I392 K, I437 T, Y438H, N439D, Y441H, N442D, N448D, N451D, D457A, D457Y, I458 L and N460 K) were predicted to have deleterious effects on the protein GPI8p. Additionally, I-Mutant 2.0 and MuPro both showed a decrease in stability after mutation as a result of these nsSNPs, suggesting the destabilization of protein. ConSurf findings suggest that most of the regions were highly conserved. In addition, COACH server was used to predict the ligand binding sites. It was found that no mutation was present at the predicted ligand binding site. The results of the STRING database showed that the protein GPI8p interacts with those proteins which either involve the biosynthetic process of attaching GPI anchor to protein or GPI anchor. The present study suggested that the GPI8p could be a novel target for anti-malarial drugs, which provides significant details for further experimentation.     

A study on the structural transition of a single polymer chain by parallel tempering molecular dynamics simulation

The structural transition of a single polymer chain with chain length of 100,200 and 300 beads was investigated by parallel tempering MD simulation.Our simulation results can capture the structural change from random coil to orientationally ordered structure with decreasing temperature.The clear transition was observed on the curves of radius of gyration and global orientational order parameter P as the function of temperature,which demonstrated structural formation of a single polymer chain.The linear relationships between three components of square radius of gyration R_gx²,R_gx²,R_gz² and global orientational order P can be obtained under the structurally transformational process.The slope of the linear relationship between x（or y-axis） component R_gx²(or R_gy²) and P is negative,while that of RL as the function of P is positive.The absolute value of slope is proportional to the chain length.Once the single polymer chain takes the random coil or ordered configuration,the linear relationship is invalid.The conformational change was also analyzed on microscopic scale.The polymer chain can be treated as the construction of rigid stems connecting by flexible loops.The deviation from exponentially decreased behavior of stem length distribution becomes prominent,indicating a stiffening of the chain arises leading to more and more segments ending up in the trans state with decreasing temperature.The stem length N_tr is about 21 bonds indicating the polymer chain is ordered with the specific fold length.So,the simulation results,which show the prototype of a liquid-crystalline polymer chain,are helpful to understand the crystallization process of crystalline polymers.     

Investigation of chiral recognition by molecular micelles with molecular dynamics simulations

Molecular dynamics simulations were used to characterize the binding of the chiral drugs chlorthalidone and lorazepam to the molecular micelle poly-(sodium undecyl-(L)-leucine-valine). The project’s goal was to characterize the nature of chiral recognition in capillary electrophoresis separations that use molecular micelles as the chiral selector. The shapes and charge distributions of the chiral molecules investigated, their orientations within the molecular micelle chiral binding pockets, and the formation of stereoselective intermolecular hydrogen bonds with the molecular micelle were all found to play key roles in determining where and how lorazepam and chlorthalidone enantiomers interacted with the molecular micelle.     

Glass transition in single poly(ethylene oxide) chain: A molecular dynamics simulation study

Local dynamics of single poly(ethylene oxide) chain in various environments (bulk, film, and isolated systems) has been characterized by the reorientation functions of various backbone bond vectors. Within any observation time, the variations of these reorientation functions with the temperature can be well described by the Kohlrausch?Williams?Watts (KWW) like equation, in which the fitted temperature parameter is identified as the glass transition temperature (T _g). The so‐obtained T _g for that polymer faithfully reveals the effects of the observation time, chain flexibility and vector range on the local dynamics. Furthermore, it is found that the KWW like relation is also applicable to the temperature‐dependence of the fraction of frozen atoms or torsions defined by the trajectory radii of gyration or the conformational transitions. Consequently, different motions lead to different values of T _g for the same system. Despite all, the consistent trend can be yielded, namely, T _g (bulk) > T _g (film) > T _g (isolated), which captures the effects of free surfaces on enhanced dynamics. In addition, dynamics heterogeneity in the systems can be quantitatively revealed. The newly proposed method holds a bright promise to predict the T _g values of complex polymers especially for comparisons. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2017 , 55 , 178–188     

Nonlinear shear of entangled polymers from nonequilibrium molecular dynamics

This study aims to use molecular dynamics (MD) simulations of Kremer–Grest (KG) chains to inform future developments of models of entangled polymer dynamics. We perform nonequilibrium MD simulations, under shear flow, for well‐entangled KG chains. We study chains of 512 and 1000 KG beads, corresponding to 8 and 15 entanglements, respectively. We compute the linear rheological properties from equilibrium simulations of the stress autocorrelation and obtain from these data the tube model parameters. Under nonlinear shear flow, we compute the shear viscosity, the first and second normal stress differences, and chain contour length. For chains of 512 monomers, we obtain agreement with the results of Cao and Likhtman (ACS Macro. Lett. 2015, 4, 1376). We also compare our nonlinear results with the Graham, Likhtman and Milner‐McLeish (GLaMM) model. We identify some systematic disagreement that becomes larger for the longer chains. We made a comparison of the transient shear stress maximum from our simulations, two nonlinear models and experiments on a wide range of melts and solutions, including polystyrene (PS), polybutadiene, and styrene–butadiene rubber. This comparison establishes that the PS melt data show markedly different behavior to all other melts and solutions and KG simulations reproduce the PS data more closely than either the GLaMM or Xie and Schweizer models. We discuss the performance of these models against the data and simulations. Finally, by imposing a rapid reversing flow, we produce a method to extract the recoverable strain from MD simulations, valid for sufficiently entangled monodisperse polymers. We explore how the resulting data can probe the melt state just before the reversing flow. © 2019 The Authors. Journal of Polymer Science Part B: Polymer Physics published by Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2019 , 57, 1692–1704     

Adsorption of a single polymer chain on a surface: A molecular dynamics simulation study

We performed simulations of the physical adsorption of a single globular chain on a surface of hemispherical shape by means of molecular dynamics simulations. For the chain, we took advantage of a united atom model. Interactions within the chain were limited to stretching, bending, and torsional as well as nonbonded interactions between the nonadjacent atoms. The interaction between each chain element and the surface formation are reigned by a Lennard–Jones potential. In this article, we focused on differences in the behavior of the adsorbed globule to the free unadsorbed one particularly in two different zones of the immediate vicinity of the surface. There were strong indications for a localized acceleration of the dynamics as compared with the bulk that appears in an increase of trans–gauche switches. For explanation we came up with an adsorption scenario. Special attention was given to the shift of the percentage of trans and gauche conformations within the globule in dependence on the strength of the adsorption potential that might be related to crystallization or glass transition. © 2001 John Wiley & Sons, Inc. J Polym Sci Part B: Polym Phys 39: 2333–2339, 2001     

The absorption and diffusion of polyethylene chains on the carbon nanotube: The molecular dynamics study

Classical molecular dynamics simulations have been used to investigate the absorption and diffusion behavior of polyethylene (PE) chains on the surface of the side‐wall of the carbon nanotube (CNT). Different degrees of polymerization from 50 to 80 at separate temperatures of 300, 400, 500, and 600 K are considered. Through the simulation, it is examined that the PE chains are absorbed on the surface of CNT and form stable composites with the nanotube as capsules. It is found that the most probable distance between the CNT and the C atoms in backbone of PE molecules only attribute to the temperature, and at T = 300 K, this distance is about 3.8 Å. Furthermore, the pattern of the composites mainly depends on the temperature and the length matching of the chains and the CNT. In particular, the PE chains keep approximately linear conformation, and extend along the axis of the CNT at the room temperature. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 46: 272–280, 2008     

The flying ice cube: Velocity rescaling in molecular dynamics leads to violation of energy equipartition

This article describes an unexpected phenomenon encountered during MD simulations: velocity rescaling using standard protocols can systematically change the proportion of total kinetic energy (KE) found in motions associated with the various degrees of freedom. Under these conditions, the simulation violates the principle of equipartition of energy, which requires a mean kinetic energy of RT/2 in each degree of freedom. A particularly pathological form of this problem occurs if one does not periodically remove the net translation of (and rotation about) the center of mass. In this case, almost all of the kinetic energy is converted into these two kinds of motion, producing a system with almost no kinetic energy associated with the internal degrees of freedom. We call this phenomenon “the flying ice cube.” We present a mathematical analysis of a simple diatomic system with two degrees of freedom, to document the origin of the problem. We then present examples from three kinds of MD simulations, one being an in vacuo simulation on a diatomic system, one involving a low resolution model of DNA in vacuo, and the third using a traditional all-atom DNA model with full solvation, periodic boundary conditions, and the particle mesh Ewald method for treating long-range electrostatics. Finally, we discuss methods for avoiding the problem. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 726–740, 1998     

Force-field parametrization and molecular dynamics simulations of Congo red

Congo red, a diazo dye widely used in medical diagnosis, is known to form supramolecular systems in solution. Such a supramolecular system may interact with various proteins. In order to examine the nature of such complexes empirical force field parameters for the Congo red molecule were developed. The parametrization of bonding terms closely followed the methodology used in the development of the charmm22 force field, except for the calculation of charges. Point charges were calculated from a fit to a quantum mechanically derived electrostatic potential using the CHELP-BOW method. Obtained parameters were tested in a series of molecular dynamics simulations of both a single molecule and a micelle composed of Congo red molecules. It is shown that newly developed parameters define a stable minimum on the hypersurface of the potential energy and crystal and ab initio geometries and rotational barriers are well reproduced. Furthermore, rotations around C-N bonds are similar to torsional vibrations observed in crystals of diphenyl-diazene, which confirms that the flexibility of the molecule is correct. Comparison of results obtained from micelles molecular dynamics simulations with experimental data shows that the thermal dependence of micelle creation is well reproduced.     

A comprehensive in Silico analysis of the functional and structural impact of single nucleotide polymorphisms (SNPs) in the human IL-33 gene

Interleukin 33 (IL-33) is the latest member of the IL-1 cytokine family, which plays both pro - and anti-inflammatory functions. Numerous Single-nucleotide polymorphisms (SNPs) in the IL-33 gene have been recognized to be associated with a vast variety of inflammatory disorders. SNPs associated studies have become a crucial approach in uncovering the genetic background of human diseases. However, distinguishing the functional SNPs in a disease-related gene from a pool of both functional and neutral SNPs is a major challenge and needs multiple experiments of hundreds or thousands of SNPs in candidate genes. This study aimed to identify the possible deleterious SNPs in the IL-33 gene using bioinformatics predictive tools. The nonsynonymous SNPs (nsSNPs) were analyzed by SIFT, PolyPhen, PROVEAN, SNP&GO, MutPred, SNAP, PhD SNP, and I-Mutant tools. The Non-coding SNPs (ncSNPs) were also analyzed by SNPinfo and RegulomeDB tools. In conclusion, our in-silico analysis predicted 5 nsSNPs and 22 ncSNPs as potential candidates in the IL-33 gene for future genetic association studies.     

Understanding the adsorption of branched polyamine on surface of gold nanoparticles by molecular dynamics simulations

Molecular dynamics simulations were conducted to characterize the adsorption behavior of branched polyethylenimine (br‐PEI) on the surface of gold nanoparticles (AuNPs). We observed the preferential adsorption of br‐PEI on the [111] surface of AuNPs. Furthermore, br‐PEI maintained a flat arrangement on the surface and wrapped the AuNPs in a stable manner, thereby blocking the adsorption of H₂O molecules and other free br‐PEI molecules. The model and computational results provide theoretical support for relevant experiments. Copyright © 2016 John Wiley & Sons, Ltd.     

Molecular dynamics study of crystalline water ices

The behavior of structures of H₂O crystalline ices Ih, Ic, XI, VII, VIII, VI is studied in molecular dynamics experiment using the potential offered by Poltev and Malenkov. The behavior of the system consisting of one of the two identical interpenetrating, but without any common hydrogen bonds, water frameworks comprising the ice VI structure is also simulated. As a result of simulations, the ice VII structure has collapsed, whereas other systems proved to be stable. The reasons of instability of the ice VII and previously studied ice IV structures in molecular dynamics experiments are discussed. Based on the simulation results of the above-mentioned ices and previous simulation of ices II, III, IX, IV, and XII, the general regularities of dynamic properties of water molecules in crystalline water ices are formulated. Unreliability of results obtained by molecular dynamics in the investigation of self-organizing processes in aqueous systems is shown.     

Structure and dynamics of liquid water with different long-range interaction truncation and temperature control methods in molecular dynamics simulations

We have used molecular dynamics simulations to study the physical properties of modified TIP3P water model included in the CHARMM program, using four different methods-the Ewald summation technique, and three different spherical truncation methods-for the treatment of the long-range interactions. Both the structure and dynamics of the liquid water model were affected by the methods used to truncate the long-range interactions. For some of the methods artificial structuring of the model liquid was observed around the cutoff radius. The model liquid properties were also affected by the commonly applied temperature control methods. Four different methods for controlling the temperature of the system were studied, and the effects of these methods on the bulk properties for liquid water were analyzed. The system size was also found to change the dynamics of the model liquid water. Two control simulations with the SPC/E water model were carried out. The self-diffusion coefficient (D), the radial distribution function (g(OO)), the distance dependent Kirkwood G-factor [G(k)(r)] and the intermolecular potential energy (E(pot)) were determined from the different trajectories and compared with the experimental data.     

Free energy calculations using flexible-constrained, hard-constrained and non-constrained molecular dynamics simulations.

A comparison of different treatments of bond-stretching interactions in molecular dynamics simulation is presented. Relative free energies from simulations using rigid bonds maintained with the SHAKE algorithm, using partially rigid bonds maintained with a recently introduced flexible constraints algorithm, and using fully flexible bonds are compared in a multi-configurational thermodynamic integration calculation of changing liquid water into liquid methanol. The formula for the free energy change due to a changing flexible constraint in a flexible constraint simulation is derived. To allow for a more direct comparison between these three methods, three different pairs of models for water and methanol were used: a flexible model (simulated without constraints and with flexible constraints), a rigid model (simulated with standard hard constraints), and an alternative flexible model (simulated with flexible constraints and standard hard constraints) in which the ideal or constrained bond lengths correspond to the average bond lengths obtained from a short simulation of the unconstrained flexible model. The particular treatment of the bonds induces differences of up to 2 % in the liquid densities, whereas (excess) free energy differences of up to 5.7 (4.3) kJ mol(-1) are observed. These values are smaller than the differences observed between the three different pairs of methanol/water models: up to 5 % in density and up to 8.5 kJ mol(-1) in (excess) free energy.     

Collective motion artifacts arising in long‐duration molecular dynamics simulations

We tested a variety of molecular dynamics simulation strategies in long‐duration (up to several nanoseconds) constant‐temperature simulations of liquid water under periodic boundary conditions. Such long durations are necessary to achieve adequate conformational sampling in simulations of membrane assemblies and other large biomolecular systems. Under a variety of circumstances, serious artifacts arise in the form of spurious collective behavior that becomes obvious only after the simulation has gone at least several hundred picoseconds. The potential energy of the system drops and the system changes from a liquid to an icy or glassy state. The underlying cause is accumulated center‐of‐mass motion of the system, coupled with velocity rescaling associated with constant‐temperature control. The velocity rescaling in the constant‐temperature algorithm reduces the thermal velocity as the net center‐of‐mass velocity grows, effectively causing the kinetic energy of the system to drain from thermal motions into coordinated motions. We found that the incidence and magnitude of the underlying artifactual motion leading to the spurious transition is mediated by: choice of method for computing electrostatic interactions; choice of ensemble; size of the simulation cell; SHAKE tolerance; frequency of nonbonded pairlist updating; and closeness of coupling to the temperature bath. The appearance of the spurious transition can be avoided by periodically subtracting net center‐of‐mass motion during the dynamics, or by improving the accuracy of the simulation by means of tightening SHAKE tolerance and updating nonbonded pairlists every timestep. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 121–131, 2000     

Molecular dynamics study of polymer crystallization in the presence of a particle

We have used molecular dynamics simulations with a coarse‐grained model to study the effect of a particle on the crystallization of polymer melt. We analyzed in particular a bond order parameter to characterize the nucleation and crystallization process. Our calculations show that the presence of a particle modifies the free energy landscape of polymer melts, locally induces the ordering of polymer melts near the particle surface, and thus enhances the polymer crystallization. Because the interaction between the particle and polymers is repulsive, our results suggest that the origin of the enhancement for polymer crystallization is entropic. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 2161–2166, 2007     

The classical molecular dynamics simulation of graphene on Ru(0001) using a fitted Tersoff interface potential

The accurate molecular dynamics simulation of weakly bound adhesive complexes, such as supported graphene (gr), is challenging because of the lack of an adequate interface potential. Instead of the widely used Lennard‐Jones potential for weak and long‐range interactions, we use a newly parameterized Tersoff potential for gr/Ru(0001) system. The new interfacial force field provides adequate moire superstructures in accordance with scanning tunneling microscopy images and with density functional theory (DFT) results. In particular, the corrugation of ξ ≈ 1.0 ± 0.2 Å is found that is somewhat smaller than found by DFT approaches (ξ ≈ 1.2 Å) and is close to scanning tunneling microscope measurements (ξ ≈ 0.8 ± 0.3 Å). The new potential could open the way toward large‐scale simulations of supported gr with adequate moire supercells in many fields of gr research. Moreover, the new interface potential might provide a new strategy in general for obtaining accurate interaction potentials for weakly bound adhesion in large‐scale systems in which atomic dynamics is inaccessible yet by accurate DFT calculations. Copyright © 2013 John Wiley & Sons, Ltd.     

Molecular dynamics study of ferrocene topology under various temperatures

We present the first quantum mechanical Atom-Centered Density-Matrix Propagation molecular dynamic (MD) study to investigate ferrocene (Fc) conformation in gas phase. The MD simulations were performed at several temperatures (7, 18, 80, 120, 180, 293, and 500 K) for a period of 10 ps. It is found that, at very low temperatures (≤18 K), ferrocene prefers eclipsed-like conformation. At higher temperatures (>18 K), the cyclopentadienyl rings (Cp) of ferrocene exhibit apparent fluxional rotations, leading to configurations with the rotational angle δ distributing within a range of 0° (eclipsed) to 18° (approximately half of 36° for the staggered conformation), accompanied by the cyclopentadienyl ring tilt up to approximately 12° at 500 K. The simulated mean inferred (IR) spectrum of ferrocene at 7 K is clearly dominant by a doublet-splitting band of eclipsed-like Fc features in the region of 400 to 600 cm⁻¹, in agreement with previous IR studies. The animation obtained from the MD simulations indicates that, at room temperature, the Fe-C distances in ferrocene are in fact not strictly congruent but 2:2:1-fold.