Cell culture and transfection
HEK293T (ATCC, CRL-3216) and Piezo1-knocked out HEK293T (HEK-P1KO; Gift from Ardem Patapoutian) cell lines were maintained in Dulbecco’s modified eagle medium (Gibco, 11–995) supplemented with 10% foetal bovine serum (sh30919.03), and 1% penicillin-streptomycin (Gibco, 15140). Both cell lines were maintained and cultured in a humidified incubator at 37°C with 5% CO2.
Cells were seeded in 35 mm dishes and cultured until approximately 50% cell confluency was reached. At a 1:3 ratio of DNA to transfection reagent ratio, 1 µg of DNA was transfected with 3 µL of FuGENE HD transfection reagent (Promega, E2311) for 15 min. TTN3 ortholog constructs, murine TTN3 single mutants and murine Piezo1 DNA constructs were subcloned into pIRES2-acGFP vectors and were utilised for transfection procedures. After 24 h, the transfected cells were collected and plated into 35 mm dishes containing glass coverslips, which were used for whole-cell patch-clamp recording the following day.
All basic chemicals were purchased from Sigma Aldrich, unless specified otherwise.
Gene preparation and mutagenesis
Tentonin 3/TMEM150C (TTN3) ortholog genes (mouse (Genscript (OMu08247)), zebrafish (Genscript (Oda15731)), frog (Genscript (OXa04470)), chicken (Genscript (OGa35266)), turtle (Genscript (Och161608)), cat (Genscript (OFb05814)), and human (Genscript (OHu05487))) were subcloned into pIRES2-acGFP vectors 5,6. TTN3 target constructs were amplified with Dreamtaq Green Polymerase (Thermo Fisher Scientific) through PCR. DNA fragments separated in agarose gels were extracted and purified with a FavorPrep™ Gel/PCR Purification Kit (Favorgen). The fragments underwent digestion with restriction enzymes (Thermo Fisher Scientific) and were ligated with T4 ligase (Thermo Fisher Scientific) using different molar insert-to-vector ratios; 1:1, 3:1, and 5:1. Transformation of the ligated products was performed with DH5α-chemically competent E. Coli (Enzynomics).
The primers used for subcloning are as follows:
zebrafish TTN3: forward primer 5′-CGCCCTTCCGCGGGAGTTGGTCAGAC-3′, reverse primer 5′-GTCTGACCAACTCCCGCGGAAGGGCG-3′,
frog TTN3: forward primer 5′-TAATACGACTCACTATAGGG-3′, reverse primer 5′-TAGAAGGCACAGTCGAGG-3′,
chicken TTN3: forward primer 5′-CTCGAGGCCACCATGGACGGGAA-3′, reverse primer 5′-GTCGACCTACACCTGATCCGTCTGGTA-3′,
turtle TTN3: forward primer 5′-GCTAGCGCCACCATGGATGGGAAGAAAT-3′, reverse primer 5′-AGTATCAAACAGATCAGGTGTAGCTCGAG-3′,
cat TTN3: forward primer 5′-CTCGAGGCCACCATGGATGGGAAG-3′, and reverse primer 5′-GTCGACTTACACCTGGTCGGTTTGATA-3′,
human TTN3: forward primer 5`-TAGGTGAACTCGAGGCCACCATGGATGGGAAGA-3′, and reverse primer 5′-CGATAACCGCGGTACCttACACCTGGTCAG-3′
Mutants of mTTN3 were generated using either a Muta-Direct™ Site-directed Mutagenesis Kit (Intron, 15071) or EZchange™ Site-Directed Mutagenesis Kit (Enzynomics, EZ004S). Primers used for site-directed mutagenesis were specifically generated for each mutant type. The mutation products underwent transformation in DH5α-chemically competent E. coli cells and later cultured on Luria-Bertani (LB) agar plates. The TTN3 mutations in the constructs were confirmed through sequencing.
Whole cell current recording
Whole-cell currents of HEK293T or N2A cells were recorded through the voltage-clamp technique. After forming a gigaseal with a glass pipette, the membrane patch under the glass pipette was ruptured with gentle suction to make a whole cell. Tip resistance of the glass pipettes was 2–3 MΩ. After forming a whole cell, capacitive transients were cancelled. The holding potential was set at -60 mV. Whole cell currents were recorded with an Axopatch 200B amplifier (Molecular Devices). The amplifier output was filtered at 1 kHz and digitised with Digidata 1550B (Molecular Devices) at a sampling rate of 5 kHz and stored on a PC for later analysis. Whole cell currents were analysed with Clampfit 10.0 software (Molecular Devices). The inactivation curves were fitted to two exponentials. The intracellular (pipette) solution (mM) contained 130 CsCl, 2 MgCl2, 4 Mg-ATP, 0.4 Na-GTP, 25 D-Mannitol and 10 HEPES adjusted to pH 7.2 with CsOH. The extracellular (bath) solution (mM) contained 130 NaCl, 5 KCl, 1 CaCl2, 2 MgCl2, 10 D-Mannitol and 10 HEPES adjusted to pH 7.2 with NaOH. All solution osmolarities were further adjusted to 295 mOsm with D-Mannitol. For the pharmacological experiments, GdCl3 (Sigma), Yoda1 (Tocris Bioscience), or NMB-1 (Anygen, Korea) were added to the bath solution.
Mechanical stimulation
Whole cells of HEK293T cells transfected with various TTN3 orthologs or Piezo1 were physically stimulated with a glass pipette as previously described 5,47. Briefly, after a glass pipette was pulled, the tip was fire-polished to produce a blunt tip 2 ~ 3 µm in diameter. The glass probe was fixed to the end of a Nano-controller NC4 micromanipulator (MM3A-LS, Kleindiek Nanotechnik, Reutlingen, Germany). The movement of the glass probe was controlled with a joystick or via a computer program provided by Kleindiek. The nano-controller was a Piezo-driven small motor that had two modes of operation; fine and coarse movement mode. We calibrated the actual movements of the glass probe. In the coarse mode, the glass probe moved on average of 0.266 µm per step. Distances of the mechanical step stimuli were set at 24, 28 and 32 steps for 6.4, 7.4, and 8.5 µm displacement, respectively. We generally used the coarse 8 modes for mechanical stimulation, where the piezoelectric device vibrated at 2.7 kHz corresponding to the speed of 718 µm/s except for the speed-current relationship experiments. For the speed-current relationship experiments, the frequencies of the nano-controller motion were set at 0.3, 0.6, 1.2, 2.4, 4.8, and 9.6 kHz, which corresponded to the velocity being approximately 80, 160, 319, 638, 1277, and 2550 µm/s, respectively. The displacement of the probe was fixed at 8 µm for a 600 ms duration. Before physical indentation, the tip of the mechanical probe was placed on top of the cell surface at a 50° angle and withdrawn by 1 µm, which was regarded as an initial point for the mechanical step.
Expression and purification of hTTN3 in E. coli
The P9 sequence was attached to the N-terminus of hTTN3 (or drTTN3) to facilitate the expression of hTTN3 in the membrane fraction of E. coli 24. For the mass production of the P9-TTN3 fusion protein, BL21 (DE3) star-pRARE cells harbouring the P9* expression vector, were grown in 5 mL LB containing 0.05 mg/mL carbenicillin and 0.05 mg/mL chloramphenicol for 16 h at 37°C. One millilitre of the pre-culture was inoculated into 50 mL of YTN (Yeast extract 1%, Bactotryptone 2%, NaCl 2%) medium containing 0.05 mg/mL carbenicillin and 0.05 mg/mL chloramphenicol and cultured at 37°C until the OD600nm reached 1.0. One-third of the culture was inoculated into 150 mL of the YTN medium containing 0.05 mg/mL carbenicillin and 0.05 mg/mL chloramphenicol and cultured at 37°C until the OD600nm reached between 1.0 and 2.0. The culture was then inoculated into a fermenter (Marado-PDA, Korea) containing 2 L culture medium (K2HPO4 0.3%, KH2PO4 0.5%, Yeast extract 2%, glucose 2%, MgSO4·7H2O 0.06%) containing 0.05 mg/mL carbenicillin and 0.05 mg/mL chloramphenicol. The cells were grown at 37°C, pH 6.7–6.8 and PO2 of 30–40%. The expression of P9-TTN3 was induced by the addition of 1 mM IPTG for 3 h at 37°C when the OD600nm of the culture was 46. The pH of the culture in the fermenter was maintained at 6.8 by the addition of a small volume of the feeding medium (yeast extract 27.4%, (NH4)2SO4 0.15%, glucose 21.1%, MgSO4·7H2O 0.1%). The cells were harvested by centrifugation for 20 min at 3,000 × g and stored at -80°C until further use.
Approximately 30 g of E. coli cells (wet weight) were resuspended in buffer A (50 mM Tris-HCl, pH 8.5, 1 mM EDTA) containing 1 mM phenylmethylsulphonyl fluoride (Merck, Germany) and 100 µg protease inhibitor cocktail (genDEPOT), and lysed using a microfluidizer M-110P (Microfluidics). The lysate was centrifuged at 15,000 × g for 30 min, and the resulting supernatant was further centrifuged at 100,000 × g for 1 h. The membrane fraction in the pellet was solubilised in treatment with buffer A containing 0.5% sarkosyl (Sigma Aldrich) for 2 h at 4°C. After removing the insoluble materials by centrifugation at 30,000 × g for 30 min, the recombinant hTTN3 in the solubilised membrane fraction was purified using a Ni-NTA agarose (Qiagen, Germany) column. The fractions containing hTTN3 were mixed with 4 mg APG for 4 h at 4°C. The hTTN3 complexed with APG was further purified by size exclusion chromatography using a Superdex™ 200 Increase 10/300 GL column (GE Healthcare).
Preparation of proteoliposome
Each 40 mg/mL of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) stock solution was dissolved in chloroform at a 1:1 ratio. The mixture was dried on a glass vial surface spread out as wide as possible under gentle N2 gas blowing for 30 min. The dried lipid film was rehydrated over 12 h in 140 mM KCl (Junsei, Japan) and 10 mM HEPES buffer to a final concentration of 10 mg/mL. The liposome solution was sonicated for 10 min with a 40% amplitude using a probe tip sonicator (Q125, Qsonica). Impurities in the liposomal solution were removed by centrifugation at 10,000 × g for 1 min. Clear milky liposome solution was mixed with the solution containing purified hTTN3, drTTN3 or MscL proteins at a 2:1 ratio, and incubated for more than 12 h at 22oC. The final lipid-to-protein ratio was 1:20,000. To remove residual APG, the proteoliposome solution was filtered using 10 k molecular weight cut-off Amicon centrifugal filters (Merck Millipore) at 15,400 × g for 15 min. All prepared proteoliposomes were used within 12 h.
Formation of the lipid bilayer incorporated with hTTN3, drTTN3 or MscL
The lipid bilayer membrane was formed with 3% DPhPC dissolved in n-decane (MP Biomedicals) using the painting method. An aperture of 80–100 µm was created using a spark generator (EM-09A, DAEDALON) on a 10 µm thick polytetrafluoroethylene film. The polytetrafluoroethylene film was clamped between two 300-µL chambers covered with transparent glass on one side. Both chambers were then filled with physiological buffer (140 mM KCl, 10 mM HEPES) and the Ag/AgCl electrodes were placed in both chambers separately. Then 0.2 µL of the DPhPC solution was painted on the aperture using a gel loading pipette (Corning). The formation of the bilayer was confirmed by microscope and capacitance measurement. After the formation of the lipid bilayer, 10–20 µL of the prepared proteoliposome was added to the cis chamber. Purified hTTN3, drTTN3 or MscL proteins in the proteoliposome were successfully reconstituted into the spanned bilayer lipid membrane by fusion.
Single-channel currents of the lipid bilayers were recorded using a patch-clamp amplifier (Axopatch 200B; molecular device). The acquired signals were processed with Digidata 1440A (molecular device). Signal recording conditions included a 10 kHz software filter, 5 kHz lowpass filter and 50 kHz sampling rate with the appropriate voltage bias in all experiments.
Protein structure prediction
The monomeric structure of mTTN3 (UniProt ID: Q8C8S3) was predicted using the AlphaFold2 program (v2.1.2) 20. The standard algorithm that was used consisted of two transformer-based modules for the sequence-residue and residue-residue pairwise features. The sequence-residue graph was built using the multiple sequence alignments (MSA) embedding with UniProt Reference Clusters 90 (UniRef90), Big Fantastic Database, Uniclust30, and MGnify databases (updated on 27 Nov 2021). The residue-residue graph contains information regarding the relationship between the residues. The MSA-based representation updates the residue-residue pair information through an element-wise outer product and sum over the MSA sequence dimension in the Evoformer block. Five initial models were created using different random seeds—while two models were generated using template structures, the others were generated without templates 20. The quality of the predicted monomer structures was evaluated using the pLDDT metric 20 applied to the Cα atoms’ distance differences between the predicted and target structures 33 (Supplementary Fig. 4).
The monomer model (i.e., the one with the highest pLDDT value) was subsequently used as input into the multi-body interface implemented in the HADDOCK program (version 2.4) 34,35,48 to generate the tetrameric structure. The computations were performed using the standard protocol where the centre of mass and non-crystallographic symmetry restraints were applied. The multi-body protocol consists of two main steps. In the first step, a set of 10,000 tetrameric models were generated via rigid-body docking followed by a refinement iteration of the 400 best-scored models introducing semi-flexible conditions. In the second step, only the 200 best-scored models obtained in the semi-flexible docking iteration were refined in the presence of an explicit water model. In this step the molecular mechanics potential energy as the sum of the individual energy terms such as van der Waals, electrostatic, desolvation, and restraints violation energies were used to evaluate the optimal orientation/conformation of the tetrameric structures. The obtained tetrameric structures were grouped into four types (types 1–4) according to the conformational similarity and two most relevant energy terms; van der Waals and electrostatic energies (Fig. 4b,c).
The model of the Type 1 group was further assessed with the ProQM program 49. ProQM analysis showed that the highest accuracy for the model was the residues located at the six transmembrane helices (averaged ProQM score/residue: 0.695 ± 0.10) whereas the lowest accuracy was shown for the residues located at the C-terminus (0.282 ± 0.18). The score for the transmembrane α-helices in the tetrameric complex was significantly higher than those of the monomer (0.572 ± 0.10) indicating that the tetrameric configuration might be preferred over the monomeric.
MD simulations
The mTTN3 tetramer incorporated into the explicit bilayer model was simulated at the atomic level to obtain the molecular conformations of mTTN3 in the lipid membrane environment. The protonation states of the amino-acid residues of the mTTN3 tetramer were titrated at pH = 7.4 using PROPKA 3.1 50, and the hydrogen bond network of the tetrameric structure was optimised via PDB2PQR 2.1 51. The protein was then embedded into the 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) bilayer using CHARMM-GUI 52. Each leaflet of the bilayer contained 350 DOPC. The system was fully hydrated with 130 TIP3P water per lipid 53. Neutralising cations (four K+ ions) for the negative charge of the protein and additional KCl salt were added to achieve a 150 mM KCl bulk concentration. A molecular snapshot of the MD simulation system is shown in the Supplementary Fig. 5a. Each system was simulated for three replicas (R1, R2, R3) with different random seeds for the initial velocities.
Simulations were conducted using the OpenMM (version 7.6.0) 54 program. The coordinates of all the atoms in the MD system were adjusted to minimise the energy and heated to 293.15 K over 125 ps followed by system equilibration for 1.75 ns. The production simulation was run using a Langevin thermostat at 298.15 K with a collision frequency of 1/ps. A Monte Carlo barostat 55 specifically designed for membrane simulations was used. The pressure, normal to the membrane plane, was 1 atm, corresponding to the normal laboratory conditions, and the surface tension applied parallel to the plane at the water-lipid headgroup interface was set to 5 dyn/cm. This surface tension corresponded to the negative lateral pressure that is capable of activating mechanosensitive ion channels, such as MscL, by laterally stretching the bilayer 56–58. The van der Waals term used a standard 6–12 LJ form with a force-switching function between 8 and 12 Å. Particle mesh Ewald was used for long-range electrostatic measurements beyond 12 Å. The integration time step equals 2 fs. Coordinate sets were saved for every 5 ps. The CHARMM C36 lipid 59 and protein 60–62 force fields were used. Three replicas were simulated with different random seeds for the initial velocities to ensure the convergence of the simulations.
Statistical analysis
All data presented in the figures represent the mean ± standard error of the mean. Student’s t-test was used to compare two means. One-way analysis of variance (ANOVA) followed by Tukey test was used to compare multiple means. P < 0.05 was considered as significant.