Coarse grained force field for the molecular simulation of natural gases and condensates

AbstractThe atomistically-detailed molecular modelling of petroleum fluids is challenging, amongst other aspects, due to the very diverse multicomponent and asymmetric nature of the mixtures in question. Complicating matters further, the time scales for many important processes can be much larger than the current and foreseeable capacity of modern computers running fully-atomistic models. To overcome these limitations, a coarse grained (CG) model is proposed where some of the less-important degrees of freedom are safely integrated out, leaving as key parameters the average energy levels, the molecular conformations and the range of the Mie intermolecular potentials employed as the basis of the model. The parametrization is performed by using an analytical equation of state of the statistical associating fluid theory (SAFT) family to link the potential parameters to macroscopically observed thermophysical properties. The parameters found through this top-down approach are used directly in molecular dynamics simulations of multi-component multi-phase systems. The procedure is exemplified by calculating the phase envelope of the methane–decane binary and of two synthetic light condensate mixtures. A methodology based on the discrete expansion of a mixture is used to determine the bubble points of these latter mixtures, with an excellent agreement to experimental data. The model presented is entirely predictive and an abridged table of parameters for some fluids of interest is provided

Functional dynamics of the folded ankyrin repeats of I kappa B alpha revealed by nuclear magnetic resonance.

Inhibition of nuclear factor kappaB (NF-kappaB) is mainly accomplished by IkappaB alpha, which consists of a signal response sequence at the N-terminus, a six-ankyrin repeat domain (ARD) that binds NF-kappaB, and a C-terminal PEST sequence. Previous studies with the ARD revealed that the fifth and sixth repeats are only partially folded in the absence of NF-kappaB. Here we report NMR studies of a truncated version of IkappaB alpha, containing only the first four ankyrin repeats, IkappaB alpha(67-206). This four-repeat segment is well-structured in the free state, enabling full resonance assignments to be made. H-D exchange, backbone dynamics, and residual dipolar coupling (RDC) experiments reveal regions of flexibility. In addition, regions consistent with the presence of micro- to millisecond motions occur periodically throughout the repeat structure. Comparison of the RDCs with the crystal structure gave only moderate agreement, but an ensemble of structures generated by accelerated molecular dynamics gave much better agreement with the measured RDCs. The regions showing flexibility correspond to those implicated in entropic compensation for the loss of flexibility in ankyrin repeats 5 and 6 upon binding to NF-kappaB. The regions showing micro- to millisecond motions in the free protein are the ends of the beta-hairpins that directly interact with NF-kappaB in the complex

Aqueous Amino Acids and Proteins Near the Surface of Gold in Hydrophilic and Hydrophobic Force Fields

We calculate potentials of the mean force for twenty amino acids in the vicinity of the (111) surface of gold, for several dipeptides, and for some analogs of the side chains, using molecular dynamics simulations and the umbrella sampling method. We compare results obtained within three different force fields: one hydrophobic (for a contaminated surface) and two hydrophilic. All of these fields lead to good binding with very different specificities and different patterns in the density and polarization of water. The covalent bond with the sulfur atom on cysteine is modeled by the Morse potential. We demonstrate that binding energies of dipeptides are different than the combined binding energies of their amino-acidic components. For the hydrophobic gold, adsorption events of a small protein are driven by attraction to the strongest binding amino acids. This is not so in the hydrophilic cases - a result of smaller specificities combined with the difficulty for proteins, but not for single amino acids, to penetrate the first layer of water. The properties of water near the surface sensitively depend on the force field