Design of a specific peptide against phenolic glycolipid-1 from Mycobacterium leprae and its implications in leprosy bacilli entry

Abstract

BACKGROUND  Mycobacterium leprae, the causative agent of Hansen’s disease, causes neural damage through the specific interaction between the external phenolic glycolipid-1 (PGL-1) and laminin subunit alpha-2 (LAMA2) from Schwann cells.

OBJECTIVE  To design a LAMA2-based peptide that targets PGL-1 from M. leprae.

METHODS  We retrieved the protein sequence of human LAMA2 and designed a specific peptide using the Antimicrobial Peptide Database and physicochemical parameters for antimycobacterial peptide-lipid interactions. We used the AlphaFold2 server to predict its three-dimensional structure, AUTODOCK-VINA for docking, and GROMACS programs for molecular dynamics simulations.

FINDINGS  We analysed 52 candidate peptides from LAMA2, and subsequent screening resulted in a single 60-mer peptide. The mapped peptide comprises four β-sheets and a random coiled region. This peptide exhibits a 45% hydrophobic ratio, in which one-third covers the same surface. Molecular dynamics simulations show that our predicted peptide is stable in aqueous solution and remains stable upon interaction with PGL-1 binding. In addition, we found that PGL-1 has a preference for one of the two faces of the predicted peptide, which could act as the preferential binding site of PGL-1.

MAIN CONCLUSIONS  Our LAMA2-based peptide targeting PGL-1 might have the potential to specifically block this key molecule, suggesting that the preferential region of the peptide is involved in the initial contact during the attachment of leprosy bacilli to Schwann cells.

Key words: leprosy; Mycobacterium leprae; PGL-1; antimicrobial peptide; drug design

Hansen’s disease (leprosy) is an ancient infection that remains a significant health impairment in susceptible populations and is still endemic in several countries, such as Brazil, India, and Colombia.¹ The long-term vision of the World Health Organization is to eradicate leprosy by 2030. The strategy of controlling leprosy still must be reinforced with new diagnostic tools in combination with improved therapeutic regimens.² Moreover, the risk of drug resistance remains a latent threat; thus, cautious surveillance is necessary for preventing the spread of drug-resistant strains.³ Innovative therapies have been proposed as strategies to combat infection and antibiotic resistance by targeting pivotal bacterial processes, such as adhesion, cell wall permeability, quorum sensing, virulence regulons, and toxin production.⁴ In mycobacteria, this approach has been explored only in Mycobacterium tuberculosis and Mycobacterium marinum, and the approach targeted the PhoPR regulon, SapM, and ESX-1 secretion system.^4,5,6 Some promising candidates are in the preclinical stages and are being tested in animal models.

Hansen’s disease is characterised by loss of sensitivity at the peripheral nerve level due to irreversible tissue damage and subsequent weakening by the infection chronicity.⁷ This process begins with the invasion of the causative agent, Mycobacterium leprae, through a specific interaction of phenolic glycolipid-1 (PGL-1) with human laminin subunit alpha-2 (LAMA2) to promote the attachment of mycobacterial to the basal lamina of Schwann cells and pathogen internalisation.⁸

The PGL-1 molecule is surface exposed in the mycobacterial cell wall and capsule, and its structure is composed of trisaccharide units, which are defined as methyl-rhamnose derivatives bound to a phenyl group, a mycocerosic acid, and a phthiocerol region.⁹ Since PGL-1 is a well-known diagnostic marker for Hansen’s disease, this molecule has been useful for the specific differentiation of M. leprae from other mycobacteria or even for the quantification of bacterial loads to monitor the treatment outcome during multibacillary infection.¹⁰ Furthermore, the PGL-1 molecule has been reported to induce a proinflammatory response and nerve damage in patients by inducing the activation of nitric oxide synthase in infected macrophages.¹¹

PGL-1 binds specifically to the laminin multiprotein complex of the axon and is among the first steps during the Schwann cell interaction.¹² LAMA2 is involved in Schwann cell differentiation and is a key component that mediates cell-surface interaction, migration, and assembly into tissues through the promotion of laminin connections with other extracellular matrix components.¹³ We hypothesised that dissecting the LAMA2 subunits into peptides could target the specific region that binds PGL-1 and provide a further application for therapeutic or diagnostic purposes. In this study, we addressed this strategy to design in silico a LAMA2-specific peptide that targets the PGL-1 molecule from M. leprae.

MATERIALS AND METHODS

Peptide design and parameters - We retrieved the LAMA2 sequence from the UniProt database (accession code P24043) and evaluated peptide properties by using the antimicrobial peptide calculator implemented in the Antimicrobial Peptide Database (APD, https://aps.unmc.edu/home).¹⁴

We defined screening parameters based on peptides that were 60 residues long and included properties that were expected to support lipid binding, such as the following: hydrophobic ratio percentage, total net charge, GRAVY (grand average hydropathy value of the peptide), Wimley-White whole-residue hydrophobicity of the peptide, protein-binding potential (Boman index) and the total hydrophobic residues on the same surface. The 60-residue peptide was designed to preserve the functional regions of LAMA2, ensuring full PGL-1 coating within an exposed protein area; preferentially, the peptide was without disulfide bonds or any posttranslational modification and had a long peptide size to reduce the occurrence of alternative biological activity due to its length. We preferred to avoid bulky carbohydrate modifications since they might hinder contact with extracellular ligand molecules from the cell surface.

Since the expected PGL-1 molecule displayed a low solubility in aqueous solutions, we selected the peptide based on the best hydrophobic scores, and the key criteria was that peptides with positive values, compared to those with negative values, are more hydrophobic and thereby less soluble.¹⁵ Another parameter was the Wimley-White whole residue hydrophobicity; more negative values for peptides indicate a higher hydrophobicity.¹⁶ Other properties, such as the highest hydrophobic ratio percentage and the maximum number of hydrophobic residues on the same surface, were considered key for peptide selection.

LAMA2 was scanned for domain and functional motifs in SMART and visualised in the DOG program.^17,18 The DISULFIND server (http://disulfind.disi.unitn.it/) was used to predict the disulfide bridges between cysteines and their connectivity pattern.¹⁹ A prediction for N-glycosylation was performed in NetNGlyc 1.0 (https://services.healthtech.dtu.dk/service.php?NetNGlyc-1.0) based on the consensus sequence Asn-Xaa-Ser/Thr.²⁰

Prediction of possible biological properties - The peptide with the best score was checked for allergenicity using AlgPred 2.0 (https://webs.iiitd.edu.in/raghava/algpred2/index.html),²¹ toxicity in ToxinPred (https://webs.iiitd.edu.in/raghava/toxinpred/algo.php),²² and hemolytic activity by HemoPred (http://codes.bio/hemopred/).²³

Three-dimensional structure of the peptide - The peptide was mapped in the LAMA2 protein and modeled by using I-TASSER²⁴ and AlphaFold2 (https://colab.research.google.com/github/sokrypton/ColabFold/blob/main/AlphaFold2.ipynb).²⁵ The AlphaFold2 method as implemented in Google ColabFold was used as suggested by Mirdita et al.,²⁶ which differs from the original implementation from Deepmind²⁷ by replacing the homology detection of AlphaFold2 with MMseqs2 (many-against-many sequence searching).²⁸ This 3D model was subjected to a minimisation stage using the GROMACS package²⁹ for 50000 steps using a steep descent algorithm, with a maximal force tolerance of 1000 kJ mol^-1 nm^-1. The peptide was optimised using the amber99sb-ildn force field³⁰ and solvated in a dodecahedron box using the TIP3P water model.³¹ Na⁺ and Cl^- ions were added to neutralise the system’s charges and to reach a NaCl physiological concentration of 0.15 The stereochemical quality of the model before and after the minimisation stage was inspected using the Ramachandran plot, which was obtained using the Molprobity web tool (http://molprobity.biochem.duke.edu/).³²

Molecular dynamics simulation of PGL-1 binding to the peptide - We carried out molecular dynamics simulations with the aim of inspecting the modes of interaction between the proposed peptide and the M. leprae trisaccharide PGL-1. The peptide-PGL-1 complex was obtained from molecular docking calculations using AUTODOCK VINA software;³³ for this, the structure of PGL-1 was retrieved from the PubChem database (CID: 45480571). The search box was configured in such a way that it covered the entire surface of the peptide. The built complex with the best pose predicted by AUTODOCK VINA was subjected to MD simulations with the GROMACS 2019 package.³⁴ The amber ff99sb-ILDN force field and the TIP3P model were used to represent the behavior of protein in water as a solvent. PGL-1 was parameterised using the ACPYPE web server (https://www.bio2byte.be/acpype) to obtain ligand parameters for GROMACS.³⁵ The complex was neutralised with Na⁺ and Cl^- ions, brought to a concentration of 0.15 M NaCl and then subjected to a potential energy minimisation step for 50,000 steps (similar to that used for peptide minimisation), followed by two equilibration steps, including one NVT (constant volume and temperature) and a series of NPT (constant pressure and temperature) equilibrations, which were carried out for 250 ps using position restrictions on all heavy atoms. Finally, a production stage of 1000 ns (1 μs) was carried out, with a temperature of 310 K, which was controlled with the V-rescale thermostat, and a 1 bar pressure, which was controlled with the Parrinello-Rahman barostat. A time step of 2 fs was used. As a control, the peptide in water was also simulated following the same procedure. All visualisations were created with Chimera UCSF.³⁶

Prediction of dimer/PGL-1 interactions - Using the minimised structure of the peptide, peptide association and oligomerisation were calculated by an ab initio strategy with the Galaxy-Homomer server (http://galaxy.seoklab.org/index.html).³⁷ This server calculates the interface area (Ų) between predetermined chains (by user) Molprobity score³³ and a docking score, in which high values determine a greater probability of peptide interaction and the model quality, respectively. As a preliminary method in which the predicted homodimer structure was used, we carried out molecular docking of PGL-1 following the same procedure used for the single peptide. Peptide oligomer interactions and homodimer/PGL-1 interactions were assessed with Ligplot+ software.³⁸

RESULTS

Structural features of the laminin subunit alpha 2 - To understand the structural features of the LAMA2 protein, domain mapping and functional motif analysis were performed along with analysis of the whole protein. Our prediction found 28 sites of N-glycosylation in LAMA2, which agrees with the functional annotation in the UniProt database. LAMA2 analysis in the SMART tool allowed us to identify the modular composition of four domains, including a single LamNT domain and the modular arrangement of EGF-Lam, laminin B (LamB), and laminin G (LamG) domains (Fig. 1). The LamNT domain is located between residues 33-285 through the N-terminus for protein insertion in the cell membrane. In contrast, LAMA2 contains 16 EGF-laminin domains, which are characterised by the presence of many cysteine residues that form disulfide bonds. We predicted 66 disulfide bonds that were distributed across the whole protein (data not shown). Furthermore, our analysis showed two LamB domains located between residues 578-710 and 1229-1364 that are interspaced by a set of EGF laminin domains. The C-termini exhibit an arrangement of five LamG domains at the C-terminus of the LAMA2 protein.

Fig. 1:
modular structure of laminin subunit alpha-2 from Homo sapiens (Human). N-linked glycosylation (GlcNAc) labeled asparagine (Asn) positions 55, 89, 303, 363, 380, 470, 746, 1061, 1597, 1614, 1700, 1810, 1901, 1916, 1920, 2017, 2028, 2045, 2126, 2240, 2360, 2435, 2478, 2551, 2558, 2648, 2868, and 2893 as predicted in NetNGlyc 1.0. The peptide position is shown in the red line above the first LamG domain.

Peptide selection and pharmacological properties - LAMA2 contained 3,122 amino acids, and our computational screening resulted in 52 candidate peptides. However, based on our computer-aided peptide design strategy and physicochemical criteria, a 60-mer peptide was designed with a predominance of hydrophobic residues that are solvent-exposed to ensure a potential hydrophobic interaction with PGL-1 (Table). Thus, the calculated percentages of the most frequent amino acids were Val ratio = 13%, Lys = 10%, Gly = 10%, Leu = 8%, and 7% for Ile, Phe, Ala, Ser, and Asn residues, respectively. This designed peptide showed a hydrophobic ratio of 45%, in which 22 hydrophobic residues (37%) were located on the same surface. Other calculated parameters included 0.18 as a GRAVY value, a Wimley-White whole-residue hydrophobicity of the peptide of 7.94, and a protein-binding potential of 0.78 (Table). The peptide mapping indicated that our peptide was located between the 2161-2220 positions of the LAMA2 protein, specifically in the first LamG domain (positioned between 2166-2311), and seemed to be surface exposed in the three-dimensional structure.

Other pharmaceutical properties of our peptide indicated that it was nonallergenic (score: 0.31), nontoxic (-1.00), and nonhemolytic.

TABLE
List of predicted peptides from the human LAMA2 protein

Region	Hydrophobic ratio (%)	Total net charge	GRAVY	Wimley-White whole-residue hydrophobicity	Protein binding Potential (Kcal/mol)	Total hydrophobic residues on the same surface
1-60	43	+2.2	0.09	6.84	0.73	14
61-120	27	+3.5	-1.09	16.95	2.45	7
121-180	40	-2.5	-0.16	6.19	1.18	15
181-240	28	-3.5	-0.58	15.76	1.95	8
241-300	42	+4.5	-0.22	9.70	2.11	9
301-360	33	-3.2	-0.93	20.63	2.57	0
361-420	32	+4.2	-0.47	9.78	1.93	0
421-480	32	+2.5	-0.66	16.97	2.31	ND
481-540	35	-2.0	-0.60	10.36	2.01	7
541-600	33	-0.7	-0.51	7.97	1.44	2
601-660	35	-13.5	-0.37	24.8	1.90	13
661-720	45	+2.5	0.26	5.8	1.06	21
721-780	35	-3.25	0.40	14.49	1.87	ND
781-840	32	-2.7	-0.28	8.34	1.08	ND
841-900	32	-2.0	-0.34	10.81	1.57	2
901-960	37	+0.25	0.55	14.78	2.17	9
961-1020	37	-0.5	0.45	10.53	1.73	2
1021-1080	33	+2.2	-0.49	10.91	1.77	ND
1081-1140	33	+1.5	-0.63	12.98	2.23	2
1141-1200	38	+0.25	-0.21	12.42	1.48	ND
1201-1260	40	0	-0.35	10.71	1.00	9
1261-1320	32	+0.75	-0.66	17.98	2.27	7
1321-1380	37	-1.5	-0.32	15.36	2.04	13
1381-1440	32	-0.75	-0.29	8.23	1.40	4
1441-1500	38	-1.5	-0.22	9.06	1.55	ND
1501-1560	27	-3	-0.60	14.06	1.66	ND
1561-1620	43	-2.5	-0.15	9.86	0.92	18
1621-1680	32	-1.75	-0.75	24.7	2.62	14
1681-1740	37	-2	-0.91	33.88	3.02	17
1741-1800	37	-5	-1.04	36.57	3.1	18
1801-1860	35	-4	-0.83	30.02	2.81	16
1861-1920	35	-5.75	-0.70	27.67	2.41	17
1921-1980	40	+1.25	-0.55	23.59	1.92	20
1981-2040	35	+1.25	-0.85	25.72	2.69	17
2041-2100	38	+2.25	-0.64	23.44	2.2	18
2101-2160	32	+2	-0.85	27.87	2.65	14
2161-2220	45	+2	0.18	7.94	0.78	22
2221-2280	35	-0.5	-0.17	7.69	1.64	11
2281-2340	33	-1	-0.44	17.35	1.64	9
2341-2400	43	+1.25	0.07	3.61	1.4	15
2401-2460	30	-1.75	-0.62	12.57	2.16	12
2461-2520	32	+3	-0.45	13.76	1.74	6
2521-2580	32	+2	-0.23	13.95	1.49	7
2581-2640	35	+4	-0.42	15.12	2.42	12
2641-2700	40	0	-0.25	14.14	1.37	13
2701-2760	38	-4.5	-0.31	22.51	1.29	ND
2761-2820	42	+4.5	-0.21	13.71	2.02	15
2821-2880	30	+1.5	-0.58	11.66	1.77	11
2881-2940	40	+0.5	0.04	7.53	1.13	ND
2941-3000	42	-2.7	0.24	7.55	0.56	14
3001-3060	33	-1.2	-0.58	16.92	1.92	8
3061-3120	39	+3	-0.07	11.44	1.12	ND

• ND: nondetermined.

Peptide 3D structure - The peptide modeling was based on artificial intelligence through Alphafold. The predicted structure exhibited a high structural quality according to its Ramachandra plot, which indicates that 96.6% of the residues have a favorable stereochemistry (Fig. 2A). In addition, AlphaFold provided two confidence indicators to determine the reliability of the results. The first indicator, called pLDDT (predicted lDDT-Cα), provided a measurement of the local confidence (for each residue) on a scale from 0 to 100. The results presented in Fig. 2B show that most of the residues had a plDDT > 80, which corresponded to confidence ranging between high and very high. The second indicator, called PAE (Predicted Aligned Error), represented the expected error associated with the relative positions of the different domains of the protein. The error values calculated for the peptide were consistently low except for the first two and last two residues (data not shown), thus indicating good confidence in the positions of the beta-strands. Even so, after a minimisation step, a structure with improved structural quality was obtained based on its Ramachandran plot, which indicates that all (100%) of the amino acids have favorable stereochemistry (Fig. 2C). An overlay of the peptide structure before and after minimisation is shown in Fig. 2D. In general, the peptide exhibited a random coiled region in the N-termini followed by four antiparallel β-strands. Our peptide lacked Cys residues, and its three-dimensional structure resembled a defensin-like beta structure without N-glycosylation sites. Additionally, the sequence of our peptide has 89.92% sequential identity, which corresponds to the crystal structure of the LG1-3 region of LAMA2 (PDB 1QU0) from Mus musculus.

Fig. 2:
results for the Ramachandran plot for the structure of the selected peptide before (top-left side) and after energy minimisation (top-right side). On the bottom side, the superposition of the structures before (green) and after (blue) the minimisation are shown. The N- and C-termini are shown schematically. The peptide is shown as ribbons with all its residues as sticks.

Docking and molecular dynamics analysis - The peptide-ligand complex was obtained by molecular docking calculations. Given the stochastic nature of AutoDock Vina’s search algorithm, we performed a triplicate run. All 27 conformations (nine for each run) were located in the same hydrophobic region of the chosen peptide, as shown on the left side of Fig. 3A. The foregoing suggests that the PGL-1 ligand would have a greater affinity or a greater preference to interact with this face of the peptide (from now on, this face will be called the front face of the peptide) than with the back face, which has a more hydrophilic character (Fig. 3A-right). Fig. 3B shows the best pose obtained for each run, and the best pose had a binding affinity value of -5.1 kcal/mol ± 0.0. One of these conformations was used as the initial coordinates of the peptide/PGL-1 complex for the MD simulations.

Fig. 3:
complete results of molecular docking triplicate (A-left) and the back side of the peptide (A-right). The first pose of each molecular docking run is shown in blue, pink and green (B). The protein is shown as a surface, and the residues are coloured according to their hydrophobicity, with blue being very hydrophilic and red being very hydrophobic. PGL-1 is shown in sticks.

In our MD simulations, we first inspected the structural stability of the peptide (Fig. 4). The RMSD results (Fig. 4A) suggest that the peptide without PGL-1 was more flexible during the simulation, or the peptide interaction with PGL-1 favors conformational changes. These greater changes in the RMSD of the peptide without PGL-1 are associated with greater flexibility in the N-termini coil region, as shown in Fig. 4B. The fluctuations (RMSF) of the loop regions in the peptide without PGL-1 were considerably greater than those in the peptide/PGL-1 complex (Fig. 4B-C). Additionally, we tracked the secondary structure of the peptide (Fig. 4D), suggesting that the main structure of four antiparallel beta sheets is highly stable since after 1000 ns of simulation, the beta sheets were preserved in the presence and absence of PGL-1. It was also observed that after approximately 650 ns, a new beta sheet was formed at the N-terminal end of the peptide by PGL-1 binding (Fig. 4D), which may be related to the fact that the RMSD of the peptide in complex with PGL-1 stabilises at ~0.8 nm.

To map the modes of interaction between PGL-1 and the selected peptide, we calculated the contact frequencies between them by defining each contact within a range of 3 Å in the interaction.

Fig. 4:
(A) RMSD of the backbone and (B) RMSF of the protein residues for the simulations with (green line) and without PGL-1 (purple line). (C) Peptide B-factor ratio. Blue and red represent the increase in movement for the residues in the peptide with and without PGL-1, respectively. Similarly, a thicker ribbon indicates greater flexibility. (D) Secondary structure as a function of time of the protein with and without PGL-1. α-helix, beta sheets and disordered regions are shown in orange, green and white, respectively.

As shown in Fig. 5, the highest interaction frequencies (from 40% to 50%), which define a hydrophobic pocket that is delineated by residues Tyr7, Val11, Phe24, Ile37 and Phe46. Tyr7 and Val11, are located in the loop of the N-terminal end, and the other three residues Phe24, Ile37 and Phe46 are part of the peptide front face, and each is located in a different beta sheet. Additionally, several residues, such as Leu22, Tyr25, Ala36, Glu38, Ser45 and Leu47, present interactions of less than 20%. Those residues are located on the rear face of the peptide in the preferential face of PGL-1 to interact with the peptide front face.

Fig. 5:
interaction frequency in the peptide/PGL-1 complex (left). Schematically, the last conformation of the simulation is shown, highlighting the main interactions between the peptide and PGL-1 (right). The N- and C-termini are shown schematically. The peptide is illustrated with ribbons, while the main residues and PGL-1 are represented by sticks.

Homodimer interactions of LAMA2-derived peptide with PGL-1 - Oligomer modeling results were used to calculate the potential peptide-peptide interactions. Oligomer calculations were based on the structure of a laminin G-like module of LAMA2 (peptide coverage between 2157-2216), leading to a homodimer of 7410.55 Ų (Fig. 6C). Peptide interacting chains were analysed in Ligplot+ software to calculate dimer interactions. Interacting chains might be associated with 26 nonbonded contacts that involve neutral, aliphatic, aromatic, and positively charged residues (Fig. 6A). The molecular docking result for PGL-1 in the homodimer docks at a site with similar hydrophobicity to that of the peptide alone brought the mycobacterial ligand together (Fig. 6B-C). We inspected the interactions, which included hydrophobic interactions with residues such as Asn13, Ile37, Met39, Val44, Phe46, and Tyr59, all of which were in chain A, and hydrogen bond interactions with Val12 and Gly42 in the A chain and with Ser52 in chain B (Fig. 6D).

Fig. 6:
(A) 2D interactions between peptide chains in the homodimer. Structure of the complex formed by the homodimer and PGL-1 using ribbons (B) for the homodimer and hydrophobic surface (in which blue indicates the most hydrophilic residues, red indicates the most hydrophobic residues) (C); in both cases, sticks are used for PGL-1. (D) 2D interaction between the homodimer and PGL-1. (A and B) Chains A and B are shown in red and pink, respectively. (A and D) Hydrogen bonds are shown as green dashed lines, and hydrophobic interactions are shown as half circles with lines through them.

DISCUSSION

Since ancient times, Hansen’s disease has been a public health problem worldwide. Despite the many efforts to eradicate this disease, concerns over concomitants, such as poverty, delayed diagnosis, and drug resistance, have emerged.^2,39,40 The PGL-1 molecule has significantly contributed to the serodiagnosis of Hansen’s disease, providing a specific target to identify M. leprae.¹⁰ Currently, the PGL-1-based ELISA test still contributes to determining whether a patient is free from leprosy bacilli.⁴¹ In addition to its importance in the process of infection, the PGL-1 saccharide fraction is species-specific for M. leprae, constituting a highly immunogenic molecule that is quickly recognised by the immune system.⁴² Thus, the conjugate composition of PGL-1 has been associated with its neural tropism from M. leprae.⁸

The domains present in the LAMA2 protein, including a set of EGF-laminin domains, are essential for signal transduction and protein dimerization.¹³ The LN domain is involved in basement membrane assembly, and its role might be a cooperative process in which laminins polymerise through their N-terminal domain (VI) and anchor to the cell surface. Finally, LamG was associated with cellular adhesion to laminins and was mediated by a repetitive region of five laminin G-like (LG) domains.^43,44 Consistently, our peptide was designed on the basis of the C-terminal of the LAMA2 chain, which contains five LamG domains, in which PGL-1 seems to contact LAMA2 during the mycobacterial invasion of Schwann cells.^8,45

Our peptide contained a portion of the LamG domain, and further analysis showed that it is compatible with the steroid-binding site of related proteins containing laminin G-like domains, including those with sexual hormone-binding sites that resemble lipid-related ligands.⁴⁶ In this study, we proposed that this site might anchor PGL-1 to the cell surface through LAMA2, representing the initial bacilli-specific interaction. Consistently, a study has shown that M. leprae is strongly bound to the LAMA2 C-terminal but not to the N-terminal region in the proximal G1-G3 subdomains.^12,47

In this paper, we attempted to identify the region from human LAMA2 that was necessary for triggering the attachment of M. leprae to Schwann cells. Our computational peptide design supported the 60-residue simulation that might be sufficient to capture peptide folding.⁴⁸ The peptide encompassing residues 2161-2220 of LAMA2 might be able to bind flexibly to PGL-1, allowing bacterial attachment and subsequent pathogenesis. In our MD simulations, the peptide conformation was more stable in the presence of PGL-1 than in the absence, which hypothetically seems to be an early step in bacilli membrane attachment. Additionally, our simulations suggest that PGL-1 prefers to interact with one of the two faces of the peptide, i.e., the face called the front face in this study. For example, PGL-1 might be positioned by a network of van der Waals and hydrophobic interactions in the loop region and residues located in three of the four main beta-sheets. The docking and simulation results also supported that the binding of PGL-1 to the LAMA2-derived peptide is flexible on its N-terminal region (ΔG -5.1 Kcal/mol) and, thereby, might increase the affinity between the ligand and the peptide or even homodimers. Our peptide-aided design in the selection of peptide candidates promotes hydrophobic interactions that might define the PGL-1 binding site and proper adjustment to the basal laminin of Schwann cells. Similarly, hydrophobic interactions contribute strongly to steroid pocket binding and fine-tuned interactions with hydrophobic ligands in proteins containing laminin-like domains.^46,49 Our molecular dynamic results suggested that a flexible loop region assisted as the gate for PGL-1 (through interactions with residues, such as Tyr7 and Val11), as reported similarly in sex steroid hormones, resembling the capability to bind lipidic ligands with their loop segment for ligand-specific rearrangement.⁴⁹

Protein interactions with cellular membranes have been thoroughly studied as computational models in antimicrobial peptides.^50,51 However, we performed peptide-glycolipid molecular modeling predictions as a novel approach for antimicrobial peptide design against a key conjugate molecule from leprosy bacilli. Our calculations included the possibility of peptide homodimerization and homodimer interaction with PGL-1. Although the predicted affinity for the homodimer was lower (4.8 kcal/mol) than that for the peptide alone (5.1 kcal/mol), it is necessary to mention that this bond is flexible; that is, this affinity can vary over time and, as expected, depends on the conformational changes that occur in the receptor. Interaction percentages less than 40% and a variety of structural conformations were reported for the organic compounds that interacted with highly flexible proteins.⁵² Given that we found higher percentages of interaction and a high conformational stability, our results are promising.

The LAMA2-based peptide might interact with membranes and has a chance to be an antimicrobial peptide in which aligning is performed to find the most similar peptides in the APD database. This peptide showed a 32% similarity with halocin-like peptides (halocin S8 and halocin R1), which have activity against Gram-positive and Gram-negative bacteria and a similar hydrophobic ratio.⁵³ The results showed that peptides with a low similarity (below 30%) presented a shared hydrophobicity percentage with PGL-1 based on the sequence comparison, which was performed to find antimicrobial peptides that most resembled our input peptide sequence in the APD. Consistently, previous approaches with host-based peptides were tested successfully to control mycobacterial growth, and hydrophobicity was determined to be a key parameter for enhancing mycobactericidal activity and selectivity.^54,55,56,57

Our findings might contribute to deciphering the first step of how M. leprae establishes initial host contact for manipulating signaling pathways leading to axonal damage and hindering myelin maintenance in the basal lamina.⁵⁸ Similar mechanisms have been reported involving other host receptors that contain laminin-like domains for the entry of viral pathogens.^59,60 For example, the virus that causes Lassa fever targets Schwann cells, selectively interfering with the myelination process through its viral receptor dystroglycan, leading to neurological disorders.⁶¹

Our study might provide insights into how mycobacterial glycolipids interact with host laminin and suggests a new strategy for exploring the development of new diagnostic or therapeutic options based on the druggable proteome from the leprosy bacillus.^62,63,64

In conclusion - We identified N-glycosylation sites, disulfide bridges, and domains along human LAMA2. We dissected the LAMA2 sequence into 52 peptides, including the potential PGL-1 binding site from M. leprae. Our peptide targeting PGL-1 is located between residues 2161-2220, and its structure exhibits a combination of β-sheets and random coiled region that might flexibly bind PGL-1. Thus, our approach with a specific peptide could block the interaction of M. leprae with the host cell, thereby preventing long therapeutic regimens, disease chronicity, and possibly nerve damage in Hansen’s disease patients.

REFERENCES

• ¹ WHO - World Health Organization. Towards zero leprosy. Global leprosy (Hansen's disease) strategy 2021-2030. WHO team, control of neglected tropical diseases. 2017. Available from: https://www.who.int/publications/i/item/9789290228509
  » https://www.who.int/publications/i/item/9789290228509
• ² Chaves LL, Patriota Y, Soares-Sobrinho JL, Vieira ACC, Lima SAC, Reis S. Drug delivery systems on leprosy therapy: moving towards eradication? Pharmaceutics. 2020; 12: 1202.
• ³ Smith CS, Aerts A, Saunderson P, Kawuma J, Kita E, Virmond M. Multidrug therapy for leprosy: a game changer on the path to elimination. Lancet Infect Dis. 2017; 17: e293-7.
• ⁴ Dickey SW, Cheung GYC, Otto M. Different drugs for bad bugs: antivirulence strategies in the age of antibiotic resistance. Nat Rev Drug Discov. 2017; 16: 457-71.
• ⁵ Chen JM, Pojer F, Blasco B, Cole ST. Towards anti-virulence drugs targeting ESX-1 mediated pathogenesis of Mycobacterium tuberculosis. Drug Discov Today Dis Mech. 2010; 7: e25-e31.
• ⁶ Gries R, Sala C, Rybniker J. Host-directed therapies and anti-virulence compounds to address anti-microbial resistant tuberculosis infection. Appl Sci. 2020; 10: 2688.
• ⁷ Maymone MBC, Venkatesh S, Laughter M, Abdat R, Hugh J, Dacso MM, et al. Leprosy: treatment and management of complications. J Am Acad Dermatol. 2020; 83: 17-30.
• ⁸ Ng V, Zanazzi G, Timpl R, Talts JF, Salzer JL, Brennan PJ, et al. Role of the cell wall phenolic glycolipid-1 in the peripheral nerve predilection of Mycobacterium leprae. Cell. 2000; 103: 511-24.
• ⁹ Gautam S, Sharma D, Goel A, Patil SA, Bisht D. Insights into Mycobacterium leprae proteomics and biomarkers-an overview. Proteomes. 2021; 9(1): 7.
• ¹⁰ Nath I, Saini C, Valluri VL. Immunology of leprosy and diagnostic challenges. Clin Dermatol. 2015; 33: 90-8.
• ¹¹ Oldenburg R, Demangel C. Pathogenic and immunosuppressive properties of mycobacterial phenolic glycolipids. Biochimie. 2017; 141: 3-8.
• ¹² Rambukkana A, Salzer JL, Yurchenco PD, Tuomanen EI. Neural targeting of Mycobacterium leprae mediated by the G domain of the laminin-alpha2 chain. Cell. 1997; 88: 811-21.
• ¹³ Sasaki T, Fässler R, Hohenester E. Laminin: the crux of basement membrane assembly. J Cell Biol. 2004; 164: 959-63.
• ¹⁴ Wang G, Li X, Wang Z. APD3: the antimicrobial peptide database as a tool for research and education. Nucleic Acids Res. 2016; 44: 1087-93.
• ¹⁵ Kyte J, Doolittle RF. A simple method for displaying the hydropathic character of a protein. J Mol Biol. 1982; 157: 105-32.
• ¹⁶ Wimley WC, White SH. Experimentally determined hydrophobicity scale for proteins at membrane interfaces. Nat Struct Biol. 1996; 3: 842-8.
• ¹⁷ Letunic I, Bork P. 20 years of the SMART protein domain annotation resource. Nucleic Acids Res. 2018; 46: 493-6.
• ¹⁸ Ren J, Wen L, Gao X, Jin C, Xue Y, Yao X. DOG 1.0: illustrator of protein domain structures. Cell Res. 2009; 19: 271-3.
• ¹⁹ Ceroni A, Passerini A, Vullo A, Frasconi P. DISULFIND: a disulfide bonding state and cysteine connectivity prediction server. Nucleic Acids Res. 2006; 34: W177-81.
• ²⁰ Gupta R, Brunak S. Prediction of glycosylation across the human proteome and the correlation to protein function. Pac Symp Biocomput. 2002; 310-22.
• ²¹ Sharma N, Patiyal S, Dhall A, Pande A, Arora C, Raghava GPS. AlgPred 2.0: an improved method for predicting allergenic proteins and mapping of IgE epitopes. Brief Bioinform. 2021; 22(4): bbaa294.
• ²² Gupta S, Kapoor P, Chaudhary K, Gautam A, Kumar R, Raghava GP. In silico approach for predicting toxicity of peptides and proteins. PLoS One. 2013; 8: e73957.
• ²³ Win TS, Malik AA, Prachayasittikul V, WikbergJES, Nantasenamat C, Shoombuatong W. HemoPred: a web server for predicting the hemolytic activity of peptides. Future Med Chem. 2017; 9(3): 275-91.
• ²⁴ Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008; 9: 40.
• ²⁵ Skolnick J, Gao M, Zhou H, Singh S. AlphaFold 2: Why it works and its implications for understanding the relationships of protein sequence, structure, and function. J Chem Inf Model. 2021; 61(10): 4827-31.
• ²⁶ Mirdita M, Schütze K, Moriwaki Y, Heo L, Ovchinnikov S, Steinegger M. ColabFold: making protein folding accessible to all. Nat Methods. 2022; 19(6): 679-82.
• ²⁷ Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, et al. Highly accurate protein structure prediction with AlphaFold. Nature. 2021; 596: 583-9.
• ²⁸ Mirdita M, Steinegger M, Söding J. MMseqs2 desktop and local web server app for fast, interactive sequence searches. Bioinformatics. 2019; 35: 2856-8.
• ²⁹ Rakhshani H, Dehghanian E, Rahati A. Enhanced GROMACS: toward a better numerical simulation framework. J Mol Model. 2019; 25: 355.
• ³⁰ Lindorff-Larsen K, Piana S, Palmo K, Maragakis P, Klepeis JL, Dror RO, et al. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins. 2010; 78: 1950-8.
• ³¹ Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. Comparison of simple potential functions for simulating liquid water. J Chem Phys. 1983; 79(2): 926-35.
• ³² Williams CJ, Headd JJ, Moriarty NW, Prisant MG, Videau LL, Deis LN, et al. MolProbity: more and better reference data for improved all-atom structure validation. Protein Sci. 2018; 27: 293-315.
• ³³ Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010; 31: 455-61.
• ³⁴ Abraham MJ, Murtola T, Schulz R, Páll S, Smith JC, Hess B, et al. GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX. 2015; 1: 19-25.
• ³⁵ Sousa da Silva AW, Vranken WF. ACPYPE-Antechamber python parser interface. BMC Res Notes. 2012; 5(1): 1-8.
• ³⁶ Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera - a visualization system for exploratory research and analysis. J Comput Chem. 2004; 25: 1605-12.
• ³⁷ Heo L, Lee H, Seok C. GalaxyRefineComplex: refinement of protein-protein complex model structures driven by interface repacking. Sci Rep. 2016; 6: 32153.
• ³⁸ Laskowski RA, Swindells MB. LigPlot+: multiple ligand-protein interaction diagrams for drug discovery. J Chem Inf Model. 2011; 10: 2778-86.
• ³⁹ Guerrero MI, Colorado CL, Torres JF, León CI. Is drug-resistant Mycobacterium leprae a real cause for concern?: First approach to molecular monitoring of multibacillary Colombian patients with and without previous leprosy treatment. Biomedica. 2014; 34(Suppl. 1): 137-47.
• ⁴⁰ Chavarro-Portillo B, Soto CY, Guerrero MI. Mycobacterium leprae's evolution and environmental adaptation. Acta Trop. 2019; 197: 105041.
• ⁴¹ Araujo S, Rezende MM, Sousa DC, Rosa MR, Santos DC, Goulart LR, et al. Risk-benefit assessment of Bacillus Calmette-Guérin vaccination, anti-phenolic glycolipid I serology, and Mitsuda test response: 10-year follow-up of household contacts of leprosy patients. Rev Soc Bras Med Trop. 2015; 48: 739-45.
• ⁴² Hunter SW, Brennan PJ. A novel phenolic glycolipid from Mycobacterium leprae possibly involved in immunogenicity and pathogenicity. J Bacteriol. 1981; 147: 728-35.
• ⁴³ Beckmann G, Hanke J, Bork P, Reich JG. Merging extracellular domains: fold prediction for laminin G-like and amino-terminal thrombospondin-like modules based on homology to pentraxins. J Mol Biol. 1998; 275: 725-30.
• ⁴⁴ Hohenester E, Tisi D, Talts JF, Timpl R. The crystal structure of a laminin G-like module reveals the molecular basis of a-dystroglycan binding to laminins, perlecan, and agrin. Mol Cell. 1999; 4: 783-92.
• ⁴⁵ Carafoli F, Clout NJ, Hohenester E. Crystal structure of the LG1-3 region of the laminin alpha2 chain. J Biol Chem. 2009; 284: 22786-92.
• ⁴⁶ Grishkovskaya I, Avvakumov GV, Hammond GL, Catalano MG, Muller YA. Steroid ligands bind human sex hormone-binding globulin in specific orientations and produce distinct changes in protein conformation. J Biol Chem. 2002; 277: 32086-93.
• ⁴⁷ Rambukkana A. Molecular basis for the peripheral nerve predilection of Mycobacterium leprae. Curr Opin Microbiol. 2001; 4: 21-7.
• ⁴⁸ Chen CH, Melo MC, Berglund N, Khan A, de la Fuente-Nunez C, Ulmschneider JP, et al. Understanding and modelling the interactions of peptides with membranes: from partitioning to self-assembly. Curr Opin Struct Biol. 2020; 61: 160-6.
• ⁴⁹ Round P, Das S, Wu TS, Wähälä K, Van Petegem F, Hammond GL. Molecular interactions between sex hormone-binding globulin and nonsteroidal ligands that enhance androgen activity. J Biol Chem. 2020; 295: 1202-11.
• ⁵⁰ Palmer N, Maasch J, Torres MDT, de la Fuente-Nunez C. Molecular dynamics for antimicrobial peptide discovery. Infect Immun. 2021; 89: e00703-20.
• ⁵¹ King GM, Kosztin I. Towards a quantitative understanding of protein-lipid bilayer interactions at the single molecule level: opportunities and challenges. J Membr Biol. 2021; 254: 17-28.
• ⁵² Jin F, Yu C, Lai L, Liu Z. Ligand clouds around protein clouds: a scenario of ligand binding with intrinsically disordered proteins. PLoS Comput Biol. 2013; 9(10): e1003249.
• ⁵³ Kumar V, Singh B, van Belkum MJ, Diep DB, Chikindas ML, Ermakov AM, et al. Halocins, natural antimicrobials of Archaea: exotic or special or both? Biotechnol Adv. 2021; 53: 107834.
• ⁵⁴ Noschka R, Wondany F, Kizilsavas G, Weil T, Weidinger G, Walther P, et al. Gran1: a granulysin-derived peptide with potent activity against intracellular Mycobacterium tuberculosis. Int J Mol Sci. 2021; 22(16): 8392.
• ⁵⁵ Raj S, Venugopal U, Pant G, Kalyan M, Arockiaraj J, Krishnan MY, et al. Anti-mycobacterial activity evaluation of designed peptides: cryptic and database filtering based approach. Arch Microbiol. 2021; 203: 4891-9.
• ⁵⁶ Khara JS, Lim FK, Wang Y, Ke XY, Voo ZX, Yang YY, et al. Designing a-helical peptides with enhanced synergism and selectivity against Mycobacterium smegmatis: Discerning the role of hydrophobicity and helicity. Acta Biomater. 2015; 28: 99-108.
• ⁵⁷ Pearson CS, Kloos Z, Murray B, Tabe E, Gupta M, Kwak JH, et al. Combined bioinformatic and rational design approach to develop antimicrobial peptides against Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2016; 60: 2757-64.
• ⁵⁸ Mietto BS, de Souza BJ, Rosa PS, Pessolani MCV, Lara FA, Sarno EN. Myelin breakdown favours Mycobacterium leprae survival in Schwann cells. Cell Microbiol. 2020; 22: e13128.
• ⁵⁹ Jae LT, Raaben M, Riemersma M, van Beusekom E, Blomen VA, Velds A, et al. Deciphering the glycosylome of dystroglycanopathies using haploid screens for lassa virus entry. Science. 2013; 340: 479-83.
• ⁶⁰ Kunz S, Borrow P, Oldstone MB. Receptor structure, binding, and cell entry of arenaviruses. Curr Top Microbiol Immunol. 2002; 262: 111-37.
• ⁶¹ Rambukkana A, Kunz S, Min J, Campbell KP, Oldstone MB. Targeting Schwann cells by nonlytic arenaviral infection selectively inhibits myelination. Proc Natl Acad Sci USA. 2003; 100(26): 16071-6.
• ⁶² Soares de Lima C, Zulianello L, Marques MA, Kim H, Portugal MI, Antunes SL, et al. Mapping the laminin-binding and adhesive domain of the cell surface-associated Hlp/LBP protein from Mycobacterium leprae. Microbes Infect. 2005; 7: 1097-109.
• ⁶³ Vedithi SC, Malhotra S, Acebrón-García-de-Eulate M, Matusevicius M, Torres PHM, Blundell TL. Structure-guided computational approaches to unravel druggable proteomic landscape of Mycobacterium leprae. Front Mol Biosci. 2021; 8: 663301.
• ⁶⁴ Van Dijk JHM, van Hooij A, Groot LM, Geboers J, Moretti R, Verhard-Seymonsbergen E, et al. Synthetic phenolic glycolipids for application in diagnostic tests for leprosy. Chembiochem. 2021; 22(8): 1487-93.

• Financial support: This work was financed by the Ministry of Science, Innovation and Technology with the grant “Determinación de polimorfismos genéticos asociados con susceptibilidad para enfermedad de Hansen en población del Departamento de Norte de Santander”, code 212084368694, MINCIENCIAS contract number 847-2019.

Publication Dates

• Publication in this collection
  18 July 2022
• Date of issue
  2022

History

• Received
  26 Jan 2022
• Accepted
  14 June 2022