929 resultados para MOLECULAR SIMULATIONS
Resumo:
Estimation of the dissociation constant, or pK(a), of weak acids continues to be a central goal in theoretical chemistry. Here we show that ab initio Car-Parrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free energy profile of the dissociation reaction can provide reasonable estimates of the successive pK(a) values of polyprotic acids. We use the distance-dependent coordination number of the protons bound to the hydroxyl oxygen of the carboxylic group as the collective variable to explore the free energy profile of the dissociation process. Water molecules, sufficient to complete three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. Two distinct minima corresponding to the dissociated and un-dissociated states of the acid are observed and the difference in their free energy values provides the estimate for pK(a), the acid dissociation constant. We show that the method predicts the pK(a) value of benzoic acid in good agreement with experiment and then show using phthalic acid (benzene dicarboxylic acid) as a test system that both the first and second pK(a) values as well, as the subtle difference in their values for different isomers can be predicted in reasonable agreement with experimental data.
Resumo:
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.
Resumo:
The de novo purine biosynthesis is one of the highly conserved pathways among all organisms and is essential for the cell viability. A clear understanding of the enzymes in this pathway would pave way for the development of antimicrobial and anticancer drugs. Phosphoribosylaminoimidazole-succinocar boxamide (SAICAR) synthetase is one of the enzymes in this pathway that catalyzes ATP dependent ligation of carboxyaminoimidazole ribotide (CAIR) with L-aspartate (ASP). Here, we describe eight crystal structures of this enzyme, in C222(1) and H3 space groups, bound to various substrates and substrate mimics from a hyperthermophilic archaea Pyrococcus horikoshii along with molecular dynamics simulations of the structures with substrates. Complexes exhibit minimal deviation from its apo structure. The CAIR binding site displays a preference for pyrimidine nucleotides. In the ADP.TMP-ASP complex, the ASP binds at a position equivalent to that found in Saccharomyces cerevisiae structure (PDB: 2CNU) and thus, clears the ambiguity regarding ASP's position. A possible mode for the inhibition of the enzyme by CTP and UTP, observed earlier in the yeast enzyme, is clearly illustrated in the structures bound to CMP and UMP. The ADP.Mg2+.PO4.CD/MP complex having a phosphate ion between the ATP and CAIR sites strengthens one of the two probable pathways (proposed in Escherichia coli study) of catalytic mechanism and suggests the possibility of a phosphorylation taking place before the ASP's attack on CAIR. Molecular dynamic simulations of this enzyme along with its substrates at 90 degrees C reveal the relative strengths of substrate binding, possible antagonism and the role of Mg2+ ions. (C) 2015 Elsevier Inc. All rights reserved.
Resumo:
We report molecular dynamics (MD) simulations to explore the influence of a counterion on the structure and dynamics of cationic and anionic solvation shells for various ions in methanol at 298 K. We show that the variation in ionic size of either the cation or the anion in an ion pair influences the solvation structure of the other ion as well as the diffusivity in an electrolyte solution of methanol. The extent of ionic association between the cation and its counteranion of different ionic sizes has been investigated by analyzing the radial distribution functions (RDFs) and the orientation of methanol molecules in the first solvation shell (FSS) of ions. It is shown that the methanol in the FSS of the anion as well the cation exhibit quite different radial and orientational structures as compared to methanol which lie in the FSS of either the anion or the cation but not both. We find that the coordination number (CN) of F-, Cr-, and I- ions decreases with increasing size of the anion which is contrary to the trend reported for the anions in H2O. The mean residence time (MRT) of methanol molecules in the FSS of ions has been calculated using the stable states picture (SSP) approach. It is seen that the ion-counterion interaction has a considerable influence on the MRT of methanol molecules in the FSS of ions. We also discuss the stability order of the ion-counterion using the potentials of mean force (PMFs) for ion pairs with ions of different sizes. The PMF plots reveal that the Li+-F- pair (small-small) is highly stable and the Li+-I- pair is least stable (small-large) in electrolyte solutions.
Resumo:
A new class of dendrimers, the poly(propyl ether imine) (PETIM) dendrimer, has been shown to be a novel hyperbranched polymer having potential applications as a drug delivery vehicle. Structure and dynamics of the amine terminated PETIM dendrimer and their changes with respect to the dendrimer generation are poorly understood. Since most drugs are hydrophobic in nature, the extent of hydrophobicity of the dendrimer core is related to its drug encapsulation and retention efficacy. In this study, we carry out fully atomistic molecular dynamics (MD) simulations to characterize the structure of PETIM (G2-G6) dendrimers in salt solution as a function of dendrimer generation at different protonation levels. Structural properties such as radius of gyration (R-g), radial density distribution, aspect ratio, and asphericity are calculated. In order to assess the hydrophilicity of the dendrimer, we compute the number of bound water molecules in the interior of dendrirner as well as the number of dendrimer-water hydrogen bonds. We conclude that PETIM dendrimers have relatively greater hydrophobicity and flexibility when compared with their extensively investigated PAMAM counterparts. Hence PETIM dendrimers are expected to have stronger interactions with lipid membranes as well as improved drug encapsulation and retention properties when compared with PAMAM dendrimers. We compute the root-mean-square fluctuation of dendrimers as well as their entropy to quantify the flexibility of the dendrimer. Finally we note that structural and solvation properties computed using force field parameters derived based on the CHARMM general purpose force field were in good quantitative agreement with those obtained using the generalized Amber force field (GAFF).
Resumo:
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.
Resumo:
Molecular dynamics simulations were employed to investigate the specimen thickness-dependent tensile behavior of a series of Cu(x)Z(100-x) (x = 20, 40, 50, 64 and 80 at%) metallic glass (MG) films, with a particular focus on the critical thickness, tc, below which non-localized plastic flow takes place. The simulation results reveal that while the transition occurs in all the alloys examined, t(c) is sensitive to the composition. We rationalize t(c) by postulating that the strain energy stored in the sample at the onset of plastic deformation has to be sufficient for the formation of shear bands. The composition-dependence of t(c) was found to correlate with the average activation energy of the atomic level plastic deformation events. (C) 2015 Elsevier Ltd. All rights reserved.
Resumo:
CucurbitacinE (CurE) has been known to bind covalently to F-actin and inhibit depolymerization. However, the mode of binding of CurE to F-actin and the consequent changes in the F-actin dynamics have not been studied. Through quantum mechanical/molecular mechanical (QM/MM) and density function theory (DFT) simulations after the molecular dynamics (MD) simulations of the docked complex of F-actin and CurE, a detailed transition state (TS) model for the Michael reaction is proposed. The TS model shows nucleophilic attack of the sulphur of Cys257 at the beta-carbon of Michael Acceptor of CurE producing an enol intermediate that forms a covalent bond with CurE. The MD results show a clear difference between the structure of the F-actin in free form and F-actin complexed with CurE. CurE affects the conformation of the nucleotide binding pocket increasing the binding affinity between F-actin and ADP, which in turn could affect the nucleotide exchange. CurE binding also limits the correlated displacement of the relatively flexible domain 1 of F-actin causing the protein to retain a flat structure and to transform into a stable ``tense'' state. This structural transition could inhibit depolymerization of F-actin. In conclusion, CurE allosterically modulates ADP and stabilizes F-actin structure, thereby affecting nucleotide exchange and depolymerization of F-actin. (C) 2015 Elsevier Inc. All rights reserved.
Resumo:
The serotonin(1A) receptor belongs to the superfamily of G protein-coupled receptors (GPCRs) and is a potential drug target in neuropsychiatric disorders. The receptor has been shown to require membrane cholesterol for its organization, dynamics and function. Although recent work suggests a close interaction of cholesterol with the receptor, the structural integrity of the serotonin(1A) receptor in the presence of cholesterol has not been explored. In this work, we have carried out all atom molecular dynamics simulations, totaling to 3s, to analyze the effect of cholesterol on the structure and dynamics of the serotonin(1A) receptor. Our results show that the presence of physiologically relevant concentration of membrane cholesterol alters conformational dynamics of the serotonin(1A) receptor and, on an average lowers conformational fluctuations. Our results show that, in general, transmembrane helix VII is most affected by the absence of membrane cholesterol. These results are in overall agreement with experimental data showing enhancement of GPCR stability in the presence of membrane cholesterol. Our results constitute a molecular level understanding of GPCR-cholesterol interaction, and represent an important step in our overall understanding of GPCR function in health and disease.
Resumo:
Molecular dynamics (MD) simulations are performed to study adhesion and peeling of a short fragment of single strand DNA (ssDNA) molecule from a graphite surface. The critical peel-off force is found to depend on both the peeling angle and the elasticity of ssDNA. For the short ssDNA strand under investigation, we show that the simulation results can be explained by a continuum model of an adhesive elastic band on substrate. The analysis suggests that it is often the peak value, rather than the mean value, of adhesion energy which determines the peeling of a nanoscale material.
Resumo:
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.
Resumo:
A correlative reference model for computer molecular dynamics simulations is proposed. Based on this model, a flexible displacement boundary scheme is introduced and the dislocations emitted from a crack tip can continuously pass through the border of the inner discrete atomic region and pile up at the outer continuum region. The effect of the emitted dislocations within the plastic zone on the inner atomistic region can be clearly demonstrated. The simulations for a molybdinum crystal show that a full dislocation in a bcc crystal is dissociated into three partial dislocations and interaction between the crack and the emitted dislocations results in gradual decrease of the local stress intensity factor.
Resumo:
Molecular dynamics simulations on diffusion bonding of Cu-Ag showed that the thickness of the interfacial region depended on the stress. The interfacial region became amorphous during diffusion bonding, and it would normally transform from amorphous into crystalline structure when the structure was cooled to the room temperature.
Resumo:
A correlative reference model for a computer simulation of molecular dynamics is proposed in this paper. Based on this model, a flexible displacement boundary scheme is naturally introduced and the dislocations emitted from a crack tip are presumed to continuously pass through the border of an inner discrete atomic region to pile up at an outer continuum region. The simulations for a Mo crystal show that the interaction between a crack and emitted dislocations results in the decrease in local stress intensity factor gradually.
Resumo:
Selectin-ligand interactions are crucial to such biological processes as inflammatory cascade or tumor metastasis. How transient formation and dissociation of selectin-ligand bonds in blood flow are coupled to molecular conformation at atomic level, however, has not been well understood. In this study, steered molecular dynamics (SMD) simulations were used to elucidate the intramolecular and intermolecular conformational evolutions involved in forced dissociation of three selectin-ligand systems: the construct consisting of P-selectin lectin (Lec) and epidermal growth factor (EGF)-like domains (P-LE) interacting with synthesized sulfoglycopeptide or SGP-3, P-LE with sialyl Lewis X (sLeX), and E-LE with sLeX. SMD simulations were based on newly built-up force field parameters including carbohydrate units and sulfated tyrosine(s) using an analogy approach. The simulations demonstrated that the complex dissociation was coupled to the molecular extension. While the intramolecular unraveling in P-LESGP-3 system mainly resulted from the destroy of the two anti-parallel sheets of EGF domain and the breakage of hydrogen-bond cluster at the Lec-EGF interface, the intermolecular dissociation was mainly determined by separation of fucose (FUC) from Ca2+ ion in all three systems. Conformational changes during forced dissociations depended on pulling velocities and forces, as well as on how the force was applied. This work provides an insight into better understanding of conformational changes and adhesive functionality of selectin-ligand interactions under external forces.