Membrane Curvature Sensing by Amphipathic Helices: Insights from Implicit Membrane Modeling

Sensing and generation of lipid membrane curvature, mediated by the binding of specific proteins onto the membrane surface, play crucial roles in cell biology. A number of mechanisms have been proposed, but the molecular understanding of these processes is incomplete. All-atom molecular dynamics simulations have offered valuable insights but are extremely demanding computationally. Implicit membrane simulations could provide a viable alternative, but current models apply only to planar membranes. In this work, the implicit membrane model 1 is extended to spherical and tubular membranes. The geometric change from planar to curved shapes is straightforward but insufficient for capturing the full curvature effect, which includes changes in lipid packing. Here, these packing effects are taken into account via the lateral pressure profile. The extended implicit membrane model 1 is tested on the wild-types and mutants of the antimicrobial peptide magainin, the ALPS motif of arfgap1, α-synuclein, and an ENTH domain. In these systems, the model is in qualitative agreement with experiments. We confirm that favorable electrostatic interactions tend to weaken curvature sensitivity in the presence of strong hydrophobic interactions but may actually have a positive effect when those are weak. We also find that binding to vesicles is more favorable than binding to tubes of the same diameter and that the long helix of α-synuclein tends to orient along the axis of tubes, whereas shorter helices tend to orient perpendicular to it. Adoption of a specific orientation could provide a mechanism for coupling protein oligomerization to tubule formation.     

The amino acid sequences of two alpha chains of hemoglobins from Komodo dragon Varanus komodoensis and phylogenetic relationships of amniotes

To elucidate phylogenetic relationships among amniotes and the evolution of alpha globins, hemoglobins were analyzed from the Komodo dragon (Komodo monitor lizard) Varanus komodoensis, the world's largest extant lizard, inhabiting Komodo Islands, Indonesia. Four unique globin chains (alpha A, alpha D, beta B, and beta C) were isolated in an equal molar ratio by high performance liquid chromatography from the hemolysate. The amino acid sequences of two alpha chains were determined. The alpha D chain has a glutamine at E7 as does an alpha chain of a snake, Liophis miliaris, but the alpha A chain has a histidine at E7 like the majority of hemoglobins. Phylogenetic analyses of 19 globins including two alpha chains of Komodo dragon and ones from representative amniotes showed the following results: (1) The a chains of squamates (snakes and lizards), which have a glutamine at E7, are clustered with the embryonic alpha globin family, which typically includes the alpha D chain from birds; (2) birds form a sister group with other reptiles but not with mammals; (3) the genes for embryonic and adult types of alpha globins were possibly produced by duplication of the ancestral alpha gene before ancestral amniotes diverged, indicating that each of the present amniotes might carry descendants of the two types of alpha globin genes; (4) squamates first split off from the ancestor of other reptiles and birds.      

TGF-β-Elicited Induction of Tissue Inhibitor of Metalloproteinases (TIMP)-3 Expression in Fibroblasts Involves Complex Interplay between Smad3, p38α, and ERK1/2

Transforming growth factor-β (TGF-β) promotes extracellular matrix deposition by down-regulating the expression of matrix degrading proteinases and upregulating their inhibitors. Tissue inhibitor of metalloproteinases (TIMP)-3 is an ECM-associated specific inhibitor of matrix degrading metalloproteinases. Here, we have characterized the signaling pathways mediating TGF-β-induced expression of TIMP-3. Basal and TGF-β-induced TIMP-3 mRNA expression was abolished in Smad4-deficient mouse embryonic fibroblasts and restoring Smad4 expression rescued the response. Inhibition of Smad signaling by expression of Smad7 and dominant negative Smad3 completely abolished TGF-β-elicited expression of TIMP-3 in human fibroblasts, whereas overexpression of Smad3 enhanced it. Inhibition of extracellular signal-regulated kinase 1/2 (ERK1/2) activation with PD98059 and p38 mitogen-activated protein kinase activity by SB203580 resulted in suppression of TGF-β-induced TIMP-3 expression, indicating that ERK1/2 and p38 MAPK mediate the effect of TGF-β on TIMP-3 expression. Specific activation of p38α and ERK1/2 by constitutively active mutants of MKK3b or MEK1, respectively, and simultaneous co-expression of Smad3 resulted in induction of TIMP-3 expression in the absence of TGF-β indicating that Smad3 co-operates with p38 and ERK1/2 in the induction of TIMP-3 expression. These results demonstrate the complex interplay between Smad3, p38α, and ERK1/2 signaling in the regulation of TIMP-3 gene expression in fibroblasts, which may play a role in inflammation, tissue repair, and fibrosis.     

Fibroblast-specific expression of a kinase-deficient type II transforming growth factor beta (TGFbeta) receptor leads to paradoxical activation of TGFbeta signaling pathways with fibrosis in transgenic mice

To better understand the role of disrupted transforming growth factor beta (TGFbeta) signaling in fibrosis, we have selectively expressed a kinase-deficient human type II TGFbeta receptor (TbetaRIIDeltak) in fibroblasts of transgenic mice, using a lineage-specific expression cassette subcloned from the pro-alpha2(I) collagen gene. Surprisingly, despite previous studies that characterized TbetaRIIDeltak as a dominant negative inhibitor of TGFbeta signaling, adult mice expressing this construct demonstrated TGFbeta overactivity and developed dermal and pulmonary fibrosis. Compared with wild type cells, transgenic fibroblasts proliferated more rapidly, produced more extracellular matrix, and showed increased expression of key markers of TGFbeta activation, including plasminogen activator inhibitor-1, connective tissue growth factor, Smad3, Smad4, and Smad7. Smad2/3 phosphorylation was increased in transgenic fibroblasts. Overall, the gene expression profile of explanted transgenic fibroblasts using cDNA microarrays was very similar to that of littermate wild type cells treated with recombinant TGFbeta1. Despite basal up-regulation of TGFbeta signaling pathways, transgenic fibroblasts were relatively refractory to further stimulation with TGFbeta1. Thus, responsiveness of endogenous genes to TGFbeta was reduced, and TGFbeta-regulated promoter-reporter constructs transiently transfected into transgenic fibroblasts showed little activation by recombinant TGFbeta1. Responsiveness was partially restored by overexpression of wild type type II TGFbeta receptors. Activation of MAPK pathways by recombinant TGFbeta1 appeared to be less perturbed than Smad-dependent signaling. Our results show that expression of TbetaRIIDeltak selectively in fibroblasts leads to paradoxical ligand-dependent activation of downstream signaling pathways and causes skin and lung fibrosis. As well as confirming the potential for nonsignaling receptors to regulate TGFbeta activity, these findings support a direct role for perturbed TGFbeta signaling in fibrosis and provide a novel genetically determined animal model of fibrotic disease.     

Investigation of pathways for the low-pH conformational transition in influenza hemagglutinin

Targeted molecular dynamics simulations were used to study the conformational transition of influenza hemagglutinin (HA) from the native conformation to putative fusogenic or postfusion conformations populated at low pH. Three pathways for this conformational change were considered. Complete dissociation of the globular domains of HA was observed in one pathway, whereas smaller rearrangements were observed in the other two. The fusion peptides became exposed and moved toward the target membrane, although occasional movement toward the viral membrane was also observed. The effective energy profiles along the paths show multiple barriers. The final low-pH structures, which are consistent with available experimental data, are comparable in effective energy to native HA. As a control, the uncleaved precursor HA0 was also forced along the same pathway. In this case both the final energy and the energy barrier were much higher than in the cleaved protein. This study suggests that 1) as proposed, the native conformation is the global minimum energy conformation for the uncleaved precursor but a metastable state for cleaved HA; 2) the spring-loaded conformational change is energetically plausible in full-length HA; and 3) complete globular domain dissociation is not necessary for extension of the coiled coil and fusion peptide exposure, but the model with complete dissociation has lower energy.     

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.     

Immune Response Gene Expression in Colorectal Cancer Carries Distinct Prognostic Implications According to Tissue,Stage and Site: A Prospective Retrospective Translational Study in the Context of a Hellenic Cooperative Oncology Group Randomised Trial

Background

Although host immune response is an emerging prognostic factor for colorectal cancer, there is no consensus on the optimal methodology, surrogate markers or tissue for study.

Patients and Methods

Tumour blocks were prospectively collected from 344 patients with stage II/III colorectal cancer (CRC) treated with adjuvant chemotherapy. Whole section lymphocytic infiltration was studied along with mRNA expression of CD3Z, CD8, CD4, CXCL9, CXCL13, IGHM, FOXP3, SNAI2 and ESR1 by qRT-qPCR in tissue microarray (TMA) cores from the centre of tumour, invasive margin and adjacent normal mucosa.

Results

Lymphocytic infiltration, deficient MMR (10.9%), KRAS (40.7%) and BRAF (4.9%) mutations or single mRNA gene expression were not prognostic. Tumour ESR1 gene expression (Hazard Ratio [HR] for relapse 2.33, 95% CI 1.35-4.02; HR for death 1.74, 95% CI 1.02-2.97) and absence of necrosis (HR for relapse 1.71, 95% CI 1.05-2.71; HR for death 1.98, 95% CI 1.14-3.43) were adverse prognostic features. We used CD3Z and CD8 expression in order to devise the mRNA-based Immune Score (mIS) and proceeded to partitioning analysis in 267 patients, with age, stage, tumour site (Right vs Left CRC), KRAS mutation and tumour mIS as input factors. Only in patients with stage III right-sided colon cancer, a low immune response was associated with inferior disease-free survival (mIS-low, HR for relapse 2.28, 95% CI 1.05-8.02). No prognostic significance was seen for tumour mIS in any other stage or site of CRC, or for a similar mIS score derived from adjacent normal mucosa. Independent adverse prognostic significance was retained in multivariable analysis for absence of necrosis, tumour ESR1 expression in all patients and low tumour mIS in stage III right-sided CRC.

Conclusions

In localised CRC, mRNA-based CD3Z/CD8 profiling of tumour immune response may have stage, site and tissue-specific prognostic significance, along with ESR1 expression.

Trial Registration

ANZCTR.org.au ACTRN12610000509066     

Implicit membrane treatment of buried charged groups: Application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

Reduced transforming growth factor-beta signaling in cartilage of old mice: role in impaired repair capacity

Osteoarthritis (OA) is a common joint disease, mainly effecting the elderly population. The cause of OA seems to be an imbalance in catabolic and anabolic factors that develops with age. IL-1 is a catabolic factor known to induce cartilage damage, and transforming growth factor (TGF)-beta is an anabolic factor that can counteract many IL-1-induced effects. In old mice, we observed reduced responsiveness to TGF-beta-induced IL-1 counteraction. We investigated whether expression of TGF-beta and its signaling molecules altered with age. To mimic the TGF-beta deprived conditions in aged mice, we assessed the functional consequence of TGF-beta blocking. We isolated knee joints of mice aged 5 months or 2 years, half of which were exposed to IL-1 by intra-articular injection 24 h prior to knee joint isolation. Immunohistochemistry was performed, staining for TGF-beta1, -2 or -3, TGF-betaRI or -RII, Smad2, -3, -4, -6 and -7 and Smad-2P. The percentage of cells staining positive was determined in tibial cartilage. To mimic the lack of TGF-beta signaling in old mice, young mice were injected with IL-1 and after 2 days Ad-LAP (TGF-beta inhibitor) or a control virus were injected. Proteoglycan (PG) synthesis (³⁵S-sulfate incorporation) and PG content of the cartilage were determined. Our experiments revealed that TGF-beta2 and -3 expression decreased with age, as did the TGF-beta receptors. Although the number of cells positive for the Smad proteins was not altered, the number of cells expressing Smad2P strongly dropped in old mice. IL-1 did not alter the expression patterns. We mimicked the lack of TGF-beta signaling in old mice by TGF-beta inhibition with LAP. This resulted in a reduced level of PG synthesis and aggravation of PG depletion. The limited response of old mice to TGF-beta induced-IL-1 counteraction is not due to a diminished level of intracellular signaling molecules or an upregulation of intracellular inhibitors, but is likely due to an intrinsic absence of sufficient TGF-beta receptor expression. Blocking TGF-beta distorted the natural repair response after IL-1 injection. In conclusion, TGF-beta appears to play an important role in repair of cartilage and a lack of TGF-beta responsiveness in old mice might be at the root of OA development.     

The contribution of C alpha-H...O hydrogen bonds to membrane protein stability depends on the position of the amide

Structural analyses of membrane proteins reveal a large number of C(alpha)-H...O contacts between transmembrane helices, presumed to be hydrogen bonds. Recent experiments produced conflicting results for the contribution of such hydrogen bonds to membrane protein stability. An FTIR study estimated an energy of -0.88 kcal/mol for the G79-C(alpha)-H...I76-O hydrogen bond in glycophorin A, whereas a mutagenesis study showed that the A51-C(alpha)-H...T24-O(gamma) hydrogen bond does not stabilize bacteriorhodopsin. Here, we reconcile these results using molecular mechanics calculations and an implicit membrane model (IMM1). With explicit hydrogen atoms, the potential energy of the G79-C(alpha)-H...I76-O interaction in GpA ranges from -0.54 to -0.9 kcal/mol and its contribution to stability (effective energy) from -0.49 to -0.83 kcal/mol, depending on the structural model used. The average values of these quantities in GpA-like motifs are similar. In bR, the corresponding numbers for the A51-C(alpha)-H...T24-O(gamma) interaction are +0.15 and +0.32 kcal/mol. The difference results from the different arrangement of the interacting groups and specifically the position of the acceptor with respect to the C(alpha) and N atoms. This conclusion likely applies to soluble proteins as well.