MD investigations for mechanical properties of copper nanowires with and without surface defects

Based on the molecular dynamics (MD) method, the single-crystalline copper nanowire with different surface defects is investigated through tension simulation. For comparison, the MD tension simulations of perfect nanowire are firstly carried out under different temperatures, strain rates, and sizes. It has concluded that the surface-volume ratio significantly affects the mechanical properties of nanowire. The surface defects on nanowires are then systematically studied in considering different defect orientation and distribution. It is found that the Young’s modulus is insensitive of surface defects. However, the yield strength and yield point show a significant decrease due to the different defects. Different defects are observed to serve as a dislocation source.

Failure to launch: The self-regulating Md -MYB10R6 gene from apple is active in flowers but not leaves of Petunia

Coloured foliage due to anthocyanin pigments (bronze/red/black) is an attractive trait that is often lacking in many bedding, ornamental and horticultural plants. Apples (Malus × domestica) containing an allelic variant of the anthocyanin regulator, Md-MYB10R6, are highly pigmented throughout the plant, due to autoregulation by MYB10 upon its own promoter. We investigated whether Md-MYB10R6 from apple is capable of functioning within the heterologous host Petunia hybrida to generate plants with novel pigmentation patterns. The Md-MYB10R6 transgene (MYB10–R6pro:MYB10:MYB10term) activated anthocyanin synthesis when transiently expressed in Antirrhinumroseadorsea petals and petunia leaf discs. Stable transgenic petunias containing Md-MYB10R6 lacked foliar pigmentation but had coloured flowers, complementing the an2 phenotype of ‘Mitchell’ petunia. The absence of foliar pigmentation was due to the failure of the Md-MYB10R6 gene to self-activate in vegetative tissues, suggesting that additional protein partners are required for Md-MYB10 to activate target genes in this heterologous system. In petunia flowers, where endogenous components including MYB-bHLH-WDR (MBW) proteins were present, expression of the Md-MYB10R6 promoter was initiated, allowing auto-regulation to occur and activating anthocyanin production. Md-MYB10 is capable of operating within the petunia MBW gene regulation network that controls the expression of the anthocyanin biosynthesis genes, AN1 (bHLH) and MYBx (R3-MYB repressor) in petals.

Allostery in tRNA Synthetases Elucidated from MD Simulations and Protein Structure Networks

tRNA synthetases (aaRS) are enzymes crucial in the translation of genetic code. The enzyme accylates the acceptor stem of tRNA by the congnate amino acid bound at the active site, when the anti-codon is recognized by the anti-codon site of aaRS. In a typical aaRS, the distance between the anti-codon region and the amino accylation site is approximately 70 Å. We have investigated this allosteric phenomenon at molecular level by MD simulations followed by the analysis of protein structure networks (PSN) of non-covalent interactions. Specifically, we have generated conformational ensembles by performing MD simulations on different liganded states of methionyl tRNA synthetase (MetRS) from Escherichia coli and tryptophenyl tRNA synthetase (TrpRS) from Human. The correlated residues during the MD simulations are identified by cross correlation maps. We have identified the amino acids connecting the correlated residues by the shortest path between the two selected members of the PSN. The frequencies of paths have been evaluated from the MD snapshots[1]. The conformational populations in different liganded states of the protein have been beautifully captured in terms of network parameters such as hubs, cliques and communities[2]. These parameters have been associated with the rigidity and plasticity of the protein conformations and can be associated with free energy landscape. A comparison of allosteric communication in MetRS and TrpRS [3] elucidated in this study highlights diverse means adopted by different enzymes to perform a similar function. The computational method described for these two enzymes can be applied to the investigation of allostery in other systems.

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.

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.

Engagement of CF3 Group in N-H center dot center dot center dot F-C Hydrogen Bond in the Solution State: NMR Spectroscopy and MD Simulation Studies

Unambiguous evidence for the engagement of CF3 group in N-H center dot center dot center dot F-C hydrogen bond in a low polarity solvent, the first observation of its kind, is reported. The presence of such weak molecular interactions in the solution state is convincingly established by one and two-dimensional H-1, F-19, and natural abundant N-15 NMR spectroscopic studies. The strong and direct evidence is derived by the observation of through-space couplings, such as, (1h)J(FH), (1h)J(FN), and (2h)J(FF), where the spin polarization is transmitted through hydrogen bond. In an interesting example of a molecule containing two CF3 groups getting simultaneously involved in hydrogen bond, where hydrogen bond mediated couplings are not reflected in the NMR spectrum, F-19-F-19 NOESY experiment yielded confirmatory evidence. Significant deviations in the strengths of (1)J(NH), variable temperature, and the solvent induced perturbations yielded additional support. The NMR results are corroborated by both DFT calculations and MD simulations, where the quantitative information on different ways of involvement of fluorine in two and three centered hydrogen bonds, their percentage of occurrences, and geometries have been obtained. The hydrogen bond interaction energies have also been calculated.

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.

New insight into the architecture of oxy-anion pocket in unliganded conformation of GAT domains: A MD-simulation study

Human Guanine Monophosphate Synthetase (hGMPS) converts XMP to GMP, and acts as a bifunctional enzyme with N-terminal ``glutaminase'' (GAT) and C-terminal ``synthetase'' domain. The enzyme is identified as a potential target for anticancer and immunosuppressive therapies. GAT domain of enzyme plays central role in metabolism, and contains conserved catalytic residues Cys104, His190, and Glu192. MD simulation studies on GAT domain suggest that position of oxyanion in unliganded conformation is occupied by one conserved water molecule (W1), which also stabilizes that pocket. This position is occupied by a negatively charged atom of the substrate or ligand in ligand bound crystal structures. In fact, MD simulation study of Ser75 to Val indicates that W1 conserved water molecule is stabilized by Ser75, while Thr152, and His190 also act as anchor residues to maintain appropriate architecture of oxyanion pocket through water mediated H-bond interactions. Possibly, four conserved water molecules stabilize oxyanion hole in unliganded state, but they vacate these positions when the enzyme (hGMPS)-substrate complex is formed. Thus this study not only reveals functionally important role of conserved water molecules in GAT domain, but also highlights essential role of other non-catalytic residues such as Ser75 and Thr152 in this enzymatic domain. The results from this computational study could be of interest to experimental community and provide a testable hypothesis for experimental validation. Conserved sites of water molecules near and at oxyanion hole highlight structural importance of water molecules and suggest a rethink of the conventional definition of chemical geometry of inhibitor binding site.

Micropropagación en piña (Ananas comosus (L.) Merr) cultivar MD-2

El presente estudio se realizó en el laboratorio de cultivo de tejidos de la Universidad Nacional Agraria de abril del 2007 a marzo del 2008, con el objetivo de mejorar la eficiencia de producción de piña (Ananas comosus (L.) Merr.) en diferentes fases de la micropropagación. En establecimiento, enraizamiento y aclimatización se evaluaron medios de cultivos con difere ntes concentraciones de reguladores de crecimiento. En la multiplicación se estudió el tipo de consistencia del medio, frasco y número de explantes por frasco. En la fase de establecimiento los medios de cultivo que contenían las sales de Murashige y Skoog, (1962) con 1 mg/l BAP y 0.1 mg/l AIA ó ANA favorecieron la separación de los primordios de hojas; el color verde se presentó en los medios con ANA y BAP. En la fase de multiplicación los mejores medios de cultivo contenían 2 y 3 mg/l BAP de consistencia líquida, con brotes de 0.96 y 1.20. Los medios de cultivo que contenían las sales MS de consistencia semi sólida fueron superiores en número de hojas y raíces. La formación de brotes fue superior al usar frascos de 220 ml con 8 ex plantes en medio semisólido con un promedio de 3.43. La longitud de planta y el número de hojas presentaron promedios superiores en medio semisólido con 5 explantes y en frascos de 220 m l y 430 ml, respectivamente. La inducción de raíces fue favorecida en los medios de cultivo que contenían 1 mg/l de BAP con 40 g/l de sacarosa. Para la variable número de brotes, hojas y longitud de planta no hubo diferencia significativa entre los tratamientos. Las vitroplantas provenientes de plantas enteras presentaron 100% de sobrevivencia. El promedio de espinas fue más acentuado en los tratamientos de plantas partidas longitudinalmente. El número de hoja resultó similar estadísticamente en los medios con diferentes niveles de sacarosa en combinación con 0.1 y 1 mg/l de AIA. La longitud de planta obtuvo el mejor promedio en el medio que contenía 0.5 mg/l de AIA con 30 g/l de sacarosa.

Embriogenesis somatica en reproduccion in vitro del cultivo de piña (Ananas comosus L. Merr) cultivar MD-2

La piña es reconocida como una de las frutas más finas de las re giones tropicales y se denomina la reina de todas las frutas. En Nicaragua se consume principalmente como frut a fresca, refresco y otros productos derivados de la industrialización de la piña. Para el estudio de embriogénesis somática en piña cultivar MD - 2 se realizaron diferentes variantes de medios de cultivo en las fases de inducción, multiplicación, formación de plantas , enraizamiento y aclimatación de plantas. En la fase de inducción de callo s con el empleo de explante s de hojas no se produjo iniciación de callos, pero con explantes de ápices meristemáticos se obtuvieron los mejores resultados en los medios qu e contenían 1 mg/l de 2,4 - D y 0 y 0.2 mg/l 6 - B encil - amino - purina . La multiplicación de callos fue mejor en los medios de cultivos que se le suministraron 2 mg/l de Ácido Indolacético y 1.8 mg/l de A cido - naftaleno - acético con agua de coco al 1 5% y 400 mg/ l de carbón activado . La formaci ón de pla n tas se favoreció en el medio de cultivo que contenía 0.38 mg/l de Ácido Indolbutírico y dos mg/l de 6 - Bencil - Amino - Purina. En l a fase de enraizamiento de las plantas resultó ser mejor el medio de cultivo que con tenía 0.5 mg/l de Ácido Indolacé tico y 40 g/l de sacarosa . En la fase de aclimat iz ación de plantas los medios que previamente se desarrollaron en las cinco variantes de medios en la fase de en raizamiento presentaron los mejores resultados con respecto al po rcen taje de sobrevivencia.

Md Simulation Of Dislocation Mobility During Cutting With Diamond Tip On Silicon

By means of Tersoff and Morse potentials, a three-dimensional molecular dynamics simulation is performed to study atomic force microscopy cutting on silicon monocrystal surface. The interatomic forces between the workpiece and the pin tool and the atoms of workpiece themselves are simulated. Two partial edge dislocations are introduced into workpiece Si, it is found that the motion of dislocations does not occur during the atomic force microscopy cutting processing. Simulation results show that the shear stress acting on dislocations is far below the yield strength of Si. (c) 2008 Elsevier Ltd. All rights reserved.

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.