Bioinformatics-driven identification and examination of candidate genes for non-alcoholic fatty liver disease

Objective

Candidate genes for non-alcoholic fatty liver disease (NAFLD) identified by a bioinformatics approach were examined for variant associations to quantitative traits of NAFLD-related phenotypes.

Research Design and Methods

By integrating public database text mining, trans-organism protein-protein interaction transferal, and information on liver protein expression a protein-protein interaction network was constructed and from this a smaller isolated interactome was identified. Five genes from this interactome were selected for genetic analysis. Twenty-one tag single-nucleotide polymorphisms (SNPs) which captured all common variation in these genes were genotyped in 10,196 Danes, and analyzed for association with NAFLD-related quantitative traits, type 2 diabetes (T2D), central obesity, and WHO-defined metabolic syndrome (MetS).

Results

273 genes were included in the protein-protein interaction analysis and EHHADH, ECHS1, HADHA, HADHB, and ACADL were selected for further examination. A total of 10 nominal statistical significant associations (P<0.05) to quantitative metabolic traits were identified. Also, the case-control study showed associations between variation in the five genes and T2D, central obesity, and MetS, respectively. Bonferroni adjustments for multiple testing negated all associations.

Conclusions

Using a bioinformatics approach we identified five candidate genes for NAFLD. However, we failed to provide evidence of associations with major effects between SNPs in these five genes and NAFLD-related quantitative traits, T2D, central obesity, and MetS.     

Validation of the ALS assay in adult patients with culture confirmed pulmonary tuberculosis

Background

We have earlier shown that Bacille Calmette-Guérin (BCG) vaccine-specific IgG Antibodies in Lymphocyte Supernatant (ALS) can be used for diagnosis of active tuberculosis (TB) in adults and children.

Methodology/Principal Findings

The ALS method was validated in a larger cohort (n = 212) of patients with suspicion of pulmonary TB using multiple antigens (BCG, LAM, TB15.3, TB51A, CFP10-ESAT6-A, CFP, CW) from Mycobacterium tuberculosis. The sensitivity and specificity of the ALS assay was calculated using non-TB patients as controls. The sensitivity and the specificity were highest with BCG vaccine (90% and 88% respectively) followed by LAM (89% and 87% respectively). Simultaneous assessment of multiple antigen-specific antibodies increased sensitivity (91%) and specificity (88%). Using higher lymphocyte count in smaller volume of culture media increased detection and reduced the assay duration to ∼30 hrs. Twenty one patients with clinical findings strongly suggestive of TB finally diagnosed as non-TB patients were positive by the ALS assay, of which 9 (43%) were positive for 7 antigens and 19 (90%) for at least 3 antigens.

Conclusions/Significance

Our findings show that simultaneous detection of antigens improves the diagnostic potential of the ALS assay; the modified method increases sensitivity and can provide results in <48 hours, and enable detection of some cases of pulmonary TB that are not detectable by standard methods.     

CCN2 promotes keratinocyte adhesion and migration via integrin α5β1

Background

CCN2, (a.k.a. connective tissue growth factor and CTGF) has emerged as a regulator of cell migration. While the importance of CCN2 for the fibrotic process in wound healing has been well studied, the effect of CCN2 on keratinocyte function is not well understood. In this study, we investigated the mechanism behind CCN2-driven keratinocyte adhesion and migration.Materials and methods: Adhesion assays were performed by coating wells with 10 μg/ml fibronectin (FN) or phosphate-buffered saline (PBS). Keratinocytes were seeded in the presence or absence of 200 ng/ml CCN2, 5 mmol/l ethylenediaminetetraacetic acid, 10 mmol/l cations, 500 μl arginine–glycine–aspartic acid (RGD), 500 μM arginine–glycine–glutamate–serine (RGES), and 10 μg/ml anti-integrin blocking antibodies. Migration studies were performed using a modified Boyden chamber assay. Quantitative PCR was used to study the effect of CCN2 on integrin subunit mRNA expression. To block intracellular pathways, keratinocytes were pretreated with 20 μM PD98059 (MEK-1 inhibitor) or 20 μM PF573228 (FAK inhibitor) for 60 min prior the addition of CCN2. Western blot was used to measure CCN2, p-ERK1/2, and ERK1/2.Results: CCN2 enhanced keratinocyte adhesion to fibronectin via integrin α5β1. The addition of anti-integrin α5β1 antibodies reduced CCN2-mediated keratinocyte migration. In addition, CCN2 regulated mRNA and protein expression of integrin subunits α5 and β1. CCN2 activated the FAK-MAPK signaling pathway, and pretreatment with MEK1-specific inhibitor PD98059 markedly reduced CCN2-induced keratinocyte migration.Conclusions: Our results demonstrate that CCN2 enhances keratinocyte adhesion and migration through integrin α5β1 and activation of the FAK-MAPK signaling cascade.     

Conformational States of Melittin at a Bilayer Interface

The distribution of peptide conformations in the membrane interface is central to partitioning energetics. Molecular-dynamics simulations enable characterization of in-membrane structural dynamics. Here, we describe melittin partitioning into dioleoylphosphatidylcholine lipids using CHARMM and OPLS force fields. Although the OPLS simulation failed to reproduce experimental results, the CHARMM simulation reported was consistent with experiments. The CHARMM simulation showed melittin to be represented by a narrow distribution of folding states in the membrane interface.Unstructured peptides fold into the membrane interface because partitioned hydrogen-bonded peptide bonds are energetically favorable compared to free peptide bonds (1–3). This folding process is central to the mechanisms of antimicrobial and cell-penetrating peptides, as well as to lipid interactions and stabilities of larger membrane proteins (4). The energetics of peptide partitioning into membrane interfaces can be described by a thermodynamic cycle (Fig. 1). State A is a theoretical state representing the fully unfolded peptide in water, B is the unfolded peptide in the membrane interface, C is the peptide in water, and D is the folded peptide in the membrane. The population of peptides in solution (State C) is best described as an ensemble of folded and unfolded conformations, whereas the population of peptides in State D generally is assumed to have a single, well-defined helicity, as shown in Fig. 1 A (5). Given that, in principle, folding in solution and in the membrane interface should follow the same basic rules, peptides in state D could reasonably be assumed to also be an ensemble. A fundamental question (5) is therefore whether peptides in state D can be correctly described as having a single helicity. Because differentiating an ensemble of conformations and a single conformation may be an impossible experimental task (5), molecular-dynamics (MD) simulations provide a unique high-resolution view of the phenomenon.Open in a separate windowFigure 1Thermodynamic cycles for peptide partitioning into a membrane interface. States A and B correspond to the fully unfolded peptide in solution and membrane interface, respectively. The folded peptide in solution is best described as an ensemble of unfolded and folded conformations (State C). State D is generally assumed to be one of peptides with a narrow range of conformations, but the state could actually be an ensemble of states as in the case of State C.Melittin is a 26-residue, amphipathic peptide that partitions strongly into membrane interfaces and therefore has become a model system for describing folding energetics (3,6–8). Here, we describe the structural dynamics of melittin in a dioleoylphosphatidylcholine (DOPC) bilayer by means of two extensive MD simulations using two different force fields.We extended a 12-ns equilibrated melittin-DOPC system (9) by 17 μs using the Anton specialized hardware (10) with the CHARMM22/36 protein/lipid force field and CMAP correction (11,12) (see Fig. S1 and Fig. S2 in the Supporting Material). To explore force-field effects, a similar system was simulated for 2 μs using the OPLS force field (13) (see Methods in the Supporting Material). In agreement with x-ray diffraction measurements on melittin in DOPC multilayers (14), melittin partitioned spontaneously into the lipid headgroups at a position below the phosphate groups at similar depth as glycerol/carbonyl groups (Fig. 2).Open in a separate windowFigure 2Melittin partitioned into the polar headgroup region of the lipid bilayer. (A) Snapshot of the simulation cell showing two melittin molecules (MLT1 and MLT2, in yellow) at the lipid-water interface. (B) Density cross-section of the simulation cell extracted from the 17-μs simulation. The peptides are typically located below the lipid phosphate (PO4) groups, in a similar depth as the glycerol/carbonyl (G/C) groups.To describe the secondary structure for each residue, we defined helicity by backbone dihedral angles (φ, ψ) within 30° from the ideal α-helical values (–57°, –47°). The per-residue helicity in the CHARMM simulation displays excellent agreement with amide exchange rates from NMR measurements that show a proline residue to separate two helical segments, which are unfolded below Ala⁵ and above Arg²² (15) (Fig. 3 A). In contrast, the OPLS simulation failed to reproduce the per-residue helicity except for a short central segment (see Fig. S3).Open in a separate windowFigure 3Helicity and conformational distribution of melittin as determined via MD simulation. (A) Helicity per residue for MLT1 and MLT2. (B) Corresponding evolution of the helicity. (C) Conformational distributions over the entire 17-μs simulation.Circular dichroism experiments typically report an average helicity of ∼70% for melittin at membrane interfaces (3,6,16,17), but other methods yield average helicities as high as 85% (15,18). Our CHARMM simulations are generally consistent with the experimental results, especially amide-exchange measurements (15); melittin helicity averaged to 78% for MLT1, whereas MLT2 transitioned from 75% to 89% helicity at t ≈ 8 μs, with an overall average helicity of 82% (Fig. 3 B). However, in the OPLS simulation, melittin steadily unfolds over the first 1.3 μs, after which the peptide remains only partly folded, with an average helicity of 33% (see Fig. S3). Similar force-field-related differences in peptide helicity were recently reported, albeit at shorter timescales (19). Although suitable NMR data are not presently available, we have computed NMR quadrupolar splittings for future reference (see Fig. S4).To answer the question asked in this article—whether the conformational space of folded melittin in the membrane interface can be described by a narrow distribution—the helicity distributions for the equilibrated trajectories are shown in Fig. 3 C. Whereas MLT1 in the CHARMM simulation produces a single, narrow distribution of the helicity, MLT2 has a bimodal distribution as a consequence of the folding event at t ≈ 8 μs (Fig. 3 C). We note that CHARMM force fields have a propensity for helix-formation and this transition might therefore be an artifact. We performed a cluster analysis to describe the structure of the peptide in the membrane interface. The four most populated conformations in the CHARMM simulation are shown in Fig. 4. The dominant conformation for both peptides was a helix kinked at G12 and unfolded at the last 5–6 residues of the C-terminus. The folding transition of MLT2 into a complete helix is visible by the 48% occupancy of a fully folded helix.Open in a separate windowFigure 4Conformational clusters of the two melittin peptides (MLT1 and MLT2) from the 17-μs CHARMM simulation in DOPC. Clustering is based on Cα-RMSD with a cutoff criterion of 2 Å.We conclude that the general assumption when calculating folding energetics holds: Folded melittin partitioned into membrane interfaces can be described by a narrow distribution of conformations. Furthermore, extended (several microsecond) simulations are needed to differentiate force-field effects. Although the CHARMM and OPLS simulations would seem to agree for the first few hundred nanoseconds, the structural conclusions differ drastically with longer trajectories, with CHARMM parameters being more consistent with experiments. However, as implied by the difference in substate distributions between MLT1 and MLT2, 17 μs might not be sufficient to observe the fully equilibrated partitioning process. The abrupt change in MLT2 might indicate that the helicity will increase to greater than experimentally observed in a sufficiently long simulation. On the other hand, it could be nothing more than a transient fluctuation. Increased sampling will provide further indicators of convergence of the helix partitioning process.     

Structure of the Homodimeric Glycine Decarboxylase P-protein from Synechocystis sp. PCC 6803 Suggests a Mechanism for Redox Regulation

Glycine decarboxylase, or P-protein, is a pyridoxal 5′-phosphate (PLP)-dependent enzyme in one-carbon metabolism of all organisms, in the glycine and serine catabolism of vertebrates, and in the photorespiratory pathway of oxygenic phototrophs. P-protein from the cyanobacterium Synechocystis sp. PCC 6803 is an α₂ homodimer with high homology to eukaryotic P-proteins. The crystal structure of the apoenzyme shows the C terminus locked in a closed conformation by a disulfide bond between Cys⁹⁷² in the C terminus and Cys³⁵³ located in the active site. The presence of the disulfide bridge isolates the active site from solvent and hinders the binding of PLP and glycine in the active site. Variants produced by substitution of Cys⁹⁷² and Cys³⁵³ by Ser using site-directed mutagenesis have distinctly lower specific activities, supporting the crucial role of these highly conserved redox-sensitive amino acid residues for P-protein activity. Reduction of the 353–972 disulfide releases the C terminus and allows access to the active site. PLP and the substrate glycine bind in the active site of this reduced enzyme and appear to cause further conformational changes involving a flexible surface loop. The observation of the disulfide bond that acts to stabilize the closed form suggests a molecular mechanism for the redox-dependent activation of glycine decarboxylase observed earlier.     

Antibodies Biotinylated Using a Synthetic Z-domain from Protein A Provide Stringent In Situ Protein Detection

Antibody-based protein profiling on a global scale using immunohistochemistry constitutes an emerging strategy for mapping of the human proteome, which is crucial for an increased understanding of biological processes in the cell. Immunohistochemistry is often performed indirectly using secondary antibodies for detection, with the benefit of signal amplification. Direct immunohistochemistry instead brings the advantage of multiplexing; however, it requires labeling of the primary antibody. Many antibody-labeling kits do not specifically target IgG and may therefore cause labeling of stabilizing proteins present in the antibody solution. A new conjugation method has been developed that utilizes a modified Z-domain of protein A (ZBPA) to specifically target the Fc part of antibodies. The aim of the present study was to compare the ZBPA conjugation method and a commercially available labeling kit, Lightning-Link, for in situ protein detection. Fourteen antibodies were biotinylated with each method and stained using immunohistochemistry. For all antibodies tested, ZBPA biotinylation resulted in distinct immunoreactivity without off-target staining, regardless of the presence of stabilizing proteins in the buffer, whereas the majority of the Lightning-Link biotinylated antibodies displayed a characteristic pattern of nonspecific staining. We conclude that biotinylated ZBPA domain provides a stringent method for antibody biotinylation, advantageous for in situ protein detection in tissues.     

Monoclonal antibodies selective for α‐synuclein oligomers/protofibrils recognize brain pathology in Lewy body disorders and α‐synuclein transgenic mice with the disease‐causing A30P mutation

Inclusions of intraneuronal alpha‐synuclein (α‐synuclein) can be detected in brains of patients with Parkinson's disease and dementia with Lewy bodies. The aggregation of α‐synuclein is a central feature of the disease pathogenesis. Among the different α‐synuclein species, large oligomers/protofibrils have particular neurotoxic properties and should therefore be suitable as both therapeutic and diagnostic targets. Two monoclonal antibodies, mAb38F and mAb38E2, with high affinity and strong selectivity for large α‐synuclein oligomers were generated. These antibodies, which do not bind amyloid‐beta or tau, recognize Lewy body pathology in brains from patients with Parkinson's disease and dementia with Lewy bodies and detect pathology earlier in α‐synuclein transgenic mice than linear epitope antibodies. An oligomer‐selective sandwich ELISA, based on mAb38F, was set up to analyze brain extracts of the transgenic mice. The overall levels of α‐synuclein oligomers/protofibrils were found to increase with age in these mice, although the levels displayed a large interindividual variation. Upon subcellular fractionation, higher levels of α‐synuclein oligomers/protofibrils could be detected in the endoplasmic reticulum around the age when behavioral disturbances develop. In summary, our novel oligomer‐selective α‐synuclein antibodies recognize relevant pathology and should be important tools to further explore the pathogenic mechanisms in Lewy body disorders. Moreover, they could be potential candidates both for immunotherapy and as reagents in an assay to assess a potential disease biomarker.