Ananya Debnath

Entropy and dynamics of water in hydration layers of a bilayer

We compute the entropy and transport properties of water in the hydration layer of dipalmitoylphosphatidylcholine bilayer by using a recently developed theoretical scheme [two-phase thermodynamic model, termed as 2PT method; S.-T. Lin et al., J. Chem. Phys. 119, 11792 (2003)] based on the translational and rotational velocity autocorrelation functions and their power spectra. The weights of translational and rotational power spectra shift from higher to lower frequency as one goes from the bilayer interface to the bulk. Water molecules near the bilayer head groups have substantially lower entropy (48.36 J/mol/K) than water molecules in the intermediate region (51.36 J/mol/K), which have again lower entropy than the molecules (60.52 J/mol/K) in bulk. Thus, the entropic contribution to the free energy change (TΔS) of transferring an interface water molecule to the bulk is 3.65 kJ/mol and of transferring intermediate water to the bulk is 2.75 kJ/mol at 300 K, which is to be compared with 6.03 kJ/mol for melting of ice at 273 K. The translational diffusion of water in the vicinity of the head groups is found to be in a subdiffusive regime and the rotational diffusion constant increases going away from the interface. This behavior is supported by the slower reorientational relaxation of the dipole vector and OH bond vector of interfacial water. The ratio of reorientational relaxation time for Legendre polynomials of order 1 and 2 is approximately 2 for interface, intermediate, and bulk water, indicating the presence of jump dynamics in these water molecules.

Structure and Dynamics of Phospholipid Nanodiscs from All-Atom and Coarse-Grained Simulations.

We investigated structural and dynamical properties of nanodiscs comprising dimyristoylphosphatidylcholine (DMPC) lipids and major scaffold protein MSP1Δ(1-22) from human apolipoprotein A-1 using combined all-atom and coarse-grained (CG) molecular dynamics (MD) simulations. The computational efficiency of the Martini-CG force field enables the spontaneous self-assembly of lipids and scaffold proteins into stable nanodisc structures on time scales up to tens of microseconds. Subsequent all-atom and CG-MD simulations reveal that the lipids in the nanodisc have lower configurational entropy and higher acyl tail order than in a lamellar bilayer phase. These altered average properties arise from rather differential behavior of lipids, depending on their location in the nanodisc. Since the scaffold proteins exert constrictive forces from the outer rim of the disc toward its center, lipids at the center of the nanodisc are highly ordered, whereas annular lipids that are in contact with the MSP proteins are remarkably disordered due to perturbed packing. Although specific differences between all-atom and CG simulations are also evident, the results obtained at both levels of resolution are in overall good agreement with each other and provide atomic level interpretations of recent experiments. Thus, the present study highlights the applicability of multiscale simulation approaches for nanodisc systems and opens the way for future applications, including the study of nanodisc-embedded membrane proteins.

Rate processes with dynamical disorder: A direct variational approach

Using path integral approach, we develop variational approximations to the calculation of survival probability for rate processes with dynamical disorder. We derive both upper and lower bounds to the survival probability using Jensens inequality. The inequalities involve the use of a trial action for which the path integrals can be evaluated exactly. Any parameter in the trial action can be varied to optimize the bounds. We have also derived a lower bound to the rate of the process. As a simple illustration, we apply the method to the problem of a particle undergoing Brownian motion in a harmonic potential well, in the presence of a delta function sink, for which one can calculate the exact survival probability numerically. The calculation confirms the two inequalities. The method should be very useful in similar but more complex problems where even numerical solution is not possible.

Polymer in a double well: dynamics of translocation of short chains over a barrier

We consider the dynamics of a short chain polymer crossing over a free energy barrier in space. Adopting the continuum version of the Rouse model, we find exact expressions for the activation energy and the rate of crossing. For this model, the analysis of barrier crossing is analogous to semiclassical treatment of quantum tunnelling. Finding the saddle point for the process requires solving a Newton-like equation of motion for a fictitious particle. The analysis shows that short chains would cross the barrier as a globule. The activation free energy for this would increase linearly with the number of units N in the polymer. The saddle point for longer chains is an extended conformation, in which the chain is stretched out. The stretching out lowers the energy and hence the activation free energy is no longer linear in N. The rates in both the cases are calculated using a multidimensional approach and analytical expressions are derived, using a new formula for evaluating the infinite products. However, due to the harmonic approximation made in the derivation, the rates are found to diverge at the point where the saddle point changes over from the globule to the stretched out conformation. The reason for this is identified to be the bifurcation of the saddle to give two new saddles, and a correction formula is derived for the rate in the vicinity of this point. Numerical results using the formulae are presented. As a function of N, it is possible for the rate to have a minimum. This is due to confinement effects in the initial state.

Solvent Assisted Tuning of Morphology of a Peptide-Perylenediimide Conjugate: Helical Fibers to Nano-Rings and their Differential Semiconductivity

Understanding the regulatory factors of self-assembly processes is a necessity in order to modulate the nano-structures and their properties. Here, the self-assembly mechanism of a peptide-perylenediimide (P-1) conjugate in mixed solvent systems of THF/water is studied and the semiconducting properties are correlated with the morphology. In THF, right handed helical fibers are formed while in 10% THF-water, the morphology changes to nano-rings along with a switch in the helicity to left-handed orientation. Experimental results combined with DFT calculations reveal the critical role of thermodynamic and kinetic factors to control these differential self-assembly processes. In THF, P-1 forms right handed helical fibers in a kinetically controlled fashion. In case of 10% THF-water, the initial nucleation of the aggregate is controlled kinetically. Due to differential solubility of the molecule in these two solvents, elongation of the nuclei into fibers is restricted after a critical length leading to the formation of nano-rings which is governed by the thermodynamics. The helical fibers show superior semi-conducting property to the nano-rings as confirmed by conducting-AFM and conventional I-V characteristics.

Derivation of coarse-grained simulation models of chlorophyll molecules in lipid bilayers for applications in light harvesting systems

The correct interplay of interactions between protein, pigment and lipid molecules is highly relevant for our understanding of the association behavior of the light harvesting complex (LHCII) of green plants. To cover the relevant time and length scales in this multicomponent system, a multi-scale simulation ansatz is employed that subsequently uses a classical all atomistic (AA) model to derive a suitable coarse grained (CG) model which can be backmapped into the AA resolution, aiming for a seamless conversion between two scales. Such an approach requires a faithful description of not only the protein and lipid components, but also the interaction functions for the indispensable pigment molecules, chlorophyll b and chlorophyll a (referred to as chl b/chl a). In this paper we develop a CG model for chl b and chl a in a dipalmitoylphosphatidyl choline (DPPC) bilayer system. The structural properties and the distribution behavior of chl within the lipid bilayer in the CG simulations are consistent with those of AA reference simulations. The non-bonded potentials are parameterized such that they fit to the thermodynamics based MARTINI force-field for the lipid bilayer and the protein. The CG simulation shows chl aggregation in the lipid bilayer which is supported by fluorescence quenching experiments. It is shown that the derived chl model is well suited for CG simulations of stable, structurally consistent, trimeric LHCII and can in the future be used to study their large scale aggregation behavior.

Trigonella seed extract ameliorates inflammation via regulation of the inflammasome adaptor protein ASC

Trigonella foenum-graecum (fenugreek) is an important medicinal plant, well known for its anti-inflammatory properties. However, the underlying cellular and molecular mechanisms of its action remain largely unknown. The apoptosis associated speck like protein containing a caspase recruitment domain (CARD) (ASC) is central to inflammatory and cell death pathways in innate and adaptive immunity. Here, we show that fenugreek seed extract provides cytoprotection to bacterial lipopolysaccharide (LPS) inflammed and nanosilica-treated fibroblasts via a reactive oxygen species independent pathway. All atom molecular dynamics simulations of ASC-ligand complex reveal that individual phytochemicals in fenugreek can bind to ASC via specific non-covalent interactions. These data show that a synergistic effect of fenugreek phytochemicals with the ASC protein alters its molecular properties resulting in altered cellular function. Such information is crucial to the development of targeted therapeutic interventions for inflammatory diseases.

Diffusion in an elastic medium: A model for macromolecule transport across the nuclear pore complex

Nuclear pore complexes (NPCs) are very selective filters that monitor the transport between the cytoplasm and the nucleoplasm. Two models have been suggested for the plug of the NPC. They are (i) it is a reversible hydrogel or (ii) it is a polymer brush. We propose a mesoscopic model for the transport of a protein through the plug, that is general enough to cover both. The protein stretches the plug and creates a local deformation. The bubble so created (prtoein+deformation) executes random walk in the plug. We find that for faster relaxation of the gel, the diffusion of the bubble is greater. Further, on using parameters appropriate for the brush, we find that the diffusion coefficient is much lower. Hence the gel model seems to be more likely explanation for the workings of the plug.

Hydration dynamics of a lipid membrane: Hydrogen bond networks and lipid-lipid associations

Dynamics of hydration layers of a dimyristoylphosphatidylcholine (DMPC) bilayer are investigated using an all atom molecular dynamics simulation. Based upon the geometric criteria, continuously residing interface water molecules which form hydrogen bonds solely among themselves and then concertedly hydrogen bonded to carbonyl, phosphate, and glycerol head groups of DMPC are identified. The interface water hydrogen bonded to lipids shows slower relaxation rates for translational and rotational dynamics compared to that of the bulk water and is found to follow sub-diffusive and non-diffusive behaviors, respectively. The mean square displacements and the reorientational auto-correlation functions are slowest for the interfacial waters hydrogen bonded to the carbonyl oxygen since these are buried deep in the hydrophobic core among all interfacial water studied. The intermittent hydrogen bond auto-correlation functions are calculated, which allows breaking and reformations of the hydrogen bonds. The auto-correlation functions for interfacial hydrogen bonded networks develop humps during a transition from cage-like motion to eventual power law behavior of t−3/2. The asymptotic t−3/2 behavior indicates translational diffusion dictated dynamics during hydrogen bond breaking and formation irrespective of the nature of the chemical confinement. Employing reactive flux correlation analysis, the forward rate constant of hydrogen bond breaking and formation is calculated which is used to obtain Gibbs energy of activation of the hydrogen bond breaking. The relaxation rates of the networks buried in the hydrophobic core are slower than the networks near the lipid-water interface which is again slower than bulk due to the higher Gibbs energy of activation. Since hydrogen bond breakage follows a translational diffusion dictated mechanism, chemically confined hydrogen bond networks need an activation energy to diffuse through water depleted hydrophobic environments. Our calculations reveal that the slow relaxation rates of interfacial waters in the vicinity of lipids are originated from the chemical confinement of concerted hydrogen bond networks. The analysis suggests that the networks in the hydration layer of membranes dynamically facilitate the water mediated lipid-lipid associations which can provide insights on the thermodynamic stability of soft interfaces relevant to biological systems in the future.Dynamics of hydration layers of a dimyristoylphosphatidylcholine (DMPC) bilayer are investigated using an all atom molecular dynamics simulation. Based upon the geometric criteria, continuously residing interface water molecules which form hydrogen bonds solely among themselves and then concertedly hydrogen bonded to carbonyl, phosphate, and glycerol head groups of DMPC are identified. The interface water hydrogen bonded to lipids shows slower relaxation rates for translational and rotational dynamics compared to that of the bulk water and is found to follow sub-diffusive and non-diffusive behaviors, respectively. The mean square displacements and the reorientational auto-correlation functions are slowest for the interfacial waters hydrogen bonded to the carbonyl oxygen since these are buried deep in the hydrophobic core among all interfacial water studied. The intermittent hydrogen bond auto-correlation functions are calculated, which allows breaking and reformations of the hydrogen bonds. The auto-corre...