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.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

Excited state structure and dynamics of p-benzoquinone and bromanil from time-resolved resonance Raman spectra and simulation

p-Benzoquinone and its halogen substituted derivatives are known to have differing reactivities in the triplet excited state. While bromanil catalyzes the reduction of octaethylporphyrin most efficiently among the halogenated p-benzoquinones, the reaction does not take place in presence of the unsubstituted p-benzoquinone (T. Nakano and Y. Mori, Bull. Chem. Soc. Jpn., 67, 2627 (1994)). Understanding of such differences requires a detailed knowledge of the triplet state structures, normal mode compositions and excited state dynamics. In this paper, we apply a recently presented scheme (M. Puranik, S. Umapathy, J. G. Snijders, and J. Chandrasekhar, J. Chem, Phys., 115, 6106 (2001)) that combines parameters from experiment and computation in a wave packet dynamics simulation to the triplet states of p-benzoquinone and bromanil. The absorption and resonance Raman spectra of both the molecules have been simulated. The normal mode compositions and mode specific excited state displacements have been presented and compared. Time-dependent evolution of the absorption and Raman overlaps for all the observed modes has been discussed in detail. In p-benzoquinone, the initial dynamics is along the C=C stretching and C-H bending modes whereas in bromanil nearly equal displacements are observed along all the stretching coordinates.

Computer simulation of liquid crystals

In this article we review the current status in the modelling of both thermotropic and lyotropic Liquid crystal. We discuss various coarse-graining schemes as well as simulation techniques such as Monte Carlo (MC) and Molecular dynamics (MD) simulations.In the area of MC simulations we discuss in detail the algorithm for simulating hard objects such as spherocylinders of various aspect ratios where excluded volume interaction enters in the simulation through overlap test. We use this technique to study the phase diagram, of a special class of thermotropic liquid crystals namely banana liquid crystals. Next we discuss a coarse-grain model of surfactant molecules and study the self-assembly of the surfactant oligomers using MD simulations. Finally we discuss an atomistically informed coarse-grained description of the lipid molecules used to study the gel to liquid crystalline phase transition in the lipid bilayer system.

A spectral multiscale method for wave propagation analysis: Atomistic-continuum coupled simulation

In this paper, we present a new multiscale method which is capable of coupling atomistic and continuum domains for high frequency wave propagation analysis. The problem of non-physical wave reflection, which occurs due to the change in system description across the interface between two scales, can be satisfactorily overcome by the proposed method. We propose an efficient spectral domain decomposition of the total fine scale displacement along with a potent macroscale equation in the Laplace domain to eliminate the spurious interfacial reflection. We use Laplace transform based spectral finite element method to model the macroscale, which provides the optimum approximations for required dynamic responses of the outer atoms of the simulated microscale region very accurately. This new method shows excellent agreement between the proposed multiscale model and the full molecular dynamics (MD) results. Numerical experiments of wave propagation in a 1D harmonic lattice, a 1D lattice with Lennard-Jones potential, a 2D square Bravais lattice, and a 2D triangular lattice with microcrack demonstrate the accuracy and the robustness of the method. In addition, under certain conditions, this method can simulate complex dynamics of crystalline solids involving different spatial and/or temporal scales with sufficient accuracy and efficiency. (C) 2014 Elsevier B.V. All rights reserved.

Investigation of selective sensing of a diamine for aldehyde by experimental and simulation studies

An organic molecule-o-phenylene diamine (OPD)-is selected as an aldehyde sensing material. It is studied for selectivity to aldehyde vapours both by experiment and simulation. A chemiresistor based sensor for detection of aldehyde vapours is fabricated. An o-phenylene diamine-carbon black composite is used as the sensing element. The amine groups in the OPD would interact with the carbonyl groups of the aldehydes. The selectivity and cross-sensitivity of the OPD-CB sensor to VOCs aldehyde, ketone and alcohol-are studied. The sensor shows good response to aldehydes compared to other VOCs. The higher response for aldehydes is attributed to the interaction of the carbonyl oxygen of aldehydes with-NH2 groups of OPD. The surface morphology of the sensing element is studied by scanning electron microscopy. The OPD-CB sensor is responsive to 10 ppm of formaldehyde. The interaction of the VOCs with the OPD-CB nanocomposite is investigated by molecular dynamics studies. The interaction energies of the analyte with the OPD-CB nanocomposite were calculated. It is observed that the interaction energies for aldehydes are higher than those for other analytes. Thus the OPD-CB sensor shows selectivity to aldehydes. The simulated radial distribution function is calculated for the O-H pair of analyte and OPD which further supports the finding that the amine groups are involved in the interaction. These results suggest that it is important and easy to identify appropriate sensing materials based on the understanding of analyte interaction properties.

Mode coupling theory analysis of electrolyte solutions: Time dependent diffusion, intermediate scattering function, and ion solvation dynamics

A self-consistent mode coupling theory (MCT) with microscopic inputs of equilibrium pair correlation functions is developed to analyze electrolyte dynamics. We apply the theory to calculate concentration dependence of (i) time dependent ion diffusion, (ii) intermediate scattering function of the constituent ions, and (iii) ion solvation dynamics in electrolyte solution. Brownian dynamics with implicit water molecules and molecular dynamics method with explicit water are used to check the theoretical predictions. The time dependence of ionic self-diffusion coefficient and the corresponding intermediate scattering function evaluated from our MCT approach show quantitative agreement with early experimental and present Brownian dynamic simulation results. With increasing concentration, the dispersion of electrolyte friction is found to occur at increasingly higher frequency, due to the faster relaxation of the ion atmosphere. The wave number dependence of intermediate scattering function, F(k, t), exhibits markedly different relaxation dynamics at different length scales. At small wave numbers, we find the emergence of a step-like relaxation, indicating the presence of both fast and slow time scales in the system. Such behavior allows an intriguing analogy with temperature dependent relaxation dynamics of supercooled liquids. We find that solvation dynamics of a tagged ion exhibits a power law decay at long times-the decay can also be fitted to a stretched exponential form. The emergence of the power law in solvation dynamics has been tested by carrying out long Brownian dynamics simulations with varying ionic concentrations. The solvation time correlation and ion-ion intermediate scattering function indeed exhibit highly interesting, non-trivial dynamical behavior at intermediate to longer times that require further experimental and theoretical studies. (c) 2015 AIP Publishing LLC.

Atomistic Simulation of Stacked Nucleosome Core Particles: Tail Bridging, the H4 Tail, and Effect of Hydrophobic Forces

We report the first atomistic simulation of two stacked nucleosome core particles (NCPs), with an aim to understand, in molecular detail, how they interact, the effect of salt concentration, and how different histone tails contribute to their interaction, with a special emphasis on the H4 tail, known to have the largest stabilizing effect on the NCP-NCP interaction. We do not observe specific K16-mediated interaction between the H4 tail and the H2A-H2B acidic patch, in contrast with the findings from crystallographic studies, but find that the stacking was stable even in the absence of this interaction. We perform simulations with the H4 tail (partially/completely) removed and find that the region between LYS-16 and LYS-20 of the H4 tail holds special importance in mediating the inter-NCP interaction. Performing similar tail-clipped simulations with the H3 tail removed, we compare the roles of the H3 and H4 tails in maintaining the stacking. We discuss the relevance of our simulation results to the bilayer and other liquid-crystalline phases exhibited by NCPs in vitro and, through an analysis of the histone-histone interface, identify the interactions that could possibly stabilize the inter-NCP interaction in these columnar mesophases. Through the mechanical disruption of the stacked nucleosome system using steered molecular dynamics, we quantify the strength of inter-NCP stacking in the presence and absence of salt. We disrupt the stacking at some specific sites of internucleosomal tail-DNA contact and perform a comparative quantification of the binding strengths of various tails in stabilizing the stacking. We also examine how hydrophobic interactions may contribute to the overall stability of the stacking and find a marked difference in the role of hydrophobic forces as compared with electrostatic forces in determining the stability of the stacked nucleosome system.

Numerical simulation of rock failure and earthquake process on mesoscopic scale

On the basis of the lattice model of MORA and PLACE, Discrete Element Method, and Molecular Dynamics approach, another kind of numerical model is developed. The model consists of a 2-D set of particles linked by three kinds of interactions and arranged into triangular lattice. After the fracture criterion and rules of changes between linking states are given, the particle positions, velocities and accelerations at every time step are calculated using a finite-difference scheme, and the configuration of particles can be gained step by step. Using this model, realistic fracture simulations of brittle solid (especially under pressure) and simulation of earthquake dynamics are made.

MD simulation of the effect of contact area and tip radius on nanoindentation

Molecular dynamics simulations of nanoindentation are performed on monocrystal copper. A new "contact atoms" method is presented for calculating the contact area. Compared with conventional methods, this method can provide the contact area more accurately not only for sink-in but also for pile-up situation. The effect of tip radius on indentation is investigated too. The results indicate that the measured hardness of the material will become higher as the tip radius increases.

Computer simulation on the collision-sticking dynamics of two colloidal particles in an optical trap

Collisions of a particle pair induced by optical tweezers have been employed to study colloidal stability. In order to deepen insights regarding the collision-sticking dynamics of a particle pair in the optical trap that were observed in experimental approaches at the particle level, the authors carry out a Brownian dynamics simulation. In the simulation, various contributing factors, including the Derjaguin-Landau-Verwey-Overbeek interaction of particles, hydrodynamic interactions, optical trapping forces on the two particles, and the Brownian motion, were all taken into account. The simulation reproduces the tendencies of the accumulated sticking probability during the trapping duration for the trapped particle pair described in our previous study and provides an explanation for why the two entangled particles in the trap experience two different statuses. (c) 2007 American Institute of Physics.

Hybrid Qm/Mm Simulation Of The Hydration Phenomena Of Dipalmitoylphosphatidylcholine Headgroup

The polar headgroup of dipalmitoylphosphatidylcholine (DPPC) molecule both in gas phase and aqueous Solution is investigated by the hybrid quantum mechanical/molecular mechanical (QM/MM) method, in which the polar head of DPPC molecule and the bound water molecules are treated with density functional theory (DFT), while the apolar hydrocarbon chain of DPPC molecule is treated with MM method. It is demonstrated that the hybrid QM/MM method is both accurate and efficient to describe the conformations of DPPC headgroup. Folded structures of headgroup are found in gas phase calculations. In this work, both monohydration and polyhydration phenomena are investigated. In monohydration, different water association sites are studied. Both the hydration energy and the quantum properties of DPPC and water molecules are calculated at the DFT level of theory after geometry optimization. The binding force of monohydration is estimated by using the scan method. In polyhydration, more extended conformations are found and hydration energies in different polyhydration styles are estimated. (C) 2008 Elsevier Inc. All rights reserved.

Dynamics of RNA with A6 bulge loop

In order to understand the behavior of RNAs with large bulges In solution, molecular dynamics was performed on the RNA molecule in water with A6 bulge. The result of simulation showed that nonstacked conformation Is the main conformation in large bulges, and the backbone of large bulge is of great conformational flexibility, but bulges-induced bends are relatively rigid. The fluctuation in bulge has little influence on the bend angle of RNAs.

Computer simulation of miscibility and self-assembly structure for polymer-containing systems with special interactions

The miscibility and structure of A-B copolymer/C homopolymer blends with special interactions were studied by a Monte Carlo simulation in two dimensions. The interaction between segment A and segment C was repulsive, whereas it was attractive between segment B and segment C. In order to study the effect of copolymer chain structure on the morphology and structure of A-B copolymer/C homopolymer blends, the alternating, random and block A-B copolymers were introduced into the blends, respectively. The simulation results indicated that the miscibility of A-B block copolymer/C homopolymer blends depended on the chain structure of the A-B copolymer. Compared with alternating or random copolymer, the block copolymer, especially the diblock copolymer, could lead to a poor miscibility of A-B copolymer/C homopolymer blends. Moreover, for diblock A-B copolymer/C homopolymer blends, obvious self-organized core-shell structure was observed in the segment B composition region from 20% to 60%. However, if diblock copolymer composition in the blends is less than 40%, obvious self-organized core-shell structure could be formed in the B-segment component region from 10 to 90%. Furthermore, computer statistical analysis for the simulation results showed that the core sizes tended to increase continuously and their distribution became wider with decreasing B-segment component.

Monte Carlo simulation of phase behavior of polymer blends with special interactions

Full Paper: The phase, behavior of A-B-random copolymer/C-homopolymer, blends with special interaction was studied by a. Monte, Carlo simulation in two dimensions. The interaction between I segment A and segment C was repulsive, whereas it was attractive between segment B and segment C. The simulation results showed that the blend became two large co-continuous phase domains at lower segment-B component compositions, indicating that the blend showed spinodal decomposition. With an increase of the segment-B component, the miscibility between the copolymer,and the polymer was gradually improved up to being miscible. In addition, it was found that segment B tended to move to the surface of the copolymer phase in the case of a lower component of segment B. On the other hand, if was observed that the average, end-to-end distances ((h) over bar) for both copolymer and polymer changed slowly with increasing segment-B component of the copolymer up to 40%, thereafter they increased considerably with increasing segment B component. Moreover, it was found that the (h) over bar of the copolymer was obviously shorter than that of the homopolymer for the segment-B composition, region from 0% to 80%. Finally, a, phase diagram showing I phase and - II phase regions under the condition of constant-temperature is presented.