Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

In addition to the biologically active monomer of the protein insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The insulin monomer contains two distinct surfaces, namely, the dimer forming surface (DFS) and the hexamer forming surface (HFS), that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behavior of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. In particular, we find that more structured water molecules, with higher residence times (similar to 300-500 ps), are present near HFS than those near DFS. A significant slowing down is observed in the decay of associated rotational auto time correlation functions of O-H bond vector of water in the vicinity of HFS. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly an outward H-atoms directed, clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inward H-atoms directed orientation owing to the flat curvature of hydrophobic surface and the interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces that reveal the significant differences between the two hydration layers.

Conserved water-mediated recognition and dynamics of NAD(+) (carboxamide group) to hIMPDH enzyme: water mimic approach toward the design of isoform-selective inhibitor

Inosine monophosphate dehydrogenase (IMPDH) enzyme involves in GMP biosynthesis pathway. Type I hIMPDH is expressed at lower levels in all cells, whereas type II is especially observed in acute myelogenous leukemia, chronic myelogenous leukemia cancer cells, and 10 ns simulation of the IMP-NAD(+) complex structures (PDB ID. 1B3O and 1JCN) have revealed the presence of a few conserved hydrophilic centers near carboxamide group of NAD(+). Three conserved water molecules (W1, W, and W1 `) in di-nucleotide binding pocket of enzyme have played a significant role in the recognition of carboxamide group (of NAD(+)) to D274 and H93 residues. Based on H-bonding interaction of conserved hydrophilic (water molecular) centers within IMP-NAD(+)-enzyme complexes and their recognition to NAD(+), some covalent modification at carboxamide group of di-nucleotide (NAD(+)) has been made by substituting the -CONH(2)group by -CONHNH2 (carboxyl hydrazide group) using water mimic inhibitor design protocol. The modeled structure of modified ligand may, though, be useful for the development of antileukemic agent or it could be act as better inhibitor for hIMPDH-II.

Correlation between local structural dynamics of proteins inferred from NMR ensembles and evolutionary dynamics of homologues of known structure

Conformational changes in proteins are extremely important for their biochemical functions. Correlation between inherent conformational variations in a protein and conformational differences in its homologues of known structure is still unclear. In this study, we have used a structural alphabet called Protein Blocks (PBs). PBs are used to perform abstraction of protein 3-D structures into a 1-D strings of 16 alphabets (a-p) based on dihedral angles of overlapping pentapeptides. We have analyzed the variations in local conformations in terms of PBs represented in the ensembles of 801 protein structures determined using NMR spectroscopy. In the analysis of concatenated data over all the residues in all the NMR ensembles, we observe that the overall nature of inherent local structural variations in NMR ensembles is similar to the nature of local structural differences in homologous proteins with a high correlation coefficient of .94. High correlation at the alignment positions corresponding to helical and beta-sheet regions is only expected. However, the correlation coefficient by considering only the loop regions is also quite high (.91). Surprisingly, segregated position-wise analysis shows that this high correlation does not hold true to loop regions at the structurally equivalent positions in NMR ensembles and their homologues of known structure. This suggests that the general nature of local structural changes is unique; however most of the local structural variations in loop regions of NMR ensembles do not correlate to their local structural differences at structurally equivalent positions in homologues.

Dynamics of a driven magneto-martensitic ribbon

We investigate the dynamics of a sinusoidally driven ferromagnetic martensitic ribbon by adopting a recently introduced model that involves strain and magnetization as order parameters. Retaining only the dominant mode of excitation we reduce the coupled set of partial differential equations for strain and magnetization to a set of coupled ordinary nonlinear equations for the strain and magnetization amplitudes. The equation for the strain amplitude takes the form of parametrically driven oscillator. Finite strain amplitude can only be induced beyond a critical value of the strength of the magnetic field. Chaotic response is seen for a range of values of all the physically interesting parameters. The nature of the bifurcations depends on the choice of temperature relative to the ordering of the Curie and the martensite transformation temperatures. We have studied the nature of response as a function of the strength and frequency of the magnetic field, and magneto-elastic coupling. In general, the bifurcation diagrams with respect to these parameters do not follow any standard route. The rich dynamics exhibited by the model is further illustrated by the presence of mixed mode oscillations seen for low frequencies. The geometric structure of the mixed mode oscillations in the phase space has an unusual deep crater structure with an outer and inner cone on which the orbits circulate. We suggest that these features should be seen in experiments on driven magneto-martensitic ribbons. (C) 2014 Elsevier B. V. All rights reserved.

Non-equilibrium segmental dynamics driven by multiwall carbon nanotubes in PS/PVME blends

The present paper discusses the effect of multiwall carbon nanotubes (MWNTs) on the structural relaxation and the intermolecular cooperativity in dynamically asymmetric blends of PS/PVME (polystyrene/poly(vinyl methyl ether)). The temperature regime where chain connectivity effects dominate the thermodynamic concentration fluctuation (T/T-g > 0.75, T-g is the glass transition temperature of the blends) was studied using dielectric spectroscopy (DS). Interestingly, in the blends with MWNTs a bimodal distribution of relaxation was obtained in the loss modulus spectra. This plausibly is due to different environments experienced by the faster component (PVME) in the presence of MWNTs. The segmental dynamics of PVME was observed to be significantly slowed down in the presence of MWNTs and an Arrhenius-type behavior, weakly dependent on temperature, is observed at higher frequencies. This non-equilibrium dynamics of PVME is presumed to be originating from interphase regions near the surface of MWNTs. The length scale of the cooperative rearranging region (xi CRR) at T-g, assessed by calorimetric measurements, was observed to be higher in the case of blends with MWNTs. An enhanced molecular level miscibility driven by MWNTs in the blends corroborates with the larger xi CRR and comparatively more number of segments in CRR (in contrast to neat blends) around T-g. The configurational entropy and length scale of the cooperative volume was mapped as a function of temperature in the temperature regime, Tg < T < T-g + 60 K. The blends phase separated by spinodal decomposition which further led to an interconnected PVME network in PS. This further led to materials with very high electrical conductivity upon demixing.

Anomalous power law decay in solvation dynamics of DNA: a mode coupling theory analysis of ion contribution

Several time dependent fluorescence Stokes shift (TDFSS) experiments have reported a slow power law decay in the hydration dynamics of a DNA molecule. Such a power law has neither been observed in computer simulations nor in some other TDFSS experiments. Here we observe that a slow decay may originate from collective ion contribution because in experiments DNA is immersed in a buffer solution, and also from groove bound water and lastly from DNA dynamics itself. In this work we first express the solvation time correlation function in terms of dynamic structure factors of the solution. We use mode coupling theory to calculate analytically the time dependence of collective ionic contribution. A power law decay in seen to originate from an interplay between long-range probe-ion direct correlation function and ion-ion dynamic structure factor. Although the power law decay is reminiscent of Debye-Falkenhagen effect, yet solvation dynamics is dominated by ion atmosphere relaxation times at longer length scales (small wave number) than in electrolyte friction. We further discuss why this power law may not originate from water motions which have been computed by molecular dynamics simulations. Finally, we propose several experiments to check the prediction of the present theoretical work.

Effect of angle of incidence on mixed convective wake dynamics and heat transfer past a square cylinder in cross flow at Re=100

In this paper, a numerical investigation is performed to study the mixed convective flow and heat transfer characteristics past a square cylinder in cross flow at incidence. Utilizing air (Pr = 0.71) as an operating fluid, computations are carried out at a representative Reynolds number (Re) of 100. Angles of incidences are varied as, 0 degrees <= alpha <= 45 degrees. Effect of superimposed positive and negative cross-flow buoyancy is brought about by varying the Richardson number (RI) in the range -1.0 <= Ri <= 1.0. The detail features of flow topology and heat transport are analyzed critically for different angles of incidences. The thermo fluidic forces acting on the cylinder during mixed convection are captured in terms of the drag (C-D), lift (C-L), and moment (C-M) coefficients. The results show that the lateral width of the cylinder wake reduces with increasing alpha and the isotherms spread out far wide. In the range 0 degrees < alpha < 45 degrees, C-D reduces with increasing Ri. The functional dependence of C-M with Ri reveals a linear relationship. Thermal boundary layer thickness reduces with increasing angle of incidences. The global rate of heat transfer from the cylinder increases with increasing alpha. (C) 2014 Elsevier Ltd. All rights reserved.

Energy funneling and macromolecular conformational dynamics: a 2D Raman correlation study of PEG melting

Cascading energy landscapes through funneling has been postulated as a mechanistic route for achieving the lowest energy configuration of a macromolecular system (such as proteins and polymers). In particular, understanding the molecular mechanism for the melting and crystallization of polymers is a challenging fundamental question. The structural modifications that lead to the melting of poly(ethylene glycol) (PEG) are investigated here. Specific Raman bands corresponding to different configurations of the PEG chain have been identified, and the molecular structural dynamics of PEG melting have been addressed using a combination of Raman spectroscopy, 2D Raman correlation and density functional theory (DFT) calculations. The melting dynamics of PEG have been unambiguously explained along the C-O bond rotation coordinate.

Length-scale dependent mechanical properties of Al-Cu eutectic alloy: Molecular dynamics based model and its experimental verification

This paper attempts to gain an understanding of the effect of lamellar length scale on the mechanical properties of two-phase metal-intermetallic eutectic structure. We first develop a molecular dynamics model for the in-situ grown eutectic interface followed by a model of deformation of Al-Al2Cu lamellar eutectic. Leveraging the insights obtained from the simulation on the behaviour of dislocations at different length scales of the eutectic, we present and explain the experimental results on Al-Al2Cu eutectic with various different lamellar spacing. The physics behind the mechanism is further quantified with help of atomic level energy model for different length scale as well as different strain. An atomic level energy partitioning of the lamellae and the interface regions reveals that the energy of the lamellae core are accumulated more due to dislocations irrespective of the length-scale. Whereas the energy of the interface is accumulated more due to dislocations when the length-scale is smaller, but the trend is reversed when the length-scale is large beyond a critical size of about 80 nm. (C) 2014 Author(s).

Dynamics of the reaction between the free end of a tethered self-avoiding polymer and a flat penetrable surface: A renormalization group study

The average time tau(r) for one end of a long, self-avoiding polymer to interact for the first time with a flat penetrable surface to which it is attached at the other end is shown here to scale essentially as the square of the chain's contour length N. This result is obtained within the framework of the Wilemski-Fixman approximation to diffusion-limited reactions, in which the reaction time is expressed as a time correlation function of a ``sink'' term. In the present work, this sink-sink correlation function is calculated using perturbation expansions in the excluded volume and the polymer-surface interactions, with renormalization group methods being used to resum the expansion into a power law form. The quadratic dependence of tau(r) on N mirrors the behavior of the average time tau(c) of a free random walk to cyclize, but contrasts with the cyclization time of a free self-avoiding walk (SAW), for which tau(r) similar to N-2.2. A simulation study by Cheng and Makarov J. Phys. Chem. B 114, 3321 (2010)] of the chain-end reaction time of an SAW on a flat impenetrable surface leads to the same N-2.2 behavior, which is surprising given the reduced conformational space a tethered polymer has to explore in order to react. (C) 2014 AIP Publishing LLC.

Allostery and Conformational Dynamics in cAMP-binding Acyltransferases

Mycobacteria harbor unique proteins that regulate protein lysine acylation in a cAMP-regulated manner. These lysine acyltransferases from Mycobacterium smegmatis (KATms) and Mycobacterium tuberculosis (KATmt) show distinctive biochemical properties in terms of cAMP binding affinity to the N-terminal cyclic nucleotide binding domain and allosteric activation of the C-terminal acyltransferase domain. Here we provide evidence for structural features in KATms that account for high affinity cAMP binding and elevated acyltransferase activity in the absence of cAMP. Structure-guided mutational analysis converted KATms from a cAMP-regulated to a cAMP-dependent acyltransferase and identified a unique asparagine residue in the acyltransferase domain of KATms that assists in the enzymatic reaction in the absence of a highly conserved glutamate residue seen in Gcn5-related N-acetyltransferase-like acyltransferases. Thus, we have identified mechanisms by which properties of similar proteins have diverged in two species of mycobacteria by modifications in amino acid sequence, which can dramatically alter the abundance of conformational states adopted by a protein.

Topological blocking in quantum quench dynamics

We study the nonequilibrium dynamics of quenching through a quantum critical point in topological systems, focusing on one of their defining features: ground-state degeneracies and associated topological sectors. We present the notion of ``topological blocking,'' experienced by the dynamics due to a mismatch in degeneracies between two phases, and we argue that the dynamic evolution of the quench depends strongly on the topological sector being probed. We demonstrate this interplay between quench and topology in models stemming from two extensively studied systems, the transverse Ising chain and the Kitaev honeycomb model. Through nonlocal maps of each of these systems, we effectively study spinless fermionic p-wave paired topological superconductors. Confining the systems to ring and toroidal geometries, respectively, enables us to cleanly address degeneracies, subtle issues of fermion occupation and parity, and mismatches between topological sectors. We show that various features of the quench, which are related to Kibble-Zurek physics, are sensitive to the topological sector being probed, in particular, the overlap between the time-evolved initial ground state and an appropriate low-energy state of the final Hamiltonian. While most of our study is confined to translationally invariant systems, where momentum is a convenient quantum number, we briefly consider the effect of disorder and illustrate how this can influence the quench in a qualitatively different way depending on the topological sector considered.

Fish distribution dynamics in the Aghanashini estuary of Uttara Kannada, west coast of India

Fish diversity (77 species) in the Aghanashini River estuary of the Indian west coast is linked to variable salinity conditions and zones I, II and III for high, medium and low salinity respectively. Zone I, the junction between Arabian Sea and the estuary, had all species in yearly succession due to freshwater conditions in monsoon to high salinity in pre-monsoon. The medium (zone II) and low (zone III) salinity mid and upstream portions had maximum of 67 and 39 fish species respectively. Maintenance of natural salinity regimes in estuary, among other ecological factors, is critical for its fish diversity.