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

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.

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

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.

Inverse Molecular Docking Study of NS3-Helicase and NS5-RNA Polymerase of Zika Virus as Possible Therapeutic Targets of Ligands Derived from Marcetia taxifolia and Its Implications to Dengue Virus.

Dengue and Zika are two mosquito-borne diseases of great impact on public health around the world in tropical and subtropical countries. DENV and ZIKV belong to the Flaviviridae family and the Flavivirus genus. Currently, there are no effective therapeutic agents to treat or prevent these pathologies. The main objective of this work was to evaluate potential inhibitors from active compounds obtained from Marcetia taxifolia by performing inverse molecular docking on ZIKV-NS3-helicase and ZIKV-NS5-RNA polymerase as targets. This computational strategy is based on renormalizing the binding scores of the compounds to these two proteins, allowing a direct comparison of the results across the proteins. The crystallographic structures of the ZIKV-NS3-helicase and ZIKV-NS5-RNA-polymerase proteins share a great similarity with DENV homologous proteins. The P-loop active site of the crystallographic structure of ZIKV-NS3-helicase presents a high percentage of homology with the four dengue serotypes. It was found that most ligands of the active compounds (5,3'-dihydroxy-3,6,7,8,4'-pentamethoxyflavone (5DP); 5-hydroxy-3,6,7,8,3',4'-hexamethoxyflavone (5HH); myricetin-3-O-rhamnoside (M3OR)) from Marcetia taxifolia had a better affinity for ZIKV-NS3-helicase than for ZIKV-NS5-RNA polymerase, as indicated by the negative multiple active site correction (MASC) score, except for M3RG that showed a higher affinity for ZIKV-NS5-RNA polymerase. On the other hand, the AutoDock Vina scores showed that M3OR had the highest score value (-9.60 kcal/mol) and the highest normalized score (1.13) against ZIKV-NS3-helicase. These results in silico demonstrated that the nonstructural proteins NS3-helicase and NS5-RNA polymerase, which share similar molecular structures between the selected viruses, could become therapeutic targets for some bioactive compounds derived from Marcetia taxifolia.

Functional Heterologous Expression of Mature Lipase LipA from Pseudomonas aeruginosa PSA01 in Escherichia coli SHuffle and BL21 (DE3): Effect of the Expression Host on Thermal Stability and Solvent Tolerance of the Enzyme Produced.

This study aimed to express heterologously the lipase LipA from Pseudomonas aeruginosa PSA01 obtained from palm fruit residues. In previous approaches, LipA was expressed in Escherichia coli fused with its signal peptide and without its disulfide bond, displaying low activity. We cloned the mature LipA with its truncated chaperone Lif in a dual plasmid and overexpressed the enzyme in two E. coli strains: the traditional BL21 (DE3) and the SHuffle® strain, engineered to produce stable cytoplasmic disulfide bonds. We evaluated the effect of the disulfide bond on LipA stability using molecular dynamics. We expressed LipA successfully under isopropyl ß-d-1-thio-galactopyranoside (IPTG) and slow autoinducing conditions. The SHuffle LipA showed higher residual activity at 45 °C and a greater hyperactivation after incubation with ethanol than the enzyme produced by E. coli BL21 (DE3). Conversely, the latter was slightly more stable in methanol 50% and 60% (t½: 49.5 min and 9 min) than the SHuffle LipA (t½: 31.5 min and 7.4 min). The molecular dynamics simulations showed that removing the disulfide bond caused some regions of LipA to become less flexible and some others to become more flexible, significantly affecting the closing lid and partially exposing the active site at all times.

Exploring the local elastic properties of bilayer membranes using molecular dynamics simulations.

Membrane mechanical elastic properties regulate a variety of cellular processes involving local membrane deformation, such as ion channel function and vesicle fusion. In this work, we used molecular dynamics simulations to estimate the local elastic properties of a membrane. For this, we calculated the energy needed to extract a DOPE lipid molecule, modified with a linker chain, from a POPC bilayer membrane using the umbrella sampling technique. Although the extraction energy entails several contributions related not only to elastic deformation but also to solvation, careful analysis of the potential of mean force (PMF) allowed us to dissect the elastic contribution. With this information, we calculated an effective linear spring constant of 44 ± 4 kJ·nm(-2)·mol(-1) for the DOPC membrane, in agreement with experimental estimates. The membrane deformation profile was determined independently during the stretching process in molecular detail, allowing us to fit this profile to a previously proposed continuum elastic model. Through this approach, we calculated an effective membrane spring constant of 42 kJ·nm(-2)·mol(-1), which is in good agreement with the PMF calculation. Furthermore, the solvation energy we derived from the data is shown to match the solvation energy estimated from critical micelle formation constants. This methodology can be used to determine how changes in lipid composition or the presence of membrane modifiers can affect the elastic properties of a membrane at a local level.

Binding ensemble profiling with photoaffinity labeling (BEProFL) approach: mapping the binding poses of HDAC8 inhibitors.

A binding ensemble profiling with (f)photoaffinity labeling (BEProFL) approach that utilizes photolabeling of HDAC8 with a probe containing a UV-activated aromatic azide, mapping of the covalent modifications by liquid chromatography-tandem mass spectrometry, and a computational method to characterize the multiple binding poses of the probe is described. By use of the BEProFL approach, two distinct binding poses of the HDAC8 probe were identified. The data also suggest that an "upside-down" pose with the surface binding group of the probe bound in an alternative pocket near the catalytic site may contribute to the binding.

Parameterization of aromatic azido groups: application as photoaffinity probes in molecular dynamics studies.

The accuracy of molecular dynamics (MD) simulations is limited by the availability of parameters for the molecular system of interest. In most force fields, parameters of common chemical groups are already present. With the development of novel small organic molecules as probes to study biological systems, more chemical groups require parameterization. An azide group is often used in studies of biological systems but computational studies are still impeded by the lack of parameters. In this paper, we present a set of molecular mechanics (MM) parameters for aromatic and aliphatic azido groups, and their application in MD simulations of a photoaffinity probe currently used in our laboratory for mapping binding modes available in the active site of histone deacetylases. The parameters were developed for the generalized Amber force field (GAFF) using density functional theory (DFT) calculations at B3LYP 6-311G(d) level. The parameters were validated by geometry optimization and MD simulations.

Sampling and convergence in free energy calculations of protein-ligand interactions: the binding of triphenoxypyridine derivatives to factor Xa and trypsin.

The binding of a set of 10 triphenoxypyridine derivatives to two serine proteases, factor Xa and trypsin, has been used to analyze factors related to sampling and convergence in free energy calculations based on molecular dynamics simulation techniques. The inhibitors investigated were initially proposed as part of the Critical Assessment of Techniques for Free Energy Evaluation (CATFEE) project for which no experimental results nor any assessment of the predictions submitted by various groups have ever been published. The inhibitors studied represent a severe challenge for explicit free energy calculations. The mutations from one compound to another involve up to 19 atoms, the creation and annihilation of net charge and several alternate binding modes. Nevertheless, we demonstrate that it is possible to obtain highly converged results (+/- 5-10 kJ/mol) even for such complex multi-atom mutations by simulating on a nanosecond time scale. This is achieved by using soft-core potentials to facilitate the creation and deletion of atoms and by a careful choice of mutation pathway. The results show that given modest computational resources, explicit free energy calculations can be successfully applied to realistic problems in drug design.