Specificity of Loxosceles α clade phospholipase D enzymes for choline-containing lipids: Role of a conserved aromatic cage.

Spider venom GDPD-like phospholipases D (SicTox) have been identified to be one of the major toxins in recluse spider venom. They are divided into two major clades: the α clade and the ß clade. Most α clade toxins present high activity against lipids with choline head groups such as sphingomyelin, while activities in ß clade toxins vary and include preference for substrates containing ethanolamine headgroups (Sicarius terrosus, St_ßIB1). A structural comparison of available structures of phospholipases D (PLDs) reveals a conserved aromatic cage in the α clade. To test the potential influence of the aromatic cage on membrane-lipid specificity we performed molecular dynamics (MD) simulations of the binding of several PLDs onto lipid bilayers containing choline headgroups; two SicTox from the α clade, Loxosceles intermedia αIA1 (Li_αIA) and Loxosceles laeta αIII1 (Ll_αIII1), and one from the ß clade, St_ßIB1. The simulation results reveal that the aromatic cage captures a choline-headgroup and suggest that the cage plays a major role in lipid specificity. We also simulated an engineered St_ßIB1, where we introduced the aromatic cage, and this led to binding with choline-containing lipids. Moreover, a multiple sequence alignment revealed the conservation of the aromatic cage among the α clade PLDs. Here, we confirmed that the i-face of α and ß clade PLDs is involved in their binding to choline and ethanolamine-containing bilayers, respectively. Furthermore, our results suggest a major role in choline lipid recognition of the aromatic cage of the α clade PLDs. The MD simulation results are supported by in vitro liposome binding assay experiments.

DOPA Homeostasis by Dopamine: A Control-Theoretic View.

Dopamine (DA) is an important signal mediator in the brain as well as in the periphery. The term "dopamine homeostasis" occasionally found in the literature refers to the fact that abnormal DA levels can be associated with a variety of neuropsychiatric disorders. An analysis of the negative feedback inhibition of tyrosine hydroxylase (TH) by DA indicates, with support from the experimental data, that the TH-DA negative feedback loop has developed to exhibit 3,4-dihydroxyphenylalanine (DOPA) homeostasis by using DA as a derepression regulator. DA levels generally decline when DOPA is removed, for example, by increased oxidative stress. Robust DOPA regulation by DA further implies that maximum vesicular DA levels are established, which appear necessary for a reliable translation of neural activity into a corresponding chemical transmitter signal. An uncontrolled continuous rise (windup) in DA occurs when Levodopa treatment exceeds a critical dose. Increased oxidative stress leads to the successive breakdown of DOPA homeostasis and to a corresponding reduction in DA levels. To keep DOPA regulation robust, the vesicular DA loading requires close to zero-order kinetics combined with a sufficiently high compensatory flux provided by TH. The protection of DOPA and DA due to a channeling complex is discussed.

Molecular determinants of the N-terminal acetyltransferase Naa60 anchoring to the Golgi membrane.

Nα-Acetyltransferase 60 (Naa60 or NatF) was recently identified as an unconventional N-terminal acetyltransferase (NAT) because it localizes to organelles, in particular the Golgi apparatus, and has a preference for acetylating N termini of the transmembrane proteins. This knowledge challenged the prevailing view of N-terminal acetylation as a co-translational ribosome-associated process and suggested a new mechanistic functioning for the enzymes responsible for this increasingly recognized protein modification. Crystallography studies on Naa60 were unable to resolve the C-terminal tail of Naa60, which is responsible for the organellar localization. Here, we combined modeling, in vitro assays, and cellular localization studies to investigate the secondary structure and membrane interacting capacity of Naa60. The results show that Naa60 is a peripheral membrane protein. Two amphipathic helices within the Naa60 C terminus bind the membrane directly in a parallel position relative to the lipid bilayer via hydrophobic and electrostatic interactions. A peptide corresponding to the C terminus was unstructured in solution and only folded into an α-helical conformation in the presence of liposomes. Computational modeling and cellular mutational analysis revealed the hydrophobic face of two α-helices to be critical for membranous localization. Furthermore, we found a strong and specific binding preference of Naa60 toward membranes containing the phosphatidylinositol PI(4)P, thus possibly explaining the primary residency of Naa60 at the PI(4)P-rich Golgi. In conclusion, we have defined the mode of cytosolic Naa60 anchoring to the Golgi apparatus, most likely occurring post-translationally and specifically facilitating post-translational N-terminal acetylation of many transmembrane proteins.

Phase transitions in coarse-grained lipid bilayers containing cholesterol by molecular dynamics simulations.

Coarse-grained simulations of model membranes containing mixtures of phospholipid and cholesterol molecules at different concentrations and temperatures have been performed. A random mixing without tendencies for segregation or formation of domains was observed on spatial scales corresponding to a few thousand lipids and timescales up to several microseconds. The gel-to-liquid crystalline phase transition is successively weakened with increasing amounts of cholesterol without disappearing completely even at a concentration of cholesterol as high as 60%. The phase transition temperature increases slightly depending on the cholesterol concentration. The gel phase system undergoes a transition with increasing amounts of cholesterol from a solid-ordered phase into a liquid-ordered one. In the solid phase, the amplitude of the oscillations in the radial distribution function decays algebraically with a prefactor that goes to zero at the solid-liquid transition.

Kinetics and mechanisms of catalyzed dual-E (antithetic) controllers.

Homeostasis plays a central role in our understanding how cells and organisms are able to oppose environmental disturbances and thereby maintain an internal stability. During the last two decades there has been an increased interest in using control engineering methods, especially integral control, in the analysis and design of homeostatic networks. Several reaction kinetic mechanisms have been discovered which lead to integral control. In two of them integral control is achieved, either by the removal of a single control species E by zero-order kinetics ("single-E controllers"), or by the removal of two control species by second-order kinetics ("antithetic or dual-E control"). In this paper we show results when the control species E1 and E2 in antithetic control are removed enzymatically by ping-pong or ternary-complex mechanisms. Our findings show that enzyme-catalyzed dual-E controllers can work in two control modes. In one mode, one of the two control species is active, but requires zero-order kinetics in its removal. In the other mode, both controller species are active and both are removed enzymatically. Conditions for the two control modes are put forward and biochemical examples with the structure of enzyme-catalyzed dual-E controllers are discussed.

An amplified derepression controller with multisite inhibition and positive feedback.

How organisms are able to maintain robust homeostasis has in recent years received increased attention by the use of combined control engineering and kinetic concepts, which led to the discovery of robust controller motifs. While these motifs employ kinetic conditions showing integral feedback and homeostasis for step-wise perturbations, the motifs' performance differ significantly when exposing them to time dependent perturbations. One type of controller motifs which are able to handle exponentially and even hyperbolically growing perturbations are based on derepression. In these controllers the compensatory reaction, which neutralizes the perturbation, is derepressed, i.e. its reaction rate is increased by the decrease of an inhibitor acting on the compensatory flux. While controllers in this category can deal well with different time-dependent perturbations they have the disadvantage that they break down once the concentration of the regulatory inhibitor becomes too low and the compensatory flux has gained its maximum value. We wondered whether it would be possible to bypass this restriction, while still keeping the advantages of derepression kinetics. In this paper we show how the inclusion of multisite inhibition and the presence of positive feedback loops lead to an amplified controller which is still based on derepression kinetics but without showing the breakdown due to low inhibitor concentrations. By searching for the amplified feedback motif in natural systems, we found it as a part of the plant circadian clock where it is highly interlocked with other feedback loops.

Undulation contributions to the area compressibility in lipid bilayer simulations.

It is here shown that there is a considerable system size-dependence in the area compressibility calculated from area fluctuations in lipid bilayers. This is caused by the contributions to the area fluctuations from undulations. This is also the case in experiments. At present, such a contribution, in most cases, is subtracted from the experimental values to obtain a true area compressibility. This should also be done with the simulation values. Here, this is done by extrapolating area compressibility versus system size, down to very small (zero) system size, where undulations no longer exist. The area compressibility moduli obtained from such simulations do not agree with experimental true area compressibility moduli as well as the uncorrected ones from contemporary or earlier simulations, but tend, instead, to be approximately 50% too large. As a byproduct, the bending modulus can be calculated from the slope of the compressibility modulus versus system-size. The values obtained in this way for the bending modulus are then in good agreement with experiment.

Interfacial Aromatics Mediating Cation-π Interactions with Choline-Containing Lipids Can Contribute as Much to Peripheral Protein Affinity for Membranes as Aromatics Inserted below the Phosphates.

Membrane-binding interfaces of peripheral proteins are restricted to a small part of their exposed surface, so the ability to engage in strong selective interactions with membrane lipids at various depths in the interface, both below and above the phosphates, is an advantage. Driven by their hydrophobicity, aromatic amino acids preferentially partition into membrane interfaces, often below the phosphates, yet enthalpically favorable interactions with the lipid headgroups, above the phosphate plane, are likely to further stabilize high interfacial positions. Using free-energy perturbation, we calculate the energetic cost of alanine substitution for 11 interfacial aromatic amino acids from 3 peripheral proteins. We show that the involvement in cation-π interactions with the headgroups (i) increases the ΔΔGtransfer as compared with insertion at the same depth without cation-π stabilization and (ii) can contribute at least as much as deeper insertion below the phosphates, highlighting the multiple roles of aromatics in peripheral membrane protein affinity.

Molecular characterization of the lipophorin receptor in the crustacean ectoparasite Lepeophtheirus salmonis.

The Salmon louse (Lepeophtheirus salmonis) is a marine ectoparasite of salmonid fish in the Northern Hemisphere and considered as a major challenge in aquaculture and a threat to wild populations of salmonids. Adult female lice produce a large number of lipid-rich eggs, however, the mechanism of maternal lipid transport into developing eggs during salmon louse reproduction has not been described. In the present study, a full-length L. salmonis lipophorin receptor (LsLpR) consisting of 16 exons was obtained by RACE and RT-PCR. The predicted ORF was 952 amino acids and structural analysis showed five functional domains that are similar to LpR of insects and decapods. Phylogenetic analysis placed the LsLpR together with LpRs from decapods and insects. Expression analysis revealed that the relative abundance of LsLpR transcripts was highest in the larvae and adult female lice. In adult females, the LsLpR transcripts and protein were found in the ovary and vitellogenic oocytes whereas, in larvae, the LsLpR transcripts were found in the neuronal somata of the brain and the intestine. Oil Red O stain results revealed that storage of neutral lipids was found in vitellogenic oocytes and ovaries of adult females, and in the yolk of larvae. Moreover, RNA interference (RNAi) was conducted to demonstrate the function of LsLpR in reproduction and lipid metabolism in L. salmonis. In larvae, the transcription of LsLpR was decreased by 44-54% while in an experiment LsLpR knockdown female lice produced 72% less offspring than control lice.

Quantum Corrections to Classical Molecular Dynamics Simulations of Water and Ice.

Classical simulations of simple water models reproduce many properties of the liquid and ice but overestimate the heat capacity by about 65% at ordinary temperatures and much more for low temperature ice. This is due to the fact that the atomic vibrations are quantum mechanical. The application of harmonic quantum corrections to the molecular motion results in good heat capacities for the liquid and for ice at low temperatures but a successively growing positive deviation from experimental results for ice above 200 K that reaches 15% just below melting. We suggest that this deviation is due to the lack of quantum corrections to the anharmonic motions. For the liquid, the anharmonicities are even larger but also softer and thus in less need of quantum correction. Therefore, harmonic quantum corrections to the classically calculated liquid heat capacities result in agreement with the experimental values. The classical model underestimates the heat of melting by 15%, while the application of quantum corrections produces fair agreement. On the other hand, the heat of vaporization is overestimated by 10% in the harmonically corrected classical model.