Selective robust optimization: A new intensity-modulated proton therapy optimization strategy.

PURPOSE: To develop a new robust optimization strategy for intensity-modulated proton therapy as an important step in translating robust proton treatment planning from research to clinical applications. METHODS: In selective robust optimization, a worst-case-based robust optimization algorithm is extended, and terms of the objective function are selectively computed from either the worst-case dose or the nominal dose. Two lung cancer cases and one head and neck cancer case were used to demonstrate the practical significance of the proposed robust planning strategy. The lung cancer cases had minimal tumor motion less than 5 mm, and, for the demonstration of the methodology, are assumed to be static. RESULTS: Selective robust optimization achieved robust clinical target volume (CTV) coverage and at the same time increased nominal planning target volume coverage to 95.8%, compared to the 84.6% coverage achieved with CTV-based robust optimization in one of the lung cases. In the other lung case, the maximum dose in selective robust optimization was lowered from a dose of 131.3% in the CTV-based robust optimization to 113.6%. Selective robust optimization provided robust CTV coverage in the head and neck case, and at the same time improved controls over isodose distribution so that clinical requirements may be readily met. CONCLUSIONS: Selective robust optimization may provide the flexibility and capability necessary for meeting various clinical requirements in addition to achieving the required plan robustness in practical proton treatment planning settings.

Association of lipidome remodeling in the adipocyte membrane with acquired obesity in humans.

Identification of early mechanisms that may lead from obesity towards complications such as metabolic syndrome is of great interest. Here we performed lipidomic analyses of adipose tissue in twin pairs discordant for obesity but still metabolically compensated. In parallel we studied more evolved states of obesity by investigating a separated set of individuals considered to be morbidly obese. Despite lower dietary polyunsaturated fatty acid intake, the obese twin individuals had increased proportions of palmitoleic and arachidonic acids in their adipose tissue, including increased levels of ethanolamine plasmalogens containing arachidonic acid. Information gathered from these experimental groups was used for molecular dynamics simulations of lipid bilayers combined with dependency network analysis of combined clinical, lipidomics, and gene expression data. The simulations suggested that the observed lipid remodeling maintains the biophysical properties of lipid membranes, at the price, however, of increasing their vulnerability to inflammation. Conversely, in morbidly obese subjects, the proportion of plasmalogens containing arachidonic acid in the adipose tissue was markedly decreased. We also show by in vitro Elovl6 knockdown that the lipid network regulating the observed remodeling may be amenable to genetic modulation. Together, our novel approach suggests a physiological mechanism by which adaptation of adipocyte membranes to adipose tissue expansion associates with positive energy balance, potentially leading to higher vulnerability to inflammation in acquired obesity. Further studies will be needed to determine the cause of this effect.

Role of lipids in spheroidal high density lipoproteins.

We study the structure and dynamics of spherical high density lipoprotein (HDL) particles through coarse-grained multi-microsecond molecular dynamics simulations. We simulate both a lipid droplet without the apolipoprotein A-I (apoA-I) and the full HDL particle including two apoA-I molecules surrounding the lipid compartment. The present models are the first ones among computational studies where the size and lipid composition of HDL are realistic, corresponding to human serum HDL. We focus on the role of lipids in HDL structure and dynamics. Particular attention is paid to the assembly of lipids and the influence of lipid-protein interactions on HDL properties. We find that the properties of lipids depend significantly on their location in the particle (core, intermediate region, surface). Unlike the hydrophobic core, the intermediate and surface regions are characterized by prominent conformational lipid order. Yet, not only the conformations but also the dynamics of lipids are found to be distinctly different in the different regions of HDL, highlighting the importance of dynamics in considering the functionalization of HDL. The structure of the lipid droplet close to the HDL-water interface is altered by the presence of apoA-Is, with most prominent changes being observed for cholesterol and polar lipids. For cholesterol, slow trafficking between the surface layer and the regimes underneath is observed. The lipid-protein interactions are strongest for cholesterol, in particular its interaction with hydrophobic residues of apoA-I. Our results reveal that not only hydrophobicity but also conformational entropy of the molecules are the driving forces in the formation of HDL structure. The results provide the first detailed structural model for HDL and its dynamics with and without apoA-I, and indicate how the interplay and competition between entropy and detailed interactions may be used in nanoparticle and drug design through self-assembly.

Membrane proteins diffuse as dynamic complexes with lipids.

We describe how membrane proteins diffuse laterally in the membrane plane together with the lipids surrounding them. We find a number of intriguing phenomena. The lateral displacements of the protein and the lipids are strongly correlated, as the protein and the neighboring lipids form a dynamical protein-lipid complex, consisting of approximately 50-100 lipids. The diffusion of the lipids in the complex is much slower compared to the rest of the lipids. We also find a strong directional correlation between the movements of the protein and the lipids in its vicinity. The results imply that in crowded membrane environments there are no "free" lipids, as they are all influenced by the protein structure and dynamics. Our results indicate that, in studies of cell membranes, protein and lipid dynamics have to be considered together.

Composition and lipid spatial distribution of HDL particles in subjects with low and high HDL-cholesterol.

A low level of high density lipoprotein cholesterol (HDL-C) is a powerful risk factor for cardiovascular disease. However, despite the reported key role of apolipo-proteins, specifically, apoA-I, in HDL metabolism, lipid molecular composition of HDL particles in subjects with high and low HDL-C levels is currently unknown. Here lipidomics was used to study HDL derived from well-characterized high and low HDL-C subjects. Low HDL-C subjects had elevated triacylglycerols and diminished lysophosphatidylcholines and sphingomyelins. Using information about the lipid composition of HDL particles in these two groups, we reconstituted HDL particles in silico by performing large-scale molecular dynamics simulations. In addition to confirming the measured change in particle size, we found that the changes in lipid composition also induced specific spatial distributions of lipids within the HDL particles, including a higher amount of triacylglycerols at the surface of HDL particles in low HDL-C subjects. Our findings have important implications for understanding HDL metabolism and function. For the first time we demonstrate the power of combining molecular profiling of lipoproteins with dynamic modeling of lipoprotein structure.

Concerted diffusion of lipids in raft-like membranes.

Currently, there is no comprehensive model for the dynamics of cellular membranes. The understanding of even the basic dynamic processes, such as lateral diffusion of lipids, is still quite limited. Recent studies of one-component membrane systems have shown that instead of single-particle motions, the lateral diffusion is driven by a more complex, concerted mechanism for lipid diffusion (E. Falck et al., J. Am. Chem. Soc., 2008, 130, 44-45), where a lipid and its neighbors move in unison in terms of loosely defined clusters. In this work, we extend the previous study by considering the concerted lipid diffusion phenomena in many-component raft-like membranes. This nature of diffusion phenomena emerge in all the cases we have considered, including both atom-scale simulations of lateral diffusion within rafts and coarse-grained MARTINI simulations of diffusion in membranes characterized by coexistence of raft and non-raft domains. The data allows us to identify characteristic time scales for the concerted lipid motions, which turn out to range from hundreds of nanoseconds to several microseconds. Further, we characterize typical length scales associated with the correlated lipid diffusion patterns and find them to be about 10 nm, or even larger if weak correlations are taken into account. Finally, the concerted nature of lipid motions is also found in dissipative particle dynamics simulations of lipid membranes, clarifying the role of hydrodynamics (local momentum conservation) in membrane diffusion phenomena.

Bioinformatics and computational methods for lipidomics.

Large amounts of lipidomics data are rapidly becoming available. However, there is a lack of tools capable of taking the full advantage of the wealth of new information. Lipid bioinformatics is thus an emerging need as well as challenge for lipid research. Lipid concentration changes in biological systems reflect regulation at multiple spatial and dynamic scales, e.g., biochemical reactions in the cells, intercellular lipid trafficking, changes in cell membrane composition, systemic lipid metabolism or lipid oxidation. In order to address the complexity of lipids and their regulation, four areas of bioinformatics need to be developed: (1) data processing and lipid identification, (2) statistical data analysis, (3) pathway analysis, and (4) lipid modeling in systems and biophysical contexts. In this paper we overview the current state of the lipid bioinformatics field as well as suggest few potential new areas of research.

Conformational changes and slow dynamics through microsecond polarized atomistic molecular simulation of an integral Kv1.2 ion channel.

Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins-not to mention numerous applications in drug design. Here, we present a full 1 micros atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120 degrees rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation ( approximately 35 degrees ) of the extracellular end of all S4 segments is present also in a reference 0.5 micros simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 3(10) helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4-lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5-1 micros). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations.

Atom-scale molecular interactions in lipid raft mixtures.

We review the relationship between molecular interactions and the properties of lipid environments. A specific focus is given on bilayers which contain sphingomyelin (SM) and sterols due to their essential role for the formation of lipid rafts. The discussion is based on recent atom-scale molecular dynamics simulations, complemented by extensive comparison to experimental data. The discussion is divided into four sections. The first part investigates the properties of one-component SM bilayers and compares them to bilayers with phosphatidylcholine (PC), the focus being on a detailed analysis of the hydrogen bonding network in the two bilayers. The second part deals with binary mixtures of sterols with either SM or PC. The results show how the membrane properties may vary substantially depending on the sterol and SM type available, the membrane order and interdigitation being just two of the many examples of this issue. The third part concentrates on the specificity of intermolecular interactions in three-component mixtures of SM, PC and cholesterol (CHOL) under conditions where the concentrations of SM and CHOL are dilute with respect to that of PC. The results show how SM and CHOL favor one another, thus acting as nucleation sites for the formation of highly ordered nanosized domains. Finally, the fourth part discusses the large-scale properties of raft-like membrane environments and compares them to the properties of non-raft membranes. The differences turn out to be substantial. As a particularly intriguing example of this, the lateral pressure profiles of raft-like and non-raft systems indicate that the lipid composition of membrane domains may have a major impact on membrane protein activation.

Structure of spheroidal HDL particles revealed by combined atomistic and coarse-grained simulations.

Spheroidal high-density lipoprotein (HDL) particles circulating in the blood are formed through an enzymatic process activated by apoA-I, leading to the esterification of cholesterol, which creates a hydrophobic core of cholesteryl ester molecules in the middle of the discoidal phospholipid bilayer. In this study, we investigated the conformation of apoA-I in model spheroidal HDL (ms-HDL) particles using both atomistic and coarse-grained molecular dynamics simulations, which are found to provide consistent results for all HDL properties we studied. The observed small contribution of cholesteryl oleate molecules to the solvent-accessible surface area of the entire ms-HDL particle indicates that palmitoyloleoylphosphatidylcholines and apoA-I molecules cover the hydrophobic core comprised of cholesteryl esters particularly well. The ms-HDL particles are found to form a prolate ellipsoidal shape, with sizes consistent with experimental results. Large rigid domains and low mobility of the protein are seen in all the simulations. Additionally, the average number of contacts of cholesteryl ester molecules with apoA-I residues indicates that cholesteryl esters interact with protein residues mainly through their cholesterol moiety. We propose that the interaction of annular cholesteryl oleate molecules contributes to apoA-I rigidity stabilizing and regulating the structure and function of the ms-HDL particle.

Assessing the nature of lipid raft membranes.

The paradigm of biological membranes has recently gone through a major update. Instead of being fluid and homogeneous, recent studies suggest that membranes are characterized by transient domains with varying fluidity. In particular, a number of experimental studies have revealed the existence of highly ordered lateral domains rich in sphingomyelin and cholesterol (CHOL). These domains, called functional lipid rafts, have been suggested to take part in a variety of dynamic cellular processes such as membrane trafficking, signal transduction, and regulation of the activity of membrane proteins. However, despite the proposed importance of these domains, their properties, and even the precise nature of the lipid phases, have remained open issues mainly because the associated short time and length scales have posed a major challenge to experiments. In this work, we employ extensive atom-scale simulations to elucidate the properties of ternary raft mixtures with CHOL, palmitoylsphingomyelin (PSM), and palmitoyloleoylphosphatidylcholine. We simulate two bilayers of 1,024 lipids for 100 ns in the liquid-ordered phase and one system of the same size in the liquid-disordered phase. The studies provide evidence that the presence of PSM and CHOL in raft-like membranes leads to strongly packed and rigid bilayers. We also find that the simulated raft bilayers are characterized by nanoscale lateral heterogeneity, though the slow lateral diffusion renders the interpretation of the observed lateral heterogeneity more difficult. The findings reveal aspects of the role of favored (specific) lipid-lipid interactions within rafts and clarify the prominent role of CHOL in altering the properties of the membrane locally in its neighborhood. Also, we show that the presence of PSM and CHOL in rafts leads to intriguing lateral pressure profiles that are distinctly different from corresponding profiles in nonraft-like membranes. The results propose that the functioning of certain classes of membrane proteins is regulated by changes in the lateral pressure profile, which can be altered by a change in lipid content.

Insight into the putative specific interactions between cholesterol, sphingomyelin, and palmitoyl-oleoyl phosphatidylcholine.

The effects of cholesterol (Chol) on phospholipid bilayers include ordering of the fatty acyl chains, condensing of the lipids in the bilayer plane, and promotion of the liquid-ordered phase. These effects depend on the type of phospholipids in the bilayer and are determined by the nature of the underlying molecular interactions. As for Chol, it has been shown to interact more favorably with sphingomyelin than with most phosphatidylcholines, which in given circumstances leads to formation of lateral domains. However, the exact origin and nature of Chol-phospholipid interactions have recently been subjects of speculation. We examine interactions between Chol, palmitoylsphingomyelin (PSM) and palmitoyl-oleoyl-phosphatidylcholine (POPC) in hydrated lipid bilayers by extensive atom-scale molecular dynamics simulations. We employ a tailored lipid configuration: Individual PSM and Chol monomers, as well as PSM-Chol dimers, are embedded in a POPC lipid bilayer in the liquid crystalline phase. Such a setup allows direct comparison of dimeric and monomeric PSMs and Chol, which ultimately shows how the small differences in PSM and POPC structure can lead to profoundly different interactions with Chol. Our analysis shows that direct hydrogen bonding between PSM and Chol does not provide an adequate explanation for their putative specific interaction. Rather, a combination of charge-pairing, hydrophobic, and van der Waals interactions leads to a lower tilt in PSM neighboring Chol than in Chol with only POPC neighbors. This implies improved Chol-induced ordering of PSM's chains over POPC's chains. These findings are discussed in the context of the hydrophobic mismatch concept suggested recently.

Tilt: major factor in sterols' ordering capability in membranes.

Using extensive atomistic simulations, we show that there is a single experimentally accessible parameter--the sterol tilt--that can be used to determine a sterol's capability to induce order, and thus to promote, e.g., formation of lipid rafts. The observations also facilitate designing new effective sterols.

Influence of chain length and unsaturation on sphingomyelin bilayers.

Sphingomyelins (SMs) are among the most common phospholipid components of plasma membranes, usually constituting a mixture of several molecular species with various fatty acyl chain moieties. In this work, we utilize atomistic molecular dynamics simulations to study the differences in structural and dynamical properties of bilayers comprised of the most common natural SM species. Keeping the sphingosine moiety unchanged, we vary the amide bonded acyl chain from 16 to 24 carbons in length and examine the effect of unsaturation by comparing lipids with saturated and monounsaturated chains. As for structural properties, we find a slight decrease in average area per lipid and a clear linear increase in bilayer thickness with increasing acyl chain length both in saturated and unsaturated systems. Increasing the acyl chain length is found to further the interdigitation across the bilayer center. This is related to the dynamics of SM molecules, as the lateral diffusion rates decrease slightly for an increasing acyl chain length. Interdigitation also plays a role in interleaflet friction, which is stronger for unsaturated chains. The effect of the cis double bond is most significant on the local order parameters and rotation rates of the chains, though unsaturation shows global effects on overall lipid packing and dynamics as well. Regarding hydrogen bonding or properties related to the lipid/water interface region, no significant effects were observed due to varying chain length or unsaturation. The significance of the findings presented is discussed.

Structure and dynamics of sphingomyelin bilayer: insight gained through systematic comparison to phosphatidylcholine.

Sphingomyelin, one of the main lipid components of biological membranes, is actively involved in various cellular processes such as protein trafficking and signal transduction. In particular, specific lateral domains enriched in sphingomyelin and cholesterol have been proposed to play an important functional role in biomembranes, although their precise characteristics have remained unclear. A thorough understanding of the functional role of membranes requires detailed knowledge of their individual lipid components. Here, we employ molecular dynamics simulations to conduct a systematic comparison of a palmitoylsphingomyelin (PSM, 16:0-SM) bilayer with a membrane that comprises dipalmitoylphosphatidylcholine (DPPC) above the main phase transition temperature. We clarify atomic-scale properties that are specific to sphingomyelin due to its sphingosine moiety, and further discuss their implications for SM-rich membranes. We find that PSM bilayers, and in particular the dynamics of PSM systems, are distinctly different from those of a DPPC bilayer. When compared with DPPC, the strong hydrogen bonding properties characteristic to PSM are observed to lead to considerable structural changes in the polar headgroup and interface regions. The strong ordering of PSM acyl chains and specific ordering effects in the vicinity of a PSM-water interface reflect this issue clearly. The sphingosine moiety and related hydrogen bonding further play a crucial role in the dynamics of PSM bilayers, as most dynamic properties, such as lateral and rotational diffusion, are strongly suppressed. This is most evident in the rotational motion characterized by spin-lattice relaxation times and the decay of hydrogen bond autocorrelation functions that are expected to be important in complexation of SM with other lipids in many-component bilayers. A thorough understanding of SM bilayers would greatly benefit from nuclear magnetic resonance experiments for acyl chain ordering and dynamics, allowing full comparison of these simulations to experiments.

Pacing-induced calcineurin activation controls cardiac Ca2+ signalling and gene expression.

Calcineurin, a Ca(2+)-calmodulin-dependent protein phosphatase (PP2B) is one of the links between Ca(2+) signals and regulation of gene transcription in cardiac muscle. We studied the Ca(2+) signal specificity of calcineurin activation experimentally and with modelling. In the rat atrial preparation, an increase in pacing frequency increased nuclear activity of the calcineurin-sensitive transcription factor, nuclear factor of activated T-cells (NFAT), 2-fold in a cyclosporin A (CsA)-sensitive manner. In line with this, modelling results predicted that the frequency of cardiac Ca(2+) transients encodes the stimulus for calcineurin activation. We further observed experimentally that calcineurin inhibition by CsA modulated Ca(2+) release in a Ca(2+)-dependent manner. CsA had no effect on [Ca(2+)](i) at a pacing frequency of 1 Hz but it significantly suppressed the amplitude of Ca(2+) transients, systolic [Ca(2+)](i) and time averaged [Ca(2+)](i) at 6 Hz. Calcineurin had a differential role in the expression of immediate-early genes B-type natriuretic peptide (BNP) and c-fos. CsA inhibited the pacing-induced BNP gene expression, whereas pacing alone had no effect on the expression of c-fos. However, in the presence of CsA, c-fos mRNA levels were significantly augmented by increased pacing frequency. These results show that frequency-dependent calcineurin activation has a specific role in [Ca(2+)](i) regulation and gene expression, constantly recruited by varying cardiac Ca(2+) signals.

Calmodulin kinase modulates Ca2+ release in mouse skeletal muscle.

Activation of the contractile machinery in skeletal muscle is initiated by the action-potential-induced release of Ca2+ from the sarcoplasmic reticulum (SR). Several proteins involved in SR Ca2+ release are affected by calmodulin kinase II (CaMKII)-induced phosphorylation in vitro, but the effect in the intact cell remains uncertain and is the focus of the present study. CaMKII inhibitory peptide or inactive control peptide was injected into single isolated fast-twitch fibres of mouse flexor digitorum brevis muscles, and the effect on free myoplasmic [Ca2+] ([Ca2+]i) and force during different patterns of stimulation was measured. Injection of the inactive control peptide had no effect on any of the parameters measured. Conversely, injection of CaMKII inhibitory peptide decreased tetanic [Ca2+]i by ~25 %, but had no significant effect on the rate of SR Ca2+ uptake or the force-[Ca2+]i relationship. Repeated tetanic stimulation resulted in increased tetanic [Ca2+]i, and this increase was smaller after CaMKII inhibition. In conclusion, CaMKII-induced phosphorylation facilitates SR Ca2+ release in the basal state and during repeated contractions, providing a positive feedback between [Ca2+]i and SR Ca2+ release.