Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Electronic properties of liquid hydrogen fluoride: A sequential quantum mechanical/Born-Oppenheimer molecular dynamics approach

The electronic properties of liquid hydrogen fluoride (HF) were investigated by carrying out sequential quantum mechanics/Born-Oppenheimer molecular dynamics. The structure of the liquid is in good agreement with recent experimental information. Emphasis was placed on the analysis of polarisation effects, dynamic polarisability and electronic excitations in liquid HF. Our results indicate an increase in liquid phase of the dipole moment (similar to 0.5 D) and isotropic polarisability (5%) relative to their gas-phase values. Our best estimate for the first vertical excitation energy in liquid HF indicates a blue-shift of 0.4 +/- 0.2 eV relative to that of the gas-phase monomer (10.4 eV). (C) 2010 Elsevier B.V. All rights reserved.

Oxygen dynamics and the effects of heat treatments in SBCO samples investigated by mechanical spectroscopy

The compound SmBa(2)Cu(3)O(7-delta) (SBCO)-obtained by substituting rare-earth Sm for Y in the well-known and most studied YBa(2)Cu(3)O(6+delta) (YBCO)-is potentially attractive to study in order to understand the superconductivity mechanism in physics and in electronic device applications. For SBCO, the possibility of variable stoichiometry and the high mobility of oxygen in CuO(x) planes give rise to a rich phase diagram. This study reports on the effect of heat treatments in an oxygen atmosphere on the anelastic properties of this oxide, in which relaxation processes were observed, attributed to oxygen atom jumps present in the Cu-O planes during the orthorhombic phase.

Dynamics of meso and thermo citrate synthases with implicit solvation

The dynamics of hydration of meso and thermo citrate synthases has been investigated using the EEF1 methodology implemented with the CHARNM program. The native enzymes are composed of two identical subunits, each divided into a small and large domain. The dynamics behavior of both enzymes at 30 degrees C and 60 degrees C has been compared. The results of simulations show that during the hydration process, each subunit follows a different pathway of hydration, in spite of the identical sequence. The hydrated structures were compared with the crystalline structure, and the root mean square deviation (RMSD) of each residue along the trajectory was calculated. The regions with larger and smaller mobility were identified. In particular, helices belonging to the small domain are more mobile than those of the large domain. In contrast, the residues that constitute the active site show a much lower displacement compared with the crystalline structure. Hydration free energy calculations point out that Thermoplasma acidophilum citrate synthase (TCS) is more stable than chicken citrate synthase (CCS), at high temperatures. Such result has been ascribed to the higher number of superficial charges in the thermophilic homologue, which stabilizes the enzyme, while the mesophilic homologue denatures. These results are in accord with the experimental found that TCS keeps activity at temperatures farther apart from the catalysis regular temperature than the CCS. (c) 2005 Wiley Periodicals, Inc.

Mechanical properties and fracture dynamics of silicene membranes

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Oxygen Interstitial Dynamics in Bi(2)Sr(2)CaCu(2)O(8+delta) Oxides Measured by Mechanical Spectroscopy

Many researchers became interested in the discovery of Bi(2)Sr(2)CaCu(2)O(8+delta) oxides with critical temperature of around 80 K. It is known that the critical temperature is related to the CuO2 planes of the material. For this reason, the study of the interstitial oxygen in these oxides is of great relevance. The samples were prepared by means of conventional solid state reactions, through the stoichiometric mixture of precursory powders. After the sinterization, the samples were submitted to measurements of density, electrical resistivity, x-ray diffraction, scanning electron microscopy and energy dispersion spectroscopy, with the objective of performing their characterization. The measurements of mechanical spectroscopy were performed by a torsion pendulum. The results show three relaxation processes in the temperature range of 200 and 700 K, with activation energy of approximately 0.9 eV, which has been attributed to the dynamics of the interstitial oxygen present in the material.

Spectroscopic properties from hybrid QM, MM molecular dynamics simulation

The central aim of this thesis work is the application and further development of a hybrid quantum mechanical/molecular mechanics (QM/MM) based approach to compute spectroscopic properties of molecules in complex chemical environments from electronic structure theory. In the framework of this thesis, an existing density functional theory implementation of the QM/MM approach is first used to calculate the nuclear magnetic resonance (NMR) solvent shifts of an adenine molecule in aqueous solution. The findings show that the aqueous solvation with its strongly fluctuating hydrogen bond network leads to specific changes in the NMR resonance lines. Besides the absolute values, also the ordering of the NMR lines changes under the influence of the solvating water molecules. Without the QM/MM scheme, a quantum chemical calculation could have led to an incorrect assignment of these lines. The second part of this thesis describes a methodological improvement of the QM/MM method that is designed for cases in which a covalent chemical bond crosses the QM/MM boundary. The development consists in an automatized protocol to optimize a so-called capping potential that saturates the electronic subsystem in the QM region. The optimization scheme is capable of tuning the parameters in such a way that the deviations of the electronic orbitals between the regular and the truncated (and "capped") molecule are minimized. This in turn results in a considerable improvement of the structural and spectroscopic parameters when computed with the new optimized capping potential within the QM/MM technique. This optimization scheme is applied and benchmarked on the example of truncated carbon-carbon bonds in a set of small test molecules. It turns out that the optimized capping potentials yield an excellent agreement of NMR chemical shifts and protonation energies with respect to the corresponding full molecules. These results are very promising, so that the application to larger biological complexes will significantly improve the reliability of the prediction of the related spectroscopic properties.

Predicting the Mechanical Properties of Carbon-Based Materials Using Molecular Dynamics

The development of innovative carbon-based materials can be greatly facilitated by molecular modeling techniques. Although molecular modeling has been used extensively to predict elastic properties of materials, modeling of more complex phenomenon such as fracture has only recently been possible with the development of new force fields such as ReaxFF, which is used in this work. It is not fully understood what molecular modeling parameters such as thermostat type, thermostat coupling, time step, system size, and strain rate are required for accurate modeling of fracture. Selection of modeling parameters to model fracture can be difficult and non-intuitive compared to modeling elastic properties using traditional force fields, and the errors generated by incorrect parameters may be non-obvious. These molecular modeling parameters are systematically investigated and their effects on the fracture of well-known carbon materials are analyzed. It is determined that for coupling coefficients of 250 fs and greater do not result in substantial differences in the stress-strain response of the materials using any thermostat type. A time step of 0.5 fs of smaller is required for accurate results. Strain rates greater than 2.2 ns-1 are sufficient to obtain repeatable results with slower strain rates for the materials studied. The results of this study indicate that further refinement of the Chenoweth parameter set is required to accurately predict the mechanical response of carbon-based systems. The ReaxFF has been used extensively to model systems in which bond breaking and formation occur. In particular ReaxFF has been used to model reactions of small molecules. Some elastic and fracture properties have been successfully modeled using ReaxFF in materials such as silicon and some metals. However, it is not clear if current parameterizations for ReaxFF are able to accurately reproduce the elastic and fracture properties of carbon materials. The stress-strain response of a new ReaxFF parameterization is compared to the previous parameterization and density functional theory results for well-known carbon materials. The new ReaxFF parameterization makes xv substantial improvements to the predicted mechanical response of carbon materials, and is found to be suitable for modeling the mechanical response of carbon materials. Finally, a new material composed of carbon nanotubes within an amorphous carbon (AC) matrix is modeled using the ReaxFF. Various parameters that may be experimentally controlled are investigated such as nanotube bundling, comparing multi-walled nanotube with single-walled nanotubes, and degree of functionalization of the nanotubes. Elastic and fracture properties are investigated for the composite systems and compared to results of pure-nanotube and pure-AC models. It is found that the arrangement of the nanotubes and degree of crosslinking may substantially affect the properties of the systems, particularly in the transverse directions.

Review of recent results using computational fluid dynamics simulations in patients receiving mechanical assist devices for end-stage heart failure

Many end-stage heart failure patients are not eligible to undergo heart transplantation due to organ shortage, and even those under consideration for transplantation might suffer long waiting periods. A better understanding of the hemodynamic impact of left ventricular assist devices (LVAD) on the cardiovascular system is therefore of great interest. Computational fluid dynamics (CFD) simulations give the opportunity to study the hemodynamics in this patient population using clinical imaging data such as computed tomographic angiography. This article reviews a recent study series involving patients with pulsatile and constant-flow LVAD devices in which CFD simulations were used to qualitatively and quantitatively assess blood flow dynamics in the thoracic aorta, demonstrating its potential to enhance the information available from medical imaging.

Simulating the remote handling of the blanket segments in DEMO fusion reactor with thermo-mechanical meshfree multibody dynamics

The present study shows a first approach to the simulation of the remote handling oper- ation which takes into account the thermal and flexible behavior of the blanket segments and its implications on the remote handling equipment, in order to validate and improve its design.

Five hundred and seven mechanical movements : embracing all those which are most important in dynamics, hydraulics, hydrostatics, pneumatics, steam engines, mill and other gearing, presses, horology, and miscellaneous machinery : and including many movements never before published, and several which have only recently come into use /

Mode of access: Internet.

Towards a better understanding of granular material processes: The coupling of a computational fluid dynamics framework with micro-mechanical models

Abstract not available

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties. (c) 2008 American Institute of Physics.

Neuromuscular-skeletal constraints upon the dynamics of unimanual and bimanual coordination

In the first of three experiments, 11 participants generated pronation and supination movements of the forearm, in time with an auditory metronome. The metronome frequency was increased in eight steps (0.25 Hz) from a base frequency of 1.75 Hz. On alternating trials, participants were required to coordinate either maximum pronation or maximum supination with each beat of the metronome. In each block of trials, the axis of rotation was either coincident with the long axis of the forearm, above this axis, or below this axis. The stability of the pronate-on-the-beat pattern, as indexed by the number of pattern changes, and the time of onset of pattern change, was greatest when the axis of rotation of the movement was below the long axis of the forearm. In contrast, the stability of the supinate-on-the-beat pattern was greatest when the axis of rotation of the movement was above the long axis of the forearm. In a second experiment, we examined how changes in the position of the axis of rotation alter the activation patterns of muscles that contribute to pronation and supination of the forearm. Variations in the relative dominance of the pronation and supination phases of the movement cycle across conditions were accounted for primarily by changes in the activation profile of flexor carpi radialis (FCR) and extensor carpi radialis longus (ECR). In the Final experiment we examined how these constraints impact upon the stability of bimanual coordination. Thirty-two participants were assigned at random to one of four conditions, each of which combined an axis of rotation configuration (bottom or top) for each limb. The participants generated both inphase (both limbs pronating simultaneously, and supinating simultaneously) and antiphase (left limb pronating and right limb supinating simultaneously, and vice versa) patterns of coordination. When the position of the axis of rotation was equivalent for the left and the right limb, transitions from antiphase to inphase patterns of coordination were Frequently observed. In marked contrast, when the position of the axis of rotation for the left and right limb was contradistinct, transitions From inphase to antiphase patterns of coordination occurred. The results demonstrated that when movements are performed in an appropriate mechanical context, inphase patterns of coordination are less stable than antiphase patterns.

Two-dimensional approximation for tidal dynamics in coastal aquifers: Capillarity correction

Most previous investigations on tide-induced watertable fluctuations in coastal aquifers have been based on one-dimensional models that describe the processes in the cross-shore direction alone, assuming negligible along-shore variability. A recent study proposed a two-dimensional approximation for tide-induced watertable fluctuations that took into account coastline variations. Here, we further develop this approximation in two ways, by extending the approximation to second order and by taking into account capillary effects. Our results demonstrate that both effects can markedly influence watertable fluctuations. In particular, with the first-order approximation, the local damping rate of the tidal signal could be subject to sizable errors.