Differential peptide binding to CD40 evokes counteractive responses.

The antigen-presenting cell-expressed CD40 is implied in the regulation of counteractive immune responses such as induction of pro-inflammatory and anti-inflammatory cytokines interleukin (IL)-12 and IL-10, respectively. The mechanism of this duality in CD40 function remains unknown. Here, we investigated whether such duality depends on ligand binding. Based on CD40 binding, we identifed two dodecameric peptides, peptide-7 and peptide-19, from the phage peptide library. Peptide-7 induces IL-10 and increases Leishmania donovani infection in macrophages, whereas peptide-19 induces IL-12 and reduces L. donovani infection. CD40-peptide interaction analyses by surface plasmon resonance and atomic force microscopy suggest that the functional differences are not associated with the studied interaction parameters. The molecular dynamic simulation of the CD40-peptides interaction suggests that these two peptides bind to two different places on CD40. Thus, we suggest for the first time that differential binding of the ligands imparts functional duality to CD40.

Heat transfer between nanoparticles: Thermal conductance for near-field interactions

We analyze the heat transfer between two nanoparticles separated by a distance lying in the near-field domain in which energy interchange is due to the Coulomb interactions. The thermal conductance is computed by assuming that the particles have charge distributions characterized by fluctuating multipole moments in equilibrium with heat baths at two different temperatures. This quantity follows from the fluctuation-dissipation theorem for the fluctuations of the multipolar moments. We compare the behavior of the conductance as a function of the distance between the particles with the result obtained by means of molecular dynamics simulations. The formalism proposed enables us to provide a comprehensive explanation of the marked growth of the conductance when decreasing the distance between the nanoparticles.

Metastable liquid-liquid phase transition in a single-component system with only one crystal phase and no density anomaly

We investigate the phase behavior of a single-component system in three dimensions with spherically-symmetric, pairwise-additive, soft-core interactions with an attractive well at a long distance, a repulsive soft-core shoulder at an intermediate distance, and a hard-core repulsion at a short distance, similar to potentials used to describe liquid systems such as colloids, protein solutions, or liquid metals. We showed [Nature (London) 409, 692 (2001)] that, even with no evidence of the density anomaly, the phase diagram has two first-order fluid-fluid phase transitions, one ending in a gas¿low-density-liquid (LDL) critical point, and the other in a gas¿high-density-liquid (HDL) critical point, with a LDL-HDL phase transition at low temperatures. Here we use integral equation calculations to explore the three-parameter space of the soft-core potential and perform molecular dynamics simulations in the interesting region of parameters. For the equilibrium phase diagram, we analyze the structure of the crystal phase and find that, within the considered range of densities, the structure is independent of the density. Then, we analyze in detail the fluid metastable phases and, by explicit thermodynamic calculation in the supercooled phase, we show the absence of the density anomaly. We suggest that this absence is related to the presence of only one stable crystal structure.

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.

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.

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.

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.

Mekanismien dynamiikan simuloinnissa sovellettavia numeerisia- ja mallinnusmenetelmiä

Koneet voidaan usein jakaa osajärjestelmiin, joita ovat ohjaus- ja säätöjärjestelmät, voimaa tuottavat toimilaitteet ja voiman välittävät mekanismit. Eri osajärjestelmiä on simuloitu tietokoneavusteisesti jo usean vuosikymmenen ajan. Osajärjestelmien yhdistäminen on kuitenkin uudempi ilmiö. Usein esimerkiksi mekanismien mallinnuksessa toimilaitteen tuottama voimaon kuvattu vakiona, tai ajan funktiona muuttuvana voimana. Vastaavasti toimilaitteiden analysoinnissa mekanismin toimilaitteeseen välittämä kuormitus on kuvattu vakiovoimana, tai ajan funktiona työkiertoa kuvaavana kuormituksena. Kun osajärjestelmät on erotettu toisistaan, on niiden välistenvuorovaikutuksien tarkastelu erittäin epätarkkaa. Samoin osajärjestelmän vaikutuksen huomioiminen koko järjestelmän käyttäytymissä on hankalaa. Mekanismien dynamiikan mallinnukseen on kehitetty erityisesti tietokoneille soveltuvia numeerisia mallinnusmenetelmiä. Useimmat menetelmistä perustuvat Lagrangen menetelmään, joka mahdollistaa vapaasti valittaviin koordinaattimuuttujiin perustuvan mallinnuksen. Numeerista ratkaisun mahdollistamiseksi menetelmän avulla muodostettua differentiaali-algebraaliyhtälöryhmää joudutaan muokkaamaan esim. derivoimalla rajoiteyhtälöitä kahteen kertaan. Menetelmän alkuperäisessä numeerisissa ratkaisuissa kaikki mekanismia kuvaavat yleistetyt koordinaatit integroidaan jokaisella aika-askeleella. Tästä perusmenetelmästä johdetuissa menetelmissä riippumattomat yleistetyt koordinaatit joko integroidaan ja riippuvat koordinaatit ratkaistaan rajoiteyhtälöiden perusteella tai yhtälöryhmän kokoa pienennetään esim. käyttämällä nopeus- ja kiihtyvyysanalyyseissä eri kiertymäkoordinaatteja kuin asema-analyysissä. Useimmat integrointimenetelmät on alun perin tarkoitettu differentiaaliyhtälöiden (ODE) ratkaisuunjolloin yhtälöryhmään liitetyt niveliä kuvaavat algebraaliset rajoiteyhtälöt saattavat aiheuttaa ongelmia. Nivelrajoitteiden virheiden korjaus, stabilointi, on erittäin tärkeää mekanismien dynamiikan simuloinnin onnistumisen ja tulosten oikeellisuuden kannalta. Mallinnusmenetelmien johtamisessa käytetyn virtuaalisen työn periaatteen oletuksena nimittäin on, etteivät rajoitevoimat tee työtä, eli rajoitteiden vastaista siirtymää ei tapahdu. Varsinkaan monimutkaisten järjestelmien pidemmissä analyyseissä nivelrajoitteet eivät toteudu tarkasti. Tällöin järjestelmän energiatasapainoei toteudu ja järjestelmään muodostuu virtuaalista energiaa, joka rikkoo virtuaalisen työn periaatetta, Tästä syystä tulokset eivät enää pidäpaikkaansa. Tässä raportissa tarkastellaan erityyppisiä mallinnus- ja ratkaisumenetelmiä, ja vertaillaan niiden toimivuutta yksinkertaisten mekanismien numeerisessa ratkaisussa. Menetelmien toimivuutta tarkastellaan ratkaisun tehokkuuden, nivelrajoitteiden toteutumisen ja energiatasapainon säilymisen kannalta.

Mass transfer and particleinteractions in granular systems and suspension flows

The purpose of this study was to investigate some important features of granular flows and suspension flows by computational simulation methods. Granular materials have been considered as an independent state ofmatter because of their complex behaviors. They sometimes behave like a solid, sometimes like a fluid, and sometimes can contain both phases in equilibrium. The computer simulation of dense shear granular flows of monodisperse, spherical particles shows that the collisional model of contacts yields the coexistence of solid and fluid phases while the frictional model represents a uniform flow of fluid phase. However, a comparison between the stress signals from the simulations and experiments revealed that the collisional model would result a proper match with the experimental evidences. Although the effect of gravity is found to beimportant in sedimentation of solid part, the stick-slip behavior associated with the collisional model looks more similar to that of experiments. The mathematical formulations based on the kinetic theory have been derived for the moderatesolid volume fractions with the assumption of the homogeneity of flow. In orderto make some simulations which can provide such an ideal flow, the simulation of unbounded granular shear flows was performed. Therefore, the homogeneous flow properties could be achieved in the moderate solid volume fractions. A new algorithm, namely the nonequilibrium approach was introduced to show the features of self-diffusion in the granular flows. Using this algorithm a one way flow can beextracted from the entire flow, which not only provides a straightforward calculation of self-diffusion coefficient but also can qualitatively determine the deviation of self-diffusion from the linear law at some regions nearby the wall inbounded flows. Anyhow, the average lateral self-diffusion coefficient, which was calculated by the aforementioned method, showed a desirable agreement with thepredictions of kinetic theory formulation. In the continuation of computer simulation of shear granular flows, some numerical and theoretical investigations were carried out on mass transfer and particle interactions in particulate flows. In this context, the boundary element method and its combination with the spectral method using the special capabilities of wavelets have been introduced as theefficient numerical methods to solve the governing equations of mass transfer in particulate flows. A theoretical formulation of fluid dispersivity in suspension flows revealed that the fluid dispersivity depends upon the fluid properties and particle parameters as well as the fluid-particle and particle-particle interactions.