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.     

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.     

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.     

Calculating the free energy of association of transmembrane helices

A large number of experimental studies have been devoted to quantifying the interaction between transmembrane (TM) helices in detergent micelles and, more recently, in bilayers. Theoretical calculation of association free energy of TM helices would be useful for predicting the propensity of given sequences to oligomerize and for understanding the difference between association in micelles and in bilayers. In this article, the theoretical foundation for calculating the standard association free energy of TM helices is laid out and is applied to glycophorin A in both micelles and bilayers. The standard association free energy is decomposed into the effective energy, translational, rotational, and conformational entropy terms. The effective energy of association is obtained by molecular dynamics simulations in an implicit membrane model. The translational and rotational entropy of association is calculated from the probability distribution of the translational and rotational degrees of freedom obtained from the molecular dynamics simulations. The side-chain conformational entropy of association is estimated from the probability distribution obtained by rigid rotation of all side-chain dihedral angles. The calculated standard association free energy of glycophorin A in N-dodecylphosphocholine micelles is in good agreement with the experimental value. The translational entropy cost is larger, whereas the rotational entropy cost is smaller in bilayers than in micelles. The standard association free energy in 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayers is calculated to be approximately 1.3 kcal/mol more favorable than in N-dodecylphosphocholine micelles, consistent with available experimental data.     

Charge distribution and imperfect amphipathicity affect pore formation by antimicrobial peptides

Antimicrobial peptides often permeabilize biological membranes via a pore mechanism. Two pore types have been proposed: toroidal, where the pore is partly lined by lipid, and barrel-stave, where a cylindrical pore is completely lined by peptides. What drives the preference of antimicrobial peptides for a certain pore type is not yet fully understood. According to neutron scattering and oriented circular dichroism, melittin and MG-H2 induce toroidal pores whereas alamethicin forms barrel-stave pores. In previous work we found that indeed melittin seems to favor toroidal pores whereas alamethicin favors cylindrical pores. Here we designed mutants of these two peptides and the magainin analog MG-H2, aimed to probe how the distribution of charges along the helix and its imperfectly amphipathic structure influence pore formation. Molecular dynamics (MD) simulations of the peptides in a pre-formed cylindrical pore have been performed. The duration of the simulations was 136ns to 216ns. We found that a melittin mutant with lysine 7 neutralized favors cylindrical pores whereas a MG-H2 mutant with lysines in the N-terminal half of these peptides neutralized and an alamethicin mutant with a positive charge at the position 7 form semitoroidal pores. These results suggest that charged residues within the N-terminal half are important for toroidal pore formation. Toroidal pores produced by MG-H2 are more disordered than the melittin pores, likely because of the charged residues located in the middle of the MG-H2 helix (K11 and K14). Imperfect amphipathicity of melittin seems to play a role in its preference for toroidal pores since the substitutions of charged residues located within the nonpolar face by hydrophobic residues suppress evolution of a toroidal pore. The mutations change the position of lysine 7 near the N-terminus, relative to the lower leaflet headgroups. The MD simulations also show that the melittin P14A mutant forms a toroidal pore, but its configuration diverges from that of melittin and it is probably metastable.