Crack propagation in graphene

The crack initiation and growth mechanisms in an 2D graphene lattice structure are studied based on molecular dynamics simulations. Crack growth in an initial edge crack model in the arm-chair and the zig-zag lattice configurations of graphene are considered. Influence of the time steps on the post yielding behaviour of graphene is studied. Based on the results, a time step of 0.1 fs is recommended for consistent and accurate simulation of crack propagation. Effect of temperature on the crack propagation in graphene is also studied, considering adiabatic and isothermal conditions. Total energy and stress fields are analyzed. A systematic study of the bond stretching and bond reorientation phenomena is performed, which shows that the crack propagates after significant bond elongation and rotation in graphene. Variation of the crack speed with the change in crack length is estimated. (C) 2015 AIP Publishing LLC.

DNA Elasticity from Short DNA to Nucleosomal DNA

Active biological processes like transcription, replication, recombination, DNA repair, and DNA packaging encounter bent DNA. Machineries associated with these processes interact with the DNA at short length (<100 base pair) scale. Thus, the study of elasticity of DNA at such length scale is very important. We use fully atomistic molecular dynamics (MD) simulations along with various theoretical methods to determine elastic properties of dsDNA of different lengths and base sequences. We also study DNA elasticity in nucleosome core particle (NCP) both in the presence and the absence of salt. We determine stretch modulus and persistence length of short dsDNA and nucleosomal DNA from contour length distribution and bend angle distribution, respectively. For short dsDNA, we find that stretch modulus increases with ionic strength while persistence length decreases. Calculated values of stretch modulus and persistence length for DNA are in quantitative agreement with available experimental data. The trend is opposite for NCP DNA. We find that the presence of histone core makes the DNA stiffer and thus making the persistence length 3-4 times higher than the bare DNA. Similarly, we also find an increase in the stretch modulus for the NCP DNA. Our study for the first time reports the elastic properties of DNA when it is wrapped around the histone core in NCP. We further show that the WLC model is inadequate to describe DNA elasticity at short length scale. Our results provide a deeper understanding of DNA mechanics and the methods are applicable to most protein-DNA complexes.

Suspensions of polymer-grafted nanoparticles with added polymers-Structure and effective pair-interactions

We present the results of combined experimental and theoretical (molecular dynamics simulations and integral equation theory) studies of the structure and effective interactions of suspensions of polymer grafted nanoparticles (PGNPs) in the presence of linear polymers. Due to the absence of systematic experimental and theoretical studies of PGNPs, it is widely believed that the structure and effective interactions in such binary mixtures would be very similar to those of an analogous soft colloidal material-star polymers. In our study, polystyrene-grafted gold nanoparticles with functionality f = 70 were mixed with linear polystyrene (PS) of two different molecular weights for obtaining two PGNP: PS size ratios, xi = 0.14 and 2.76 (where, xi = M-g/M-m, M-g and M-m being the molecular weights of grafting and matrix polymers, respectively). The experimental structure factor of PGNPs could be modeled with an effective potential (Model-X), which has been found to be widely applicable for star polymers. Similarly, the structure factor of the blends with xi = 0.14 could be modeled reasonably well, while the structure of blends with xi = 2.76 could not be captured, especially for high density of added polymers. A model (Model-Y) for effective interactions between PGNPs in a melt of matrix polymers also failed to provide good agreement with the experimental data for samples with xi = 2.76 and high density of added polymers. We tentatively attribute this anomaly in modeling the structure factor of blends with xi = 2.76 to the questionable assumption of Model-X in describing the added polymers as star polymers with functionality 2, which gets manifested in both polymer-polymer and polymer-PGNP interactions especially at higher fractions of added polymers. The failure of Model-Y may be due to the neglect of possible many-body interactions among PGNPs mediated by matrix polymers when the fraction of added polymers is high. These observations point to the need for a new framework to understand not only the structural behavior of PGNPs but also possibly their dynamics and thermo-mechanical properties as well. (C) 2015 AIP Publishing LLC.

Mono layer to Interdigitated Partial Bilayer Smectic C Transition in Thiophene-Based Spacer Mesogens: X-ray Diffraction and C-13 Nuclear Magnetic Resonance Studies

Mesophase organization of molecules built with thiophene at the center and linked via flexible spacers to rigid side arm core units and terminal alkoxy chains has been investigated. Thirty homologues realized by varying the span of the spacers as well as the length of the terminal chains have been studied. In addition to the enantiotropic nematic phase observed for all the mesogens, the increase of the spacer as well as the terminal chain lengths resulted in the smectic C phase. The molecular organization in the smectic phase as investigated by temperature dependent X-ray diffraction measurements revealed an interesting behavior that depended on the length of the spacer vis-a-vis the length of the terminal chain. Thus, a tilted interdigitated partial bilayer organization was observed for molecules with a shorter spacer length, while a tilted monolayer arrangement was observed for those with a longer spacer length. High-resolution solid state C-13 NMR studies carried out for representative mesogens indicated a U-shape for all the molecules, indicating that intermolecular interactions and molecular dynamics rather than molecular shape are responsible for the observed behavior. Models for the mesophase organization have been considered and the results understood in terms of segregation of incompatible parts of the mesogens combined with steric frustration leading to the observed lamellar order.

Stabilization of the high-temperature and high-pressure cubic phase of ZnO by temperature-controlled milling

The reversible transition of wurtzite to rock salt phase under pressure is well reported in literature. The cubic phase is unstable under ambient conditions both in the bulk and in nanoparticles. This paper reports defect-induced stabilization of cubic ZnO phase in sub 20 nm ZnO particles and explores their optical properties. The size reduction was achieved by ball milling in a specially designed mill which allows a control of the milling temperature. The process of synthesis involved both variation of milling temperature (including low temperature similar to 150 K) and impact pressure. We show that these have profound influence in the introduction of defects and stabilization of the cubic phase. A molecular dynamics simulation is presented to explain the observed results. The measured optical properties have further supported the observations of defect-induced stabilization of cubic ZnO and reduction in particle size.

The putative role of some conserved water molecules in the structure and function of human transthyretin

Human transthyretin (hTTR) is a multifunctional protein that is involved in several neurodegenerative diseases. Besides the transportation of thyroxin and vitamin A, it is also involved in the proteolysis of apolipoprotein A1 and A beta peptide. Extensive analyses of 32 high-resolution X-ray and neutron diffraction structures of hTTR followed by molecular-dynamics simulation studies using a set of 15 selected structures affirmed the presence of 44 conserved water molecules in its dimeric structure. They are found to play several important roles in the structure and function of the protein. Eight water molecules stabilize the dimeric structure through an extensive hydrogen-bonding network. The absence of some of these water molecules in highly acidic conditions (pH <= 4.0) severely affects the interfacial hydrogen-bond network, which may destabilize the native tetrameric structure, leading to its dissociation. Three pairs of conserved water molecules contribute to maintaining the geometry of the ligand-binding cavities. Some other water molecules control the orientation and dynamics of different structural elements of hTTR. This systematic study of the location, absence, networking and interactions of the conserved water molecules may shed some light on various structural and functional aspects of the protein. The present study may also provide some rational clues about the conserved water-mediated architecture and stability of hTTR.

The SPL7013 dendrimer destabilizes the HIV-1 gp120-CD4 complex

The poly (l-lysine)-based SPL7013 dendrimer with naphthalene disulphonate surface groups blocks the entry of HIV-1 into target cells and is in clinical trials for development as a topical microbicide. Its mechanism of action against R5 HIV-1, the HIV-1 variant implicated in transmission across individuals, remains poorly understood. Using docking and fully atomistic MD simulations, we find that SPL7013 binds tightly to R5 gp120 in the gp120-CD4 complex but weakly to gp120 alone. Further, the binding, although to multiple regions of gp120, does not occlude the CD4 binding site on gp120, suggesting that SPL7013 does not prevent the binding of R5 gp120 to CD4. Using MD simulations to compute binding energies of several docked structures, we find that SPL7013 binding to gp120 significantly weakens the gp120-CD4 complex. Finally, we use steered molecular dynamics (SMD) to study the kinetics of the dissociation of the gp120-CD4 complex in the absence of the dendrimer and with the dendrimer bound in each of the several stable configurations to gp120. We find that SPL7013 significantly lowers the force required to rupture the gp120-CD4 complex and accelerates its dissociation. Taken together, our findings suggest that SPL7013 compromises the stability of the R5 gp120-CD4 complex, potentially preventing the accrual of the requisite number of gp120-CD4 complexes across the virus-cell interface, thereby blocking virus entry.

Ab Initio MD Simulations of the Bronsted Acidity of Glutathione in Aqueous Solutions: Predicting pK(a) Shifts of the Cysteine Residue

The tripeptide glutathione (GSH) is one of the most abundant peptides and the major repository for nonprotein sulfur in both animal and plant cells. It plays a critical role in intracellular oxidative stress management by the reversible formation of glutathione disulfide with the thiol-disulfide pair acting as a redox buffer. The state of charge of the ionizable groups of GSH can influence the redox couple, and hence the pK(a) value of the cysteine residue of GSH is critical to its functioning. Here we report ab initio Car-Parrinello molecular dynamics simulations of glutathione solvated by 200 water molecules, all of which are considered in the simulation. We show that the free-energy landscape for the protonation-deprotonation reaction of the cysteine residue of GSH computed using metadynamics sampling provides shift in the dissociation constant values as compared with the isolated accurate estimates of the pK(a) and correctly predicts the cysteine amino acid.

Vesicular Nanostructure Formation by Self-Assembly of Anisotropic Penta-phenol-Substituted Fullerene in Water

A study on self-assembly of anisotropically substituted penta-aryl fullerenes in water has been reported. The penta-phenol-substituted amphiphilic fullerene derivative C60Ph5(OH)(5)],exhibited self-assembled vesicular nanostructures in water under the experimental conditions. The size of the vesicles Was observed to depend upon the kinetics of self-assembly and could be varied from similar to 300 to similar to 70 nm. Our mechanistic study indicated that the self-assembly of C60Ph5(OH)(5) is driven by extensive intermolecular as well as water-mediated hydrogen :bonding along with fullerene-fullerene hydrophobic interaction in water. The cumulative effect of these interactions is responsible for the stability of vesicular structures even on the removal of solvent. The substitution of phenol with anisole resulted in different packing and interaction of the fullerene derivative, as Indicated in the molecular dynamics studies, thus resulting in different self-assembled nanostructures. The hollow vesicles were further encapsulated with a hydrophobic conjugated polymer and water-soluble dye as guest molecules. Such confinement of pi-conjugated polymers in fullerene has significance in bulk heterojunction devices for efficient exciton diffusion.

Nanoparticle induced miscibility in LCST polymer blends: critically assessing the enthalpic and entropic effects

The use of copolymer and polymer blends widened the possibility of creating materials with multilayered architectures. Hierarchical polymer systems with a wide array of micro and nanostructures are generated by thermally induced phase separation (TIPS) in partially miscible polymer blends. Various parameters like the interaction between the polymers, concentration, solvent/non-solvent ratio, and quenching temperature have to be optimized to obtain these micro/nanophase structures. Alternatively, the addition of nanoparticles is another strategy to design materials with desired hetero-phase structures. The dynamics of the polymer nanocomposite depends on the statistical ordering of polymers around the nanoparticle, which is dependent on the shape of the nanoparticle. The entropic loss due to deformation of polymer chains, like the repulsive interactions due to coiling and the attractive interactions in the case of swelling has been highlighted in this perspective article. The dissipative particle dynamics has been discussed and is correlated with the molecular dynamics simulation in the case of polymer blends. The Cahn Hillard Cook model on variedly shaped immobile fillers has shown difference in the propagation of the composition wave. The nanoparticle shape has a contributing effect on the polymer particle interaction, which can change the miscibility window in the case of these phase separating polymer blends. Quantitative information on the effect of spherical particles on the demixing temperature is well established and further modified to explain the percolation of rod shaped particles in the polymer blends. These models correlate well with the experimental observations in context to the dynamics induced by the nanoparticle in the demixing behavior of the polymer blend. The miscibility of the LCST polymer blend depends on the enthalpic factors like the specific interaction between the components, and the solubility product and the entropic losses occurring due to the formation of any favorable interactions. Hence, it is essential to assess the entropic and enthalpic interactions induced by the nanoparticles independently. The addition of nanoparticles creates heterogeneity in the polymer phase it is localized. This can be observed as an alteration in the relaxation behavior of the polymer. This changes the demixing behavior and the interaction parameter between the polymers. The compositional changes induced due to the incorporation of nanoparticles are also attributed as a reason for the altered demixing temperature. The particle shape anisotropy causes a direction dependent depletion, which changes the phase behavior of the blend. The polymer-grafted nanoparticles with varying grafting density show tremendous variation in the miscibility of the blend. The stretching of the polymer chains grafted on the nanoparticles causes an entropy penalty in the polymer blend. A comparative study on the different shaped particles is not available up to date for understanding these aspects. Hence, we have juxtaposed the various computational studies on nanoparticle dynamics, the shape effect of NPs on homopolymers and also the cases of various polymer blends without nanoparticles to sketch a complete picture on the effect of various particles on the miscibility of LCST blends.

Potency and pharmacokinetics of broad spectrum and isoform-specific p110 gamma and delta inhibitors in cancers

Emerging data on cancer suggesting that target-based therapy is promising strategy in cancer treatment. PI3K-AKT pathway is extensively studied in many cancers; several inhibitors target this pathway in different levels. Recent finding on this pathway uncovered the therapeutic applications of PI3K-specific inhibitors; PI3K, AKT, and mTORC broad spectrum inhibitors. Noticeably, class I PI3K isoforms, p110 and p110 catalytic subunits have rational therapeutic application than other isoforms. Therefore, three classes of inhibitors: isoform-specific, dual-specific and broad spectrum were selected for molecular docking and dynamics. First, p110 structure was modelled; active site was analyzed. Then, molecular docking of each class of inhibitors were studied; the docked complexes were further used in 1.2ns molecular dynamics simulation to report the potency of each class of inhibitor. Remarkably, both the studies retained the similar kind of protein ligand interactions. GDC-0941, XL-147 (broad spectrum); TG100-115 (dual-specific); and AS-252424, PIK-294 (isoform-specific) were found to be potential inhibitors of p110 and p110, respectively. In addition to that pharmacokinetic properties are within recommended ranges. Finally, molecular phylogeny revealed that p110 and p110 are evolutionarily divergent; they probably need separate strategies for drug development.

Short-Time Beta Relaxation in Glass-Forming Liquids Is Cooperative in Nature

Temporal relaxation of density fluctuations in supercooled liquids near the glass transition occurs in multiple steps. Using molecular dynamics simulations for three model glass-forming liquids, we show that the short-time beta relaxation is cooperative in nature. Using finite-size scaling analysis, we extract a growing length scale associated with beta relaxation from the observed dependence of the beta relaxation time on the system size. We find, in qualitative agreement with the prediction of the inhomogeneous mode coupling theory, that the temperature dependence of this length scale is the same as that of the length scale that describes the spatial heterogeneity of local dynamics in the long-time alpha-relaxation regime.

Folding of Protein L with Implications for Collapse in the Denatured State Ensemble

A fundamental question in protein folding is whether the coil to globule collapse transition occurs during the initial stages of folding (burst phase) or simultaneously with the protein folding transition. Single molecule fluorescence resonance energy transfer (FRET) and small-angle X-ray scattering (SAXS) experiments disagree on whether Protein L collapse transition occurs during the burst phase of folding. We study Protein L folding using a coarse-grained model and molecular dynamics simulations. The collapse transition in Protein L is found to be concomitant with the folding transition. In the burst phase of folding, we find that FRET experiments overestimate radius of gyration, R-g, of the protein due to the application of Gaussian polymer chain end-to-end distribution to extract R-g from the FRET efficiency. FRET experiments estimate approximate to 6 angstrom decrease in R-g when the actual decrease is approximate to 3 angstrom on guanidinium chloride denaturant dilution from 7.5 to 1 M, thereby suggesting pronounced compaction in the protein dimensions in the burst phase. The approximate to 3 angstrom decrease is close to the statistical uncertainties of the R-g data measured from SAXS experiments, which suggest no compaction, leading to a disagreement with the FRET experiments. The transition-state ensemble (TSE) structures in Protein L folding are globular and extensive in agreement with the Psi-analysis experiments. The results support the hypothesis that the TSE of single domain proteins depends on protein topology and is not stabilized by local interactions alone.

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.

Elasticity, stability, and ideal strength of beta-SiC in plane-wave-based ab initio calculations

On the basis of the pseudopotential plane-wave method and the local-density-functional theory, this paper studies energetics, stress-strain relation, stability, and ideal strength of beta-SiC under various loading modes, where uniform uniaxial extension and tension and biaxial proportional extension are considered along directions [001] and [111]. The lattice constant, elastic constants, and moduli of equilibrium state are calculated and the results agree well with the experimental data. As the four SI-C bonds along directions [111], [(1) over bar 11], [11(1) over bar] and [111] are not the same under the loading along [111], internal relaxation and the corresponding internal displacements must be considered. We find that, at the beginning of loading, the effect of internal displacement through the shuffle and glide plane diminishes the difference among the four Si-C bonds lengths, but will increase the difference at the subsequent loading, which will result in a crack nucleated on the {111} shuffle plane and a subsequently cleavage fracture. Thus the corresponding theoretical strength is 50.8 GPa, which agrees well with the recent experiment value, 53.4 GPa. However, with the loading along [001], internal relaxation is not important for tetragonal symmetry. Elastic constants during the uniaxial tension along [001] are calculated. Based on the stability analysis with stiffness coefficients, we find that the spinodal and Born instabilities are triggered almost at the same strain, which agrees with the previous molecular-dynamics simulation. During biaxial proportional extension, stress and strength vary proportionally with the biaxial loading ratio at the same longitudinal strain.