Liquid-liquid phase transition for soft-core attractive potential

Using event-driven molecular dynamics simulations, we study a three-dimensional one-component system of spherical particles interacting via a discontinuous potential combining a repulsive square soft core and an attractive square well. In the case of a narrow attractive well, it has been shown that this potential has two metastable gas-liquid critical points. Here we systematically investigate how the changes of the parameters of this potential affect the phase diagram of the system. We find a broad range of potential parameters for which the system has both a gas-liquid critical point C1 and a liquid-liquid critical point C2. For the liquid-gas critical point we find that the derivatives of the critical temperature and pressure, with respect to the parameters of the potential, have the same signs: they are positive for increasing width of the attractive well and negative for increasing width and repulsive energy of the soft core. This result resembles the behavior of the liquid-gas critical point for standard liquids. In contrast, for the liquid-liquid critical point the critical pressure decreases as the critical temperature increases. As a consequence, the liquid-liquid critical point exists at positive pressures only in a finite range of parameters. We present a modified van der Waals equation which qualitatively reproduces the behavior of both critical points within some range of parameters, and gives us insight on the mechanisms ruling the dependence of the two critical points on the potential¿s parameters. The soft-core potential studied here resembles model potentials used for colloids, proteins, and potentials that have been related to liquid metals, raising an interesting possibility that a liquid-liquid phase transition may be present in some systems where it has not yet been observed.

Dynamic properties of Lennard Jones fluids and liquid metals

The dependence of the dynamic properties of liquid metals and Lennard-Jones fluids on the characteristics of the interaction potentials is analyzed. Molecular-dynamics simulations of liquids in analogous conditions but assuming that their particles interact either through a Lennard-Jones or a liquid-metal potential were carried out. The Lennard-Jones potentials were chosen so that both the effective size of the particles and the depth of the potential well were very close to those of the liquid-metal potentials. In order to investigate the extent to which the dynamic properties of liquids depend on the short-range attractive interactions as well as on the softness of the potential cores, molecular-dynamics simulations of the same systems but assuming purely repulsive interactions with the same potential cores were also performed. The study includes both singleparticle dynamic properties, such as the velocity autocorrelation functions, and collective dynamic properties, such as the intermediate scattering funcfunctions, and collective dynamic properties, such as the intermediate scattering functions, the dynamic structure factors, the longitudinal and transverse current correlations, and the transport coefficients.

Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores

Time correlation functions between the velocity of a tagged particle and velocities of particles within specified ranges of initial separations have been obtained by molecular dynamics simulation. These correlation functions have allowed us to analyze the momentum transfer between particles in different coordination shells. Two simple liquids at very different densities and two purely repulsive potentials with very different softnesses have been considered. The longitudinal correlations, which are the velocity cross-correlations along the initial direction defined by the centers of two given particles, have been calculated separately. It has been proven that these correlations should be attributed to particles both in front of and behind the central one. As with propagating longitudinal modes, they are strongly dependent on the softness of the potential core. Some characteristic features of the velocity correlation functions after the initial rise should be related to nonlongitudinal correlations. It has been shown that velocity cross-correlations between distinct particles cannot only be attributed to the direct interactions among particles, but also to the motions induced by the movement of a tagged particle on their neighbors.

Lack of Enantioselectivity in the SULT1A3-catalyzed Sulfoconjugation of Normetanephrine Enantiomers: An In Vitro and Computational Study.

(1R)-Normetanephrine is the natural stereoisomeric substrate for sulfotransferase 1A3 (SULT1A3)-catalyzed sulfonation. Nothing appears known on the enantioselectivity of the reaction despite its potential significance in the metabolism of adrenergic amines and in clinical biochemistry. We confronted the kinetic parameters of the sulfoconjugation of synthetic (1R)-normetanephrine and (1S)-normetanephrine by recombinant human SULT1A3 to a docking model of each normetanephrine enantiomer with SULT1A3 and the 3'-phosphoadenosine-5'-phosphosulfate cofactor on the basis of molecular modeling and molecular dynamics simulations of the stability of the complexes. The K(M) , V(max) , and k(cat) values for the sulfonation of (1R)-normetanephrine, (1S)-normetanephrine, and racemic normetanephrine were similar. In silico models were consistent with these findings as they showed that the binding modes of the two enantiomers were almost identical. In conclusion, SULT1A3 is not substrate-enantioselective toward normetanephrine, an unexpected finding explainable by a mutual adaptability between the ligands and SULT1A3 through an "induced-fit model" in the catalytic pocket. Chirality, 00:000-000, 2012.© 2012 Wiley Periodicals, Inc.

Diffusive transport properties in monovalent and divalent metal-ion halide melts: A computer simulation study

Self- and cross-velocity correlation functions and related transport coefficients of molten salts are studied by molecular-dynamics simulation. Six representative systems are considered, i.e., NaCl and KCl alkali halides, CuCl and CuBr noble-metal halides, and SrCl2 and ZnCl2 divalent metal-ion halides. Computer simulation results are compared with experimental self-diffusion coefficients and electrical conductivities. Special attention is paid to dynamic cross correlations and their dependence on the Coulomb interactions as well as on the size and mass differences between anions and cations.

Polymorphisms in fibronectin binding protein A of Staphylococcus aureus are associated with infection of cardiovascular devices.

Medical implants, like cardiovascular devices, improve the quality of life for countless individuals but may become infected with bacteria like Staphylococcus aureus. Such infections take the form of a biofilm, a structured community of bacterial cells adherent to the surface of a solid substrate. Every biofilm begins with an attractive force or bond between bacterium and substratum. We used atomic force microscopy to probe experimentally forces between a fibronectin-coated surface (i.e., proxy for an implanted cardiac device) and fibronectin-binding receptors on the surface of individual living bacteria from each of 80 clinical isolates of S. aureus. These isolates originated from humans with infected cardiac devices (CDI; n = 26), uninfected cardiac devices (n = 20), and the anterior nares of asymptomatic subjects (n = 34). CDI isolates exhibited a distinct binding-force signature and had specific single amino acid polymorphisms in fibronectin-binding protein A corresponding to E652D, H782Q, and K786N. In silico molecular dynamics simulations demonstrate that residues D652, Q782, and N786 in fibronectin-binding protein A form extra hydrogen bonds with fibronectin, complementing the higher binding force and energy measured by atomic force microscopy for the CDI isolates. This study is significant, because it links pathogenic bacteria biofilms from the length scale of bonds acting across a nanometer-scale space to the clinical presentation of disease at the human dimension.

Somatostatin Subtype‑2 Receptor-Targeted Metal-Based Anticancer Complexes

Conjugates of a dicarba analogue of octreotide, a potent somatostatin agonist whose receptors are overexpressed on tumor cells, with [PtCl2(dap)] (dap = 1-(carboxylic acid)-1,2-diaminoethane) (3), [(η6-bip)Os(4-CO2-pico)Cl] (bip = biphenyl, pico = picolinate) (4), [(η6-p-cym)RuCl(dap)]+ (p-cym = p-cymene) (5), and [(η6-p-cym)RuCl(imidazole-CO2H)(PPh3)]+ (6), were synthesized by using a solid-phase approach. Conjugates 35 readily underwent hydrolysis and DNA binding, whereas conjugate 6 was inert to ligand substitution. NMR spectroscopy and molecular dynamics calculations showed that conjugate formation does not perturb the overall peptide structure. Only 6 exhibited antiproliferative activity in human tumor cells (IC50 = 63 ± 2 μM in MCF-7 cells and IC50 = 26 ± 3 μM in DU-145 cells) with active participation of somatostatin receptors in cellular uptake. Similar cytotoxic activity was found in a normal cell line (IC50 = 45 ± 2.6 μM in CHO cells), which can be attributed to a similar level of expression of somatostatin subtype-2 receptor. These studies provide new insights into the effect of receptor-binding peptide conjugation on the activity of metal-based anticancer drugs, and demonstrate the potential of such hybrid compounds to target tumor cells specifically.

The role of Cysteine 6.47 in class A GPCRs

Abstract Background: The CWxP motif of transmembrane helix 6 (x: any residue) is highly conserved in class A GPCRs. Within this motif, W6.48 is a big star in the theory of the global “toggle switch” because of its key role in the activation mechanism of GPCRs upon ligand binding. With all footlights focused on W6.48, the reason why the preceding residue, C6.47, is largely conserved is still unknown. The present study is aimed to fill up this lack of knowledge by characterizing the role of C6.47 of the CWxP motif. Results: A complete analysis of available crystal structures has been made alongside with molecular dynamics simulations of model peptides to explore a possible structural role for C6.47. Conclusions: We conclude that C6.47 does not modulate the conformation of the TM6 proline kink and propose that C6.47 participates in the rearrangement of the TM6 and TM7 interface accompanying activation.

Amplification of 16S rRNA gene sequences to differentiate two highly related bradyrhizobia species

A 16S rRNA gene PCR-based assay was developed aiming at a fast molecular diagnostic method to differentiate the two phylogenetically closely related species Bradyrhizobium japonicum and B. elkanii, isolated from soybean nodules, in order to identify those more competitive and comprising greater nitrogen fixation ability for use in the formulation of commercial inoculants. The assay used was able to discriminate ten reference strains belonging to these two Bradyrhizobium species, as well as to efficiently identify 37 strains isolated from fields cultivated with soybean.

Matching stages of heavy-ion collision models

Heavy-ion reactions and other collective dynamical processes are frequently described by different theoretical approaches for the different stages of the process, like initial equilibration stage, intermediate locally equilibrated fluid dynamical stage, and final freeze-out stage. For the last stage, the best known is the Cooper-Frye description used to generate the phase space distribution of emitted, noninteracting particles from a fluid dynamical expansion or explosion, assuming a final ideal gas distribution, or (less frequently) an out-of-equilibrium distribution. In this work we do not want to replace the Cooper-Frye description, but rather clarify the ways of using it and how to choose the parameters of the distribution and, eventually, how to choose the form of the phase space distribution used in the Cooper-Frye formula. Moreover, the Cooper-Frye formula is used in connection with the freeze-out problem, while the discussion of transition between different stages of the collision is applicable to other transitions also. More recently, hadronization and molecular dynamics models have been matched to the end of a fluid dynamical stage to describe hadronization and freeze-out. The stages of the model description can be matched to each other on space-time hypersurfaces (just like through the frequently used freeze-out hypersurface). This work presents a generalized description of how to match the stages of the description of a reaction to each other, extending the methodology used at freeze-out, in simple covariant form which is easily applicable in its simplest version for most applications.

Unraveling the "pressure effect" in nucleation

The influence of the pressure of a chemically inert carrier gas on the nucleation rate is one of the biggest puzzles in the research of gas-liquid nucleation. Experiments can show a positive effect, a negative effect, or no effect at all. The same experiment may show both trends for the same substance depending on temperature, or for different substances at the same temperature. We show how this ambiguous effect naturally arises from the competition of two contributions: nonisothermal effects and pressure-volume work. Our model clarifies seemingly contradictory experimental results and quantifies the variation of the nucleation ability of a substance in the presence of an ambient gas. Our findings are corroborated by molecular dynamics simulations and might have important implications since nucleation in experiments, technical applications, and nature practically always occurs in the presence of an ambient gas.

Influence of the pseudopotential on the properties of liquid rubidium at 315 K

The influence of the pseudopotential on both the structure and the self-diffusion of liquid rubidium at the melting point has been investigated by means of molecular-dynamics calculations. The model potential considered has been computed from the pseudopotential of Ashcroft, the dielectric function of Geldart and Vosko, and a Born-Mayer term. Four different values for the core radius which enters as input in the pseudopotential have been considered. In this way we have been able to observe and interpret the effect of this contribution on the properties of the liquid.

Investigation of classical epidemiological links between patients harbouring identical, non-predominant meticillin-resistant Staphylococcus aureus genotypes and lessons for epidemiological tracking.

According to molecular epidemiology theory, two isolates belong to the same chain of transmission if they are similar according to a highly discriminatory molecular typing method. This has been demonstrated in outbreaks, but is rarely studied in endemic situations. Person-to-person transmission cannot be established when isolates of meticillin-resistant Staphylococcus aureus (MRSA) belong to endemically predominant genotypes. By contrast, isolates of infrequent genotypes might be more suitable for epidemiological tracking. The objective of the present study was to determine, in newly identified patients harbouring non-predominant MRSA genotypes, whether putative epidemiological links inferred from molecular typing could replace classical epidemiology in the context of a regional surveillance programme. MRSA genotypes were defined using double-locus sequence typing (DLST) combining clfB and spa genes. A total of 1,268 non-repetitive MRSA isolates recovered between 2005 and 2006 in Western Switzerland were typed: 897 isolates (71%) belonged to four predominant genotypes, 231 (18%) to 55 non-predominant genotypes, and 140 (11%) were unique. Obvious epidemiological links were found in only 106/231 (46%) patients carrying isolates with non-predominant genotypes suggesting that molecular surveillance identified twice as many clusters as those that may have been suspected with classical epidemiological links. However, not all of these molecular clusters represented person-to-person transmission. Thus, molecular typing cannot replace classical epidemiology but is complementary. A prospective surveillance of MRSA genotypes could help to target epidemiological tracking in order to recognise new risk factors in hospital and community settings, or emergence of new epidemic clones.

0 AND 02 ADSORPTION ON Ag/Cu[100]

In this work the adsorption mechanisms of atomic and molecular oxygen on Cu(100) surface are studied using ab initio simulation methods. Through the atomistic scale under-standing of the elementary oxidation processes we can further understand the large-scale oxidation. Copper is a material widely used in industry which makes it an interesting subject, and also understanding the oxidation of copper helps us understand the oxidation mechanism of other metals. First we have a look on some theory on surface alloys in general and behaviour of Ag on Cu(100) surface. After that the physical background there is behind the methods of density functional calculations are discussed, and some methods, namely potential energy surfaces and molecular dynamics, are introduced. Then there is a brief look on the numerical details used in the calculations, and after that, the results of the simulations are exhibited.

Relation between sequence and structure of HIV-1 protease inhibitor complexes: a model system for the analysis of protein flexibility.

The flexibility of different regions of HIV-1 protease was examined by using a database consisting of 73 X-ray structures that differ in terms of sequence, ligands or both. The root-mean-square differences of the backbone for the set of structures were shown to have the same variation with residue number as those obtained from molecular dynamics simulations, normal mode analyses and X-ray B-factors. This supports the idea that observed structural changes provide a measure of the inherent flexibility of the protein, although specific interactions between the protease and the ligand play a secondary role. The results suggest that the potential energy surface of the HIV-1 protease is characterized by many local minima with small energetic differences, some of which are sampled by the different X-ray structures of the HIV-1 protease complexes. Interdomain correlated motions were calculated from the structural fluctuations and the results were also in agreement with molecular dynamics simulations and normal mode analyses. Implications of the results for the drug-resistance engendered by mutations are discussed briefly.