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.

Influence of grain boundary on melting

The temperature behaviour of an Al bicrystal with surfaces consisting of (110) and (111) crystals is simulated using molecular dynamics. The result shows that the (110) crystal losses its crystalline order at 820K, whereas the disorder does not propagate through the (111) crystal at this temperature. Instead, some disordered atoms are recrystallized into the (111) crystal and the initial grain boundary changes into a stable order-disorder interface. Thus, it was discovered that at a temperature near its melting point, the (111) crystal grew and obstructed the propagation of disorder. Such an obstruction is helpful for understanding melting.

The Uniaxial Tensile Deformation of Ni Nanowire: Atomic-Scale Computer Simulations

The mechanical deformations of nickel nanowire subjected to uniaxial tensile strain at 300 K are simulated by using molecular dynamics with the quantum corrected Sutten-Chen many-body force field. We have used common neighbor analysis method to investigate the structural evolution of Ni nanowire during the elongation process. For the strain rate of 0.1%/ps, the elastic limit is up to about 11% strain with the yield stress of 8.6 GPa. At the elastic stage, the deformation is carried mainly through the uniform elongation of the distances between the layers (perpendicular to the Z-axis) while the atomic structure remains basically unchanged. With further strain, the slips in the {111} planes start to take place in order to accommodate the applied strain to carry the deformation partially, and subsequently the neck forms. The atomic rearrangements in the neck region result in a zigzag change in the stress-strain curve; the atomic structures beyond the region, however, have no significant changes. With the strain close to the point of the breaking, we observe the formation of a one-atom thick necklace in Ni nanowire. The strain rates have no significant effect on the deformation mechanism, but have some influence on the yield stress, the elastic limit, and the fracture strain of the nanowire.

分子动力学模拟的主要技术

综述了分子动力学模拟技术的发展，介绍了分子动力学的基本原理，有关的有限差分技术，势函数的发展，初始条件和边界条件的选取，平衡态系综及其调控，感兴趣量的提取，分子动力学的特别用途以及与其他计算方法的结合，最后还指出了分子动力学模拟方法本身进一步研究的方向。

Deformation twinning and its effect on crack extension

Deformation twinning near a crack tip is observed in b.c.c. metal Mo based on molecular dynamics simulation at temperature T = 50 K and loading rate (K) over dot(II) = 0.0706 MPa m(1/2)/ps. The defor mation twinning is closely controlled by both the crystal geometry orientation and the stress distribution. The width of the deformation twin band is affected by the distance between the upper and lower crack surfaces. The twin plane and twin direction are (<1(1)over bar>2) and [(1) over bar 11], respectively. The initial crack extension occurs in the deformation twin region near the crack tip. The simulation shows that the extension direction of the crack is changed as the crack propagates over the twinning boundary.

Effects Analysis on the Mechanical Properties of Nanocrystalline Materials

Nanocrystalline materials are characterized by a typical grain size from 1 to 100nm. In order to study the nanocrystalline properties of nanocrystalline materials, we chose nanocrystalline coppers as the research object. The uniaxial tensile deformation of computer produced nanocrystalline coppers is simulated by using molecular dynamics with Finnis-Sinclair potential. The mean grain size of simulated nanocrystalline coppers is varied within the 5.38 to 1.79 nm range. The strength, Young's modulus and stress-strain are strongly depended on the grain size and nanocrystalline structure. The simulated nanocrystalline coppers show a reverse Hall-Petch effect.