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.

Cluster glasses of ultrasoft particles

We present molecular dynamics (MD) simulations results for dense fluids of ultrasoft, fully penetrable particles. These are a binary mixture and a polydisperse system of particles interacting via the generalized exponential model, which is known to yield cluster crystal phases for the corresponding monodisperse systems. Because of the dispersity in the particle size, the systems investigated in this work do not crystallize and form disordered cluster phases. The clusteringtransition appears as a smooth crossover to a regime in which particles are mostly located in clusters, isolated particles being infrequent. The analysis of the internal cluster structure reveals microsegregation of the big and small particles, with a strong homo-coordination in the binary mixture. Upon further lowering the temperature below the clusteringtransition, the motion of the clusters" centers-of-mass slows down dramatically, giving way to a cluster glass transition. In the cluster glass, the diffusivities remain finite and display an activated temperature dependence, indicating that relaxation in the cluster glass occurs via particle hopping in a nearly arrested matrix of clusters. Finally we discuss the influence of the microscopic dynamics on the transport properties by comparing the MD results with Monte Carlo simulations.

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.

Analysis of Industrial Granular Flow Applications by Using Advanced Collision Models

Granular flow phenomena are frequently encountered in the design of process and industrial plants in the traditional fields of the chemical, nuclear and oil industries as well as in other activities such as food and materials handling. Multi-phase flow is one important branch of the granular flow. Granular materials have unusual kinds of behavior compared to normal materials, either solids or fluids. Although some of the characteristics are still not well-known yet, one thing is confirmed: the particle-particle interaction plays a key role in the dynamics of granular materials, especially for dense granular materials. At the beginning of this thesis, detailed illustration of developing two models for describing the interaction based on the results of finite-element simulation, dimension analysis and numerical simulation is presented. The first model is used to describing the normal collision of viscoelastic particles. Based on some existent models, more parameters are added to this model, which make the model predict the experimental results more accurately. The second model is used for oblique collision, which include the effects from tangential velocity, angular velocity and surface friction based on Coulomb's law. The theoretical predictions of this model are in agreement with those by finite-element simulation. I n the latter chapters of this thesis, the models are used to predict industrial granular flow and the agreement between the simulations and experiments also shows the validation of the new model. The first case presents the simulation of granular flow passing over a circular obstacle. The simulations successfully predict the existence of a parabolic steady layer and show how the characteristics of the particles, such as coefficients of restitution and surface friction affect the separation results. The second case is a spinning container filled with granular material. Employing the previous models, the simulation could also reproduce experimentally observed phenomena, such as a depression in the center of a high frequency rotation. The third application is about gas-solid mixed flow in a vertically vibrated device. Gas phase motion is added to coherence with the particle motion. The governing equations of the gas phase are solved by using the Large eddy simulation (LES) and particle motion is predicted by using the Lagrangian method. The simulation predicted some pattern formation reported by experiment.

Poisoning effect of s on pd surfaces

In this study we observe the poisoning effect of S to the adsorption and dissociation of 02 molecule on Pd surfaces. To perform this study we used Viennaab initio Simulation Package (VASP) and Spanish Initiative for Electronic structure with thousands of Atoms (SIESTA) ab initio softwares. To describe all Pd surfaces we selected the (100), and (211) surfaces, because we need very reactive and simple surfaces. Before studying the poison¬ing effect of S we had to study the dissociation of 02 on the surfaces. We discovered that on the (100) surface the hollow site is the most reactive site, but at room temperature the steric hinderace effect occurs very easily. If the molecule has enough vibrational energyit will dissociate. On the (211) surface the (100) micro facet's hollow site is the most reactive site and the molecule dissociates in the site without any barrier, and the molecule drifts from the terrace to this site. An S atom sticks on the Pd (100) surface in the hollow site. It affects the d-band density of states of the nearests Pd atoms; It moves the center of the d-band downin energy, when the bond between the Pd atom and the 0 atom is more antibonding. In the hollow site the S atom also blocks the dissociation site of the molecule. On the Pd(211) surface the energetically favourable site of the S atom is the(100) microfacet's hollow site. There it blocks the most reactive site, but its effect to the Pd atoms next to it is not significant.

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.