Роль липидного окружения в процессе димеризации трансмембранных сегментов гликофорина А
An efficient computational approach is developed to quantify the free energy of a spontaneous association of the α-helices of proteins in the membrane environment. The approach is based on the numerical decomposition of the free energy profiles of the transmembrane (TM) helices into components corresponding to protein-protein, protein-lipid, and protein-water interactions. The method was tested for the TM segments of human glycophorin A (GpA) and two mutant forms, Gly83Ala and Thr87Val. It was shown that lipids make a significant negative contribution to the free energy of dimerization, while amino acid residues forming the interface of the helix-helix contact may be unfavorable in terms of free energy. The detailed balance between different energy contributions is highly dependent on the amino acid sequence of the TM protein segment. The results show the dominant role of the environment in the interaction of membrane proteins that is changing our notion of the driving force behind the spontaneous association of TM α-helices. Adequate quantification of the contribution of the water-lipid environment thus becomes an extremely urgent task for a rational design of new molecules targeting bitopic membrane proteins, including receptor tyrosine kinases.
Plasmatic membranes contain high amount of membrane proteins. They perform vital functions of life, so any disruptions in their structure result in pathologies and diseases. Studies of these proteins with experimental methods are very complicated and expensive, as they require the membrane environment. Despite considerable progress achieved so far in methods of structure determination and property analysis, many computational methods are developing to predict the structural and dynamical parameters of proteins in membranes. Among the algorithms of modeling are the homology analysis, de novo structure prediction, molecular dynamics simulations and other. With growing computational capabilities, sophisticated techniques are developed taking into account more environmental factors. Combined approaches with different levels of approximation of intermolecular interactions are widely used. The major interest in studies of membrane proteins is focused on their transmembrane domains that are fundamental structural elements and are constituted by α-helices or helical bundles incorporated into lipid bilayer in most cases. Therefore, the fundamental problem of interaction of a pair of helices in membrane arises: the exact mechanism of this process is still not so clear. In place of the prevailing concept of dimerization motifs that states the importance of protein-protein contacts, a new model of the membrane as an adaptable lipid matrix is proposed. It states that biological membrane can adjust its properties around proteins and also modulates their activity. This mechanism of the mutual influence of two components is challenging modern computational methods of membrane model- ing because these systems are quite large and include many components to be treated accurately. Nowadays, investigations of the complex multi-component model systems become possible with modern methods of computational experiment.
The method ofWave Packet Molecular Dynamics Method (WPMD) is a promising replacement of the classical molecular dynamics for the simulations of many-electron systems including nonideal plasmas. In this contribution we report on a packet splitting technique where an electron is represented by multiple Gaussians, with mixing coefficients playing the role of additional dynamic variables. It provides larger flexibility and better accuracy than the original WPMD with a single Gaussian per electron. As a test case we consider ionization of hydrogen atom in a short laser pulse, where the split packets provide a basis for quantum branching.
The wave packet molecular dynamics (WPMD) method provides a variational approximation to the solution of the time-dependent Schr¨odinger equation. Its application in the field of high-temperature dense plasmas has yielded diverging electron width (spreading), which results in diminishing electron-nuclear interactions. Electron spreading has previously been ascribed to a shortcoming of the WPMD method and has been counteracted by various heuristic additions to the models used. We employ more accurate methods to determine if spreading continues to be predicted by them and how WPMD can be improved. A scattering process involving a single dynamic electron interacting with a periodic array of statically screened protons is used as a model problem for comparison. We compare the numerically exact split operator Fourier transform method, the Wigner trajectory method, and the time-dependent variational principle (TDVP). Within the framework of the TDVP, we use the standard variational form of WPMD, the single Gaussian wave packet (WP), as well as a sum of Gaussian WPs, as in the split WP method. Wave packet spreading is predicted by all methods, so it is not the source of the unphysical diminishing of electron-nuclear interactions in WPMD at high temperatures. Instead, the Gaussian WP’s inability to correctly reproduce breakup of the electron’s probability density into localized density near the protons is responsible for the deviation from more accurate predictions. Extensions of WPMD must include a mechanism for breakup to occur in order to yield dynamics that lead to accurate electron densities.
This paper describes the surface environment of the dense plasma arcs that damage rf accelerators, tokamaks, and other high gradient structures. We simulate the dense, nonideal plasma sheath near a metallic surface using molecular dynamics (MD) to evaluate sheaths in the non-Debye region for high density, low temperature plasmas. We use direct two-component MD simulations where the interactions between all electrons and ions are computed explicitly. We find that the non-Debye sheath can be extrapolated from the Debye sheath parameters with small corrections. We find that these parameters are roughly consistent with previous particle-in-cell code estimates, pointing to densities in the range 10^24–10^25 m^3. The high surface fields implied by these results could produce field emission that would short the sheath and cause an instability in the time evolution of the arc, and this mechanism could limit the maximum density and surface field in the arc. These results also provide a way of understanding how the properties of the arc depend on the properties (sublimation energy, for example) of the metal. Using these results, and equating surface tension and plasma pressure, it is possible to infer a range of plasma densities and sheath potentials from scanning electron microscope images of arc damage. We find that the high density plasma these results imply and the level of plasma pressure they would produce is consistent with arc damage on a scale 100 nm or less, in examples where the liquid metal would cool before this structure would be lost. We find that the submicron component of arc damage, the burn voltage, and fluctuations in the visible light production of arcs may be the most direct indicators of the parameters of the dense plasma arc, and the most useful diagnostics of the mechanisms limiting gradients in accelerators.
The dynamics of a two-component Davydov-Scott (DS) soliton with a small mismatch of the initial location or velocity of the high-frequency (HF) component was investigated within the framework of the Zakharov-type system of two coupled equations for the HF and low-frequency (LF) fields. In this system, the HF field is described by the linear Schrödinger equation with the potential generated by the LF component varying in time and space. The LF component in this system is described by the Korteweg-de Vries equation with a term of quadratic influence of the HF field on the LF field. The frequency of the DS soliton`s component oscillation was found analytically using the balance equation. The perturbed DS soliton was shown to be stable. The analytical results were confirmed by numerical simulations.
Radiation conditions are described for various space regions, radiation-induced effects in spacecraft materials and equipment components are considered and information on theoretical, computational, and experimental methods for studying radiation effects are presented. The peculiarities of radiation effects on nanostructures and some problems related to modeling and radiation testing of such structures are considered.
Let k be a field of characteristic zero, let G be a connected reductive algebraic group over k and let g be its Lie algebra. Let k(G), respectively, k(g), be the field of k- rational functions on G, respectively, g. The conjugation action of G on itself induces the adjoint action of G on g. We investigate the question whether or not the field extensions k(G)/k(G)^G and k(g)/k(g)^G are purely transcendental. We show that the answer is the same for k(G)/k(G)^G and k(g)/k(g)^G, and reduce the problem to the case where G is simple. For simple groups we show that the answer is positive if G is split of type A_n or C_n, and negative for groups of other types, except possibly G_2. A key ingredient in the proof of the negative result is a recent formula for the unramified Brauer group of a homogeneous space with connected stabilizers. As a byproduct of our investigation we give an affirmative answer to a question of Grothendieck about the existence of a rational section of the categorical quotient morphism for the conjugating action of G on itself.
This proceedings publication is a compilation of selected contributions from the “Third International Conference on the Dynamics of Information Systems” which took place at the University of Florida, Gainesville, February 16–18, 2011. The purpose of this conference was to bring together scientists and engineers from industry, government, and academia in order to exchange new discoveries and results in a broad range of topics relevant to the theory and practice of dynamics of information systems. Dynamics of Information Systems: Mathematical Foundation presents state-of-the art research and is intended for graduate students and researchers interested in some of the most recent discoveries in information theory and dynamical systems. Scientists in other disciplines may also benefit from the applications of new developments to their own area of study.
Let G be a connected semisimple algebraic group over an algebraically closed field k. In 1965 Steinberg proved that if G is simply connected, then in G there exists a closed irreducible cross-section of the set of closures of regular conjugacy classes. We prove that in arbitrary G such a cross-section exists if and only if the universal covering isogeny Ĝ → G is bijective; this answers Grothendieck's question cited in the epigraph. In particular, for char k = 0, the converse to Steinberg's theorem holds. The existence of a cross-section in G implies, at least for char k = 0, that the algebra k[G]G of class functions on G is generated by rk G elements. We describe, for arbitrary G, a minimal generating set of k[G]G and that of the representation ring of G and answer two Grothendieck's questions on constructing generating sets of k[G]G. We prove the existence of a rational (i.e., local) section of the quotient morphism for arbitrary G and the existence of a rational cross-section in G (for char k = 0, this has been proved earlier); this answers the other question cited in the epigraph. We also prove that the existence of a rational section is equivalent to the existence of a rational W-equivariant map T- - - >G/T where T is a maximal torus of G and W the Weyl group.