Amanda Fernandes de Medeiros1, Beatriz Blenda Pinheiro de Souza2,3, Lucas Pinheiro Coutinho4, Aline Melro Murad3, Paula Ivani Medeiros Dos Santos5, Norberto de Kássio Vieira Monteiro6, Elizeu Antunes Dos Santos1,7,8, Bruna Leal Lima Maciel9,10, Ana Heloneida de Araújo Morais1,9,10. 1. Postgraduate Biochemistry Program, Biosciences Center, Federal University of Rio Grande do Norte, Natal, Brazil. 2. Postgraduate Biological Molecular, Institute of Biological Sciences, University of Brasília, Brasília, Brazil. 3. Laboratory of Mass Spectometry-LEM, Embrapa Genetic Resources and Biotechnology, Brasília, Brazil. 4. Chemistry Course, Science Center, Federal University of Ceará, Fortaleza, Brazil. 5. Federal Institute of Education, Science and Technology of Rio Grande do Norte, Brazil. 6. Analytical Chemistry and physical Chemistry Department, Science Center, Federal University of Ceará, Fortaleza, Brazil. 7. Department of Biochemistry, Biosciences Center, Federal University of Rio Grande, Natal, Brazil. 8. Tropical Medicine Institute, Federal University of Rio Grande do Norte, Natal, Brazil. 9. Postgraduate Nutrition Program, Center for Health Sciences, Federal University of Rio Grande do Norte, Natal, Brazil. 10. Department of Nutrition, Center for Health Sciences, Federal University of Rio Grande do Norte, Natal, Brazil.
Abstract
Trypsin inhibitors from tamarind seed have been studied in vitro and in preclinical studies for the treatment of obesity, its complications and associated comorbidities. It is still necessary to fully understand the structure and behaviour of these molecules. We purifed this inhibitor, sequenced de novo by MALDI-TOF/TOF, performed its homology modelling, and assessed the interaction with the trypsin enzyme through molecular dynamics (MD) simulation under physiological conditions. We identified additional 75 amino acid residues, reaching approximately 72% of total coverage. The four best conformations of the best homology modelling were submitted to the MD. The conformation n°287 was selected considering the RMSD analysis and interaction energy (-301.0128 kcal.mol-1). Residues Ile (54), Pro (57), Arg (59), Arg (63), and Glu (78) of pTTI presented the highest interactions with trypsin, and arginine residues were mainly involved in its binding mechanism. The results favour bioprospecting of this protein for pharmaceutical health applications.
Trypsin inhibitors from tamarind seed have been studied in vitro and in preclinical studies for the treatment of obesity, its complications and associated comorbidities. It is still necessary to fully understand the structure and behaviour of these molecules. We purifed this inhibitor, sequenced de novo by MALDI-TOF/TOF, performed its homology modelling, and assessed the interaction with the trypsin enzyme through molecular dynamics (MD) simulation under physiological conditions. We identified additional 75 amino acid residues, reaching approximately 72% of total coverage. The four best conformations of the best homology modelling were submitted to the MD. The conformation n°287 was selected considering the RMSD analysis and interaction energy (-301.0128 kcal.mol-1). Residues Ile (54), Pro (57), Arg (59), Arg (63), and Glu (78) of pTTI presented the highest interactions with trypsin, and arginine residues were mainly involved in its binding mechanism. The results favour bioprospecting of this protein for pharmaceutical health applications.
Proteolytic enzymes are critical elements for several biological events, such as digestion, healing, viral replication and the blood clotting cascade, and these processes must be regulated with great precision,. Within proteolytic enzymes, more than a third are serine proteases (EC 3.4.21), which are subdivided into families according to catalytic mechanisms and specificity, especially the S1 peptidase family, with trypsin (EC 3.4.21.11) as the most representative member being mainly related to the digestive process.The proteolytic enzymes are regulated by protease inhibitors (PIs) and these natural regulators are present in multiple forms in the animals, plants and microorganisms,. The purification and characterisation of PIs biomolecules may lead to a molecule with structure and functions and used as a potential biomolecule for herbal medicine nutraceutical or even a medicine.The trypsin inhibitors from tamarind seed have already been studied from the perspective of biotechnology and health applications in vitro and in preclinical studies. In our previous studies, eutrophic Wistar rats fed with partially purified trypsin inhibitor from Tamarind seeds (TTI) increased serum levels of the cholecystokinin (CCK) hormone, reduced food consumption and weight gain without altering true digestibility. On the other hand, Costa et al., evaluating the effect of TTI on plasma CCK and leptin in Wistar rats with diet-induced obesity, found no increases in CCK with TTI treatment, despite decreased leptin concentrations. The treatment with TTI in Wistar rats with diet-induced obesity and metabolic syndrome (MS) reduced food consumption and the plasma concentration of tumour necrosis factor-alpha (TNF-α) regardless of weight loss, characterising TTI as a molecule with anti-inflammatory characteristics.To increase its potential effect, TTI was nanoencapsulated in chitosan and whey-protein (ECW), and results revealed an improvement in the function and stability of TTI. In another study, ECW induced a significant reduction in fasting blood glucose in Wistar rats fed with a high glycemic index and glycemic load diet (HGLI). In addition, nanoencapsulation protected TTI and promoted the controlled release in vitro simulations of digestion under physiological conditions. Also, TTI nanoparticles revealed no cytotoxic effect in intestinal cell lines (Caco-2 and CDD18-Co), nor significant changes in hematological parameters, liver and kidney functions monitored by blood markers in Wistar rats under a HGLI diet, representing a promising formulation for the clinical application.Another study purified the trypsin inhibitor from Tamarind seeds (pTTI), verifying it was a competitive inhibitor with 19.578 kDa, resistant to high temperatures and extreme pH. Medeiros et al. also evaluated the effect of pTTI on plasma CCK and leptin in Wistar rats with diet-induced obesity, and results similar to the partially purified TTI were found, without altering CCK concentrations, but reducing leptin, suggesting an improvement in the leptin resistance profile, a characteristic present in obesity and SM. Further, Carvalho et al. indicated pTTI, was able to reduce molecular and plasm levels of TNF-α (tumour necrosis factor-α), VLDL-C (very-low-density lipoprotein-cholesterol), and TG (triglycerides) in Wistar rats with diet-induced dyslipidemia. Thus, we have a protein with beneficial biological effects in experimental studies applied in the context of obesity and associated complications.However, information about a better understanding of the structure and physical-chemical properties of pTTI are important for pTTI function’s enlightenment. Once the variety of bioactive compounds, such as PIs, in nature is limited, synthetic peptides presenting similar biological activity may be an alternative to overcome this barrier,.Comprehension of the antitrypsin activity mechanism and its structure, through molecular dynamics (MD) under physiological conditions, may contribute to future research, such as in the synthesis of a bioactive peptide from pTTI. It is of great scientific interest to understand the interactions between biomolecules that are involved in biological reactions, especially those that are related to diseases at a global level, such as obesity, metabolic syndrome and its associated comorbidities; and to achieve this study objective, MD appear as an important computational tool.According to several discoveries in preclinical studies, biological activities attributed to TTI and considering this molecule is limited in nature, it is important to expand the protein identification to provide technical and scientific advances, such as the generation of bioactive peptides derived from the pTTI that attract the interest of the pharmaceutical industry. In this present work, the protein sequence of pTTI was extended, homology modelling was performed, and its interaction with the enzyme trypsin was validated using MD simulation.
Materials and methods
Plant material
The tamarind (Tamarindus indica L.) fruit, was obtained from the seed bank in Natal-RN, Brazil, and registered in the National System for the Management of Genetic Heritage and Associated Traditional Knowledge under the number AF6CE9C. The pulp was removed, the seeds were peeled and ground, and used for the procedures.
Purification and sequencing of the tamarind seed trypsin inhibitor
pTTI purification and preparation for MSMS analysis
The purification process was performed according to Medeiros et al.. The disulphide bridges from cys were reduced, and the free thiol groups alkylated to obtain the protein’s secondary structure. The reduction and alkylation steps were carried out according to Medeiros et al. and the reduced and alkylated pTTI fractions were subjected to chromatographic analysis. The procedure was performed on reverse-phase HPLC on an analytical column (Pharmacia Biotech µRCP C2/C18 ST 4.6/100 mm, 120 Å, cod. No.175057–01). Solvent A (analytical water grade + 0.1% of trifluoroacetic acid (TFA)) and B (acetonitrile (ACN) + 0.1% TFA) were used. The pTTI purification was performed using the following gradient steps: three minutes with 5% of solvent B, to desalinate and remove the surfactant from the sample, a linear gradient of 5–95% of solvent B in 22 min (4.09% B.min−1), and a final step, of 5 min at 95% of solvent B, with a flow of 1 mL.min−1, monitored by UV detection at 216 nm for peptide detection and 280 nm for aromatic ring detection, in a run of 30 min.
Proteolytic digestion of pTTI and MALDI-TOF/TOF MS/MS mass spectrometry
The cleavage of the reduced and alkylated pTTI was performed with two sequencing grade proteolytic enzymes, according to the manufacturer’s instructions: trypsin (Sequencing Grade Modified Trypsin, Promega trypsin, v511A), with the use of surfactant (RapiGest SF, Water, Part 186001861) and endoproteinase GluC (Staphylococcus aureus Protease V8–Typical GluC Digest, P8100).The products generated by the enzymatic cleavages with trypsin and GluC were centrifuged at 24,500 × g and 6 °C for 30 min, and the supernatants were eluted in a linear gradient of 5–95% solvent B in a C18 analytical column (Pharmacia Biotech µRCP C2/C18 ST 4.6/100 mm, 120 Å, cod. No.175057–01) with a flow of 1 mL.min−1. Each fraction was collected manually and analysed by mass spectrometry using MALDI-TOF/TOFMS/MS in positive mode.Then, they were analysed in an ultrafleXtreme mass spectrometer (Bruker Daltonics, Bremen–Germany) after mixing 1:3 (v/v) with saturated solution of α-cyano-4-hydroxycinnamic acid (α-CHCA) (5 mg CHCA, 250 μL of ACN, 50 μL TFA 3% and 200 μL of deionisedH2O) and application in MALDI plate. The MS and spectra were obtained in positive reflected mode in the range of m/z between 700 and 4500 and MS/MS spectra were acquired in the same mass range in LIFTTM method after external calibration using calibration standards (Protein Calibration Standard I and II) (Bruker Daltonics, Bremen–Germany). The new sequencing of the pTTI triptych fragments was performed by signalling and manual interpretation of the spectra, using the FlexAnalysis 3.4 software (Bruker Daltonics, Bremen–Germany).
Computational methods
Multiple alignments
The pTTI protein sequence obtained by MS/MS was compared with other sequences present in the protein database of NCBI (National Centre for Biotechnology Information) (https://www.ncbi.nlm.nih.gov/) and PDB (Protein Data Bank proteins) files sing the BLASTp algorithm (Basic Local Alignment Search Tool–protein) (https://blast.ncbi.nlm.nih.gov/Blast.cgi), with multiple alignments of protein sequences on the Clustal Omega server (https://www.ebi.ac.uk/Tools/msa/clustalo/).
Homology modelling and validation of predicted models
The three-dimensional structure of the pTTI was determined based on multiple alignments with the four proteins, of vegetable origin, with the highest percentage identity (% ID) deposited on the NCBI/PDB database server: 4AN7_B, 4AN6_A, 4J2K_A and 1TIE_A (https://www.rcsb.org/). Amino acids from the sequence with higher % ID protein sequence were used (4AN7_B) to fill the gaps. Homology modelling was performed in the Modeller 9.21 program (https://salilab.org/modeller/) to create 150 models, and models with the highest negative modules Discrete Optimised Protein Energy (DOPE) values were selected and validated by MolProbity server (http://molprobity.biochem.duke.edu/), according to the following parameters: residues with bad bonds, residues with bad angles, MolProbity score and Ramachandran plot. The model showing better stability and less steric impediment was chosen as the three-dimensional structure of the pTTI.
Conformations
The best homology model of the pTTI, was used to generate 500 different conformations through CONCOORD (from CONstraints to COORDinates) version 2.1.2 (https://www3.mpibpc.mpg.de/groups/de_groot/concoord/), according to geometric restrictions, which are based on the strength of the interaction and with a set of upper and lower geometric limits, for all pairs of atoms in interaction. From all conformations generated by the program, four were selected based on the sum of violations, numerical value related to the correction of the distances of the positions of atoms that are involved in interatomic distances that violate the upper or lower limit of distance. The correction process is interrupted when the sum of violations is zero.
Molecular dynamics (MD) simulation
From the four best pTTI conformations generated by CONCOORD, MD simulations were performed between each conformation and trypsin enzyme (PDB ID 2PTN). MD simulations were performed using GROningen MAchine for Chemical Simulations (GROMACS) version 2018.4 package implemented with the CHARMM36 force field. The pTTI conformers and trypsin were placed 10 Å away before the simulations started. The transferable intramolecular potential with 3 points (TIP3P) water molecules were used to solvate the simulated systems. The systems neutralisation was achieved through the addition of counter ions.The systems’ geometry was minimised by the steepest descent algorithm for 50,000 steps with a tolerance of 100 kJ mol−1 nm−1 followed by conjugate gradient algorithm for 10,000 steps with a tolerance of 100 kJ mol−1 nm−1.The Leap-Frog algorithm was applied to integrate the motion equation with time step of 2.0 fs. The long-range interactions were modelled using particle-mesh Ewald sum (PME) with a cut-off of 1.2 nm. The van der Waals interactions were also calculated using the same threshold. Bonds involving hydrogen atoms were restrained using LINCS algorithm. The Nosé–Hoover thermostat, was used to fix the system temperature (310 K) in all production simulations, while the system pressure was controlled using a Parrinello–Rahman barostat in the NPT simulations. Four shorts 200 ps equilibrium dynamics with NVT and NPT ensembles were performed. Finally, 200 ns production MD simulation using NVT ensemble was carried out for each system to determine its interaction with trypsin.After MD simulations, the most stable interaction complex model between the pTTI conformers and trypsin was analysed through the Interaction Potential Energy (IPE) (in kcal.mol−1), which can be defined as the total interaction energy between two groups (the sum of van der Waals and electrostatic contributions), calculated according to the Equation (1):
where IPE is the interaction energy between a group of atoms i and a group of atoms j, and N and N are the total number of atoms on groups i and j, V and V are the terms corresponding to electrostatic and van der Waals contribution, respectively. This parameter is often used to evaluate interaction energies in protein–ligand and protein–protein systems.
Results
Purification and determination of molecular weight of the trypsin inhibitor of purified tamarind seed
The chromatogram demonstrates the protein profile of the partially purified inhibitor (Figure 1(A)). The major peak was collected and named pTTI, eluted in 35 min in 45% of solvent B (ACN/TFA 0.1%). The pTTI presented different ions with m/z values [M + H]+ = 19,586 Da; [M + 2H]+ = 9792 Da e [M + 3H]+ = 6528 Da (Figure 1(B)).
Figure 1.
Identification of native pTTI and postreduction and alkylation by reverse-phase high-performance liquid chromatography and mass spectrometry. (A) Chromatographic profile of the TTI in the C18 Vydac reverse-phase analytical column. (B) Mass spectrum with the pTTI MALDI-TOF ionisation source. (C) Chromatographic profile of the reduced and alkylated pTTI in a C2/C18 Pharmacia Biotech µRCP analytical column. (D) Mass spectrum with reduced and alkylated pTTI MALDI-TOF ionisation source. TTI: Partially purified trypsin inhibitor from tamarind seeds. pTTI: Purified trypsin inhibitor from tamarind seeds.
Identification of native pTTI and postreduction and alkylation by reverse-phase high-performance liquid chromatography and mass spectrometry. (A) Chromatographic profile of the TTI in the C18 Vydac reverse-phase analytical column. (B) Mass spectrum with the pTTI MALDI-TOF ionisation source. (C) Chromatographic profile of the reduced and alkylated pTTI in a C2/C18 Pharmacia Biotech µRCP analytical column. (D) Mass spectrum with reduced and alkylated pTTI MALDI-TOF ionisation source. TTI: Partially purified trypsin inhibitor from tamarind seeds. pTTI: Purified trypsin inhibitor from tamarind seeds.The pTTI, when reduced and alkylated, showed a protein peak with greater absorbance with a retention time (RT) of 51 min and solvent concentration B of approximately 41% (Figure 1(C)). Consequently, pTTI was analysed by MALDI-TOF, verifying the reduction and alkylation reaction. The pTTI presented two predominant ions with m/z values [M + 3H]+ = 6602,789 Da and m/z [M + 2H]+ = 9,906,420 Da (Figure 1(D)), which indicates pTTI. It can be inferred that the pTTI is reduced and alkylated, since the difference in mass observed shows the presence of the acetamide group in the reduced thiol groups.
Proteolytic digestion of pTTI
The chromatograms of the digestion profiles of both enzymes are shown in Figure 2. Within the initial 3 min of chromatography after enzymatic digestion, the digestion salts are eluted (Figure 2(A,B)), and in the case of trypsin digestion, the surfactant is also eluted (Figure 2(A)), followed by fractions resulting from digestion.
Figure 2.
Chromatograms of the reduced and alkylated pTTI by reverse-phase high-performance liquid chromatography treated with digestive enzymes with degree sequencing for mass spectrometry. (A) Chromatographic profile of pTTI digested with trypsin in a C2/C18 Pharmacia Biotech µRCP analytical column. (B) Chromatographic profile of pTTI digested with GluC in a C2/C18 Pharmacia Biotech µRCP analytical column. pTTI: Purified trypsin inhibitor from tamarind seeds.
Chromatograms of the reduced and alkylated pTTI by reverse-phase high-performance liquid chromatography treated with digestive enzymes with degree sequencing for mass spectrometry. (A) Chromatographic profile of pTTI digested with trypsin in a C2/C18 Pharmacia Biotech µRCP analytical column. (B) Chromatographic profile of pTTI digested with GluC in a C2/C18 Pharmacia Biotech µRCP analytical column. pTTI: Purified trypsin inhibitor from tamarind seeds.
Sequencing and analysis by mass spectrometry by MALDI-TOF/TOF MS/MS
The peptide sequences from reduced and alkylated pTTI determined were show in Table 1, with peptides number 1, 2, 4 and 6 originating from digestion with trypsin and peptides number 3 and 5 coming from digestion with GluC. The ions acquired by MALDI-TOF/TOF, of each proteolytic digestion, with sequential signals are highlighted in Supplementary Fig. S1. The mass spectra of the peptides with the y and b-ions series assigned are shown in Supplementary Fig. S2.
Table 1.
Sequences of peptide fragments of the reduced and alkylated pTTI determined by MALDI TOF/TOF in LIFT method.
Peptide
Sequence
[M + H]+Theoretica
[M + H]+ (MALDI)
1
DTVHDTDGQVP(I/L)NNAQGYY(I/L) (I/L)PAQQGK
2942.42 Da
2942.46 Da
2
(I/L)FDEQSSEKGYTPVK
1727.84 Da
1727.84 Da
3
QGYT(PV)K(I/L)SDDFSSAAPFK (I/L)KQFEE
2832.37 Da
2832.30 Da
4
QGFEEDYK(I/L)VYC*SK
1765.77 Da
1765.96 Da
5
DYK(I/L)VYC*SKSE
1391.65 Da
1391.62 Da
6
(I/L)V(I/L)KEGDPFKVKFKKVDEES
2335.29 Da
2335.25 Da
C*: Carbamidomethylcysteine; pTTI: Purified trypsin inhibitor from tamarind seeds.
Sequences of peptide fragments of the reduced and alkylated pTTI determined by MALDI TOF/TOF in LIFT method.C*: Carbamidomethylcysteine; pTTI: Purified trypsin inhibitor from tamarind seeds.The sequencing identified additional 75 amino acid residues, expanding the coverage percentage of the sequencing from 28% (with the sequencing of 53 previously published residues) to 72% in total, comprising 128 amino acid residues elucidated from the pTTI. Peptide 4 showed a glycine residue, which is not found in peptide number 3 that contains part of the same fragment, this may occur due to the presence of pTTI isoform fragment. It is not possible to distinguish leucine and isoleucine residues and glutamine and lysine, by mass spectrometry because their molecular masses are identical, with 113 Da and 128 Da, respectively.Table 2 shows the proteins with the highest identity with the pTTI. According to Figure 3, it is possible to observe the alignment of the primary structure of the pTTI with the sequences obtained in the BLASTp database, which were used to perform the homology modelling. Among the acquired sequences, it was possible to fill in the gaps remaining in the sequencing of the pTTI with the sequence of the highest percentage of identity, 4AN7_B, which comprises: 28 residues between numbers 66 and 93, 21 residues from number 138 to 158 and 7 residues in C’-terminal portion 178 to 184; totalling the addition of 56 amino acid residues.
Table 2.
Values obtained from the protein sequences selected in the BLASTp search.
Accession
Max Score
Total Score
Query Cover
E value
%ID
Chain
4AN7_B
191.0
191.0
99%
7e–63
66.48%
β
4AN6_A
188.0
188.0
99%
7e–62
65.34%
α
4J2K_A
64.3
64.3
93%
2e–13
30.77%
α
1TIE_A
59.7
59.7
78%
1e–11
33.08%
α
BLASTp: Basic Local Alignment Search Tool–protein; %ID: percentage identity.
Figure 3.
Multiple alignment of the primary partial sequence of pTTI with other sequences deposited in the tool Protein BLAST/NCB tooI. pTTI: Purified trypsin inhibitor from tamarind seeds.
Multiple alignment of the primary partial sequence of pTTI with other sequences deposited in the tool Protein BLAST/NCB tooI. pTTI: Purified trypsin inhibitor from tamarind seeds.Values obtained from the protein sequences selected in the BLASTp search.BLASTp: Basic Local Alignment Search Tool–protein; %ID: percentage identity.The proteins mentioned above were used as a basic model for homology construction, and the five models with the highest DOPE score are shown in Table 3.
Table 3.
The five best three-dimensional models generated by Modeller and their respective DOPE values.
Model Number
DOPE
37
–13092.01855
56
–13057.27832
63
–13081.18359
94
–13165.39648
148
–13107.49219
DOPE: Discrete Optimised Protein Energy.
The five best three-dimensional models generated by Modeller and their respective DOPE values.DOPE: Discrete Optimised Protein Energy.The results of the validation in MolProbity to validate the predicted models are shown in Table 4.
Table 4.
MolProbity selection criteria for the five best models generated by Modeller.
Model Number
Ramachandran favoured
Residues with bad angles
Residues with bad bonds
MolProbity score
37
86.81%
0.00%
1.56%
3.77
56
87.36%
0.00%
1.11%
3.44
63
86.26%
0.00%
1.81%
3.63
94
89.01%
0.00%
1.86%
3.70
148
87.91%
0.00%
1.36%
3.67
MolProbity selection criteria for the five best models generated by Modeller.Due to a set of factors, such as residues with bad bonds and MolProbity score, the best result was achieved for model number 56, hence this model was used to generate 500 possible conformations for pTTI structure. among these conformations, four were selected according to the sum of the violations and the geometric restrictions (Figure 4).
Figure 4.
Visualisation of the four best conformations generated by CONCOORD of model n. 56 of the pTTI and their respective values of the sum of violations for conformations number 177 (A), 200 (B), 287 (C), and 323 (D). The N′-terminal portion which is a residue of Asp is highlighted in red and the C’-terminal portion, which is a residue of Val in tan, in all images. (The colours: yellow, β-sheet; tan, β-bridge; white, coil; cyan, turn; red, N’-terminal portion; tan, C’-terminal portion). pTTI: Purified trypsin inhibitor from tamarind seeds.
Visualisation of the four best conformations generated by CONCOORD of model n. 56 of the pTTI and their respective values of the sum of violations for conformations number 177 (A), 200 (B), 287 (C), and 323 (D). The N′-terminal portion which is a residue of Asp is highlighted in red and the C’-terminal portion, which is a residue of Val in tan, in all images. (The colours: yellow, β-sheet; tan, β-bridge; white, coil; cyan, turn; red, N’-terminal portion; tan, C’-terminal portion). pTTI: Purified trypsin inhibitor from tamarind seeds.The best conformation of the pTTI obtained from model number 56 was subjected to MD simulations to analyse pTTI interactions with trypsin enzyme. The MD revealed a Root Mean Square Deviation (RMSD) plot as a function of time and in an aqueous system (Supplementary Fig. S2). Regarding the total energy of interaction, taking as a minimum and maximum point the interval between sudden changes, the results demonstrated the interaction conformation number 287 has the lowest interaction energy value (Table 5), suggesting, the conformation n. 287 was the best model and, therefore, may represent the putative pTTI structure.
Table 5.
Interaction potential energy of the molecular dynamics between the four best conformations of the pTTI with the enzyme trypsin.
Conformer number
IPE (kcal.mol–1)
177
–172.8038
200
–235.1159
287
–301.0128
323
–96.5838
IPE: interaction potential energy.
Interaction potential energy of the molecular dynamics between the four best conformations of the pTTI with the enzyme trypsin.IPE: interaction potential energy.Figure 5 shows the 287 pTTI conformer and tripsin before, during and after the MD simulation. The amino acid residues of pTTI that most interacted with trypsin were Arg (59), Arg (63), Ile (54), Pro (57) and Glu (78) (Figure 6), which showed lower IPE. Among the residues, Arg (59) and (63) are those that have the greatest interaction with trypsin, and both are highlighted in grey (above Arg (63) and below Arg (59)) when MD occurs (Figure 5(B)).
Figure 5.
Visualisation of the stages of molecular dynamics simulation between model number 56 in conformation number 287 of pTTI with the enzyme trypsin. (A) Immediately before the MD, with the molecules positioned 10 A apart. (B) During MD, arginine residues are discarded, above Arg 63 and below Arg 59 (in grey, both). (C) Immediately after MD. (The colours: yellow, β-sheet; tan, β-bridge; white, coil; cyan, turn; blue, 310-helix; purple, α-helix). pTTI: Purified trypsin inhibitor from tamarind seeds.
Figure 6.
Representative graph of the interaction energy of each amino acid residue of the pTTI and enzyme trypsin. The highlighted residues have the lowest interaction energy, which represents those that have the greatest interaction with trypsin, among them Arg (59), Arg (63), Ile (54), Pro (57) and Glu (78). pTTI: Purified trypsin inhibitor from tamarind seeds.
Visualisation of the stages of molecular dynamics simulation between model number 56 in conformation number 287 of pTTI with the enzyme trypsin. (A) Immediately before the MD, with the molecules positioned 10 A apart. (B) During MD, arginine residues are discarded, above Arg 63 and below Arg 59 (in grey, both). (C) Immediately after MD. (The colours: yellow, β-sheet; tan, β-bridge; white, coil; cyan, turn; blue, 310-helix; purple, α-helix). pTTI: Purified trypsin inhibitor from tamarind seeds.Representative graph of the interaction energy of each amino acid residue of the pTTI and enzyme trypsin. The highlighted residues have the lowest interaction energy, which represents those that have the greatest interaction with trypsin, among them Arg (59), Arg (63), Ile (54), Pro (57) and Glu (78). pTTI: Purified trypsin inhibitor from tamarind seeds.Gathering all information obtained through CONCOORD and MD, the model number 56 in its conformation number 287 may represent the pTTI structure and its is interaction with trypsin, in a system under physiological conditions.
Discussion
Countless biological processes involving ligand-protein interactions and substances with potential bioactivity are investigated concerning healthcare applications. Considering previously studies carried out with TTI and its bioactive properties in preclinical studies,,,,,, it is crucial to understand acutely this biomolecule.In the current work, we pursue more information regarding a trypsin inhibitor’s identity extracted from a vegetal source of the tamarind fruit seeds’,. The protein sequence obtained through mass spectrometry considerably increased the coverage percentage compared to the previously known protein sequence. The former coverage sequence was nearby 52 amino acid residues and our efforts expanded the sequence to a total of 128 amino acid residues, representing 72% of the total pTTI molecular mass (19,586 Da).Analysis regarding the protein identification pTTI exhibited high identity to tamarind Kunitz inhibitor–TKI (66% of identity–represents 127 hits of the 128 sequenced pTTI residues). The pTTI protein sequence’s missing gaps were manually completed considering the TKI protein sequence identity deposited in BLASTp database. Nevertheless, even with high similarity between both sequences, a difference of 989 Da in molecular masses between pTTI and TKI remains. This difference may be because they are isoforms.The pTTI model was built up from other proteins recognised as trypsin inhibitors and also extracted from seeds, including two proteins from Tamarindus indica, one from Enterolobium contortisiliquum, and lastly one from Erythrina caffra. All proteins described were extracted from plants belonging to the Fabaceae family, which impacts the relation of protein identity between them. In silico modelling demonstrated to be an extremely useful tool for hypotheses application, especially in molecular biology/biochemistry providing clearance regarding putative active sites and protein-ligand interactions and specificity of purified or designed biomolecules for drug discovery. The three-dimensional structure model of the pTTI provided constructive information concerning pTTI interaction with trypsin, reinforcing the need for future studies with other molecules that may interact with pTTI, with applications in the healthcare system.The MD performed in this study was between two large proteins under physiological conditions, unusual computational parameters, nonetheless with great relevance and application in the healthcare area. The best model fitting the MD parameters revealed arginine residue in two positions, which is positively charged; the isoleucine and proline residues, which have nonpolar and aliphatic characteristics; and glutamate, which has a negative charge. Considering the two residues that presented the most significant interaction, Arg (59) and Arg (63), it is possible to evaluate the biomolecule dynamics and infer that the active site of trypsin for interaction with pTTI has a negative character. Thus, electrostatic interactions are observed between pTTI and trypsin, and these residues should receive attention as the target of further studies about pTTI.A recently published study proposed a model system that measures the binding mechanism between trypsin and its inhibitor bovine pancreatic trypsin inhibitor (BPTI)–an inhibitor type Kunitz, known as aprotinin – using MD approach. The electrostatic bonding is mainly performed by Van der Waals force, with the trypsin binding cavity being negative and BPTI positively charged. BPTIArg (17) is a crucial residue that requires conformational rearrangement to bind into trypsin’s active site perfectly and is responsible for the protein complex maintainance. Therefore, the residues Arg (59) and Arg (63) of pTTI may also be critical for the pTTI inhibition mechanism.Furthermore, the importance of the aqueous medium during the interaction process’ between proteins is also discussed. In the trypsin and BPTI interaction, the solvation layer is an agent that enhances the electrostatic force between the molecules, bringing the molecules closer together and displacing the water molecules throughout protein-protein interactions. Consequently, an MD under physiological conditions is a crucial factor for the effective study of understanding the interaction, considering water in the environment through in vivo/in vitro reactions. These insights are essential to study the mechanism of interaction between proteins, which is suggested to be similar between several protein-protein complexes and may be extrapolated to the interaction between pTTI and trypsin.Previous studies have shown the activities and health application perspectives of Tamarind enzyme inhibitors, and with pTTI, findings and insights are gathered in Figure 7. The exploration of computational technologies is a handling tool for unravelling the pTTI structure and innovation towards the bioprospecting of active molecules from plant origin, useful for the biotechnology industry. Designed bioactive peptides, MD studies with other proteins related to energy metabolism (membrane receptor), and MD with other compounds support new application perspectives, whether in silico, in vitro, or in vivo studies. Therefore, the study of these molecules directly influences future proposals for application in healthcare-related to obesity, its complications, and associated comorbidities, thus linked to a possibility of interest to the pharmaceutical industry.
Figure 7.
Summary of claims in pre-clinical studies already documented, findings and insights related to the biochemical characterisation of pTTI presented in this paper. pTTI: Purified trypsin inhibitor from tamarind seeds.
Summary of claims in pre-clinical studies already documented, findings and insights related to the biochemical characterisation of pTTI presented in this paper. pTTI: Purified trypsin inhibitor from tamarind seeds.Click here for additional data file.
Authors: Izael S Costa; Amanda F Medeiros; Fabiana M C Carvalho; Vanessa C O Lima; Raphael P Serquiz; Alexandre C Serquiz; Vivian N Silbiger; Raul H Bortolin; Bruna L L Maciel; Elizeu A Santos; Ana H A Morais Journal: Obes Facts Date: 2018-12-11 Impact factor: 3.942
Authors: Christopher J Williams; Jeffrey J Headd; Nigel W Moriarty; Michael G Prisant; Lizbeth L Videau; Lindsay N Deis; Vishal Verma; Daniel A Keedy; Bradley J Hintze; Vincent B Chen; Swati Jain; Steven M Lewis; W Bryan Arendall; Jack Snoeyink; Paul D Adams; Simon C Lovell; Jane S Richardson; David C Richardson Journal: Protein Sci Date: 2017-11-27 Impact factor: 6.725
Authors: Dongwen Zhou; Yara A Lobo; Isabel F C Batista; Rafael Marques-Porto; Alla Gustchina; Maria L V Oliva; Alexander Wlodawer Journal: PLoS One Date: 2013-04-23 Impact factor: 3.240
Authors: Paula Judith Perez Espitia; Nilda de Fátima Ferreira Soares; Jane Sélia Dos Reis Coimbra; Nélio José de Andrade; Renato Souza Cruz; Eber Antonio Alves Medeiros Journal: Compr Rev Food Sci Food Saf Date: 2012-02-29 Impact factor: 12.811
Authors: Fabiana M C Carvalho; Vanessa C O Lima; Izael S Costa; Amanda F Medeiros; Alexandre C Serquiz; Maíra C J S Lima; Raphael P Serquiz; Bruna L L Maciel; Adriana F Uchôa; Elizeu A Santos; Ana H A Morais Journal: Nutrients Date: 2016-09-27 Impact factor: 5.717
Authors: Amanda Maria de Souza Nascimento; Victor Hugo de Oliveira Segundo; Ana Júlia Felipe Camelo Aguiar; Grasiela Piuvezam; Thaís Souza Passos; Karla Suzanne Florentino da Silva Florentino da Silva Chaves Damasceno; Ana Heloneida de Araújo Morais Journal: J Enzyme Inhib Med Chem Date: 2022-12 Impact factor: 5.051