The fitness trade-off between growth and stress resistance determines the phenotypic landscape.

BACKGROUND: A central challenge in biology is to discover a principle that determines individual phenotypic differences within a species. The growth rate is particularly important for a unicellular organism, and the growth rate under a certain condition is negatively associated with that of another condition, termed fitness trade-off. Therefore, there should exist a common molecular mechanism that regulates multiple growth rates under various conditions, but most studies so far have focused on discovering those genes associated with growth rates under a specific condition. RESULTS: In this study, we found that there exists a recurrent gene expression signature whose expression levels are related to the fitness trade-off between growth preference and stress resistance across various yeast strains and multiple conditions. We further found that the genomic variation of stress-response, ribosomal, and cell cycle regulators are potential causal genes that determine the sensitivity between growth and survival. Intriguingly, we further observed that the same principle holds for human cells using anticancer drug sensitivities across multiple cancer cell lines. CONCLUSIONS: Together, we suggest that the fitness trade-off is an evolutionary trait that determines individual growth phenotype within a species. By using this trait, we can possibly overcome anticancer drug resistance in cancer cells.

Global stabilizing control of large-scale biomolecular regulatory networks.

MOTIVATION: Cellular behavior is determined by complex non-linear interactions between numerous intracellular molecules that are often represented by Boolean network models. To achieve a desired cellular behavior with minimal intervention, we need to identify optimal control targets that can drive heterogeneous cellular states to the desired phenotypic cellular state with minimal node intervention. Previous attempts to realize such global stabilization were based solely on either network structure information or simple linear dynamics. Other attempts based on non-linear dynamics are not scalable. RESULTS: Here, we investigate the underlying relationship between structurally identified control targets and optimal global stabilizing control targets based on non-linear dynamics. We discovered that optimal global stabilizing control targets can be identified by analyzing the dynamics between structurally identified control targets. Utilizing these findings, we developed a scalable global stabilizing control framework using both structural and dynamic information. Our framework narrows down the search space based on strongly connected components and feedback vertex sets then identifies global stabilizing control targets based on the canalization of Boolean network dynamics. We find that the proposed global stabilizing control is superior with respect to the number of control target nodes, scalability, and computational complexity. AVAILABILITY AND IMPLEMENTATION: We provide a GitHub repository that contains the DCGS framework written in Python as well as biological random Boolean network datasets (https://github.com/sugyun/DCGS). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dynamic embedding of effective harmonic normal mode vibrations in all-atomistic energy gap fluctuations: Case study of light harvesting 2 complex.

Environmental effects in excitation energy transfer have mostly been modeled by baths of harmonic oscillators, but to what extent such modeling provides a reliable description of actual interactions between molecular systems and environments remains an open issue. We address this issue by investigating fluctuations in the excitation energies of the light harvesting 2 complex using a realistic all-atomistic simulation of the potential energy surface. Our analyses reveal that molecular motions exhibit significant anharmonic features, even for underdamped intramolecular vibrations. In particular, we find that the anharmonicity contributes to the broadening of spectral densities and substantial overlaps between neighboring peaks, which complicates the meaning of mode frequencies constituting a bath model. Thus, we develop a strategy to construct a minimally underdamped harmonic bath that has a clear connection to all-atomistic dynamics by utilizing actual normal modes of molecules but optimizing their frequencies such that the resulting bath model can best reproduce the all-atomistic simulation results. By subtracting the underdamped contribution from the entire fluctuations, we also show that identifying a residual spectral density representing all other contributions with overdamped behavior is possible. We find that this can be fitted well with a well-established analytic form of a spectral density function or, alternatively, modeled as explicit time dependent fluctuations with muti-exponential or power law type correlation functions. We provide an assessment and the implications of these possibilities. The approach presented here can also serve as a general strategy to construct a simplified bath model that can effectively represent the underlying all-atomistic bath dynamics.

Electron-Binding Dynamics of the Dipole-Bound State: Correlation Effect on the Autodetachment Dynamics.

The nature of the electron-binding forces in the dipole-bound states (DBS) of anions is interrogated through experimental and theoretical means by investigating the autodetachment dynamics from DBS Feshbach resonances of ortho-, meta-, and para-bromophenoxide (BrPhO-). Though the charge-dipole electrostatic potential has been widely regarded to be mainly responsible for the electron binding in DBS, the effect of nonclassical electron correlation has been conceived to be quite significant in terms of its static and/or dynamic contributions toward the binding of the excess electron to the neutral core. State-specific real-time autodetachment dynamics observed by picosecond time-resolved photoelectron velocity-map imaging spectroscopy reveal that the autodetachment processes from the DBS Feshbach resonances of BrPhO- anions cannot indeed be rationalized by the conventional charge-dipole potential. Specifically, the autodetachment lifetime is drastically lengthened depending on differently positioned Br-substitution, and this rate change cannot be explained within the framework of Fermi's golden rule based on the charge-dipole assumption. High-level ab initio quantum chemical calculations with EOM-EA-CCSD, which intrinsically takes into account electron correlations, generate more reasonable predictions on the binding energies than density functional theory (DFT) calculations, and semiclassical quantum dynamics simulations based on the EOM-EA-CCSD data excellently predict the trend in the autodetachment rates. These findings illustrate that static and dynamic properties of the excess electron in the DBS are strongly influenced by correlation interactions among electrons in the nonvalence orbital of the dipole-bound electron and highly polarizable valence orbitals of the bromine atom, which, in turn, dictate the interesting chemical fate of exotic anion species.

Inhibition of 3-phosphoinositide-dependent protein kinase 1 (PDK1) can revert cellular senescence in human dermal fibroblasts.

Cellular senescence is defined as a stable, persistent arrest of cell proliferation. Here, we examine whether senescent cells can lose senescence hallmarks and reenter a reversible state of cell-cycle arrest (quiescence). We constructed a molecular regulatory network of cellular senescence based on previous experimental evidence. To infer the regulatory logic of the network, we performed phosphoprotein array experiments with normal human dermal fibroblasts and used the data to optimize the regulatory relationships between molecules with an evolutionary algorithm. From ensemble analysis of network models, we identified 3-phosphoinositide-dependent protein kinase 1 (PDK1) as a promising target for inhibitors to convert the senescent state to the quiescent state. We showed that inhibition of PDK1 in senescent human dermal fibroblasts eradicates senescence hallmarks and restores entry into the cell cycle by suppressing both nuclear factor κB and mTOR signaling, resulting in restored skin regeneration capacity. Our findings provide insight into a potential therapeutic strategy to treat age-related diseases associated with the accumulation of senescent cells.

Protein folding from heterogeneous unfolded state revealed by time-resolved X-ray solution scattering.

One of the most challenging tasks in biological science is to understand how a protein folds. In theoretical studies, the hypothesis adopting a funnel-like free-energy landscape has been recognized as a prominent scheme for explaining protein folding in views of both internal energy and conformational heterogeneity of a protein. Despite numerous experimental efforts, however, comprehensively studying protein folding with respect to its global conformational changes in conjunction with the heterogeneity has been elusive. Here we investigate the redox-coupled folding dynamics of equine heart cytochrome c (cyt-c) induced by external electron injection by using time-resolved X-ray solution scattering. A systematic kinetic analysis unveils a kinetic model for its folding with a stretched exponential behavior during the transition toward the folded state. With the aid of the ensemble optimization method combined with molecular dynamics simulations, we found that during the folding the heterogeneously populated ensemble of the unfolded state is converted to a narrowly populated ensemble of folded conformations. These observations obtained from the kinetic and the structural analyses of X-ray scattering data reveal that the folding dynamics of cyt-c accompanies many parallel pathways associated with the heterogeneously populated ensemble of unfolded conformations, resulting in the stretched exponential kinetics at room temperature. This finding provides direct evidence with a view to microscopic protein conformations that the cyt-c folding initiates from a highly heterogeneous unfolded state, passes through still diverse intermediate structures, and reaches structural homogeneity by arriving at the folded state.

Mechanism of Ir(ppy)3 Guest Exciton Formation with the Exciplex-Forming TCTA:TPBI Cohost within a Phosphorescent Organic Light-Emitting Diode Environment.

Cohosts based on hole transporting and electron transporting materials often act as exciplexes in the form of intermolecular charge transfer complexes. Indeed, exciplex-forming cohosts have been widely developed as the host materials for efficient phosphorescent organic light-emitting diodes (OLEDs). In host-guest systems of OLEDs, the guest can be excited by two competing mechanisms, namely, excitation energy transfer (EET) and charge transfer (CT). Experimentally, it has been reported that the EET mechanism is dominant and the excitons are primarily formed in the host first and then transferred to the guest in phosphorescent OLEDs based on exciplex-forming cohosts. With this, exciplex-forming cohosts are widely employed for avoiding the formation of trapped charge carriers in the phosphorescent guest. However, theoretical studies are still lacking toward elucidating the relative importance between EET and CT processes in exciting the guest molecules in such systems. Here, we obtain the kinetics of guest excitation processes in a few trimer model systems consisting of an exciplex-forming cohost pair and a phosphorescent guest. We adopt the Förster resonance energy transfer (FRET) rate constants for the electronic transitions between excited states toward solving kinetic master equations. The input parameters for calculating the FRET rate constants are obtained from density functional theory (DFT) and time-dependent DFT. The results show that while the EET mechanism is important, the CT mechanism may still play a significant role in guest excitations. In fact, the relative importance of CT over EET depends strongly on the location of the guest molecule relative to the cohost pair. This is understandable as both the coupling for EET and the interaction energy for CT are strongly influenced by the geometric constraints. Understanding the energy transfer pathways from the exciplex state of cohost to the emissive state of guest may provide insights for improving exciplex-forming materials adopted in OLEDs.

Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes.

Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.

Cotargeting BET proteins overcomes resistance arising from PI3K/mTOR blockade-induced protumorigenic senescence in colorectal cancer.

Therapeutics targeting the phosphatidylinositol 3-kinase/mammalian target of rapamycin (PI3K/mTOR) pathway initially produce potent antitumor effects, but resistance frequently occurs. Using a phosphoproteome analysis, we found that colorectal cancer (CRC) cells exhibit resistance against PI3K/mTOR inhibition through feedback activation of multiple receptor tyrosine kinases, and their downstream focal adhesion kinase, Src and extracellular signal-regulated kinases signaling. Unexpectedly, PI3K/mTOR blockade causes senescence, mediated by the activation of the stress kinase p38. The senescent cancer cells induce the secretion of various cytokines and this senescence-associated secretome increases migration and invasion capabilities of CRC cells. We found that cotargeting PI3K/mTOR and bromodomain and extra-terminal domain can suppress activation of many oncogenic kinases involved in resistance to the PI3K/mTOR inhibition, induce cell death in vitro and tumor regression in vivo, and further prolong the survival of xenograft models. Our findings provide a rationale for a novel therapeutic strategy to overcome resistance to the PI3K/mTOR inhibitors in CRC.

Cooperative activation of PI3K by Ras and Rho family small GTPases.

Phosphoinositide 3-kinases (PI3Ks) and Ras and Rho family small GTPases are key regulators of cell polarization, motility, and chemotaxis. They influence each other's activities by direct and indirect feedback processes that are only partially understood. Here, we show that 21 small GTPase homologs activate PI3K. Using a microscopy-based binding assay, we show that K-Ras, H-Ras, and five homologous Ras family small GTPases function upstream of PI3K by directly binding the PI3K catalytic subunit, p110. In contrast, several Rho family small GTPases activated PI3K by an indirect cooperative positive feedback that required a combination of Rac, CDC42, and RhoG small GTPase activities. Thus, a distributed network of Ras and Rho family small GTPases induces and reinforces PI3K activity, explaining past challenges to elucidate the specific relevance of different small GTPases in regulating PI3K and controlling cell polarization and chemotaxis.

Layer-specific serotonergic induction of long-term depression in the prefrontal cortex of rats.

Layer 2/3 pyramidal neurons (L2/3 PyNs) of the cortex extend their basal dendrites near the soma and as apical dendritic tufts in layer 1, which mainly receive feedforward and feedback inputs, respectively. It is suggested that neuromodulators such as serotonin and acetylcholine may regulate the information flow between brain structures depending on the brain state. However, little is known about the dendritic compartment-specific induction of synaptic transmission in single PyNs. Here, we studied layer-specific serotonergic and cholinergic induction of long-term synaptic plasticity in L2/3 PyNs of the agranular insular cortex, a lateral component of the orbitofrontal cortex. Using FM1-43 dye unloading, we verified that local electrical stimulation to layers 1 (L1) and 3 (L3) activated axon terminals mostly located in L1 and perisomatic area (L2/3). Independent and AMPA receptor-mediated excitatory postsynaptic potential was evoked by local electrical stimulation of either L1 or L3. Application of serotonin (5-HT, 10 µM) induced activity-dependent longterm depression (LTD) in L2/3 but not in L1 inputs. LTD induced by 5-HT was blocked by the 5-HT2 receptor antagonist ketanserin, an NMDA receptor antagonist and by intracellular Ca2+ chelation. The 5-HT2 receptor agonist α-me-5-HT mimicked the LTD induced by 5-HT. However, the application of carbachol induced muscarinic receptor- dependent LTD in both inputs. The differential layer-specific induction of LTD by neuromodulators might play an important role in information processing mechanism of the prefrontal cortex.

Open channel block of Kv1.4 potassium channels by aripiprazole.

Aripiprazole is a quinolinone derivative approved as an atypical antipsychotic drug for the treatment of schizophrenia and bipolar disorder. It acts as with partial agonist activities at the dopamine D2 receptors. Although it is known to be relatively safe for patients with cardiac ailments, less is known about the effect of aripiprazole on voltage-gated ion channels such as transient A-type K+ channels, which are important for the repolarization of cardiac and neuronal action potentials. Here, we investigated the effects of aripiprazole on Kv1.4 currents expressed in HEK293 cells using a whole-cell patch-clamp technique. Aripiprazole blocked Kv1.4 channels in a concentration-dependent manner with an IC50 value of 4.4 µM and a Hill coefficient of 2.5. Aripiprazole also accelerated the activation (time-to-peak) and inactivation kinetics. Aripiprazole induced a voltage-dependent (δ = 0.17) inhibition, which was use-dependent with successive pulses on Kv1.4 currents without altering the time course of recovery from inactivation. Dehydroaripiprazole, an active metabolite of aripiprazole, inhibited Kv1.4 with an IC50 value of 6.3 µM (p < 0.05 compared with aripiprazole) with a Hill coefficient of 2.0. Furthermore, aripiprazole inhibited Kv4.3 currents to a similar extent in a concentration-dependent manner with an IC50 value of 4.9 µM and a Hill coefficient of 2.3. Thus, our results indicate that aripiprazole blocked Kv1.4 by preferentially binding to the open state of the channels.

Efficiently Transplanting Potential Energy Interpolation Database between Two Systems: Bacteriochlorophyll Case with FMO and LH2 Complexes.

Generating a reliable potential energy surface (PES) is an important issue for studying the dynamics of any system with computational simulations. Interpolation mechanics/molecular mechanics (IM/MM) based on a PES interpolation scheme is a useful tool in that regard as it provides an accuracy of a quantum chemistry (QC) level while maintaining its computational cost comparable to conventional MM force fields. Despite this benefit, constructing the database for interpolation itself is still challenging and time-consuming. Here, we present a method with which we can construct the IM database of one system based on a preexisting data set for another related system. We adopt the case of constructing bacteriochlorophyll PESs for the light-harvesting 2 (LH2) complex by utilizing already available IM database for the BChls from the Fenna-Matthews-Olson (FMO) complex. In this method, the IM database from FMO is first transplanted to LH2 by considering BChl displacement vectors that take into account the geometry differences induced by the protein scaffolds. From this transplanted primitive database entries, a relatively small number of effective ones are selected by a survival process based on a genetic algorithm such that the IM energies evaluated at geometries in a conveniently collected prediction set can closely match with the reference QC energies. The selection process is expedited by using two different levels of basis sets for the QC calculations. To demonstrate the utility of the PES thus constructed, we carry out 1 ns of IM/MM dynamics simulations with the finally optimized BChl database for LH2. Indeed, the energy profiles of the snapshots are found to be closely matching with the reference QC calculation data, with only ∼0.07 eV of errors in the ground- and excited-state energies and ∼0.008 eV of errors in the transition energies. We also show that properly selecting data points is actually quite important for generating an IM PES toward performing molecular dynamics simulations.

Early detection of acral melanoma: A review of clinical, dermoscopic, histopathologic, and molecular characteristics.

Acral lentiginous melanoma is a distinct subtype of melanoma on acral skin. Patient presentation at later stages and delayed diagnosis by physicians contribute to a worse associated prognosis and survival rate. Despite our progress in understanding the key features of this disease, the diagnosis of early-stage acral melanoma is still challenging. It is essential to integrate clinical, dermoscopic, and histologic findings in the diagnosis of acral lentiginous melanoma. In addition, molecular studies can be helpful. In this review, we have summarized our current understanding of this disease entity from articles that were published between 1969 and 2018. We have outlined clinical and dermoscopic features as well as pathologic and molecular findings regarding acral melanoma and have presented an algorithm for diagnosis. Understanding and integrating these characteristics may assist clinicians in the early detection of acral melanomas.

Layer-specific cholinergic modulation of synaptic transmission in layer 2/3 pyramidal neurons of rat visual cortex.

It is known that top-down associative inputs terminate on distal apical dendrites in layer 1 while bottom-up sensory inputs terminate on perisomatic dendrites of layer 2/3 pyramidal neurons (L2/3 PyNs) in primary sensory cortex. Since studies on synaptic transmission in layer 1 are sparse, we investigated the basic properties and cholinergic modulation of synaptic transmission in layer 1 and compared them to those in perisomatic dendrites of L2/3 PyNs of rat primary visual cortex. Using extracellular stimulations of layer 1 and layer 4, we evoked excitatory postsynaptic current/potential in synapses in distal apical dendrites (L1-EPSC/L1-EPSP) and those in perisomatic dendrites (L4-EPSC/L4-EPSP), respectively. Kinetics of L1-EPSC was slower than that of L4-EPSC. L1-EPSC showed presynaptic depression while L4-EPSC was facilitating. In contrast, inhibitory postsynaptic currents showed similar paired-pulse ratio between layer 1 and layer 4 stimulations with depression only at 100 Hz. Cholinergic stimulation induced presynaptic depression by activating muscarinic receptors in excitatory and inhibitory synapses to similar extents in both inputs. However, nicotinic stimulation enhanced excitatory synaptic transmission by ~20% in L4-EPSC. Rectification index of AMPA receptors and AMPA/NMDA ratio were similar between synapses in distal apical and perisomatic dendrites. These results provide basic properties and cholinergic modulation of synaptic transmission between distal apical and perisomatic dendrites in L2/3 PyNs of the visual cortex, which might be important for controlling information processing balance depending on attentional state.

Systems biological approaches to the cardiac signaling network.

Recent systems biological studies of cardiac systems have greatly advanced our understanding of cardiac physiology with a particular focus on the excitation-contraction coupling. With these advancements, there is a growing interest in systems analysis of the cardiac signaling network because its dynamical property is closely associated with cardiac diseases. In this article, we review recent attempts at computational modeling of the cardiac signaling network and provide a system-level perspective on the analysis of the large-scale cardiac signaling network. We discuss why the systems biological approach is useful and what novel insights it can provide for the development of personalized therapeutic strategies for cardiac diseases in the post-genomic era.

Protein disulfide isomerase inhibition synergistically enhances the efficacy of sorafenib for hepatocellular carcinoma.

Sorafenib is the only approved targeted drug for hepatocellular carcinoma (HCC), but its effect on patients' survival gain is limited and varies over a wide range depending on pathogenetic conditions. Thus, enhancing the efficacy of sorafenib and finding a reliable predictive biomarker are crucial to achieve efficient control of HCCs. In this study, we utilized a systems approach by combining transcriptome analysis of the mRNA changes in HCC cell lines in response to sorafenib with network analysis to investigate the action and resistance mechanism of sorafenib. Gene list functional enrichment analysis and gene set enrichment analysis revealed that proteotoxic stress and apoptosis modules are activated in the presence of sorafenib. Further analysis of the endoplasmic reticulum stress network model, combined with in vitro experiments, showed that introducing an additional stress by treating the orally active protein disulfide isomerase (PDI) inhibitor (PACMA 31) can synergistically increase the efficacy of sorafenib in vitro and in vivo, which was confirmed using a mouse xenograft model. We also found that HCC patients with high PDI expression show resistance to sorafenib and poor clinical outcomes, compared to the low-PDI-expression group. CONCLUSION: These results suggest that PDI is a promising therapeutic target for enhancing the efficacy of sorafenib and can also be a biomarker for predicting sorafenib responsiveness. (Hepatology 2017;66:855-868).

Tannic acid inhibits EGFR/STAT1/3 and enhances p38/STAT1 signalling axis in breast cancer cells.

Tannic acid (TA), a naturally occurring polyphenol, is a potent anti-oxidant with anti-proliferative effects on multiple cancers. However, its ability to modulate gene-specific expression of tumour suppressor genes and oncogenes has not been assessed. This work investigates the mechanism of TA to regulate canonical and non-canonical STAT pathways to impose the gene-specific induction of G1-arrest and apoptosis. Regardless of the p53 status and membrane receptors, TA induced G1-arrest and apoptosis in breast cancer cells. Tannic acid distinctly modulated both canonical and non-canonical STAT pathways, each with a specific role in TA-induced anti-cancer effects. Tannic acid enhanced STAT1 ser727 phosphorylation via upstream serine kinase p38. This STAT1 ser727 phosphorylation enhanced the DNA-binding activity of STAT1 and in turn enhanced expression of p21Waf1/Cip1 . However, TA binds to EGF-R and inhibits the tyrosine phosphorylation of both STAT1 and STAT3. This inhibition leads to the inhibition of STAT3/BCL-2 DNA-binding activity. As a result, the expression and mitochondrial localization of BCl-2 are declined. This altered expression and localization of mitochondrial anti-pore factors resulted in the release of cytochrome c and the activation of intrinsic apoptosis cascade involving caspases. Taken together, our results suggest that TA modulates EGF-R/Jak2/STAT1/3 and P38/STAT1/p21Waf1/Cip1 pathways and induce G1-arrest and intrinsic apoptosis in breast carcinomas.

Momilactone B Inhibits Ketosis In Vitro by Regulating the ANGPTL3-LPL Pathway and Inhibiting HMGCS2.

Ketogenesis is the production of ketone bodies, which provide energy when the body lacks glucose. Under ketogenic conditions, the body switches from primarily carbohydrate to fat metabolism to maintain energy balance. However, accumulation of high levels of ketone bodies in the blood results in ketosis. Treating ketosis with natural substances is preferable, because they are unlikely to cause side-effects. Momilactone B is an active compound isolated from Korean rice. Based on previous studies, we hypothesized that momilactone B could inhibit ketosis. We constructed an in vitro ketosis model by glucose starvation. We used this model to test the anti-ketosis effects of momilactone B. A primary target for treating ketosis is angiopoietin-like-3 (ANGPTL3), which modulates lipoprotein metabolism by inhibiting lipoprotein lipase (LPL), a multifunctional enzyme that breaks down stored fat to produce triglycerides. We showed that momilactone B could regulate the ANGPTL3-LPL pathway. However, a strong anti-ketosis candidate drug should also inhibit ketogenesis. Ketogenesis can be suppressed by inhibiting the expression of 3-hydroxy-3-methylglutaryl-CoA synthase-2 (HMGCS2), a mitochondrial enzyme that converts acetyl-CoA to ketone bodies. We found that momilactone B suppressed the expression of HMGCS2 through the increased expression of STAT5b. We also elucidated the relationship of STAT5b to ANGPTL3 and LPL expression.