40 resultados para Molecular-dynamics simulations
Resumo:
In this work, we have used a combined of atomistic simulation methods to explore the effects of confinement of water molecules between silica surfaces. Firstly, the mechanical properties of water severe confined (~3A) between two silica alpha-quartz was determined based on first principles calculations within the density functional theory (DFT). Simulated annealing methods were employed due to the complex potential energry surface, and the difficulties to avoid local minima. Our results suggest that much of the stiffness of the material (46%) remains, even after the insertion of a water monolayer in the silica. Secondly, in order to access typical time scales for confined systems, classical molecular dynamics was used to determine the dynamical properties of water confined in silica cylindrical pores, with diameters varying from 10 to 40A. in this case we have varied the passivation of the silica surface, from 13% to 100% of SiOH, and the other terminations being SiOH2 and SiOH3, the distribution of the different terminations was obtained with a Monte Carlo simulation. The simulations indicates a lowering of the diffusion coefficientes as the diameter decreases, due to the structuration of hydrogen bonds of water molecules; we have also obtained the density profiles of the confined water and the interfacial tension.
Resumo:
This paper reports results for directed flow v(1) and elliptic flow v(2) of charged particles in Cu + Cu collisions at root s(NN) = 22.4 GeV at the Relativistic Heavy Ion Collider. The measurements are for the 0-60% most central collisions, using charged particles observed in the STAR detector. Our measurements extend to 22.4-GeV Cu + Cu collisions the prior observation that v1 is independent of the system size at 62.4 and 200 GeV and also extend the scaling of v(1) with eta/y(beam) to this system. The measured v(2)(p(T)) in Cu + Cu collisions is similar for root s(NN) throughout the range 22.4 to 200 GeV. We also report a comparison with results from transport model (ultrarelativistic quantum molecular dynamics and multiphase transport model) calculations. The model results do not agree quantitatively with the measured v(1)(eta), v(2)(p(T)), and v(2)(eta).
Resumo:
Background: This paper addresses the prediction of the free energy of binding of a drug candidate with enzyme InhA associated with Mycobacterium tuberculosis. This problem is found within rational drug design, where interactions between drug candidates and target proteins are verified through molecular docking simulations. In this application, it is important not only to correctly predict the free energy of binding, but also to provide a comprehensible model that could be validated by a domain specialist. Decision-tree induction algorithms have been successfully used in drug-design related applications, specially considering that decision trees are simple to understand, interpret, and validate. There are several decision-tree induction algorithms available for general-use, but each one has a bias that makes it more suitable for a particular data distribution. In this article, we propose and investigate the automatic design of decision-tree induction algorithms tailored to particular drug-enzyme binding data sets. We investigate the performance of our new method for evaluating binding conformations of different drug candidates to InhA, and we analyze our findings with respect to decision tree accuracy, comprehensibility, and biological relevance. Results: The empirical analysis indicates that our method is capable of automatically generating decision-tree induction algorithms that significantly outperform the traditional C4.5 algorithm with respect to both accuracy and comprehensibility. In addition, we provide the biological interpretation of the rules generated by our approach, reinforcing the importance of comprehensible predictive models in this particular bioinformatics application. Conclusions: We conclude that automatically designing a decision-tree algorithm tailored to molecular docking data is a promising alternative for the prediction of the free energy from the binding of a drug candidate with a flexible-receptor.
Resumo:
The structure of gold-platinum nanoparticles is heavily debated as theoretical calculations predict core-shell particles, whereas x-ray diffraction experiments frequently detect randomly mixed alloys. By calculating the structure of gold-platinum nanoparticles with diameters of up to approximate to 3.5 nm and simulating their x-ray diffraction patterns, we show that these seemingly opposing findings need not be in contradiction: Shells of gold are hardly visible in usual x-ray scattering, and the interpretation of Vegard's law is ambiguous on the nanoscale. DOI: 10.1103/PhysRevB.86.241403
Resumo:
By combining first-principles electronic structure calculations and existing time-differential gamma-gamma perturbed-angular-correlation experiments we studied the site localization, the local environment, and the electronic structure of Cd impurities in sapphire (alpha-Al2O3) single crystals in different charged states. The ab initio calculations were performed with the full-potential augmented plane wave plus local orbitals method and the projector augmented wave method. Comparing the calculated electric-field-gradient tensor at the Cd nuclei in the alpha-Al2O3 host lattice and the corresponding available experimental values, we have seen that it is equally possible for Cd to replace an Al atom (in a negative charge state) or to be placed in an interstitial site (in a neutral charge state). To finally address the issue of the Cd impurity localization, we performed formation energy calculations. These results have shown that Cd placed in the substitutional Al site, in the negatively charged state, is the most probable configuration.
Resumo:
Two new peptidic proteasome inhibitors were isolated as trace components from a Curacao collection of the marine cyanobacterium Symploca sp. Carmaphycin A (1) and carmaphycin B (2) feature a leucine-derived a,beta-epoxyketone warhead directly connected to either methionine sulfoxide or methionine sulfone. Their structures were elucidated on the basis of extensive NMR and MS analyses and confirmed by total synthesis, which in turn provided more material for further biological evaluations. Pure carmaphycins A and B were found to inhibit the beta 5 subunit (chymotrypsin-like activity) of the S. cerevisiae 20S proteasome in the low nanomolar range. Additionally, they exhibited strong cytotoxicity to lung and colon cancer cell lines, as well as exquisite antiproliferative effects in the NCI60 cell-line panel. These assay results as well as initial structural biology studies suggest a distinctive binding mode for these new inhibitors.
Resumo:
Lewy bodies and Lewy neurites, neuropathological hallmarks of several neurological diseases, are mainly made of filamentous assemblies of alpha-synuclein. However, other macromolecules including Tau, ubiquitin, glyceraldehyde-3-phosphate dehydrogenase, and glycosaminoglycans are routinely found associated with these amyloid deposits. Glyceraldehyde-3-phosphate dehydrogenase is a glycolytic enzyme that can form fibrillar aggregates in the presence of acidic membranes, but its role in Parkinson disease is still unknown. In this work, the ability of heparin to trigger the amyloid aggregation of this protein at physiological conditions of pH and temperature is demonstrated by infrared and fluorescence spectroscopy, dynamic light scattering, small angle x-ray scattering, circular dichroism, and fluorescence microscopy. Aggregation proceeds through the formation of short rod-like oligomers, which elongates in one dimension. Heparan sulfate was also capable of inducing glyceraldehyde-3-phosphate dehydrogenase aggregation, but chondroitin sulfates A, B, and C together with dextran sulfate had a negligible effect. Aided with molecular docking simulations, a putative binding site on the protein is proposed providing a rational explanation for the structural specificity of heparin and heparan sulfate. Finally, it is demonstrated that in vitro the early oligomers present in the glyceraldehyde-3-phosphate dehydrogenase fibrillation pathway promote alpha-synuclein aggregation. Taking into account the toxicity of alpha-synuclein prefibrillar species, the heparin-induced glyceraldehyde-3-phosphate dehydrogenase early oligomers might come in useful as a novel therapeutic strategy in Parkinson disease and other synucleinopathies.
The boundedness of penalty parameters in an augmented Lagrangian method with constrained subproblems
Resumo:
Augmented Lagrangian methods are effective tools for solving large-scale nonlinear programming problems. At each outer iteration, a minimization subproblem with simple constraints, whose objective function depends on updated Lagrange multipliers and penalty parameters, is approximately solved. When the penalty parameter becomes very large, solving the subproblem becomes difficult; therefore, the effectiveness of this approach is associated with the boundedness of the penalty parameters. In this paper, it is proved that under more natural assumptions than the ones employed until now, penalty parameters are bounded. For proving the new boundedness result, the original algorithm has been slightly modified. Numerical consequences of the modifications are discussed and computational experiments are presented.
Resumo:
Structural and electronic properties of the PtnTM55-n (TM = Co, Rh, Au) nanoalloys are investigated using density functional theory within the generalized gradient approximation and employing the all-electron projected augmented wave method. For TM = Co and Rh, the excess energy, which measures the relative energy stability of the nanoalloys, is negative for all Pt compositions. We found that the excess energy has similar values for a wide range of Pt compositions, i.e., n = 20-42 and n = 28-42 for Co and Rh, respectively, with the core shell icosahedron-like configuration (n = 42) being slightly more stable for both Co and Rh systems because of the larger release of the strain energy due to the smaller atomic size of the Co and Rh atoms. For TM = Au, the excess energy is positive for all compositions, except for n = 13, which is energetically favorable due to the formation of the core-shell structure (Pt in the core and Au atoms at the surface). Thus, our calculations confirm that the formation of core-shell structures plays an important role to increase the stability of nanoalloys. The center of gravity of the occupied d-states changes almost linearly as a function of the Pt composition, and hence, based on the d-band model, the magnitude of the adsorption energy of an adsorbate can be tuned by changing the Pt composition. The magnetic moments of PtnCo55-n decrease almost linearly as a function of the Pt composition; however, the same does not hold for PtRh and PtAu. We found an enhancement of the magnetic moments of PtRh by a few times by increasing Pt composition, which we explain by the compression effects induced by the large size of the Pt atoms compared with the Rh atoms.
Resumo:
The topographical character of conical intersections (CIs)-either sloped or peaked-has played a fundamental and important role in the discussion of the efficiency of CIs as photochemical "funnels." Here this perspective is employed in connection with a recent study of a model protonated Schiff base (PSB) cis to trans photoisomerization in solution [Malhado et al., J. Phys. Chem. A 115, 3720 (2011)]. In that study, the calculated reduced photochemical quantum yield for the successful production of trans product versus cis reactant in acetonitrile solvent compared to water was interpreted in terms of a dynamical solvent effect related to the dominance, for the acetonitrile case, of S-1 to S-0 nonadiabatic transitions prior to the reaching the seam of CIs. The solvent influence on the quantum yield is here re-examined in the sloped/peaked CI topographical perspective via conversion of the model's two PSB internal coordinates and a nonequilibrium solvent coordinate into an effective branching space description, which is then used to re-analyze the generalized Langevin equation/surface hopping results. The present study supports the original interpretation and enriches it in terms of topographical detail. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4754505]