Next Article in Journal
Advances in Imaging of Inflammation, Fibrosis, and Cancer in the Gastrointestinal Tract
Previous Article in Journal
Urine Molecular Biomarkers for Detection and Follow-Up of Small Renal Masses
Previous Article in Special Issue
Visualizing Intramolecular Dynamics of Membrane Proteins
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Silico Study of the Interactions of Anle138b Isomer, an Inhibitor of Amyloid Aggregation, with Partner Proteins

by
Maxim S. Kondratyev
1,
Vladimir R. Rudnev
2,
Kirill S. Nikolsky
2,
Denis V. Petrovsky
2,
Liudmila I. Kulikova
2,
Kristina A. Malsagova
2,*,
Alexander A. Stepanov
2,
Arthur T. Kopylov
2 and
Anna L. Kaysheva
2
1
Institute of Cell Biophysics, Russian Academy of Sciences, 142290 Pushchino, Russia
2
Biobanking Group, Branch of Institute of Biomedical Chemistry “Scientific and Education Center”, 109028 Moscow, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(24), 16096; https://doi.org/10.3390/ijms232416096
Submission received: 1 November 2022 / Revised: 2 December 2022 / Accepted: 15 December 2022 / Published: 17 December 2022
(This article belongs to the Special Issue Advances in Membrane Protein Research)

Abstract

:
Herein, we aimed to highlight current “gaps” in the understanding of the potential interactions between the Anle138b isomer ligand, a promising agent for clinical research, and the intrinsically disordered alpha-synuclein protein. The presence of extensive unstructured areas in alpha-synuclein determines its existence in the cell of partner proteins, including the cyclophilin A chaperone, which prevents the aggregation of alpha-synuclein molecules that are destructive to cell life. Using flexible and cascaded molecular docking techniques, we aimed to expand our understanding of the molecular architecture of the protein complex between alpha-synuclein, cyclophilin A and the Anle138b isomer ligand. We demonstrated the possibility of intricate complex formation under cellular conditions and revealed that the main interactions that stabilize the complex are hydrophobic and involve hydrogen.

1. Introduction

Amyloidosis is a heterogeneous group of diseases characterized by the deposition of amyloid fibrils in the extracellular matrix of tissues and organs [1,2]. In clinical practice, nosologies accompanied by amyloid formation are often named after fibrillar proteins (atrial amyloidosis, prolactinomas, synucleinopathies, etc.) [1,2]. Approximately 95% of amyloids consist of fibrillar proteins, and the remaining 5% comprise the P-component of serum amyloid glycoproteins [3]. Insoluble forms of fibrillar proteins are formed via the misfolding of the soluble form of the precursor proteins. To date, more than 30 proteins have been discovered to form human amyloid fibrils [4], 18 of which cause systemic amyloidosis and 22 are localized forms. Systemic amyloidosis develops due to the aggregation of light (AL) and heavy (AH) chains of immunoglobulin procalcitonin (ACal), which are α-synuclein (AαSyn), islet amyloid polypeptide C (AIAPP), prolactin (APro), Aβ protein precursor (Aβ), and lung surfactant protein (ASPCd) [4].
Aggregates of amyloid proteins are assumed to have a wide range of negative effects on the human body, including cytotoxic effects and tissue damage. Diseases characterized by the formation of amyloid or other aggregates of amyloid proteins (cytotoxic oligomers, protofibrils) can be divided into a group of neurodegenerative diseases, which include Alzheimer’s disease (amyloid protein: Aβ), Parkinson’s disease (PD; amyloid protein: AαSyn), and type 2 diabetes (amyloid AIAPP protein).
PD is a complex neurodegenerative disease characterized by bradykinesia, resting tremors, muscle rigidity, and postural instability. The pathogenesis of PD is accompanied by the loss of dopaminergic neurons in the substantia nigra and the presence of Lewy bodies (intraneuronal accumulation of aggregated AαSyn molecules) in neurons in various regions of the brain [5,6,7,8]. Owing to the key role of AαSyn aggregation in PD pathogenesis, AαSyn is an attractive target for potential disease-modifying therapies. However, the “many-sidedness” [9] of the protein due to (a) the disordered nature of AαSyn (random coil in aqueous solution, α-helical structure upon interaction with lipid membranes) [8,10,11]; (b) a complex aggregation mechanism: the nucleation-polymerization and nucleation-conversion-polymerization model [12,13,14,15]; and (c) the lack of structural information on AαSyn aggregation complicates the discovery of α-synamyloid blockers [9,16].
Turning to evidence-based data, there are low-molecular-weight compounds that block AαSyn misfolding and aggregation and, in a sense, promote the dissociation of aggregates [9,16]. Aggregation inhibitors of AαSyn molecules include chemical chaperones, ligands Anle138b [17,18,19], SynuClean-D [20,21], and peptidomimetics (NPT100-18a and NPT200-11) [22,23,24], which are at different stages of clinical trials [8,16,25].
In the present study, we aimed to focus on the insufficient knowledge about the structure and, consequently, diversity of the AαSyn protein. Ultimately, we expanded our understanding of the molecular architecture for the potential formation of protein complexes under cellular conditions involving the intrinsically disordered protein (IDP), AαSyn, and the globular protein, cyclophylin A (CypA). These proteins were selected as targets for flexible docking owing to a lack of knowledge regarding the binding process of the Anle138b isomer ligand to partner proteins. First, we evaluated the interactions that stabilized the Anle138b isomer ligand complex with the alpha-synuclein fragment. Thereafter, we analyzed the ternary complex, AαSyn, with the cellular partner protein CypA, and the target ligand as the model closest to physiological (cellular) conditions and assessed the potential binding of the Anle138b isomer ligand to the CypA partner protein in the monomeric and dimeric forms by all-atom flexible docking. Finally, we analyzed the results of cascade docking (in the presence of an excess of ligand molecules) for AαSyn and CypA protein complexes of different compositions. Molecular dynamics experiments and evaluation of binding energy components showed good agreement with the results of molecular docking for the most significant amino acids involved in the stabilization of the studied complexes.

2. Results

2.1. Characterization of the AαSyn–Anle138b Isomer Complex

The structure of the AαSyn target was solved using electron microscopy with a resolution of 3.30 Å. The length of the amino acid sequence was 61 residues (Val37–Asp98), and the intact form of AαSyn was represented by 140 residues. The binding (affinity: −5.8 kcal/mol) of the Anle138b isomer molecule to a fragment of the alpha-synuclein monomer (PDB ID: 7V4C) is illustrated in Figure 1a.
The Val37–Asp98 region of AαSyn is rich in disorder-promoting amino acid residues [26,27] (Figure 1). The protein itself is a prominent representative of internally disordered proteins, which are characterized by extensive disordered regions. Based on an analysis of the obtained interactions (LigPlot), the geometry of the complexes was identified to be mainly stabilized by hydrophobic interactions with valine residues Val48, Val52, Val66, and Val71, and hydrogen bonding with Val52 of the main chain (Figure 1b).
If the target for the Anle138b isomer ligand is the complete AαSyn structure in the amino acid sequence (140 amino acid residues, PDB ID 1XQ8), the amino acids responsible for this binding were also identified (affinity: −6.7 kcal/mol). The results are shown in Figure 2.
The AαSyn molecule (PDB ID 1XQ8) is characterized by two extensively structured regions and two alpha helices, Met1–Val37 and Lys45–Thr92 (Figure 2). Alpha helices form a supersecondary α-hairpin structure which is stabilized by hydrogen and hydrophobic interactions [28]. A short Leu38–Thr44 constriction exists between the alpha helices. According to the molecular docking data, the amino acid residues of this unstructured region are involved in the stabilization of the ligand/target complex. In this case, the condensed hydrophobic rings of the ligand are oriented inside the hydrophobic core of the α-hairpin. Of note, the AαSyn structure retained a large C-terminal unstructured tail in the Gly93-Ala140 region.
Owing to an analysis of the microenvironment of the bound ligand molecule in the Leu38–Thr44 region of the target protein, Val40, Lys43, Lys32, and Tyr39 were identified as the main partner residues for stabilizing the spatial position of Anle138b isomer. Of note, several different interactions were detected for the Lys43 residue: hydrophobic, hydrogen bonding, and π-cationic interactions, enabling the identification of the Lys43 amino acid residue as a key partner for binding to the Anle138b isomer ligand (Figure 2d). Attention should be paid to the rather small distances of bonds between the ligand and target for the identified types of interactions, which indicates the fundamental possibility of the formation of such complexes in the cell. However, data from the literature indicate the formation of the Anle138b isomer complex with different sizes of AαSyn aggregates, but not with a monomeric form [16].

2.2. Characterization of the AαSyn–CypA and Anle138b Isomer Complex

In cells, AαSyn molecules form complexes with CypA [29,30]. CypA belongs to the family of immunophilins or chaperones containing the peptidyl prolyl cis/trans isomerase domain (PPIase). CypA is involved in basic neuronal functions, such as axonal transport and synaptic vesicle assembly; however, its most interesting role is neuroprotection. CypA prevents protein aggregation and binds unstructured or irregularly organized proteins [30]. As mentioned before, AαSyn is an unstructured protein, and its cellular partners can be different proteins and ligands. To study the binding features of the Anle138b isomer ligand and its partner protein, the experimentally known complex, AαSyn (region: Val48-Lys60), with the globular partner protein, CypA (PDB ID 6I42), as a target must be considered. The structure was resolved using X-ray crystallography with a resolution of 1.38 Å [29]. The results for the molecular docking of the ternary ligand-protein–protein complex are shown in Figure 3 (affinity: −7.1 kcal/mol).
As shown in Figure 3, the ligand was located on the mutual surface of AαSyn–CypA partner proteins, forming non-covalent interactions with each of the partner proteins. The binding site of the ligand with AαSyn also shifted toward the C-terminus of the protein, namely the Leu38–Thr44 region for the ligand–AαSyn complex compared to the Val48–Val52 region for the ligand–AαSyn–CypA complex. This phenomenon may be due to the experimental AαSyn–CypA structure, which was obtained for the AαSyn fragment. When interactions that stabilize the Anle138b isomer ternary complex and the AαSyn fragment with CypA are characterized, the hydrophobic interactions (His50, Asn106, Thr107) and hydrogen bonding (Glu81, Lys82) should be noted. The histidine residue belongs to alpha-synuclein, and the chart of interactions reveals its participation in π-stacking (Figure 3c). Based on our data, AαSyn in the ternary complex plays an important role in ligand binding. Currently, the map of interactions of the AαSyn–CypA protein complex in the absence of a ligand is well studied [29] and is described in more detail in the Section 3.

2.3. Molecular Docking: Characterization of the CypA and Anle138b Isomer Complex

During AαSyn–CypA complex formation under physiological conditions, the potential role of CypA in the binding of the studied ligand must be examined. To determine whether the evaluated Anle138b isomer ligand can compete with the binding sites of the globular cellular partner protein (CypA) as a target (PDB ID 6I42) instead of binding to alpha-synuclein, all-atom flexible docking was performed (affinity: −7.3 kcal/mol). The results are shown in Figure 4.
Owing to the chart of this complex, the amino acid residues responsible for its stabilization could be determined (Figure 4d). The main residues that determine the binding of the Anle138b isomer ligand to cyclophilin include Glu81 and Lys82, which are involved in both hydrophobic interactions (part of the side chains) and hydrogen bonding (end groups). Notably, distances of the revealed probable hydrophobic and hydrogen interactions are quite small. Accordingly, CypA can also act as a target for the ligand by stabilizing the ligand molecule on the surface of the protein globule. Similar to the case of the Anle138b isomer–AαSyn–CypA ternary complex, the same amino acid residues were involved in the stabilization of the Anle138b isomer–CypA complex. Further, the binding site remained unchanged.
Recent studies by Lutomski and Zhang have provided information on the formation of CypA dimers and trimers [31,32]. However, the biological significance of this phenomenon has not been evaluated. Aligning with this data, the binding of Anle138b isomer to the CypA dimer proved interesting. This dimer was obtained using ZDOCK (https://zdock.umassmed.edu/), (accessed on 25 November 2022). The results are shown in Figure 5 (affinity: −7.5 kcal/mol).
The location of ligand binding was found to differ from that of the monomer (Figure 5a,b). In the dimer, the ligand binds to the cleft surface formed by Arg55 (hydrogen bonding and π-cation interaction), Phe113 (π-stacking), and hydrophobic interactions through the Gln111 and Trp121 residues (Figure 5d). The chart of the bonds that stabilize the complex is illustrated in Figure 6c. Thus, for the formation of a CypA dimer, the ligand was located in the hollow formed by the interaction of the two CypA molecules.
To further predict the possible interaction of the Anle138b isomer ligand with cellular proteins, a molecular docking experiment was conducted for a more intricate complex by employing the dimeric form of the known AαSyn/CypA (PDB ID 6I42) complex. The dimer was obtained using ZDOCK, and the results obtained for the binding to Anle138b isomer are shown in Figure 6 (affinity: −7.7 kcal/mol).
In this version of the analysis, a different position was identified for the most favorable binding site compared to the dimer without the AαSyn fragment. It was found that the Anle138b isomer molecule binds on the surface of the globule in the groove under the chains of the AαSyn residues. Based on an analysis of the chart of the complex, the residues responsible for such binding could be identified (Figure 6b). The binding of the ligand occurs in a hollow formed by the surfaces of two cyclophilin globules, covered from above by the alpha-synuclein chains; the main residues stabilizing the complex are the Lys76 chains A and C (hydrophobic interactions through long chains and π-cationic interactions), and Thr68 residues (both chains, A and C). Further, Gly75A and Arg69C are involved in the formation of hydrogen bonds (Figure 6b). Simultaneously, the AαSyn fragment is directly involved in the formation of the binding site.
When the binding of a ligand to a molecule is assessed, the interactions that occur with an excess ligand must be analyzed. Such analysis is provided by cascade docking, in which a ligand is added in a step-by-step manner, molecule by molecule, after the optimal position of the previous ligand is incorporated into the target (see the Section 4). This approach allows full-atomic modeling of the binding of many ligand molecules and the detection of amino acid residues in the protein globule that can stabilize the formed complexes. Such calculations were performed for experimentally known molecules: series 1, the AαSyn fragment (PDB ID 7V4C); series 2, the full form of AαSyn (140 amino acid residues, PDB ID 1XQ8); and series 3, the CypA complex–AαSyn fragment (PDB ID 6I42). Each structure was surrounded by five Anle138b isomer ligand molecules.
Series 1. Results of cascade docking for a short AαSyn fragment are shown in Figure 7.
The close localization of binding sites in the target protein should be noted (region: His50 andThr72). The ligand molecules were oriented close to parallel stacking and were characterized by similar energy characteristics, −5.7 kcal/mol to −6.9 kcal/mol (Table 1).
The sweeps for this type of interaction are shown in Table 1. It is demonstrated that structures are mainly stabilized via hydrophobic interactions with the residues, Val66 (Ligands 1 and 4), Val71 (Ligands 1 and 2), and Thr72 (Ligands 3 and 4), and hydrogen bonds with the residues, His50 (Ligands 1 and 2) and Val52 (Ligands 1 and 2).
Series 2. Results of cascade docking for AαSyn (140 amino acid residues) are shown in Figure 8.
The analysis of the stabilization features of complexes in the presence of five molecules of the Anle138b isomer ligand with a full-length AαSyn, as in series 1, revealed an extremely close localization of the binding sites, which is especially evident against the background of the extended structure of the molecule and its large molecular surface. The energy characteristics of the positions of the ligand calculated herein (and for the variant of cascade docking of the ligand molecules and the AαSyn fragment (series 1) range from −6.7 kcal/mol to −5.8 kcal/mol (Table 2).
The sweeps for this type of interaction are shown in Figure 8. Consequently, the structures were stabilized by hydrophobic interactions, hydrogen bonds, π-stacking, and π-cation interactions. The most significant hydrophobic interactions in the ligand–target complex include Lys43 (ligands 1 and 3), Val48 (ligands 3 and 4), and hydrogen interactions with the Lys32 residue (ligands 1–3). Of note, the fourth ligand was only bound to ligands 2 and 3, without contact with the protein chain, forming a second (stacking) layer of Anle138b isomer ligands (Table 2).
Series 3. At the final stage of the work, cascade docking was carried out for the experimentally known (PDB ID 6I42) complex of the AαSyn fragment and its partner protein, CypA, as shown in Figure 9.
As this AαSyn–CypA target largely consists of a partner protein globule (and only a short peptide strand of AαSyn), the distribution of binding sites for all five Anle138b isomer ligands over the surface of the globule was expected. Only ligands 1 and 2 were found at the peptide interface. The energy characteristics of the five positions of the ligand calculated herein range from −7.1 kcal/mol to −6.6 kcal/mol (Table 3).
The sweeps for these interactions are shown in Figure 9. It was found that structures were stabilized by hydrophobic interactions, hydrogen bonds, π-stacking, and π-cationic interactions. For ligand 2, a halogen bond was found due to the presence of a bromine atom in the structure of the Anle138b isomer molecule.
In summary, the generally selective nature of the binding of Anle138b isomer ligand to alpha-synuclein molecules in different forms could be observed. For AαSyn complexes with the CypA partner protein, the preferred binding site is the interface between the two proteins. This study sought to assess the interactions of the ligand, whose interaction with AαSyn was postulated. Further, other ligands with such activity should be examined.

2.4. Molecular Dinamics: Characterization of the CypA and Anle138b Isomer Complex

Molecular dynamics experiments (GROMACS) were performed to determine the binding energies for the studied target–ligand complexes. The cumulative evaluation of the binding energies and the distribution of its components: van der Waal energy, electrostatic energy, polar solvation energy, and solvent-accessible surface energy (SASA energy) are presented in Table 4.
All energy components contribute significantly to the stabilization of the investigated variants of sets, with the exception of polar solvation energy. The calculated binding energy for ligand Anle138b isomer and AαSyn–CypA complex was −48.44 ± 26.11 kJ/mol, ligand Anle138b isomer and CypA dimer −98.53 ± 22.74 kJ/mol, ligand Anle138b isomer and dimer of the AαSyn–CypA complex −46.11 ± 24.69 kJ/mol, and ligand Anle138b isomer and AαSyn filament (homotrimer) −91± 27.29 kJ/mol. The time bases of the calculated binding energies are shown in Figure 10. The most stable in terms of cumulative binding energies are the ligand complexes with CypA dimer and AαSyn (filament).
Next, we calculated the contribution of each amino acid to the calculated values of the binding energy of each of the studied complex (Figure 10).
It was found that for the Anle138b isomer and AαSyn–CypA complex, the amino acid residues Ile57, Phe60, Met61, Trp121, Leu122, and Phe113 contribute most to the calculated binding energies (Figure 10a). At the same time, amino acids that were previously identified as significant in the stabilization of this complex in molecular docking experiments did not make a significant contribution to the calculated values of the binding energy.
For the second complex of Anle138b isomer and CypA dimer, we found that the most significant amino acid residues for the stabilization of the complex are the same as have been identified by the methods of molecular docking and molecular dynamics (Figure 10b). Such amino acid residues may include Arg55, Phe60, Gln63, Thr73, Gly74, Ala101, Gln111, Phe113 and Trp121. For the third Anle138b isomer complex and the dimer of the AαSyn–CypA, there is also agreement between the results of assessing the contribution of amino acid residues to complex formation: Thr68, Arg69, and Gly75 (Figure 10c). Additionally, for the fourth variant of Anle138b isomer and the AαSyn filament, we also observed the coincidence of the results of molecular docking and molecular dynamics for amino acid residues: Val48, Val52, Val66, Val71. For the studied complexes, the van der Waal energy, electrostatic energy, and polar solvation energy make the greatest contribution to the cumulative estimation of binding energies.

3. Discussion

In this study, we aimed to highlight some “gaps” in understanding the potential interactions of the Anle138b isomer ligand and the IDP, AαSyn. IDPs act as hubs in networks of protein interactions [33,34] and are important participants in cellular signaling mechanisms, including the regulation of transcription, translation, and the cell cycle [35,36]. The functional “many-sidedness” of IDPs, which is mediated by the formation of different types of protein complexes (protein–lipid raft, protein–protein, protein–ligand, protein–nucleic acid), is extremely important for studying the molecular mechanisms implemented in the cell in the state of “health” and “disease.”
The importance of studying the most probable variants of IDP interactions under physiological conditions was highlighted herein. These variants are not only closely associated with the implementation of normal cellular ontogenesis, but also with several nosologies. Our study focused on neurodegenerative diseases mediated by the formation of protein fibrils and aggregates (PD and AαSyn). Aligning with the above, the content of IDPs (25% of the total proteome) in the cell is strictly regulated to ensure accurate signaling in the cell cycle, while mutations in the protein-coding genes of IDPs or changes in the content of IDPs in the cell are associated with the development of diseases in humans [37,38,39]. Point mutation variants that cause PD have been annotated in the literature with high penetrance: Ala30Pro [40], Glu46Lys [41], Gly51Asp [42], His50Gln [43], Ala53Val [44], Ala53Glu [45], and Ala30Gly [46]. In some cases, carriers of point mutations in the AαSyn-coding gene are characterized by a relatively early onset of the disease (30–40 years), rapid progression, and severe cognitive impairment [47]. Of note, we revealed the participation (involvement) of these sites of amino acid substitutions in the formation of different variants of complexes between the ligand and alpha-synuclein molecule.
The presence of extensive unstructured regions in the AαSyn protein determines the presence of partner proteins in the cell, primarily chaperones, which are great representatives of agents that prevent the aggregation process that is destructive to cell life. Cellular globular chaperone peptidylprolyl isomerases (PPIases) catalyze the cis/trans proline isomerization of IDPs [48]. Cyclophilins are prominent representatives of PPIases and are highly represented in cellular proteins. CypA accounts for 0.1–0.6% of all cytosolic proteins [29]. Further, CypA plays an important role in the development of neurodegenerative diseases in humans. The role of CypA in neurodegenerative disorders is associated with the pathogenic aggregation of IDP [29], including the key pathogen in PD, AαSyn [49]. In a study by Favretto, two interaction sites between AαSyn and CypA were identified: Glu46–Gln62 in the central part of AαSyn and Val118–Glu131 in the proline-rich C-terminal region of AαSyn [29].
The structural plasticity and functional diversity of IDPs are “encoded” in the amino acid sequence. Francois-Xavier Theillet (2013) aptly defined such code as the “alphabet of intrinsic disorder” [26]. The sequences of IDPs are rich in amino acid residues that promote disorders, such as proline, aspartic acid, methionine, lysine, arginine, serine, glutamine, and glutamic acid, and contain few residues that promote “order,” namely cysteine, tryptophan, tyrosine, isoleucine, phenylalanine, valine, leucine, histidine, threonine, and asparagine [27]. This alphabet was applied in our study to describe the results of different docking options. Indeed, the AαSyn protein contains approximately 70% disorder-promoting residues (Figure 11), which are organized into at least three distinct clusters (regions: Met1–Val37, Ala56–Ala91, and Gly93–Ala140).
Herein, the possible interactions of the Anle138B isomer with AαSyn and CypA proteins in different versions of the complex was described to enrich the current understanding of processes that may occur in living cells. Using molecular docking methods and molecular dynamics, we examined the most important events, in our opinion, that will expand the current understanding of the pathogenesis of PD. Comparison of the results of molecular docking and molecular dynamics showed good convergence for the most significant amino acids involved in the stabilization of the studied complexes.

4. Materials and Methods

4.1. Objects of Study

The low-molecular-weight drug, Anle138b isomer (Figure 12), is an oligomer modulator and aggregation inhibitor developed based on a systematic high-throughput screening campaign combined with medicinal chemistry optimization [25].
The value of the gas-phase heat of formation of the Anle138b isomer molecule is positive and amounts to 16.6 kcal/mol, which is typical for aromatic compounds (Figure 12). The Anle138b isomer ligand was detected after two consecutive screening rounds. Of note, the ligand does not interact with the AαSyn monomer, but binds to the hydrophobic pocket in the oligomeric states of the target protein, preventing the formation of the β-layer and AαSyn aggregates [16]. Anle138b isomer specifically inhibits the formation of AαSyn in oligomers in vitro and in the human cell line, HEK293 [50]. The physical and chemical properties of Anle138b isomer are outlined in Table 5.
AαSyn and CypA were selected as target proteins for the Anle138b isomer ligand (Table 6). AαSyn belongs to a large group of internally disordered proteins, while CypA is a cellular chaperone with a globular and compact structure.

4.2. Molecular Docking

The protein structures (as targets) were prepared for docking according to the standard scheme for the Autodock Vina package, and the atoms (together with their coordinates) of the solvent, buffer, and ligand molecules were removed from the PDB input file. Before the numerical calculations were performed, a charge was placed on the surface of the proteins using MGLTools. The center of the molecule and the parameters of the box («cells») were set manually, ensuring that the entire protease molecule was placed within the computational area of space.
We assembled the AαSyn–CypA dimer complex by rigid protein–protein docking using the Z-Dock service. This procedure has been repeatedly tested by various teams of researchers and is described in detail in the article [52].
The Anle138b isomer (ligand) structure model was drawn in the HyperChem molecular designer; this structure was consistently optimized in the AMBER force field, and then quantum-chemically in the PM3 (пaкeт MOPAC) [53]. The optimization consisted of calculating the complete electronic structure of the molecule, characterizing the geometry of all chemical bonds, and valence and torsion angles of the molecule. In the same calculation, the distribution of partial charges on the surface of the ligand and its dipole moment were determined. The ligands in the docking calculations had maximum conformational freedom. Of note, rotation of the functional groups around all single bonds was allowed. The charge arrangement on the Anle138b isomer molecule and its protonation–deprotonation was performed automatically in the MGLTools 1.5.6 package.
To obtain the in silico results, sequential (cascade, multiple) docking was applied. For Anle138b isomer, the five binding sites were successively modeled as described below. First, docking was performed for the «protein–ligand» complex. Thereafter, the model was augmented with another ligand (i.e., three-monomer chain) and so on. Thus, «blind» or cascade docking was employed. First, the optimal position of the ligand was calculated; the structure of the ligand (fixed at the docking site) became an integral part of the target. Therefore, this binding site (with a high affinity for the ligand) was blocked. The position of the second ligand, whose molecule (in its optimal position on the globule) also became a part of the target, was subsequently calculated. Ultimately, both bound ligands become part of the receptor. These iterative searches and filling of the optimal binding sites on the target surface were repeated until the positions of all five ligands were calculated. Of note, this cascade docking is reasonable for modeling ligand–receptor interactions for excess ligand molecules. Using this procedure, we successively excluded several possible binding sites of ligands via stepwise filling of the available areas on the surface of the protein–ligand complex. In the course of the analysis of the obtained complexes, the terms, «bonds» (understanding drawn interatomic chemical bonds, hydrogen bonds between their donors and acceptors) and «interactions» (mainly understanding hydrophobic bonds), were used. In general, cascade docking enables modeling and analysis of the binding between the protein and ligands in excess of ligands at the atomic level (Anle138b isomer).
The analysis of the docking results was nontrivial. To obtain the «sweeps» of the interactions, the LIGPLOT+ package [54] was employed. From the three-dimensional coordinates of the atoms of the protein–ligand complex, interaction diagrams depicting the schemes of hydrogen bonds and hydrophobic contacts between the ligand and elements of the main or side chain of the protein were generated. For in-depth analysis of the interactions between the constituent parts of the complexes, the PLIP (Protein–Ligand Interaction Profiler) package [55] was employed. They are distinguished by accuracy and content. The PLIP server can detect and display the following interactions: hydrogen bonds, hydrophobic interactions, water bridges, parallel stacking, perpendicular stacking, π-cation bonds, halogen bonds, salt bridges, and metal complexes. In addition to the figure, the package displays tables with distances, angles, and types of detected interactions, which are important for the analysis of previously obtained experimental data.

4.3. Molecular Dynamics

Molecular dynamics simulations were performed using GROMACS software packet. Ligand topology structure were obtained from the Automated Topology Builder (ATB) resource [56]. Calculations were performed using GROMOS547 forcefield modified by ATB. This forcefield was chosen because it supports bromine (contained in Anle138b isomer) and there was already a prepared topology for the ligand in the ATB service (ATB provides topology for GROMOS54a7 forcefield only). It is not as popular as AMBER Force Field but still used in MD experiments [57].
The molecular dynamics process was based on one described in GROMACS protein–ligand tutorial [58]. Simulations were run in cubic boxes. The box size was automatically calculated by GROMACS based on complex size and adding 1 nm from the complex size to the box border. The water SPC model was used as a solvent as it is supported by the chosen forcefield and grants significant simulation performance improvement with acceptable losses in accuracy. Systems were neutralized automatically with Na+ or Cl ions (one at a time). Then, energy minimization and equilibration steps were performed using the configuration found in the same tutorial [58]. The molecular dynamics themselves were performed for 100 ns in three different runs (using the same equilibrated complex topology for each run). The simulation system temperature was set to 300 K.
Analysis was performed for the binding energies of the protein–ligand complex trajectory following the MM/PBSA approach [59,60,61]. We performed analyses of the binding energy of protein–ligand complexes and binding energy distributions for every residue. This approach to trajectory analysis was based on one found in the study [62].

5. Conclusions

Pathophysiological conditions of neurodegenerative diseases are undoubtedly related to protein misfolding. The organization of such proteins into relatively ordered structures, such as fibrillar intracellular and extracellular amyloids, leads to tissue and brain damage. In these pathologies, the appearance of protein aggregates indicates inefficient cellular reactions involving molecular chaperones that contribute to the correct folding of partner proteins [30]. This study aimed to expand our understanding of the possible mechanisms of interaction of the Anle138b isomer ligand, which is promising for clinical use and prevents the aggregation of cellular amyloidogenic proteins. This field of study reflects the urgent clinical need to develop neuroprotective molecules targeting potentially novel molecular targets. Our understanding of the interaction of the Anle138b isomer ligand with variants of the alpha-synuclein molecule was expanded, and the possibility of forming more intricate complexes with CypA, a partner protein of alpha-synuclein, was revealed.
Our study paves the way for experimental confirmation, including an analysis of mutant forms of alpha-synuclein, which cause high disease penetrance, early onset, and rapid progression of PD.

Author Contributions

Conceptualization, A.L.K., M.S.K. and V.R.R.; methodology, M.S.K.; software, K.S.N., D.V.P. and A.A.S.; validation, V.R.R. and L.I.K.; formal analysis, A.T.K., and K.A.M.; investigation, A.L.K., and V.R.R.; data curation, A.L.K.; writing—original draft preparation, M.S.K. and A.L.K.; writing—review and editing, K.A.M. All authors have read and agreed to the published version of the manuscript.

Funding

The work was conducted in the framework of the Russian Federation Fundamental Research Program for the long-term period of 2021–2030 (№122092200056-9).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rambaran, R.N.; Serpell, L.C. Amyloid fibrils. Prion 2008, 2, 112–117. [Google Scholar] [CrossRef] [PubMed]
  2. Khoor, A.; Colby, T.V. Amyloidosis of the Lung. Arch. Pathol. Lab. Med. 2017, 141, 247–254. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Snyder, P.W. Chapter 5—Diseases of Immunity1. In Pathologic Basis of Veterinary Disease, 6th ed.; Zachary, J.F., Ed.; Mosby: Maryland Heights, MO, USA, 2017; pp. 242–285.e5. ISBN 978-0-323-35775-3. [Google Scholar]
  4. Benson, M.D.; Buxbaum, J.N.; Eisenberg, D.S.; Merlini, G.; Saraiva, M.J.M.; Sekijima, Y.; Sipe, J.D.; Westermark, P. Amyloid nomenclature 2020: Update and recommendations by the International Society of Amyloidosis (ISA) nomenclature committee. Amyloid Int. J. Exp. Clin. Investig. Off. J. Int. Soc. Amyloidosis 2020, 27, 217–222. [Google Scholar] [CrossRef] [PubMed]
  5. Stefanis, L. α-Synuclein in Parkinson’s Disease. Cold Spring Harb. Perspect. Med. 2012, 2, a009399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Gómez-Benito, M.; Granado, N.; García-Sanz, P.; Michel, A.; Dumoulin, M.; Moratalla, R. Modeling Parkinson’s Disease with the Alpha-Synuclein Protein. Front. Pharmacol. 2020, 11, 356. [Google Scholar] [CrossRef] [PubMed]
  7. Srinivasan, E.; Chandrasekhar, G.; Chandrasekar, P.; Anbarasu, K.; Vickram, A.S.; Karunakaran, R.; Rajasekaran, R.; Srikumar, P.S. Alpha-Synuclein Aggregation in Parkinson’s Disease. Front. Med. 2021, 8, 736978. [Google Scholar] [CrossRef]
  8. Vidović, M.; Rikalovic, M.G. Alpha-Synuclein Aggregation Pathway in Parkinson’s Disease: Current Status and Novel Therapeutic Approaches. Cells 2022, 11, 1732. [Google Scholar] [CrossRef]
  9. Lashuel, H.A.; Overk, C.R.; Oueslati, A.; Masliah, E. The many faces of α-synuclein: From structure and toxicity to therapeutic target. Nat. Rev. Neurosci. 2013, 14, 38–48. [Google Scholar] [CrossRef] [Green Version]
  10. Ulmer, T.S.; Bax, A.; Cole, N.B.; Nussbaum, R.L. Structure and dynamics of micelle-bound human alpha-synuclein. J. Biol. Chem. 2005, 280, 9595–9603. [Google Scholar] [CrossRef] [Green Version]
  11. Georgieva, E.R.; Ramlall, T.F.; Borbat, P.P.; Freed, J.H.; Eliezer, D. Membrane-Bound Alpha-Synuclein Forms an Extended Helix: Long-Distance Pulsed ESR Measurements Using Vesicles, Bicelles, and Rod-Like Micelles. J. Am. Chem. Soc. 2008, 130, 12856–12857. [Google Scholar] [CrossRef]
  12. de Oliveira, G.A.P.; Silva, J.L. Alpha-synuclein stepwise aggregation reveals features of an early onset mutation in Parkinson’s disease. Commun. Biol. 2019, 2, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Cremades, N.; Chen, S.W.; Dobson, C.M. Structural Characteristics of α-Synuclein Oligomers. Int. Rev. Cell Mol. Biol. 2017, 329, 79–143. [Google Scholar]
  14. Cremades, N.; Cohen, S.I.A.; Deas, E.; Abramov, A.Y.; Chen, A.Y.; Orte, A.; Sandal, M.; Clarke, R.W.; Dunne, P.; Aprile, F.A.; et al. Direct observation of the interconversion of normal and toxic forms of α-synuclein. Cell 2012, 149, 1048–1059. [Google Scholar] [CrossRef] [Green Version]
  15. Iljina, M.; Garcia, G.A.; Horrocks, M.H.; Tosatto, L.; Choi, M.L.; Ganzinger, K.A.; Abramov, A.Y.; Gandhi, S.; Wood, N.W.; Cremades, N.; et al. Kinetic model of the aggregation of alpha-synuclein provides insights into prion-like spreading. Proc. Natl. Acad. Sci. USA 2016, 113, E1206–E1215. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Pujols, J.; Peña-Díaz, S.; Pallarès, I.; Ventura, S. Chemical Chaperones as Novel Drugs for Parkinson’s Disease. Trends Mol. Med. 2020, 26, 408–421. [Google Scholar] [CrossRef]
  17. Levin, J.; Sing, N.; Melbourne, S.; Morgan, A.; Mariner, C.; Spillantini, M.G.; Wegrzynowicz, M.; Dalley, J.W.; Langer, S.; Ryazanov, S.; et al. Safety, tolerability and pharmacokinetics of the oligomer modulator anle138b with exposure levels sufficient for therapeutic efficacy in a murine Parkinson model: A randomised, double-blind, placebo-controlled phase 1a trial. eBioMedicine 2022, 80, 104021. [Google Scholar] [CrossRef] [PubMed]
  18. Fellner, L.; Kuzdas-Wood, D.; Levin, J.; Ryazanov, S.; Leonov, A.; Griesinger, C.; Giese, A.; Wenning, G.K.; Stefanova, N. Anle138b Partly Ameliorates Motor Deficits Despite Failure of Neuroprotection in a Model of Advanced Multiple System Atrophy. Front. Neurosci. 2016, 10, 99. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Antonschmidt, L.; Matthes, D.; Dervişoğlu, R.; Frieg, B.; Dienemann, C.; Leonov, A.; Nimerovsky, E.; Sant, V.; Ryazanov, S.; Giese, A.; et al. The clinical drug candidate anle138b binds in a cavity of lipidic α-synuclein fibrils. Nat. Commun. 2022, 13, 5385. [Google Scholar] [CrossRef]
  20. Peña-Díaz, S.; Pujols, J.; Vasili, E.; Pinheiro, F.; Santos, J.; Manglano-Artuñedo, Z.; Outeiro, T.F.; Ventura, S. The small aromatic compound SynuClean-D inhibits the aggregation and seeded polymerization of multiple α-synuclein strains. J. Biol. Chem. 2022, 298, 101902. [Google Scholar] [CrossRef]
  21. Pujols, J.; Peña-Díaz, S.; Lázaro, D.F.; Peccati, F.; Pinheiro, F.; González, D.; Carija, A.; Navarro, S.; Conde-Giménez, M.; García, J.; et al. Small molecule inhibits α-synuclein aggregation, disrupts amyloid fibrils, and prevents degeneration of dopaminergic neurons. Proc. Natl. Acad. Sci. USA 2018, 115, 10481–10486. [Google Scholar] [CrossRef] [Green Version]
  22. Wrasidlo, W.; Tsigelny, I.F.; Price, D.L.; Dutta, G.; Rockenstein, E.; Schwarz, T.C.; Ledolter, K.; Bonhaus, D.; Paulino, A.; Eleuteri, S.; et al. A de novo compound targeting α-synuclein improves deficits in models of Parkinson’s disease. Brain 2016, 139, 3217–3236. [Google Scholar] [CrossRef] [PubMed]
  23. Prots, I.; Grosch, J.; Brazdis, R.-M.; Simmnacher, K.; Veber, V.; Havlicek, S.; Hannappel, C.; Krach, F.; Krumbiegel, M.; Schütz, O.; et al. α-Synuclein oligomers induce early axonal dysfunction in human iPSC-based models of synucleinopathies. Proc. Natl. Acad. Sci. USA 2018, 115, 7813–7818. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Price, D.L.; Koike, M.A.; Khan, A.; Wrasidlo, W.; Rockenstein, E.; Masliah, E.; Bonhaus, D. The small molecule alpha-synuclein misfolding inhibitor, NPT200-11, produces multiple benefits in an animal model of Parkinson’s disease. Sci. Rep. 2018, 8, 16165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wagner, J.; Ryazanov, S.; Leonov, A.; Levin, J.; Shi, S.; Schmidt, F.; Prix, C.; Pan-Montojo, F.; Bertsch, U.; Mitteregger-Kretzschmar, G.; et al. Anle138b: A novel oligomer modulator for disease-modifying therapy of neurodegenerative diseases such as prion and Parkinson’s disease. Acta Neuropathol. 2013, 125, 795–813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Theillet, F.-X.; Kalmar, L.; Tompa, P.; Han, K.-H.; Selenko, P.; Dunker, A.K.; Daughdrill, G.W.; Uversky, V.N. The alphabet of intrinsic disorder: I. Act like a Pro: On the abundance and roles of proline residues in intrinsically disordered proteins. Intrinsically Disord. Proteins 2013, 1, e24360. [Google Scholar] [CrossRef] [Green Version]
  27. Williams, R.M.; Obradovi, Z.; Mathura, V.; Braun, W.; Garner, E.C.; Young, J.; Takayama, S.; Brown, C.J.; Dunker, A.K. The protein non-folding problem: Amino acid determinants of intrinsic order and disorder. Pac. Symp. Biocomput. 2001, 89–100. [Google Scholar] [CrossRef] [Green Version]
  28. Rudnev, V.R.; Kulikova, L.I.; Nikolsky, K.S.; Malsagova, K.A.; Kopylov, A.T.; Kaysheva, A.L. Current Approaches in Supersecondary Structures Investigation. Int. J. Mol. Sci. 2021, 22, 11879. [Google Scholar] [CrossRef]
  29. Favretto, F.; Baker, J.D.; Strohäker, T.; Andreas, L.B.; Blair, L.J.; Becker, S.; Zweckstetter, M. The Molecular Basis of the Interaction of Cyclophilin A with α-Synuclein. Angew. Chem. Int. Ed Engl. 2020, 59, 5643–5646. [Google Scholar] [CrossRef]
  30. Daneri-Becerra, C.; Ciucci, S.M.; Mazaira, G.; Galigniana, M.D. Role of Mitochondrial Heat-shock Proteins and Immunophilins in Neuro Degenerative Diseases. Curr. Drug Targets 2021, 22, 1596–1617. [Google Scholar] [CrossRef]
  31. Zhang, X.-C.; Wang, W.-D.; Wang, J.-S.; Pan, J.-C.; Zou, G.-L. Evidences of monomer, dimer and trimer of recombinant human cyclophilin A. Protein Pept. Lett. 2011, 18, 1188–1193. [Google Scholar] [CrossRef]
  32. Lutomski, C.A.; El-Baba, T.J.; Bolla, J.R.; Robinson, C.V. Multiple Roles of SARS-CoV-2 N Protein Facilitated by Proteoform-Specific Interactions with RNA, Host Proteins, and Convalescent Antibodies. JACS Au 2021, 1, 1147–1157. [Google Scholar] [CrossRef] [PubMed]
  33. Dunker, A.K.; Cortese, M.S.; Romero, P.; Iakoucheva, L.M.; Uversky, V.N. Flexible nets. The roles of intrinsic disorder in protein interaction networks. FEBS J. 2005, 272, 5129–5148. [Google Scholar] [CrossRef] [PubMed]
  34. Kim, P.M.; Sboner, A.; Xia, Y.; Gerstein, M. The role of disorder in interaction networks: A structural analysis. Mol. Syst. Biol. 2008, 4, 179. [Google Scholar] [CrossRef] [PubMed]
  35. Gandass, N.; Kajal; Salvi, P. Intrinsically disordered protein, DNA binding with one finger transcription factor (OsDOF27) implicates thermotolerance in yeast and rice. Front. Plant Sci. 2022, 13, 956299. [Google Scholar] [CrossRef] [PubMed]
  36. Wallin, S. Intrinsically disordered proteins: Structural and functional dynamics. Res. Rep. Biol. 2017, 8, 7–16. [Google Scholar] [CrossRef] [Green Version]
  37. Wright, P.E.; Dyson, H.J. Intrinsically disordered proteins in cellular signalling and regulation. Nat. Rev. Mol. Cell Biol. 2015, 16, 18–29. [Google Scholar] [CrossRef]
  38. Kulkarni, P.; Uversky, V.N. Intrinsically Disordered Proteins in Chronic Diseases. Biomolecules 2019, 9, 147. [Google Scholar] [CrossRef] [Green Version]
  39. Uversky, V.N.; Oldfield, C.J.; Dunker, A.K. Intrinsically disordered proteins in human diseases: Introducing the D2 concept. Annu. Rev. Biophys. 2008, 37, 215–246. [Google Scholar] [CrossRef]
  40. Krüger, R.; Kuhn, W.; Müller, T.; Woitalla, D.; Graeber, M.; Kösel, S.; Przuntek, H.; Epplen, J.T.; Schöls, L.; Riess, O. Ala30Pro mutation in the gene encoding alpha-synuclein in Parkinson’s disease. Nat. Genet. 1998, 18, 106–108. [Google Scholar] [CrossRef]
  41. Zarranz, J.J.; Alegre, J.; Gómez-Esteban, J.C.; Lezcano, E.; Ros, R.; Ampuero, I.; Vidal, L.; Hoenicka, J.; Rodriguez, O.; Atarés, B.; et al. The new mutation, E46K, of alpha-synuclein causes Parkinson and Lewy body dementia. Ann. Neurol. 2004, 55, 164–173. [Google Scholar] [CrossRef]
  42. Lesage, S.; Anheim, M.; Letournel, F.; Bousset, L.; Honoré, A.; Rozas, N.; Pieri, L.; Madiona, K.; Dürr, A.; Melki, R.; et al. G51D α-synuclein mutation causes a novel parkinsonian-pyramidal syndrome. Ann. Neurol. 2013, 73, 459–471. [Google Scholar] [CrossRef] [PubMed]
  43. Appel-Cresswell, S.; Vilarino-Guell, C.; Encarnacion, M.; Sherman, H.; Yu, I.; Shah, B.; Weir, D.; Thompson, C.; Szu-Tu, C.; Trinh, J.; et al. Alpha-synuclein p.H50Q, a novel pathogenic mutation for Parkinson’s disease. Mov. Disord. Off. J. Mov. Disord. Soc. 2013, 28, 811–813. [Google Scholar] [CrossRef] [PubMed]
  44. Yoshino, H.; Hirano, M.; Stoessl, A.J.; Imamichi, Y.; Ikeda, A.; Li, Y.; Funayama, M.; Yamada, I.; Nakamura, Y.; Sossi, V.; et al. Homozygous alpha-synuclein p.A53V in familial Parkinson’s disease. Neurobiol. Aging 2017, 57, 248.e7–248.e12. [Google Scholar] [CrossRef] [PubMed]
  45. Pasanen, P.; Myllykangas, L.; Siitonen, M.; Raunio, A.; Kaakkola, S.; Lyytinen, J.; Tienari, P.J.; Pöyhönen, M.; Paetau, A. A novel α-synuclein mutation A53E associated with atypical multiple system atrophy and Parkinson’s disease-type pathology. Neurobiol. Aging 2014, 35, 2180.e1–2180.e5. [Google Scholar] [CrossRef] [PubMed]
  46. Liu, H.; Koros, C.; Strohäker, T.; Schulte, C.; Bozi, M.; Varvaresos, S.; Ibáñez de Opakua, A.; Simitsi, A.M.; Bougea, A.; Voumvourakis, K.; et al. A Novel SNCA A30G Mutation Causes Familial Parkinson’s Disease. Mov. Disord. 2021, 36, 1624–1633. [Google Scholar] [CrossRef]
  47. Oliveira, L.M.A.; Gasser, T.; Edwards, R.; Zweckstetter, M.; Melki, R.; Stefanis, L.; Lashuel, H.A.; Sulzer, D.; Vekrellis, K.; Halliday, G.M.; et al. Alpha-synuclein research: Defining strategic moves in the battle against Parkinson’s disease. Npj Park. Dis. 2021, 7, 1–23. [Google Scholar] [CrossRef]
  48. Lin, W.; Bonin, M.; Boden, A.; Wieduwild, R.; Murawala, P.; Wermke, M.; Andrade, H.; Bornhäuser, M.; Zhang, Y. Peptidyl prolyl cis/trans isomerase activity on the cell surface correlates with extracellular matrix development. Commun. Biol. 2019, 2, 1–11. [Google Scholar] [CrossRef] [Green Version]
  49. Braak, H.; Ghebremedhin, E.; Rüb, U.; Bratzke, H.; Del Tredici, K. Stages in the development of Parkinson’s disease-related pathology. Cell Tissue Res. 2004, 318, 121–134. [Google Scholar] [CrossRef]
  50. Dominguez-Meijide, A.; Vasili, E.; König, A.; Cima-Omori, M.-S.; Ibáñez de Opakua, A.; Leonov, A.; Ryazanov, S.; Zweckstetter, M.; Griesinger, C.; Outeiro, T.F. Effects of pharmacological modulators of α-synuclein and tau aggregation and internalization. Sci. Rep. 2020, 10, 12827. [Google Scholar] [CrossRef]
  51. Braun, A.R.; Liao, E.E.; Horvath, M.; Kalra, P.; Acosta, K.; Young, M.C.; Kochen, N.N.; Lo, C.H.; Brown, R.; Evans, M.D.; et al. Potent inhibitors of toxic alpha-synuclein identified via cellular time-resolved FRET biosensors. Npj Park. Dis. 2021, 7, 1–17. [Google Scholar] [CrossRef]
  52. Pierce, B.G.; Wiehe, K.; Hwang, H.; Kim, B.-H.; Vreven, T.; Weng, Z. ZDOCK server: Interactive docking prediction of protein-protein complexes and symmetric multimers. Bioinf. Oxf. Engl. 2014, 30, 1771–1773. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Stewart, J.J.P. Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements. J. Mol. Model. 2007, 13, 1173–1213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Laskowski, R.A.; Swindells, M.B. LigPlot+: Multiple ligand-protein interaction diagrams for drug discovery. J. Chem. Inf. Model. 2011, 51, 2778–2786. [Google Scholar] [CrossRef]
  55. Adasme, M.F.; Linnemann, K.L.; Bolz, S.N.; Kaiser, F.; Salentin, S.; Haupt, V.J.; Schroeder, M. PLIP 2021: Expanding the scope of the protein-ligand interaction profiler to DNA and RNA. Nucleic Acids Res. 2021, 49, W530–W534. [Google Scholar] [CrossRef] [PubMed]
  56. Malde, A.K.; Zuo, L.; Breeze, M.; Stroet, M.; Poger, D.; Nair, P.C.; Oostenbrink, C.; Mark, A.E. An Automated Force Field Topology Builder (ATB) and Repository: Version 1.0. J. Chem. Theory Comput. 2011, 7, 4026–4037. [Google Scholar] [CrossRef]
  57. In Silico Study of 3-Hydroxypyrimidine-2,4-diones as Inhibitors of HIV RT-Associated RNase H Using Molecular Docking, Molecular Dynamics, 3D-QSAR, and Pharmacophore Models—New Journal of Chemistry (RSC Publishing). Available online: https://pubs.rsc.org/en/content/articlelanding/2019/nj/c9nj03353j (accessed on 1 December 2022).
  58. Protein-Ligand Complex. Available online: http://www.mdtutorials.com/gmx/complex/index.html (accessed on 1 December 2022).
  59. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef]
  60. Kumari, R.; Kumar, R.; Open Source Drug Discovery Consortium; Lynn, A. g_mmpbsa—A GROMACS tool for high-throughput MM-PBSA calculations. J. Chem. Inf. Model. 2014, 54, 1951–1962. [Google Scholar] [CrossRef]
  61. Baker, N.A.; Sept, D.; Joseph, S.; Holst, M.J.; McCammon, J.A. Electrostatics of nanosystems: Application to microtubules and the ribosome. Proc. Natl. Acad. Sci. USA 2001, 98, 10037–10041. [Google Scholar] [CrossRef] [Green Version]
  62. Sajid Jamal, Q.M.; Alharbi, A.H.; Ahmad, V. Identification of doxorubicin as a potential therapeutic against SARS-CoV-2 (COVID-19) protease: A molecular docking and dynamics simulation studies. J. Biomol. Struct. Dyn. 2022, 40, 7960–7974. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) Binding of the Anle138b isomer molecule to a fragment of the alpha-synuclein monomer (PDB ID 7V4C); (b) chart of interactions when Anle138b isomer binds to a fragment of the alpha-synuclein monomer (region: 37–98); (c) the amino acid sequence of AαSyn, the green box represents the region of the Val37–Asp98 sequence for which the molecular docking experiment was performed (PDB ID 7V4C), amino acid residues that belong to the disorder-promoting set are marked in red.
Figure 1. (a) Binding of the Anle138b isomer molecule to a fragment of the alpha-synuclein monomer (PDB ID 7V4C); (b) chart of interactions when Anle138b isomer binds to a fragment of the alpha-synuclein monomer (region: 37–98); (c) the amino acid sequence of AαSyn, the green box represents the region of the Val37–Asp98 sequence for which the molecular docking experiment was performed (PDB ID 7V4C), amino acid residues that belong to the disorder-promoting set are marked in red.
Ijms 23 16096 g001aIjms 23 16096 g001b
Figure 2. (a,b) Anle138b isomer molecule binding to full-length alpha-synuclein molecule (PDB ID 1XQ8). The 3D structure of the Anle138b isomer ligand is indicated in green. (c) Interactions diagram of Anle138b isomer with a fragment of the alpha-synuclein monomer. (a) full image size, (b) enlargement of the binding site image. The figure above shows the amino acid sequence of AαSyn, the blocks reflect the organization of secondary structure elements: yellow blocks correspond to alpha helices, gray blocks correspond to unstructured regions (coil). (d) Interaction options for Anle138b isomer binding to the alpha-synuclein monomer (PDB ID 1XQ8).
Figure 2. (a,b) Anle138b isomer molecule binding to full-length alpha-synuclein molecule (PDB ID 1XQ8). The 3D structure of the Anle138b isomer ligand is indicated in green. (c) Interactions diagram of Anle138b isomer with a fragment of the alpha-synuclein monomer. (a) full image size, (b) enlargement of the binding site image. The figure above shows the amino acid sequence of AαSyn, the blocks reflect the organization of secondary structure elements: yellow blocks correspond to alpha helices, gray blocks correspond to unstructured regions (coil). (d) Interaction options for Anle138b isomer binding to the alpha-synuclein monomer (PDB ID 1XQ8).
Ijms 23 16096 g002aIjms 23 16096 g002b
Figure 3. (a) Binding of the Anle138b isomer molecule to the AαSyn–CypA complex (PDB ID 6I42), a fragment of the surface of the AαSyn molecule is indicated in green, the surface of the CypA globular protein is indicated in gray, the molecular surface of the Anle138b isomer ligand is indicated in blue; (b) chart of interactions upon binding of Anle138b isomer to the AαSyn–CypA complex; (c) interactions diagram of Anle138b isomer with AαSyn–CypA complex.
Figure 3. (a) Binding of the Anle138b isomer molecule to the AαSyn–CypA complex (PDB ID 6I42), a fragment of the surface of the AαSyn molecule is indicated in green, the surface of the CypA globular protein is indicated in gray, the molecular surface of the Anle138b isomer ligand is indicated in blue; (b) chart of interactions upon binding of Anle138b isomer to the AαSyn–CypA complex; (c) interactions diagram of Anle138b isomer with AαSyn–CypA complex.
Ijms 23 16096 g003
Figure 4. Binding of the Anle138b isomer molecule to CypA (a cellular partner protein of alpha-synuclein) (a) full size image; (b) an enlarged image of a possible CypA binding site; (c) interactions diagram of Anle138b isomer with cyclophilin; (d) interaction options for Anle138b isomer binding to CypA.
Figure 4. Binding of the Anle138b isomer molecule to CypA (a cellular partner protein of alpha-synuclein) (a) full size image; (b) an enlarged image of a possible CypA binding site; (c) interactions diagram of Anle138b isomer with cyclophilin; (d) interaction options for Anle138b isomer binding to CypA.
Ijms 23 16096 g004
Figure 5. (a,b) binding of the Anle138b isomer molecule to the cyclophilin dimer (two views); (c) interactions diagram of Anle138b isomer with cyclophilin dimer; (d) interaction options for Anle138b isomer binding to the CypA dimer.
Figure 5. (a,b) binding of the Anle138b isomer molecule to the cyclophilin dimer (two views); (c) interactions diagram of Anle138b isomer with cyclophilin dimer; (d) interaction options for Anle138b isomer binding to the CypA dimer.
Ijms 23 16096 g005
Figure 6. (a) Binding of the Anle138b isomer molecule to the dimer of the CypA complex (PDB ID 6I42). CypA subunits are shown in purple and orange, the AαSyn fragment is shown in green and blue, the ligand structure is shown in red; (b) chart of interactions upon binding of Anle138b isomer to the dimer of the CypA and AαSyn complex; (c) interactions diagram of Anle138b isomer binding with CypA dimer.
Figure 6. (a) Binding of the Anle138b isomer molecule to the dimer of the CypA complex (PDB ID 6I42). CypA subunits are shown in purple and orange, the AαSyn fragment is shown in green and blue, the ligand structure is shown in red; (b) chart of interactions upon binding of Anle138b isomer to the dimer of the CypA and AαSyn complex; (c) interactions diagram of Anle138b isomer binding with CypA dimer.
Ijms 23 16096 g006
Figure 7. Five Anle138b isomer ligand binding sites on a filament of AαSyn (PDB ID 7V4C). Individual Anle138b isomer ligands are numbered: 1—red, 2—yellow, 3—green, 4—blue, 5—violet.
Figure 7. Five Anle138b isomer ligand binding sites on a filament of AαSyn (PDB ID 7V4C). Individual Anle138b isomer ligands are numbered: 1—red, 2—yellow, 3—green, 4—blue, 5—violet.
Ijms 23 16096 g007
Figure 8. Five Anle138b isomer ligand binding sites on micelle-bound AαSyn (PDB ID 1XQ8). Individual Anle138b isomer ligands are as follows: 1—red, 2—yellow, 3—green, 4—blue, 5—violet.
Figure 8. Five Anle138b isomer ligand binding sites on micelle-bound AαSyn (PDB ID 1XQ8). Individual Anle138b isomer ligands are as follows: 1—red, 2—yellow, 3—green, 4—blue, 5—violet.
Ijms 23 16096 g008
Figure 9. Five binding sites of the Anle138b isomer ligand of the complex of the AαSyn fragment and cyclophilin (PDB ID 6I42). Individual Anle138b isomer ligands have the following: 1—red, 2—yellow, 3—green, 4—blue, 5—violet.
Figure 9. Five binding sites of the Anle138b isomer ligand of the complex of the AαSyn fragment and cyclophilin (PDB ID 6I42). Individual Anle138b isomer ligands have the following: 1—red, 2—yellow, 3—green, 4—blue, 5—violet.
Ijms 23 16096 g009
Figure 10. Contribution of target protein amino acids to the calculated binding energies (kJ/mol) with the Anle138b isomer ligand: (a) AαSyn–CypA; (b) CypA dimer; (c) dimer of the AαSyn–CypA complex; (d) AαSyn filament. In the figures, amino acids marked in yellow make a significant contribution to the stabilization of the complex according to the results of molecular docking experiments.
Figure 10. Contribution of target protein amino acids to the calculated binding energies (kJ/mol) with the Anle138b isomer ligand: (a) AαSyn–CypA; (b) CypA dimer; (c) dimer of the AαSyn–CypA complex; (d) AαSyn filament. In the figures, amino acids marked in yellow make a significant contribution to the stabilization of the complex according to the results of molecular docking experiments.
Ijms 23 16096 g010
Figure 11. Amino acid sequence of AαSyn. Disorder-promoting amino acid residues are marked in red; pink blocks reflect clusters rich in disorder-promoting residues. Rectangles mark amino acid residues involved in the stabilization of the Anle138b isomer–AαSyn and Anle138b isomer–AαSyn–CypA complexes (according to molecular docking data).
Figure 11. Amino acid sequence of AαSyn. Disorder-promoting amino acid residues are marked in red; pink blocks reflect clusters rich in disorder-promoting residues. Rectangles mark amino acid residues involved in the stabilization of the Anle138b isomer–AαSyn and Anle138b isomer–AαSyn–CypA complexes (according to molecular docking data).
Ijms 23 16096 g011
Figure 12. Structural formulas 2D (a) and 3D (b) of the Anle138b isomer; (c) electrostatic potential distribution for the Anle138b isomer molecule (according to the PM3 calculation, MOPAC 2021).
Figure 12. Structural formulas 2D (a) and 3D (b) of the Anle138b isomer; (c) electrostatic potential distribution for the Anle138b isomer molecule (according to the PM3 calculation, MOPAC 2021).
Ijms 23 16096 g012
Table 1. Interactions of five Anle138b isomer ligand molecules with a filament of AαSyn (PDB ID 7V4C).
Table 1. Interactions of five Anle138b isomer ligand molecules with a filament of AαSyn (PDB ID 7V4C).
InteractionsResidueAADistance, ÅAngleLigand AtomAffinity (kcal/mol)
First ligand
Hydrophobic66AVal3.8695−5.8
66AVal3.6095
69 AAla3.7795
71 AVal3.61115
71 AVal3.87108
Hydrogen bonds50 AHis3.35121.47
52 AVal2.18156.21
52 AVal1.84157.17
Second ligand
Hydrophobic53 AAla3.5884−6.9
55 AVal3.9491
71 AVal3.6375
71 AVal3.7776
Hydrogen bonds50 AHis2.50113.40
52 AVal2.30130.98
54 AThr3.00136.63
Third ligand
Hydrophobic47 AGly3.40122.20−6.3
72 AThr2.41114.07
72 AThr2.91105.00
Fourth ligand
Hydrophobic52 AVal3.4845−6.3
66 AVal3.6045
72 AThr3.9230
Hydrogen bonds70 AVal2.13149.27
Fifth ligand
Hydrophobic46 AGlu3.533−5.7
49 AVal3.6423
Hydrogen bonds48 AVal2.11163.86
π-cation50 AHis4.1617–21
Table 2. Interactions of five Anle138b isomer ligand molecules with micelle-bound AαSyn (PDB ID 1XQ8).
Table 2. Interactions of five Anle138b isomer ligand molecules with micelle-bound AαSyn (PDB ID 1XQ8).
InteractionsResidueAADistance, ÅAngleLigand AtomAffinity (kcal/mol)
First ligand
Hydrophobic40 AVal3.94108−6.7
43 ALys3.84114
Hydrogen bonds32 ALys3.40111.29
40 AVal3.36149.31
40 AVal2.61108.87
43 ALys2.56112.70
π-stacking39 ATyr4.7178.5394–99
π-cation43 ALys5.48107–109, 111, 114, 115
43 ALys5.43109–113
Second ligand
Hydrophobic33 AThr3.7991−5.8
Hydrogen bonds32 ALys2.47138.24
36 AGly2.51128.51
Third ligand
Hydrophobic43 ALys3.9662−6.2
48 AVal3.7968
48 AVal3.9662
52 AVal3.9068
Hydrogen bonds32 ALys2.37129.57
52 AVal3.29112.57
Fourth ligand
Hydrophobic45 ALys3.766−5.5
48 AVal3.556
49 AVal3.997
Table 3. Interactions of five Anle138b isomer ligand molecules with complex AαSyn–CypA (PDB ID 6I42).
Table 3. Interactions of five Anle138b isomer ligand molecules with complex AαSyn–CypA (PDB ID 6I42).
InteractionsResidueAADistance, ÅAngleLigand AtomAffinity (kcal/mol)
First ligand
Hydrophobic81 AGlu3.58107−7.1
82 ALys3.63108
107 AThr3.6896
Hydrogen bonds81 AGlu2.07171.41
82 ALys3.13146.9
82 ALys2.43119.68
π-stacking50 BHis4.0015.64107–111, 114, 115
Second ligand
Hydrophobic76 ALys3.9273−7.0
76 ALys3.7472
Hydrogen bonds68 AThr2.33128.37
75 AGly3.49108.72
75 AGly2.54109.73
Halogen bonds80 AGly3.72150.45
Third ligand
Hydrophobic46 APhe3.7952−7.0
49 ALys3.2462
67 APhe3.4752
Hydrogen bonds47 AGly2.94115.96
47 AGly2.73111.32
49 ALys2.52129.52
49 ALys2.94108.80
π-cation49 ALys4.7463–97
Fourth ligand
Hydrophobic26 AAla3.7445−6.7
31 ALys3.44625
79 ATyr3.5729
π-cation37 AArg5.3640–44
Fifth ligand
Hydrophobic88 APhe3.724−6.6
105 APro3.482
125 ALys3.9415
125 ALys3.6223
125 ALys3.5914
Hydrogen bonds85 AAsp3.18121.95
102 AAsn3.26127.71
104 AGly3.06154.03
Table 4. Free interaction energies for studied ligand–target variant complexes (MM-PBSA analysis).
Table 4. Free interaction energies for studied ligand–target variant complexes (MM-PBSA analysis).
Ligand Anle138b and TargetVan der Waal Energy (kJ/mol)Electrostatic Energy (kJ/mol)Polar Solvation Energy (kJ/mol)SASA Energy (kJ/mol)Binding Energy (kJ/mol)
AαSyn–CypA complex (PDB ID 6I42)−67.96 ± 34.5−26.05 ± 29.653.2 ± 34.47−7.7 ± 3.85−48.44 ± 26.11
CypA dimer−151.56 ± 18.9−34.54 ± 13.6102.59 ± 18.1−15 ± 1.33−98.53 ± 22.74
Dimer of the AαSyn–CypA complex−91.36 ± 33.8−16.58 ± 2.371.63 ± 10.2−9.87 ± 3.36−46.11 ± 24.69
AαSyn (filament, PDB ID 7V4C)−140.93 ± 47.68−14.66 ± 578.87 ± 28.95−14.37 ± 3.32−91 ± 27.29
Table 5. Physicochemical properties of the Anle138b isomer.
Table 5. Physicochemical properties of the Anle138b isomer.
LigandIUPAC (Formula)ClassMw, DaDonor HB *Acceptor HBSASA **, Å2Dipole Moment, DebaySizesClinical PhaseIC50, nM ***Ref.
Anle138b isomer(3-(1,3-benzodioxol-5-yl)-5-(3-bromophenyl)-1H-pyrazolidinepyrazolidine347.2124313.741.6113.514 Å 6.402 Å
4.381 Å
Preclinical900[51]
* hydrogen bond; ** solvent accessible surface area; *** H-SY5Y neuroblastoma cell model of alpha-synucleinopathy.
Table 6. Structures of AαSyn and CypA used for molecular docking with Anle138b isomer.
Table 6. Structures of AαSyn and CypA used for molecular docking with Anle138b isomer.
TargetPDB ID
(WebLink)
Uniprot IDShape of ProteinResolution, ÅLigand (Locus, a.a.)MethodLength, a.a.
AαSyn (micelle-bounded)1XQ8
(https://www.rcsb.org/structure/1XQ8, accessed on 25 November 2022)
P37840IDP *NMR **140
AαSyn (filament)7V4C
(https://www.rcsb.org/structure/7V4C, accessed on 25 November 2022)
P378403.30EM ***140 ****
CypA6I42
(https://www.rcsb.org/structure/6I42, accessed on 25 November 2022)
P62937Globular1.38AαSyn (1–13)X-ray164
IDP * intrinsically disordered protein; NMR ** nuclear magnetic resonance; EM *** electron microscopy; 140 **** homotrimer.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kondratyev, M.S.; Rudnev, V.R.; Nikolsky, K.S.; Petrovsky, D.V.; Kulikova, L.I.; Malsagova, K.A.; Stepanov, A.A.; Kopylov, A.T.; Kaysheva, A.L. In Silico Study of the Interactions of Anle138b Isomer, an Inhibitor of Amyloid Aggregation, with Partner Proteins. Int. J. Mol. Sci. 2022, 23, 16096. https://doi.org/10.3390/ijms232416096

AMA Style

Kondratyev MS, Rudnev VR, Nikolsky KS, Petrovsky DV, Kulikova LI, Malsagova KA, Stepanov AA, Kopylov AT, Kaysheva AL. In Silico Study of the Interactions of Anle138b Isomer, an Inhibitor of Amyloid Aggregation, with Partner Proteins. International Journal of Molecular Sciences. 2022; 23(24):16096. https://doi.org/10.3390/ijms232416096

Chicago/Turabian Style

Kondratyev, Maxim S., Vladimir R. Rudnev, Kirill S. Nikolsky, Denis V. Petrovsky, Liudmila I. Kulikova, Kristina A. Malsagova, Alexander A. Stepanov, Arthur T. Kopylov, and Anna L. Kaysheva. 2022. "In Silico Study of the Interactions of Anle138b Isomer, an Inhibitor of Amyloid Aggregation, with Partner Proteins" International Journal of Molecular Sciences 23, no. 24: 16096. https://doi.org/10.3390/ijms232416096

APA Style

Kondratyev, M. S., Rudnev, V. R., Nikolsky, K. S., Petrovsky, D. V., Kulikova, L. I., Malsagova, K. A., Stepanov, A. A., Kopylov, A. T., & Kaysheva, A. L. (2022). In Silico Study of the Interactions of Anle138b Isomer, an Inhibitor of Amyloid Aggregation, with Partner Proteins. International Journal of Molecular Sciences, 23(24), 16096. https://doi.org/10.3390/ijms232416096

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop