Solvent Sensitivity of Protein Unfolding: Dynamical Study of Chicken Villin Headpiece Subdomain in Water Ethanol Binary Mixture

We carry out a series of long atomistic molecular dynamics simulations to study the unfolding of a small protein, chicken villin headpiece (HP-36), in water-ethanol (EtOH) binary mixture. The prime objective of this work is to explore the sensitivity of protein unfolding dynamics toward increasing concentration of the cosolvent and unravel essential features of intermediates formed in search of a dynamical pathway toward unfolding. In water ethanol binary mixtures, HP-36 is found to unfold partially, under ambient conditions, that otherwise requires temperature as high as similar to 600 K to denature in pure aqueous solvent. However, an interesting course of pathway is observed to be followed in the process, guided by the formation of unique intermediates. The first step of unfolding is essentially the separation of the cluster formed by three hydrophobic (phenylalanine) residues, namely, Phe-7, Phe-11, and Phe-18, which constitute the hydrophobic core, thereby initiating melting of helix-2 of the protein. The initial steps are similar to temperature-induced unfolding as well as chemical unfolding using DMSO as cosolvent. Subsequent unfolding steps follow a unique path. As water-ethanol shows composition-dependent anomalies, so do the details of unfolding dynamics. With an increase in cosolvent concentration, different partially unfolded intermediates are found to be formed. This is reflected in a remarkable nonmonotonic composition dependence of several order parameters, including fraction of native contacts and protein-solvent interaction energy. The emergence of such partially unfolded states can be attributed to the preferential solvation of the hydrophobic residues by the ethyl groups of ethanol. We further quantify the local dynamics of unfolding by using a Marcus-type theory.

A new microscopic insight into membrane penetration and reorganization by PETIM dendrimers

Dendrimers are highly branched polymeric nanoparticles whose structure and topology, largely, have determined their efficacy in a wide range of studies performed so far. An area of immense interest is their potential as drug and gene delivery vectors. Realizing this potential, depending on the nature of cell surface-dendrimer interactions, here we report controlled model membrane penetration and reorganization, using a model supported lipid bilayer and poly(ether imine) (PETIM) dendrimers of two generations. By systematically varying the areal density of the lipid bilayers, we provide a microscopic insight, through a combination of high resolution scattering, atomic force microscopy and atomistic molecular dynamics simulations, into the mechanism of PETIM dendrimer membrane penetration, pore formation and membrane re-organization induced by such interactions. Our work represents the first systematic observation of a regular barrel-like membrane spanning pore formation by dendrimers, tunable through lipid bilayer packing, without membrane disruption.

Structure, stability and elasticity of DNA nanotubes

DNA nanotubes are tubular structures composed of DNA crossover molecules. We present a bottom up approach for the construction and characterization of these structures. Various possible topologies of nanotubes are constructed such as 6-helix, 8-helix and tri-tubes with different sequences and lengths. We have used fully atomistic molecular dynamics simulations to study the structure, stability and elasticity of these structures. Several nanosecond long MD simulations give the microscopic details about DNA nanotubes. Based on the structural analysis of simulation data, we show that 6-helix nanotubes are stable and maintain their tubular structure; while 8-helix nanotubes are flattened to stabilize themselves. We also comment on the sequence dependence and the effect of overhangs. These structures are approximately four times more rigid having a stretch modulus of similar to 4000 pN compared to the stretch modulus of 1000 pN of a DNA double helix molecule of the same length and sequence. The stretch moduli of these nanotubes are also three times larger than those of PX/JX crossover DNA molecules which have stretch moduli in the range of 1500-2000 pN. The calculated persistence length is in the range of a few microns which is close to the reported experimental results on certain classes of DNA nanotubes.

Composition dependent multiple structural transformations of myoglobin in aqueous ethanol solution: A combined experimental and theoretical study

Experimental studies (circular dichroism and ultra-violet (UV) absorption spectra) and large scale atomistic molecular dynamics simulations (accompanied by order parameter analyses) are combined to establish a number of remarkable (and unforeseen) structural transformations of protein myoglobin in aqueous ethanol mixture at various ethanol concentrations. The following results are particularly striking. (1) Two well-defined structural regimes, one at x(EtOH) similar to 0.05 and the other at x(EtOH) similar to 0.25, characterized by formation of distinct partially folded conformations and separated by a unique partially unfolded intermediate state at x(EtOH) similar to 0.15, are identified. (2) Existence of non-monotonic composition dependence of (i) radius of gyration, (ii) long range contact order, (iii) residue specific solvent accessible surface area of tryptophan, and (iv) circular dichroism spectra and UV-absorption peaks are observed. Interestingly at x(EtOH) similar to 0.15, time averaged value of the contact order parameter of the protein reaches a minimum, implying that this conformational state can be identified as a molten globule state. Multiple structural transformations well known in water-ethanol binary mixture appear to have considerably stronger effects on conformation and dynamics of the protein. We compare the present results with studies in water-dimethyl sulfoxide mixture where also distinct structural transformations are observed along with variation of co-solvent composition. (C) 2015 AIP Publishing LLC.

Effect of Alkyl Chain Arrangement on Conformation and Dynamics in a Surfactant Intercalated Layered Double Hydroxide: Spectroscopic Measurements and MD Simulations

The anionic surfactant dodecyl sulfate (DDS) has been intercalated in an Mg-Al layered double hydroxide (LDH). Monolayer and bilayer arrangements of the alkyl chains of the intercalated surfactant can be engineered by tuning the Al/Mg ratio of the LDH. In both arrangements the anionic headgroup of the surfactant is tethered to the LDH sheets, and consequently translational mobility of the chains is absent. The degrees of freedom of the confined alkyl chains are restricted to changes in conformation. The effects of the arrangement of the intercalated surfactant chains on conformational order and dynamics have been,investigated by spectroscopic measurements and molecular dynamics simulations. Infrared, Raman, and C-13 NMR spectroscopies were used to investigate conformation of the alkyl chains in the monolayer and bilayer arrangements and variable contact time cross-polarization magic angle spinning (VCT CPMAS) NMR measurements to probe molecular motion. The alkyl chains in the monolayer arrangement of the intercalated DDS chains showed considerably greater conformational disorder and faster dynamics as compared to chains in the bilayer arrangement, in spite of the fact that the volume available per chain in the monolayer is smaller than that in the bilayer. Atomistic MD simulations of the two arrangements of the intercalated surfactant were carried out using an isothermal-isobaric ensemble. The simulations are able to reproduce the essential results of the experiment-greater conformational disorder and faster dynamics for the alkyl chains in the monolayer arrangement of the intercalated surfactant. The MD simulations show that these results are a consequence of the fact that the nature of conformational disorder in the two arrangements is different. In the monolayer arrangement the alkyl chains can sustain isolated gauche defects, whereas in the bilayer arrangement gauche conformers occur only as part of a kink a gauche(+) trans gauche(-) sequence.

Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient. (c) 2015 AIP Publishing LLC.

Conformational Dynamics in Alkyl Chains of an Anchored Bilayer: A Molecular Dynamics Study

Molecular dynamics (MD) simulations are reported for an anchored bilayer formed by the intercalation of cetyl trimethyl ammonium (CTA) and CH3(CH2)15N+(CH3) ions in a layered solid, CdPS3. The intercalated CTA ions are organized with the cationic headgroups tethered to the inorganic sheet and the hydrocarbon tails arranged as bilayers. Simulations were performed at three temperatures, 65, 180, and 298 K, using an isothermal−isobaric ensemble that was subsequently switched once macroscopic parameters had converged to a canonical isothermal−isochoric ensemble. The simulations are able to reproduce the experimental features of this system, including the formation of the bilayer and layer-to-layer separation distance. An analysis of the conformation of the chains showed that at all three temperatures a fraction of the alkyl chains retained a planar all-trans conformation, and that gauche bonds occurred as part of a “kink” (gauche+−trans−gauche−) sequence and not as isolated gauche bonds. Trans−gauche isomerization rates for the alkyl chains in the anchored bilayer are slower than those in lipid bilayers at the same temperature and show a progressive increase as the torsion numbers approach the tail. A two-dimensional periodic Voronoi tessellation analysis was performed to obtain the single-molecular area of an alkyl chain in the bilayer. The single-molecular area relaxation times are an order of magnitude longer than the trans−gauche isomerization times. The results indicate that the trans−gauche isomerization is associated with the creation and annihilation of a kink defect sequence. The results of the present MD simulation explain the apparent conflicting estimates of the gauche disorder in this system as obtained from infrared and 13C nuclear magnetic resonance measurements.

Levitation effect in zeolites: Quasielastic neutron scattering and molecular dynamics study of pentane isomers in zeolite NaY

We report the quasielastic neutron scattering (QENS) and molecular dynamics (MD) investigations into diffusion of pentane isomers in zeolite NaY. The molecular cross section perpendicular to the long molecular axis varies for the three isomers while the mass and the isomer-zeolite interaction remains essentially unchanged. Both QENS and MD results show that the branched isomers neopentane and isopentane have higher self-diffusivities as compared with n-pentane at 300 K in NaY zeolite. This result provides direct experimental evidence for the existence of nonmonotonic, anomalous dependence of self-diffusivity on molecular diameter known as the levitation effect. The energetic barrier at the bottleneck derived from MD simulations exists for n-pentane which lies in the linear regime while no such barrier is seen for neopentane which is located clearly in the anomalous regime.Activation energy is in the order E-a(n-pentane)>E-a(isopentane)>E-a(neopentane) consistent with the predictions of the levitation effect. In the liquid phase, it is seen thatD(n pentane)>D(isopentane)>D(neopentane) and E-a(n-pentane)< E-a(isopentane)< E-a(neopentane). Intermediate scattering function for small wavenumbers obtained from MD follows a single exponential decay for neopentane and isopentane. For n-pentane, a single exponential fit provides a poor fit especially at short times. Cage residence time is largest for n-pentane and lowest for neopentane. For neopentane, the width of the self-part of the dynamic structure factor shows a near monotonic decrease with wavenumber. For n-pentane a minimum is seen near k=0.5 A degrees(-1) suggesting a slowing down of motion around the 12-ring window, the bottleneck for diffusion. Finally, the result that the branched isomer has a higher diffusivity as compared with the linear analog is at variation from what is normally seen.

Complexation of siRNA with Dendrimer: A Molecular Modeling Approach

This paper reports the structural behavior and thermodynamics of the complexation of siRNA with poly(amidoamine) (PAMAM) dendrimers of generation 3 (G3) and 4 (G4) through fully atomistic molecular dynamics (MD) simulations accompanied by free energy calculations and inherent structure determination. We have also done simulation with one siRNA and two dendrimers (2 x G3 or 2xG4) to get the microscopic picture of various binding modes. Our simulation results reveal the formation of stable siRNA-dendrimer complex over nanosecond time scale. With the increase in dendrimcr generation, the charge ratio increases and hence the binding energy between siRNA and dendrimer also increases in accordance with available experimental measurements. Calculated radial distribution functions of amines groups of various subgenerations in a given generation of dendrimer and phosphate in backbone of siRNA reveals that one dendrimer of generation 4 shows better binding with siRNA almost wrapping the dendrimer when compared to the binding with lower generation dendrimer like G3. In contrast, two dendrimers of generation 4 show binding without siRNA wrapping the den-rimer because of repulsion between two dendrimers. The counterion distribution around the complex and the water molecules in the hydration shell of siRNA give microscopic picture of the binding dynamics. We see a clear correlation between water. counterions motions and the complexation i.e. the water molecules and counterions which condensed around siRNA are moved away from the siRNA backbone when dendrimer start binding to the siRNA back hone. As siRNA wraps/bind to the dendrimer counterions originally condensed onto siRNA (Na-1) and dendrimer (Cl-) get released. We give a quantitative estimate of the entropy of counterions and show that there is gain in entropy due to counterions release during the complexation. Furthermore, the free energy of complexation of IG3 and IG4 at two different salt concentrations shows that increase in salt concentration leads to the weakening of the binding affinity of siRNA and dendrimer.

Analysis of Hydrogen Bonding and Stability of Protein Secondary Structures in Molecular Dynamics Simulation

Molecular Dynamics (MD) simulations provide an atomic level account of the molecular motions and have proven to be immensely useful in the investigation of the dynamical structure of proteins. Once an MD trajectory is obtained, specific interactions at the molecular level can be directly studied by setting up appropriate combinations of distance and angle monitors. However, if a study of the dynamical behavior of secondary structures in proteins becomes important, this approach can become unwieldy. We present herein a method to study the dynamical stability of secondary structures in proteins, based on a relatively simple analysis of backbone hydrogen bonds. The method was developed for studying the thermal unfolding of beta-lactamases, but can be extended to other systems and adapted to study relevant properties.

Molecular dynamics studies on nucleoside 2',3'-cyclic phosphates

2',3'-cyclic nucleotides are intermediates and substrates of Ribonuclease (RNase)-catalysed reactions. The characterization of the equilibrium conformation as well as the flexibility inherent in these molecules helps in understanding the enzymatic action of RNases. The present study explores parameters like phase angle, glycosydic torsion angle and hydrogen bond to find possible interrelationship between them through Molecular Dynamics (MD) simulations on 3'-GMP, 3'-UMP, A>p, G>p, U>p, C>p, GpA>p and UpA>p. Interesting results of the effect of cyclisation and other constraints such as hydrogen bond between certain groups on the equilibrium ribose conformation have emerged from this study.

The stability of transmembrane helices: A molecular dynamics study on the isolated helices of bacteriorhodopsin

Bacteriorhodopsin (bR) continues to be a proven testing ground for the study of integral membrane proteins (IMPs). It is important to study the stability of the individual helices of bR, as they are postulated to exist as independently stable transmembmne helices (TMHs) and also for their utility as templates for modeling other IMPs with the postulated seven-helix bundle topology. Toward this purpose, the seven helices of bR have been studied by molecular dynamics simulation in this study. The suitability of using the backbone-dependent rotamer library of side-chain conformations arrived at from the data base of globular protein structures in the case TMHs has been tested by another set of ? helix simulations with the side-chain orientations taken from this library. The influence of the residue's net charge oil the helix stability was examined by simulating the helices III, IV, and VI (from both of the above sets of helices) with zero net charge on the side chains. The results of these 20 simulations demonstrate in general the stability of the isolated helices of bR in conformity with the two-stage hypothesis of IMP folding. However, the helices I, II, V, and VII are more stable than the other three helices. The helical nature of certain regions of III, IV, and VI are influenced by factors such as the net charge and orientation of several residues. It is seen that the residues Arg, Lys, Asp, and Glu (charged residues), and Ser, Thr, Gly, and Pro, play a crucial role in the stability of the helices of bR. The backbone-dependent rotamer library for the side chains is found to be suitable for the study of TMHs in IMP. (C) 1996 John Wiley & Sons, Inc.

Monte Carlo and molecular dynamics simulation of argon clusters andn-alkanes in the confined regions of zeolites

Geometry and energy of argon clusters confined in zeolite NaCaA are compared with those of free clusters. Results indicate the possible existence of magic numbers among the confined clusters. Spectra obtained from instantaneous normal mode analysis of free and confined clusters give a larger percentage of imaginary frequencies for the latter indicating that the confined cluster atoms populate the saddle points of the potential energy surface significantly. The variation of the percentage of imaginary frequencies with temperature during melting is akin to the variation of other properties. It is shown that confined clusters might exhibit inverse surface melting, unlike medium-to-large-sized free clusters that exhibit surface melting. Configurational-bias Monte Carte (CBMC) simulations of n-alkanes in zeolites Y and A are reported. CBMC method gives reliable estimates of the properties relating to the conformation of molecules. Changes in the conformational properties of n-butane and other longer n-alkanes such as n-hexane and n-heptane when they are confined in different zeolites are presented. The changes in the conformational properties of n-butane and n-hexane with temperature and concentration is discussed. In general, in zeolite Y as well as A, there is significant enhancement of the gauche population as compared to the pure unconfined fluid.

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.