Raef Shams1,2, Yoshihiro Ito1,3, Hideyuki Miyatake2,3. 1. Emergent Bioengineering Materials Research Team, RIKEN Center for Emergent Matter Science, Wako, Saitama 351-0198, Japan. 2. Department of Life Science, Graduate School of Science and Engineering, Saitama University, Saitama City, Saitama 338-8570, Japan. 3. Nano Medical Engineering Laboratory, RIKEN Cluster for Pioneering Research, Wako, Saitama 351-0198, Japan.
Abstract
In cancer, the mechanistic/mammalian target of rapamycin complex-1 (mTORC1) is hyperactivated to promote survival under adverse conditions. The kinase activity of mTORC1 is activated by small-GTPase RHEB-GTP. Therefore, a new modality to inhibit mTORC1 activity has emerged, through intercepting RHEB. However, due to the relatively large contact area involved in the interaction between RHEB and mTORC1, facilitating this inhibition through small molecules has been challenging. Here, we report the development of a peptide that can inhibit the RHEB-mTORC1 interaction. The peptide, P1_WT, was designed based on the α-helix (aa 101-115) of the N-heat domain of mTOR to interact with switch II of RHEB. P1_WT bound to RHEB (K D = 0.14 μM) and inhibited RHEB-mTORN-heat interaction (IC50 = 0.33 μM) in vitro. Consequently, P1_WT inhibited mTORC1 activity at a sub-micromolar level (IC50 ∼ 0.3 μM). P1_WT was predicted to be cell-permeable due to the rich content of arginine (23%), enhancing the intracellular translocation. These results show that P1_WT is a potential compound to further develop inhibitors for mTORC1 by intercepting RHEB from mTORC1.
In cancer, the mechanistic/mammalian target of rapamycin complex-1 (mTORC1) is hyperactivated to promote survival under adverse conditions. The kinase activity of mTORC1 is activated by small-GTPase RHEB-GTP. Therefore, a new modality to inhibit mTORC1 activity has emerged, through intercepting RHEB. However, due to the relatively large contact area involved in the interaction between RHEB and mTORC1, facilitating this inhibition through small molecules has been challenging. Here, we report the development of a peptide that can inhibit the RHEB-mTORC1 interaction. The peptide, P1_WT, was designed based on the α-helix (aa 101-115) of the N-heat domain of mTOR to interact with switch II of RHEB. P1_WT bound to RHEB (K D = 0.14 μM) and inhibited RHEB-mTORN-heat interaction (IC50 = 0.33 μM) in vitro. Consequently, P1_WT inhibited mTORC1 activity at a sub-micromolar level (IC50 ∼ 0.3 μM). P1_WT was predicted to be cell-permeable due to the rich content of arginine (23%), enhancing the intracellular translocation. These results show that P1_WT is a potential compound to further develop inhibitors for mTORC1 by intercepting RHEB from mTORC1.
The mechanistic/mammalian
target of rapamycin complex 1 (mTORC1)
is a serine/threonine protein kinase to regulate cell growth and proliferation.[1,2] mTORC1 consists of multiple components: the mTOR kinase unit, mammalian
lethal with SEC13 protein 8 (mLST8), regulatory-associated protein
of mTOR (RAPTOR), DEP domain-containing mTOR-interacting protein (DEPTOR),
and the 40-kDa proline-rich AKT substrate (PRAS40).[2,3] These
components work together to recruit and phosphorylate substrate proteins
downstream the signal pathways such as the eukaryotic translation
initiation factor 4E-binding protein (4E-BP1) and ribosomal protein
S6 kinase 1 (S6K1).[3,4] Conversely, upstream the signaling
pathway, the kinase activity of mTORC1 is modulated by growth factors,
nutrients, and energy transfer to manage protein, lipid, and nucleotide
synthesis.[5,6] These upstream signals are shown to control
the activation process by coupling mTORC1 with the Ras homolog enriched
in brain (RHEB) protein on the lysosomal surface.[3,7]RHEB-GTPase is a small G-protein (∼20 kDa) that is farnesylated
to the lysosome membrane.[8] The GTPase activity
of RHEB is stimulated by the tuberous sclerosis complex 1/2 (TSC1/2),
the upstream negative regulator of mTORC1.[1,8] Thus,
RHEB activates mTORC1 only when it is charged with GTP but not GDP.[3] Recently, structural studies have elucidated
the activation mechanism of mTORC1 by the RHEB-GTP complex to allosterically
bind to the constituted site formed by the N-heat, M-heat, and FAT
domains of mTOR. Upon binding, conformational changes over a wide
area of the kinase domain of mTOR occur to receive ATP in the active
site for catalysis.[3] On the lysosome surface,
mTORC1 activation is cooperatively modulated via two in-parallel pathways;
the first is the activation of RHEB by the TSC1/2 complex to enhance
RHEB charging with GTP.[7] This pathway is
controlled by growth factors, including the insulin-like growth factor
(IGF), through IGF/PI3K/AKT pathways, which negatively or positively
control TSC complexation.[1,5,9] The second pathway involves the translocation of mTORC1 onto the
lysosome surface to form a complex with the RHEB-GTP. This translocation
is regulated via the capturing process of the mTORC1 subunit, RAPTOR,
by the Rags/Ragulator complex on the lysosome surface according to
amino acid levels.[7,10,11]To facilitate RHEB binding to the M-heat, FAT, and N-heat
domains
of mTOR, the switch I and II regions of RHEB should be conformationally
changed by the GTP binding.[3] The switch
I is destabilized upon the catalysis of GTP, which enables GDP/GTP
nucleotide exchange. In contrast, the switch II remains stable during
binding to govern mTORC1 activation.[8,12] Point mutation
in the specific residues (aa 67–77) of switch II prohibits
mTORC1 activation.[13] In addition, because
RHEB is farnesylated to the lysosomal membrane by farnesyltransferases,
the inhibitors of these enzymes impair the post-translational modification
of RHEB, resulting in mTORC1 inhibition.[14,15] However, these inhibitors are nonspecific to the enzymes, which
limits precise targeting to RHEB.[16,17] Recently,
it was reported that a small-molecule inhibitor, NR1 (HY124793), specifically
binds to the switch II region of RHEB, leading to the inhibition of
mTORC1 activity with a micromolar level of half-maximal inhibitory
concentration (IC50).[8] This
result suggests a new modality of mTORC1 inhibition by interfering
with protein–protein interaction (PPI).[8]Targeting the PPI through small molecules is challenging due
to
the large interaction area. Therefore, the sizes of potential PPI
inhibitors should be considered.[18] We used
the structure-based drug design approach to develop peptide inhibitors
to interfere with RHEB binding to mTORC1. Structural studies revealed
that the N-heat domain of mTOR (mTORΔN; aa 60–167)
interfaces RHEB with the α5- and α7-helices, wherein the
α5-helix (aa 101–115) binds with the switch II region
of RHEB (Figure ).[3] Therefore, in this study, we report the development
of a small peptide, P1_WT, that mimics the α5-helix
(Figure D) with a
predicted cell penetration property to interfere with RHEB-mTORC1
interaction, resulting in mTORC1 inhibition. We developed a strategy
of in silico and in vitro methods to identify the proper peptide sequences
with the aim of targeting RHEB with high affinity.[19,20]
Figure 1
Structure-based
design of RHEB-targeting peptides. (A) Molecular
structure of monomeric mTORC1 to illustrate the mode of RHEB binding
with mTOR (PDB ID: 6BCU). Raptor and mLST8 were omitted for clarity.
(B and C) Close-up view of RHEB binding with the mTORΔN domain (aa 60–167), showing α3, α4, α5,
α6, and α7 helices of mTOR colored as indicated in (C).
(D and E) Close-up of RHEB binding with α5-helix (D; aa 101–115)
and α7-helix (E; aa 139–158) of mTOR, showing the residues
involved in binding. Sequences and secondary structures are shown
(helices = cylinders).
Structure-based
design of RHEB-targeting peptides. (A) Molecular
structure of monomeric mTORC1 to illustrate the mode of RHEB binding
with mTOR (PDB ID: 6BCU). Raptor and mLST8 were omitted for clarity.
(B and C) Close-up view of RHEB binding with the mTORΔN domain (aa 60–167), showing α3, α4, α5,
α6, and α7 helices of mTOR colored as indicated in (C).
(D and E) Close-up of RHEB binding with α5-helix (D; aa 101–115)
and α7-helix (E; aa 139–158) of mTOR, showing the residues
involved in binding. Sequences and secondary structures are shown
(helices = cylinders).
Results and Discussion
In Silico
Development of RHEB-Binding Peptides
Previously,
we studied the molecular interactions of RHEB with full-length mTOR
and mTORΔN, yielding KD values of ∼2.4 and 6.4 μM, respectively.[12] In the RHEB-mTORΔN interface,
α5- and α7-helices of mTORΔN stabilize
the mTORC1-RHEB complexation followed by the kinase activations.[3] The residues (aa 63–79) of RHEB in the
switch II region lie between α5- and α7-helices, where
α5 interacts with Q72, T73, S75, I76, and D77 residues of RHEB
(Figure D), while
α7 interacts with the same residues with the extension to bind
with N79 and N-terminal resides (K5, S6, R7, and K8; Figure E).[3]This structural information guided us in applying an in silico
point-mutagenesis approach using ICM-Pro 3.9 software (Molsoft L.L.C.,
USA) to screen and mature the potential affinities of peptides for
RHEB.[21] This was conducted by calculating
the free binding energy (ΔΔGbind) for each residue of the peptide using the following equation:where the free
binding energy was calculated
by subtracting the wild type of free energy (ΔGWT) from the point-mutant free energy (ΔGmutant).[20] Each residue
of α5- or α7-helices was mutated to other possible 19
natural amino acids. As a result, we found that some mutant variants
showed better free binding energy than the wild-type (Figure A). Then, we selected the top
four point mutations for each helix: A(101)R, N(109)M, N(109)K, and
N(113)K of α5-helix (Figure B; hereinafter named P1), and T(139)R,
F(140)R, T(141)F, and L(153)K of α7-helix (Figure C; hereinafter named P2). The wild type and mutant variants were then studied to
determine their binding stability using steered molecular dynamics
(SMD) at various pulling forces.[19] Consequently, P1 variants were found to bind RHEB with more stability than P2 variants, as suggested by the averaged root-mean square
deviation values (Figures D,E, S1, and S2). These data were
collected from 5000 points representing the displacement of the peptide
variants from the initial binding site over the whole simulation time
(10 ns). Therefore, most variants were significantly displaced from
the initial poses at both pulling forces (0.025 or 0.25 Å/ns)
with different binding stability for P1 variants. These
results alongside the calculated physical properties and predicted
cell penetration probability provide strong evidence demonstrating
the efficacy of P1 for the intracellular targeting of
RHEB (Figure F and Table S1).[22]
Figure 2
In silico selection
of RHEB-targeting peptide. (A) In silico point
mutagenesis of the interactive α5 and α7 helices of the
mTORΔN domain (aa 60–167) at the mTOR–RHEB
interface extracted from the mTORC1 structure (PDB ID: 6BCU). The
free energy was calculated using the following equation: ΔΔGbind = ΔGmutant – ΔGWT, and the best results
are indicated. (B and C) Wild type and the corresponding variant sequences
of α5-helix (B; hereinafter named P1) and α7-helix
(C; hereinafter named P2) showing the point-mutated residues
in red color. (D and E) Steered molecular dynamics studies showing
the deviation of P1 variants (D) and P2 variants (E) from the initial
binding pose under two different pulling forces (0.025 and 0.25 Å/ns)
represented by the root-mean square deviation (RMSD). The results
represent the mean of the collected points along the simulation time
(n = 5000) ± SD. Two-way ANOVA was used (****P < 0.0001; ns, P > 0.05) with Tukey’s
test correction. (F) Prediction of cell penetration probability of
the indicated peptides using the MLCPP online tool (http://www.thegleelab.org/MLCPP/MLCPP.html). The results are indicated as a percentage of total probability.
CPP, Cell-penetrating peptides; Non-CPP, Non Cell-penetrating peptides.
In silico selection
of RHEB-targeting peptide. (A) In silico point
mutagenesis of the interactive α5 and α7 helices of the
mTORΔN domain (aa 60–167) at the mTOR–RHEB
interface extracted from the mTORC1 structure (PDB ID: 6BCU). The
free energy was calculated using the following equation: ΔΔGbind = ΔGmutant – ΔGWT, and the best results
are indicated. (B and C) Wild type and the corresponding variant sequences
of α5-helix (B; hereinafter named P1) and α7-helix
(C; hereinafter named P2) showing the point-mutated residues
in red color. (D and E) Steered molecular dynamics studies showing
the deviation of P1 variants (D) and P2 variants (E) from the initial
binding pose under two different pulling forces (0.025 and 0.25 Å/ns)
represented by the root-mean square deviation (RMSD). The results
represent the mean of the collected points along the simulation time
(n = 5000) ± SD. Two-way ANOVA was used (****P < 0.0001; ns, P > 0.05) with Tukey’s
test correction. (F) Prediction of cell penetration probability of
the indicated peptides using the MLCPP online tool (http://www.thegleelab.org/MLCPP/MLCPP.html). The results are indicated as a percentage of total probability.
CPP, Cell-penetrating peptides; Non-CPP, Non Cell-penetrating peptides.Owing to the small size
of P1 variants (13 aa) and
the positive net charge derived from at least three arginine residues,
enhanced CPP was observed over the negatively charged cell membrane
(Figure F). Conversely, P2 variants are rich in glutamic acid, which negatively charges
the peptides and thus decreases their CPP (non-CPP) due to repulsion
of the cell membrane (Figure F). Based on the findings of this computational analysis,
we synthesized P1_WT, P1_A(101)R, P1_N(109)M, double mutants P1_A(101)R_N(113)K, P2_WT, and P2_L(153)K.
mTORC1 Kinase
Assay
The synthetic peptides (purity
90–95%; Figures S3–S8) were
then screened for mTORC1 inhibition through the phosphorylation inhibition
of the downstream signals S6K1 and 4E-BP1.[19] Prestarved HeLa cells were treated with single doses (1 μM)
of peptides or Torin1 (an ATP-competing inhibitor of mTOR) for 3 hours
and then induced with 100 nM of insulin for 30 min to stimulate mTORC1
activation by RHEB. As predicted, P1 variants showed
stronger inhibitory effects on mTORC1 activity than did P2 variants (Figure A).
Figure 3
Inhibition of mTORC1 activity by RHEB-targeting peptides. (A) Screening
of selected peptides for S6K1 and 4E-BP1 phosphorylation inhibition.
HeLa cells were treated with a single dose (1 μM) of peptides
or Torin1 for 3 h under starvation conditions and induced by 100 nM
insulin for 30 min. The results represent the mean of two independent
experiments as a percent of the controls (n = 3)
± SD. Two-way ANOVA was used (****P < 0.0001;
***P < 0.001; ns, P > 0.05)
with
Tukey’s test correction. (B and C) Inhibition curves of T389p-S6K1 (B) and T37/46p-4E-BP1 (C) of prestarved
HeLa cells treated with increasing concentrations of P1_WT or P1_N(109)M peptides for 3 h and induced with 100
nM insulin for 30 min. The results are representative of two independent
experiments as a percent of the controls (n = 3)
± SD. The half-maximal inhibitory concentrations (IC50) are shown.
Inhibition of mTORC1 activity by RHEB-targeting peptides. (A) Screening
of selected peptides for S6K1 and 4E-BP1 phosphorylation inhibition.
HeLa cells were treated with a single dose (1 μM) of peptides
or Torin1 for 3 h under starvation conditions and induced by 100 nM
insulin for 30 min. The results represent the mean of two independent
experiments as a percent of the controls (n = 3)
± SD. Two-way ANOVA was used (****P < 0.0001;
***P < 0.001; ns, P > 0.05)
with
Tukey’s test correction. (B and C) Inhibition curves of T389p-S6K1 (B) and T37/46p-4E-BP1 (C) of prestarved
HeLa cells treated with increasing concentrations of P1_WT or P1_N(109)M peptides for 3 h and induced with 100
nM insulin for 30 min. The results are representative of two independent
experiments as a percent of the controls (n = 3)
± SD. The half-maximal inhibitory concentrations (IC50) are shown.The arginine-rich sequence of P1 variants potentially
enhanced CPP and consequent mTORC1 inhibition, while P2 variants did not. Unexpectedly, P1_WT inhibited mTORC1
activity more strongly than the mutant variants of P1. These results were confirmed by the dose-dependent elucidation
of the inhibitory activity of P1_WT and P1_N(109)M that showed a half-maximal inhibitory concentration (IC50) of ∼0.3 μM and >10 μM for p-S6K1 and p-4E-BP1,
respectively (Figure B,C). The stronger inhibitory effect of P1_WT suggests
its efficiency to penetrate the cell membrane and bind with RHEB at
higher affinity to inhibit mTORC1 activity.
Peptide Binding Kinetics
As previously reported, we
prepared RHEB using the BL21(DE3) E.coli overexpression system.[12] Then, we measured
the in vitro binding affinity (KD) of P1_WT and P1_N(109)M with RHEB by using a bio-layer
interferometry (BLI) system including streptavidin biosensors.[23] For that, we synthesized N-terminal-modified
peptides with biotin by using NHS-(PEG)24-biotin reagent
(Figures S9 and S10). To measure the binding
kinetics, the biotinylated peptides were immobilized onto the streptavidin
biosensors and RHEB was used as analyte in a PBS buffer (pH 7.0) containing
0.02% tween-20. The results showed that P1_WT binds RHEB
with higher affinity (KD = 0.14 μM; Figures A and S11) than P1_N(109)M (KD = 2.54 μM; Figures B and S12), which was associated
with an inhibitory effect on mTORC1 activity.
Figure 4
Binding kinetics of RHEB-targeting
peptides. (A and B) Fitting
curves of the BLI kinetics for P1_WT (A) and P1_N(109)M (B) with RHEB. Biotinylated peptides (100 nM) were immobilized onto
streptavidin biosensors, and RHEB traces was used as analytes as indicated.
Binding kinetics were calculated by the global fitting (1:1 binding)
mode. KD, the equilibrium dissociation
constant; Ka, the association constant; Kd, the dissociation constant. The kinetics parameters
are shown ± standard errors. See Figures S11 and S12 for BLI analysis views. (C) AlphaLISA binding of
mTORΔN domain with RHEB (n = 3). KD value is shown. (D) Inhibition of mTORΔN-RHEB protein–protein interaction by the P1_WT peptide (n = 3). The half-maximal inhibitory
concentration (IC50) is shown.
Binding kinetics of RHEB-targeting
peptides. (A and B) Fitting
curves of the BLI kinetics for P1_WT (A) and P1_N(109)M (B) with RHEB. Biotinylated peptides (100 nM) were immobilized onto
streptavidin biosensors, and RHEB traces was used as analytes as indicated.
Binding kinetics were calculated by the global fitting (1:1 binding)
mode. KD, the equilibrium dissociation
constant; Ka, the association constant; Kd, the dissociation constant. The kinetics parameters
are shown ± standard errors. See Figures S11 and S12 for BLI analysis views. (C) AlphaLISA binding of
mTORΔN domain with RHEB (n = 3). KD value is shown. (D) Inhibition of mTORΔN-RHEB protein–protein interaction by the P1_WT peptide (n = 3). The half-maximal inhibitory
concentration (IC50) is shown.Furthermore, we assayed the inhibitory effect of P1_WT for the
PPI of RHEB-mTORΔN using the AlphaLISA system (PerkinElmer,
USA).[24] We used the 6xHis tagged RHEB (∼20
kDa) and mTORΔN (∼13 kDa) prepared by the
BL21(DE3) E. coli overexpression system.[12] RHEB was first de-tagged by thrombin and then
labeled with biotin using a NHS-(PEG)24-biotin reagent.
Biotinylated RHEB was mixed with different concentrations of 6xHis-mTORΔN followed by the addition of streptavidin-coated donor
beads and anti-6xHis-coated acceptor beads. The results showed that
mTORΔN bound with RHEB by a KD value of 7.27 μM (Figure C), corresponding to the value reported by
the BLI analysis (KD ∼ 6.5 μM).[12] Then, we evaluated the effect of P1_WT to inhibit RHEB-mTORΔN interaction by preincubating
RHEB with a series of different concentration of P1_WT before adding mTORΔN. As a result, P1_WT inhibited RHEB binding to mTORΔN with an IC50 value similar to that of mTORC1 activity inhibition (∼0.3
μM; Figures B–C and 4D).
Conclusions
A new modality of cancer therapy has been emerging by blocking
PPIs involved in mTORC1 activation.[1] In
this study, we aimed to inhibit mTORC1 by interfering with the RHEB-mTOR
interaction on the lysosome. The large PPI area of RHEB-mTOR limits
the efficacy of small molecules for inhibition. Instead, we employed
small peptides to disturb the PPI based on the structural characteristics
of the α5- and α7-helices of N-heat domain of mTOR, which
directly interact with RHEB. Thus, we designed two different peptides, P1 and P2, based on the sequences of α5-
and α7-helices, respectively. We attempted to maturate the peptides
via in silico point mutagenesis; however, the wild-type variant P1_WT
remained the best binder during SMD simulations, corresponding with
the wet experiments to confirm binding affinity. Previously, we successfully
demonstrated that the in silico mutagenesis improved the immune checkpoint
interaction of PD-1/PD-L1.[20] This suggests
that in silico mutagenesis is more effective for proteins than peptides.
In addition, the small size and positive net charge of P1 variants were advantageous to higher cell-penetrating probability,
while P2 variants were not due to the negative net charge.
This resulted in the improved inhibition of mTORC1 activity by P1_WT compared to that of the other selected variants because P1_WT bound with RHEB at the sub-micromolar level, which corresponds
to its inhibitory activity on RHEB-mTORΔN interaction.
Overall, this study is the first to demonstrate that the small peptide-based
compound, P1_WT, inhibits the kinase activity of mTORC1
by disturbing the allosteric regulation of RHEB. P1_WT was designed
based on the structural information involved in the mTORC1-RHEB binding.
Similarly, we will be able to develop other peptide inhibitors based
on PPI modes.
Experimental Section
Analysis of RHEB-mTORC1
Interaction
We used the cryo-EM
structure of mTORC1 complexed with RHEB (PDB ID: 6BCU)[3] to analyze the RHEB-mTORC1 interaction by ICM-Pro 3.9 software
(Molsoft L.L.C., USA).[21] We focused on
the interaction of RHEB with the N-heat domain of mTOR (mTORΔN; aa 60–167). The α5-helix (aa 101–115) and α7-helix
(aa 139–158) of mTORΔN interact with switch
II of RHEB (aa 63–79) to stabilize the PPI. Based on these
interactions, we focused on the helical regions including amino acid
sequence (101-ATRIGRFANYLRN-113) of α5-helix (hereinafter named P1) and (139-TFTAEYVEFEVKRALEWL-156) of α7-helix (hereinafter
named P2).
In Silico Point-Mutagenesis of P1 and P2
Peptides
To
screen and mature the binding affinities of the selected peptide sequences
with RHEB, we applied the in silico point-mutagenesis approach to
each residue independently. Each residue mutated to one of other 19
natural amino acids. We used the TryMutation mode of the ICM-Pro 3.9
software[20,21] to calculate the binding free energy (ΔΔGbind) for each mutated residue according to
the following equation:where the free
binding energy was calculated
by subtracting the wild-type free energy (ΔGWT) from the point-mutant free energy (ΔGmutant). The results indicated that lower binding
free energy correlates with higher binding affinity.
SMD Simulations
We conducted SMD simulations to study
the binding stability of the selected peptides with RHEB, as previously
described.[19] For P1 variants,
we selected P1_WT (ATRIGRFANYLRN), P1_A(101)R (RTRIGRFANYLRN), P1_N(109)M (ATRIGRFAMYLRN), P1_N(109)K (ATRIGRFAKYLRN), and P1_N(113)K (ATRIGRFANYLRK).
For P2 variants, we selected P2_WT (TFTAEYVEFEVKRALEWL), P2_T(139)R (RFTAEYVEFEVKRALEWL), P2_F(140)R (TRTAEYVEFEVKRALEWL), P2_T(141)F (TFFAEYVEFEVKRALEWL), and P2_L(153)K (TFTAEYVEFEVKRAKEWL). We used the super-computing system, SHIROKANE,
of the Human Genome Center (HGC) at the University of Tokyo. We used
the scalable molecular dynamics software NAMD-2.14 acerated with V100
GPU through the visual molecular dynamics interface.[25−27] This interface supported the QwikMD plugin to automatically generate
a rectangular box buffered with 15 Å around macromolecules filled
with 0.15 M NaCl and TIP3 water molecules in a CHARMM36 force field.[28−30] For the peptides, a stream force field (.str) was generated by the
CGenFF server (https://cgenff.paramchem.org/).[31,32] For all SMD simulations, we set the spring
constant of 7 kcal/mol/Å (1 kcal = 69.48 pN·Å). Two
different pulling speeds, 0.025 and 0.25 Å/ns, were applied for
10 ns at 310 K and 1 atm, respectively. The pulling direction was
set along with −Z.
Prediction of Cell-Penetrating
Peptides (CPP)
It is
possible to predict and optimize cell penetration for peptides based
on net charges. Therefore, we used the machine learning-based prediction
of cell-penetrating peptides (MLCPP)[22] framework
to evaluate the probability of the selected peptides entering the
cell. MLCPP is an online platform (http://www.thegleelab.org/MLCPP/MLCPP.html) that employs machine learning models of two-layer prediction framework
based on the calculated properties of the peptide sequence, considering
amino acid sequence, atomic composition, and physiochemical properties.
We entered the peptide sequences using the FASTA format in the field
and submitted the job; the tabulated results appeared within a few
minutes and indicated CPPs and non-CPPs.
Peptide Property Calculation
To analyze the physical
properties of the selected peptides, we calculated the molecular weight
(MW) and net charge (NetC) using the NovoPro peptide property calculator
(https://www.novoprolabs.com/tools/calc_peptide_property) and
the octanol/water partition coefficient (LogP) and aqueous solubility
(LogS) using the ALOGPS 2.1 online program, which is provided by the
Virtual Computational Chemistry Laboratory (VCCLAB) (http://www.vcclab.org /lab/alogps/).[33] We used the peptide sequences for calculations
in NovoPro and used the SMILES format for calculations in ALOGPS 2.1.
Peptide Synthesis, Purification, and Analysis
Peptide
synthesis, purification, and analysis were performed upon order at
the RIKEN Research Resources Division (RRD; RIKEN, Wako, Japan).Synthesis. To evaluate
the inhibition activity of the peptides, we selected P1_WT, P1_A(101)R, P1_N(109)M, P1_A(101)R_N(113)K, P2_WT, and P2_L(153)K for synthesis. P1_WT, P1_A(101)R, P1_A(101)R_N(113)K, P2_WT, and P2_L(153)K were automatically
synthesized by MultiPep CF synthesizer (CEM Corporation, formerly
INTAVIS Bioanalytical Instruments AG), while P1_N(109)M was synthesized by Liberty Blue synthesizer (CEM Corporation). For
BLI binding kinetics, P1_WT and P1_N(109)M were labeled using NHS-(PEG)24-biotin reagent at the
N-terminal.Crude peptide analysis. The crude peptides were analyzed
using a high-performance liquid
chromatography (HPLC) L-2000 system (Hitachi High-Tech Science Corporation,
Japan) at 25 °C with Inertsil ODS-3 (250 × 4.6 mm I.D.)
through a linear gradient mobile phase (1–51%) composed of
0.1% TFA in acetonitrile for 50 min.Peptide purification. The peptides
were purified by HPLC D-7000 (Hitachi High-Tech Science
Corporation, Japan) at 25 °C using InertSustain C18 (250 ×
20 mm I.D.).Pure peptide analysis. The pure peptides were analyzed by
HPLC Chromaster (Hitachi High-Tech
Science Corporation, Japan) at 25 °C using InertSustain C18 (250
× 4.6 mm I.D.).Mass measurements. The mass of the peptides was
measured by the matrix-assisted laser
desorption time-of-flight mass spectrometry (MALDI-TOF MS) using Microflex
spectrometer (Bruker Daltonics, Germany).
mTORC1 Kinase Assay
To evaluate the inhibitory activity
of the peptides on mTORC1 kinase activity, we used AlphaLISA SureFire
Ultra HV p-S6K1 (T389) or p-4E-BP1 (T37/46) assay kits (PerkinElmer,
USA) for detecting mTORC1 phosphorylated products p-S6K or p-4E-BP1,
respectively, as previously described.[19] Briefly, 104 HeLa cells/well were seeded in a 96-well
plate in high-glucose D-MEM (FUJIFILM Wako Pure Chemicals Co., Japan)
supplemented with 10% FBS and 1% P/S and were then incubated overnight
(5% CO2; 37 °C). The cells were starved in Opti-MEM
reduced serum media (Thermo Fisher, USA) for 18 h and then treated
with single dose (1 μM) or increasing concentrations of peptides
for 3 h. After, cells were induced by 100 nM insulin for 30 min to
enhance mTORC1 activation by RHEB. The cells were then lysed with
50 μL of 1× lysis buffer, which was freshly prepared by
shaking for 10 min. Cell lysate (6 μL) was transferred to a
384-well OptiPlate (PerkinElmer, USA) and mixed with 3 μL of
the Acceptor mix of anti-p(T389) S6K or anti-p(T37/46)4E-BP1. They
were then top-sealed, covered, and incubated in the dark for 1 h at
room temperature. Finally, 3 μL Donor mix of anti-S6K or anti-4E-BP1
was added under subdued light, followed by top sealing and incubation
in the dark for >1 h at room temperature. The alpha signal was
measured
by the EnSpire plate reader (PerkinElmer, USA).
Protein Expression
and Purification
In this study,
we prepared RHEB (UniProt ID: Q15382) and mTORΔN (UniProt
ID: P42345), which were expressed and purified as previously described.[12] Briefly, RHEB or mTORΔN genes
were cloned into pET15b vector by In-Fusion cloning kit (Takara, Japan)
and transformed into BL21(DE3) E. coli for protein expression. The 6xHis-tagged RHEB or mTORΔN was then purified by Ni-NTA column (GE Healthcare, USA) and Superdex-200
column (GE Healthcare, USA), concentrated, and stored at −80
°C.[12]
BLI Binding Kinetics
We evaluated the binding kinetics
of RHEB with P1_WT or P1_N(109)M peptides
by the BLI method by using the BLItz instrument (FortéBio,
USA) as previously described.[12,20] For this purpose, the
peptides were N-terminally modified by NHS-(PEG)24-biotin
to immobilize them onto streptavidin biosensors. First, the biosensors
(FortéBio, USA) were hydrated for 1 h in the kinetics buffer
(PBS, pH 7.0 containing 0.02% (v/v) Tween-20). For the measurements,
100 nM biotinylated peptide in the kinetic buffer was immobilized
onto the biosensors. The measurement cycle composed of 30 s initial
baseline (kinetic buffer), 150 s peptide immobilization, 60 s baseline
(kinetic buffer), 150 s RHEB association, and 300 s dissociation phases
(kinetic buffer). A reference cycle was applied for each sensor by
introducing RHEB only in the association phase to exclude nonspecific
binding possibilities. RHEB concentrations were 0.1, 0.3, and 0.5
μM for P1_WT binding measurement and 1, 3, and
5 μM for P1_N(109)M binding measurement. All measurements
were performed at 1000 rpm shaking speed at room temperature. Finally,
we used BLItz Pro 1.2 software (FortéBio, USA) for curve fitting
by using 1:1 binding kinetics.
RHEB de-Tagging and Biotinylation
We used our prepared
tag-cut RHEB as previously described.[12] Briefly, 6xHis-RHEB was incubated with thrombin at a ratio of 1
mg protein: 10 units thrombin incubated overnight at room temperature
on a rotator and purified by His SpinTrap column (GE Healthcare, USA)
and HiTrap Benzamidine FF column (GE Healthcare, USA) to remove 6xHis
tag and thrombin, respectively.Tag-cut RHEB was then labeled
with NHS-(PEG)24-biotin using addition reaction of click
chemistry as previously described to immobilize onto streptavidin-coated
beads for AlphaLISA measurments.[20] Briefly,
protein buffer was exchanged to PBS, pH 7.0 over 10 kDa MW-CO Amikon
filter (Millipore (Merck), Germany). After protein concentration,
RHEB was mixed with NHS-(PEG)24-biotin at a 1:20 molar
ratio in PBS and incubated at room temperature for 2 h. Finally, the
mixture was washed with PBS over 10 kDa MW-CO Amikon filter several
times to remove excess reagent, concentrated, and stored at −80
°C.
AlphaLISA Measurements
To evaluate the effect of P1_WT to inhibit RHEB-mTORΔN interaction,
we used the AlphaLISA-based assay as previously described.[12,20] First, we evaluated RHEB-mTORΔN interaction by
mixing different concentrations of 6xHis-mTORΔN with
1 μM biotinylated RHEB in a 384-well OptiPlateTM (PerkinElmer,
USA) followed by adding 100 μg/mL streptavidin-coated donor
beads and 200 μg/mL anti-6xHis-coated acceptor beads. The plate
was then sealed, covered, and incubated in the dark for >1 h at
room
temperature. The alpha signal was then measured by an EnSpire plate
reader.The effect of P1_WT on RHEB-mTORΔN interaction was evaluated using the same method with some modifications.
Increasing concentrations of P1_WT were incubated with
1 μM biotinylated RHEB; then, 1 μM of 6xHis-mTORΔN was added followed by 100 μg/mL of streptavidin-coated donor
beads and 200 μg/mL of anti-6xHis-coated acceptor beads. The
plate was then sealed, covered, and incubated in the dark for >1
h
at room temperature. The alpha signal was then measured by an EnSpire
plate reader.
Data Analysis
Statistical significance
and the number
of samples are noted in figure legends where appropriate. Data are
expressed as mean ± SD. Two-way ANOVA was used as indicated;
**** for P < 0.0001 and ns for P > 0.05 with Tukey’s test correction. Statistical analyses
were performed using GraphPad Prism software v.9.3 (GraphPad, USA).
Authors: James C Phillips; Rosemary Braun; Wei Wang; James Gumbart; Emad Tajkhorshid; Elizabeth Villa; Christophe Chipot; Robert D Skeel; Laxmikant Kalé; Klaus Schulten Journal: J Comput Chem Date: 2005-12 Impact factor: 3.376
Authors: Igor V Tetko; Johann Gasteiger; Roberto Todeschini; Andrea Mauri; David Livingstone; Peter Ertl; Vladimir A Palyulin; Eugene V Radchenko; Nikolay S Zefirov; Alexander S Makarenko; Vsevolod Yu Tanchuk; Volodymyr V Prokopenko Journal: J Comput Aided Mol Des Date: 2005-06 Impact factor: 3.686
Authors: Haijuan Yang; Derek G Rudge; Joseph D Koos; Bhamini Vaidialingam; Hyo J Yang; Nikola P Pavletich Journal: Nature Date: 2013-05-01 Impact factor: 49.962