Impact of glycosylation on stability, structure and unfolding of soybean agglutinin (SBA): an insight from thermal perturbation molecular dynamics simulations

Glycosylation has been recognized as one of the most prevalent and complex post-translational modifications of proteins involving numerous enzymes and substrates. Its effect on the protein conformational transitions is not clearly understood yet. In this study, we have examined the effect of glycosylation on protein stability using molecular dynamics simulation of legume lectin soybean agglutinin (SBA). Its glycosylated moiety consists of high mannose type N-linked glycan (Man(9)GlcNAc(2)). To unveil the structural perturbations during thermal unfolding of these two forms, we have studied and compared them to the experimental results. From the perspective of dynamics, our simulations revealed that the nonglycosylated monomeric form is less stable than corresponding glycosylated form at normal and elevated temperatures. Moreover, at elevated temperature thermal destabilization is more prominent in solvent exposed loops, turns and ends of distinct beta sheets. SBA maintains it folded structure due to some important saltbridges, hydrogen bonds and hydrophobic interactions within the protein. The reducing terminal GlcNAc residues interact with the protein residues VAL161, PRO182 and SER225 via hydrophobic and via hydrogen bonding with ASN 9 and ASN 75. Our simulations also revealed that single glycosylation (ASN75) has no significant effect on corresponding cis peptide angle orientation. This atomistic description might have important implications for understanding the functionality and stability of Soybean agglutinin.

A carboxy terminal domain of the L protein of rinderpest virus possesses RNA triphosphatase activity - The first enzyme in the viral mRNA capping pathway

The large protein L of negative-sense RNA viruses is a multifunctional protein involved in transcription and replication of genomic RNA. It also possesses enzymatic activities involved in capping and methylation of viral mRNAs. The pathway for mRNA capping followed by the L protein of the viruses in the Morbillivirus genus has not been established, although it has been speculated that these viruses may follow the unconventional capping pathway as has been shown for some viruses of Rhabdoviridae family. We had earlier shown that the large protein L of Rinderpest virus expressed as recombinant L-P complex in insect cells as well as the ribonucleoprotein complex from purified virus possesses RNA triphosphatase (RTPase) and guanylyltransferase activities, in addition to RNA dependent RNA polymerase activity. In the present work, we demonstrate that RTPase as well as nucleoside triphosphatase (NTPase) activities are exhibited by a subdomain of the L protein in the C terminal region (a.a. 1640 1840). The RTPase activity depends absolutely on a divalent cation, either magnesium or manganese. Both the RTPase and NTPase activities of the protein show dual metal specificity. Two mutant proteins having alanine mutations in the glutamic acid residues in motif-A of the RTPase domain did not show RTPase activity, while exhibiting reduced NTPase activity suggesting overlapping active sites for the two enzymatic functions. The RTPase and NTPase activities of the L subdomain resemble those of the Vaccinia capping enzyme D1 and the baculovirus LEF4 proteins. (C) 2015 Elsevier Inc. All rights reserved.

Dynamics of history-dependent epidemics in temporal networks

The structural properties of temporal networks often influence the dynamical processes that occur on these networks, e.g., bursty interaction patterns have been shown to slow down epidemics. In this paper, we investigate the effect of link lifetimes on the spread of history-dependent epidemics. We formulate an analytically tractable activity-driven temporal network model that explicitly incorporates link lifetimes. For Markovian link lifetimes, we use mean-field analysis for computing the epidemic threshold, while the effect of non-Markovian link lifetimes is studied using simulations. Furthermore, we also study the effect of negative correlation between the number of links spawned by an individual and the lifetimes of those links. Such negative correlations may arise due to the finite cognitive capacity of the individuals. Our investigations reveal that heavy-tailed link lifetimes slow down the epidemic, while negative correlations can reduce epidemic prevalence. We believe that our results help shed light on the role of link lifetimes in modulating diffusion processes on temporal networks.

Time scales of the stick-slip dynamics of the peeling of an adhesive tape

The stick-slip dynamics of the peeling of an adhesive tape is characterized by bifurcations that have been experimentally well studied. In this work, we investigate the time scale in which the the stick-slips happen leading to the bifurcations. This is fundamental to understanding the triboluminescence and acoustic emissions associated with the bifurcations. We establish a relationship between the time scale of the bifurcations and the inherent mathematical structure of the peeling dynamics by studying a characteristic time quantity associated with the dynamics.

EFFECT OF PLUME DYNAMICS FOR HEAT TRANSFER ENHANCEMENT AT SOLID-LIQUID INTERFACE

The present paper analyzes the effects of plumes for heat transfer enhancement at solid-liquid interface taking both smooth and grooved surfaces. The experimental setup consists of a tank of dimensions 265 x 265 x 300 (height) containing water. The bottom surface was heated and free surface of the water was left open to the ambient. In the experiments, the bottom plate had either a smooth surface or a grooved surface. We used 90 V-grooved rough surfaces with two groove heights, 10mm and 3mm. The experiment was done with water layer depths of 90mm and 140mm, corresponding to values of aspect ratio(AR) equal to 2.9 and 1.8 respectively. Thymol blue, a pH sensitive dye, was used to visualize the flow near the heated plate. The measured heat transfer coefficients over the grooved surfaces were higher compared that over the smooth surface. The enhanced heat transport in the rough cavities cannot be ascribed to the increase in the contact area, rather it must be the local dynamics of the thermal boundary layer that changes the heat transport over the rough surface.

Influence of a Counterion on the Ion Atmosphere of an Anion: A Molecular Dynamics Study of LiX and CsX (X = F-, Cl-, I-) in Methanol

We report molecular dynamics (MD) simulations to explore the influence of a counterion on the structure and dynamics of cationic and anionic solvation shells for various ions in methanol at 298 K. We show that the variation in ionic size of either the cation or the anion in an ion pair influences the solvation structure of the other ion as well as the diffusivity in an electrolyte solution of methanol. The extent of ionic association between the cation and its counteranion of different ionic sizes has been investigated by analyzing the radial distribution functions (RDFs) and the orientation of methanol molecules in the first solvation shell (FSS) of ions. It is shown that the methanol in the FSS of the anion as well the cation exhibit quite different radial and orientational structures as compared to methanol which lie in the FSS of either the anion or the cation but not both. We find that the coordination number (CN) of F-, Cr-, and I- ions decreases with increasing size of the anion which is contrary to the trend reported for the anions in H2O. The mean residence time (MRT) of methanol molecules in the FSS of ions has been calculated using the stable states picture (SSP) approach. It is seen that the ion-counterion interaction has a considerable influence on the MRT of methanol molecules in the FSS of ions. We also discuss the stability order of the ion-counterion using the potentials of mean force (PMFs) for ion pairs with ions of different sizes. The PMF plots reveal that the Li+-F- pair (small-small) is highly stable and the Li+-I- pair is least stable (small-large) in electrolyte solutions.

Ab Initio Molecular Dynamics Simulations of Amino Acids in Aqueous Solutions: Estimating pK(a) Values from Metadynamics Sampling

Changes in the protonation and deprotonation of amino acid residues in proteins play a key role in many biological processes and pathways. Here, we report calculations of the free-energy profile for the protonation deprotonation reaction of the 20 canonical alpha amino acids in aqueous solutions using ab initio Car-Parrinello molecular dynamics simulations coupled with metad-ynamics sampling. We show here that the calculated change in free energy of the dissociation reaction provides estimates of the multiple pK(a) values of the amino acids that are in good agreement with experiment. We use the bond-length-dependent number of the protons coordinated to the hydroxyl oxygen of the carboxylic and the amine groups as the collective variables to explore the free-energy profiles of the Bronsted acid-base chemistry of amino acids in aqueous solutions. We ensure that the amino acid undergoing dissociation is solvated by at least three hydrations shells with all water molecules included in the simulations. The method works equally well for amino acids with neutral, acidic and basic side chains and provides estimates of the multiple pK(a) values with a mean relative error, with respect to experimental results, of 0.2 pK(a) units.

Molecular Dynamics Study of the Structure, Flexibility, and Hydrophilicity of PETIM Dendrimers: A Comparison with PAMAM Dendrimers

A new class of dendrimers, the poly(propyl ether imine) (PETIM) dendrimer, has been shown to be a novel hyperbranched polymer having potential applications as a drug delivery vehicle. Structure and dynamics of the amine terminated PETIM dendrimer and their changes with respect to the dendrimer generation are poorly understood. Since most drugs are hydrophobic in nature, the extent of hydrophobicity of the dendrimer core is related to its drug encapsulation and retention efficacy. In this study, we carry out fully atomistic molecular dynamics (MD) simulations to characterize the structure of PETIM (G2-G6) dendrimers in salt solution as a function of dendrimer generation at different protonation levels. Structural properties such as radius of gyration (R-g), radial density distribution, aspect ratio, and asphericity are calculated. In order to assess the hydrophilicity of the dendrimer, we compute the number of bound water molecules in the interior of dendrirner as well as the number of dendrimer-water hydrogen bonds. We conclude that PETIM dendrimers have relatively greater hydrophobicity and flexibility when compared with their extensively investigated PAMAM counterparts. Hence PETIM dendrimers are expected to have stronger interactions with lipid membranes as well as improved drug encapsulation and retention properties when compared with PAMAM dendrimers. We compute the root-mean-square fluctuation of dendrimers as well as their entropy to quantify the flexibility of the dendrimer. Finally we note that structural and solvation properties computed using force field parameters derived based on the CHARMM general purpose force field were in good quantitative agreement with those obtained using the generalized Amber force field (GAFF).

Disorder, oscillatory dynamics and state switching: the role of c-Myc

In this paper, using the intrinsically disordered oncoprotein Myc as an example, we present a mathematical model to help explain how protein oscillatory dynamics can influence state switching. Earlier studies have demonstrated that, while Myc overexpression can facilitate state switching and transform a normal cell into a cancer phenotype, its downregulation can reverse state-switching. A fundamental aspect of the model is that a Myc threshold determines cell fate in cells expressing p53. We demonstrate that a non-cooperative positive feedback loop coupled with Myc sequestration at multiple binding sites can generate bistable Myc levels. Normal quiescent cells with Myc levels below the threshold can respond to mitogenic signals to activate the cyclin/cdk oscillator for limited cell divisions but the p53/Mdm2 oscillator remains nonfunctional. In response to stress, the p53/Mdm2 oscillator is activated in pulses that are critical to DNA repair. But if stress causes Myc levels to cross the threshold, Myc inactivates the p53/Mdm2 oscillator, abrogates p53 pulses, and pushes the cyclin/cdk oscillator into overdrive sustaining unchecked proliferation seen in cancer. However, if Myc is downregulated, the cyclin/cdk oscillator is inactivated and the p53/Mdm2 oscillator is reset and the cancer phenotype is reversed. (C) 2015 Elsevier Ltd. All rights reserved.

Rational molecular dynamics scheme for predicting optimum concentration loading of nano-additive in phase change materials

The present study deals with the diffusion and phase transition behaviour of paraffin reinforced with carbon nano-additives namely graphene oxide (GO) and surface functionalized single walled carbon nanotubes (SWCNT). Bulk disordered systems of paraffin hydrocarbons impregnated with carbon nano-additives have been generated in realistic equilibrium conformations for potential application as latent heat storage systems. Ab initio molecular dynamics(MD) in conjugation with COMPASS forcefield has been implemented using periodic boundary conditions. The proposed scheme allows determination of optimum nano-additive loading for improving thermo-physical properties through analysis of mass, thermal and transport properties; and assists in determination of composite behaviour and related performance from microscopic point of view. It was observed that nanocomposites containing 7.8% surface functionalised SWCNT and 55% GO loading corresponds to best latent heat storage system. The propounded methodology could serve as a by-pass route for economically taxing and iterative experimental procedures required to attain the optimum composition for best performance. The results also hint at the large unexplored potential of ab-initio classical MD techniques for predicting performance of new nanocomposites for potential phase change material applications. (C) 2015 Author(s).

Estimation of activation energy for electroporation and pore growth rate in liquid crystalline and gel phases of lipid bilayers using molecular dynamics simulations

Molecular dynamics simulations of electroporation in POPC and DPPC lipid bilayers have been carried out at different temperatures ranging from 230 K to 350 K for varying electric fields. The dynamics of pore formation, including threshold field, pore initiation time, pore growth rate, and pore closure rate after the field is switched off, was studied in both the gel and liquid crystalline (L-alpha) phases of the bilayers. Using an Arrhenius model of pore initiation kinetics, the activation energy for pore opening was estimated to be 25.6 kJ mol(-1) and 32.6 kJ mol(-1) in the L-alpha phase of POPC and DPPC lipids respectively at a field strength of 0.32 V nm(-1). The activation energy decreases to 24.2 kJ mol(-1) and 23.7 kJ mol(-1) respectively at a higher field strength of 1.1 V nm(-1). At temperatures below the melting point, the activation energy in the gel phase of POPC and DPPC increases to 28.8 kJ mol(-1) and 34.4 kJ mol(-1) respectively at the same field of 1.1 V nm(-1). The pore closing time was found to be higher in the gel than in the L-alpha phase. The pore growth rate increases linearly with temperature and quadratically with field, consistent with viscosity limited growth models.

Kinematics and dynamics of Jansen leg mechanism: A bond graph approach

The Jansen mechanism is a one degree-of-freedom, planar, 12-link, leg mechanism that can be used in mobile robotic applications and in gait analysis. This paper presents the kinematics and dynamics of the Jansen leg mechanism. The forward kinematics, accomplished using circle intersection method, determines the trajectories of various points on the mechanism in the chassis (stationary link) reference frame. From the foot point trajectory, the step length is shown to vary linearly while step height varies non-linearly with change in crank radius. A dynamic model for the Jansen leg mechanism is proposed using bond graph approach with modulated multiport transformers. For given ground reaction force pattern and crank angular speed, this model helps determine the motor torque profile as well as the link and joint stresses. The model can therefore be used to rate the actuator torque and in design of the hardware and controller for such a system. The kinematics of the mechanism can also be obtained from this dynamic model. The proposed model is thus a useful tool for analysis and design of systems based on the Jansen leg mechanism. (C) 2015 Elsevier B.V. All rights reserved.

Ultrafast Chemical Dynamics in Time Domain Through Fluorescence Spectroscopy

While absorption and emission spectroscopy have always been used to detect and characterize molecules and molecular complexes, the availability of ultrashort laser pulses and associated computer-aided optical detection techniques allowed study of chemical processes directly in the time domain at unprecedented time scales, through appearance and disappearance of fluorescence from participating chemical species. Application of such techniques to chemical dynamics in liquids, where many processes occur with picosecond and femtosecond time scales lead to the discovery of a host of new phenomena that in turn led to the development of many new theories. Experiment and theory together provided new and valuable insight into many fundamental chemical processes, like isomerization dynamics, electron and proton transfer reactions, vibrational energy and phase relaxation, photosynthesis, to name just a few. In this article, we shall review a few of such discoveries in attempt to provide a glimpse of the fascinating research employing fluorescence spectroscopy that changed the field of chemical dynamics forever.

Nonlinear Dynamics of a Rotating Flexible Link

This paper deals with the study of the nonlinear dynamics of a rotating flexible link modeled as a one dimensional beam, undergoing large deformation and with geometric nonlinearities. The partial differential equation of motion is discretized using a finite element approach to yield four nonlinear, nonautonomous and coupled ordinary differential equations (ODEs). The equations are nondimensionalized using two characteristic velocities-the speed of sound in the material and a velocity associated with the transverse bending vibration of the beam. The method of multiple scales is used to perform a detailed study of the system. A set of four autonomous equations of the first-order are derived considering primary resonances of the external excitation and one-to-one internal resonances between the natural frequencies of the equations. Numerical simulations show that for certain ranges of values of these characteristic velocities, the slow flow equations can exhibit chaotic motions. The numerical simulations and the results are related to a rotating wind turbine blade and the approach can be used for the study of the nonlinear dynamics of a single link flexible manipulator.