Next Article in Journal
De Novo Transcriptome Sequencing of Rough Lemon Leaves (Citrus jambhiri Lush.) in Response to Plenodomus tracheiphilus Infection
Next Article in Special Issue
Carboxypeptidase U (CPU, TAFIa, CPB2) in Thromboembolic Disease: What Do We Know Three Decades after Its Discovery?
Previous Article in Journal
Real-Time Fluorescence Image-Guided Oncolytic Virotherapy for Precise Cancer Treatment
Previous Article in Special Issue
Venous Thrombosis and Thrombocyte Activity in Zebrafish Models of Quantitative and Qualitative Fibrinogen Disorders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

NMR-Based Structural Characterization of a Two-Disulfide-Bonded Analogue of the FXIIIa Inhibitor Tridegin: New Insights into Structure–Activity Relationships

by
Thomas Schmitz
1,
Ajay Abisheck Paul George
1,2,
Britta Nubbemeyer
1,
Charlotte A. Bäuml
1,
Torsten Steinmetzer
3,
Oliver Ohlenschläger
4,
Arijit Biswas
5 and
Diana Imhof
1,*
1
Pharmaceutical Biochemistry and Bioanalytics, Pharmaceutical Institute, University of Bonn, An der Immenburg 4, D-53121 Bonn, Germany
2
BioSolveIT GmbH, An der Ziegelei 79, D-53757 Sankt Augustin, Germany
3
Institute of Pharmaceutical Chemistry, Philipps University of Marburg, Marbacher Weg 6, 35032 Marburg, Germany
4
Leibniz Institute on Aging—Fritz-Lipmann-Institute, Beutenbergstr. 11, D-07745 Jena, Germany
5
Institute of Experimental Hematology and Transfusion Medicine, University Hospital Bonn, Sigmund-Freud-Str. 25, D-53127 Bonn, Germany
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(2), 880; https://doi.org/10.3390/ijms22020880
Submission received: 23 December 2020 / Revised: 10 January 2021 / Accepted: 12 January 2021 / Published: 17 January 2021
(This article belongs to the Special Issue Fibrinogen/Fibrin, Factor XIII and Fibrinolysis in Diseases)

Abstract

:
The saliva of blood-sucking leeches contains a plethora of anticoagulant substances. One of these compounds derived from Haementeria ghilianii, the 66mer three-disulfide-bonded peptide tridegin, specifically inhibits the blood coagulation factor FXIIIa. Tridegin represents a potential tool for antithrombotic and thrombolytic therapy. We recently synthesized two-disulfide-bonded tridegin variants, which retained their inhibitory potential. For further lead optimization, however, structure information is required. We thus analyzed the structure of a two-disulfide-bonded tridegin isomer by solution 2D NMR spectroscopy in a combinatory approach with subsequent MD simulations. The isomer was studied using two fragments, i.e., the disulfide-bonded N-terminal (Lys1–Cys37) and the flexible C-terminal part (Arg38–Glu66), which allowed for a simplified, label-free NMR-structure elucidation of the 66mer peptide. The structural information was subsequently used in molecular modeling and docking studies to provide insights into the structure–activity relationships. The present study will prospectively support the development of anticoagulant-therapy-relevant compounds targeting FXIIIa.

1. Introduction

Leeches, such as the medical leech Hirudo medicinalis, have been used as a biomedical tool from ancient times up to the present day. They were mainly used to remove blood, which should help to detoxify the patient [1]. Because of these bloodsucking properties and the intent to avoid clot formation during blood intake, the leeches’ saliva contains a variety of anticoagulants that can be used for the control and treatment of hemostasis. Over 20 different compounds with antithrombotic activity have already been isolated from the salivary glands of various leeches, showing a direct or indirect influence on the blood coagulation cascade [2,3]. One of these compounds is the well-known thrombin inhibitor hirudin, which is isolated from the leech Hirudo medicinalis [4,5]. Based on this natural compound, anticoagulants such as bivalirudine and argatroban have already been successfully developed for the treatment of cardiovascular and thromboembolic diseases [6,7,8]. However, anticoagulants that directly or indirectly inhibit thrombin frequently increase the risk of undesired bleeding, making them unsuitable for many therapeutic applications [9,10]. Since thrombin catalyzes the penultimate step of the blood coagulation cascade—the formation of fibrin polymers starting from fibrinogen—an influence on any of the previously required coagulation factors always indirectly affects thrombin [11]. Therefore, one way to achieve a thrombin-independent impact on the blood coagulation cascade is the inhibition/activation of the final step, the crosslinking of the fibrin polymers catalyzed by FXIIIa [12,13].
FXIIIa is a transglutaminase that forms isopeptide bonds within the fibrin polymer by linking a glutamine side chain and a lysine side chain, providing increased blood clot stability. In this regard, the isopeptide bond formation is catalyzed by a catalytic triad composed of Cys314, His373, and Asp396 and a catalytic dyad (His342, Glu401) located in the active side of FXIIIa [12,14,15,16]. Inhibiting this factor leads to unstable blood clots that are smaller in size and can be easily dissolved by fibrinolysis [17,18,19,20]. Consequently, FXIIIa serves as an interesting target for the development of anticoagulants in order to circumvent the problem of undesired bleeding. So far, several small-molecule inhibitors such as cerulenin and alutacenoic acid A have been investigated for their inhibitory potential against FXIIIa [21,22,23,24]. However, many of these small-molecule inhibitors possess a low selectivity and/or half-life in plasma. Thus, the development of these inhibitors has not been pursued further. Nowadays, the main focus is directed toward peptidic and allosteric FXIIIa inhibitors, such as the peptidomimetic inhibitor ZED3197 and glucosaminoglycan-derived inhibitors, e.g., NSGM 13 (non-saccharide glucosaminoglycan mimetics), introduced by Al-Horani et al. [23,24,25,26]. Furthermore, the peptidic inhibitor tridegin, a 66mer peptide, first isolated from the giant Amazon leech Haementeria ghilianii in 1997, has been investigated for its high inhibitory potential against FXIIIa [23,27]. Tridegin belongs to the class of cysteine-rich peptides due to its 6 cysteines within the amino acid sequence. The natural disulfide bond connectivity of tridegin could not be elucidated so far [28,29,30,31], which is due to the fact that a peptide with 6 cysteines has the ability to form 15 different three-disulfide-bonded isomers. So far, final conclusions have not been drawn concerning the folding of tridegin to the bioactive conformation and whether it is folded in a BPTI-like (Bovine pancreatic trypsin inhibitor) fashion, according to which one isomer is preferentially formed, or in a hirudin-like manner, where a maximum of 15 possible three-disulfide-bonded isomers are formed [32,33]. A folding experiment performed by Böhm et al. in 2014 provided first hints regarding the preferentially formed disulfide-bonded isomers of tridegin [30]. Therein, three of the 15 possible three-disulfide-bonded isomers were identified, which all shared the disulfide bond between Cys19–Cys25 and possessed similar activities toward FXIIIa [28,30]. In addition to the three-disulfide-bonded tridegin isomers, several variants lacking the disulfide bond between Cys19–Cys25 were synthesized in a following approach in order to determine the significance of this bond for their inhibitory potential against FXIIIa [29]. The results revealed that the activity against FXIIIa was not lost, despite the missing disulfide bond. Due to the unchanged activity and the improved synthetic outcome, the two-disulfide-bonded tridegin analogues, in particular the isomer B[C19S,C25S] with the disulfide bonds Cys5–Cys37 and Cys17–Cys31, are used as a lead structure for further research on tridegin and its influence on FXIIIa [29].
The aim of the current study is to gain further insight into tridegin’s structure, folding, and interaction with FXIIIa. In order to obtain a more precise understanding of the mode of action against FXIIIa, a structure of the tridegin isomer B[C19S,C25S] was determined by a combination of NMR (nuclear magnetic resonance) spectroscopy and molecular dynamics (MD) simulation studies. This two-disulfide-bonded isomer, B[C19S,C25S], was subsequently applied for molecular docking studies on FXIIIa. The information received from this study was used to refine structure–activity relationships of the tridegin–FXIIIa interaction. Our results may pave the way for further investigations on FXIIIa as well as the optimization/development of lead structures for antithrombotic and thrombolytic therapy.

2. Results and Discussion

2.1. Peptide Synthesis and Characterization

To elucidate the structure of the two-disulfide-bonded isomer B[C19S,C25S] (Figure 1a) [29], two segments derived from the peptide sequence were synthesized for NMR measurements. The reason for dividing the tridegin isomer into an N-terminal and a C-terminal part was the size of the 66mer peptide and the high number of proline residues within the sequence, which complicate the assignment of the NMR signals and hamper structure calculation. Therefore, a C-terminal fragment (amino acids 38–66) already described by Böhm et al. [31] and a new N-terminal fragment (amino acids 1–37) with the disulfide bonds Cys5–Cys37 and Cys17–Cys31 were produced (Supplementary Table S1, Supplementary Figure S1). For the synthesis of the latter fragment, a solid-phase peptide synthesis using the Fmoc-strategy with acetamidomethyl (Acm)- and trityl (Trt) -protecting groups (Cys5, Cys37 and Cys17, Cys31, respectively) and a subsequent stepwise oxidation in solution was conducted, similar to the strategy used for the earlier described two-disulfide-bonded tridegin analogues [29]. The stepwise oxidation process was performed in a one-pot reaction with iodine (oxidizing agent), whereby the former Trt-protected cysteines, were deprotected during the cleavage of the peptide from the solid support and initially linked, followed by the cleavage of the Acm-protected cysteines (Figure 1b). The selective disulfide bond formation during the synthesis of the two disulfide bonds was controlled by increasing the amount of iodine as described previously [29] and by adjusting the acetic acid concentration. For the disulfide bond formation of the cysteines formerly protected by Trt, the first oxidation was performed in 100% acetic acid, whereas the deprotection of the Acm-protected cysteines and the formation of the second disulfide bond was completed in 70% acetic acid. Due to the higher acid concentration in the first oxidation, the acid-stable Acm-protected cysteines are additionally prevented from unwanted cleavage, thus ensuring the stepwise formation of the disulfide bonds [34]. To confirm the correct disulfide bond connectivity, the N-terminal fragment was digested using chymotrypsin, and the resulting fragments were analyzed subsequently by MS/MS as earlier described [29]. Several disulfide-bonded fragments representative for the disulfide bonds Cys5–Cys37 and Cys17–Cys31 were detected, confirming the correct disulfide bond connectivity of the N-terminal fragment (Figure 1c, Supplementary Table S2).
In addition, activity studies of the N-terminal fragment on FXIIIa were performed using the established FXIIIa isopeptidase activity assay [28,29,30]. As already observed in the case of a three-disulfide-bonded N-terminal fragment of tridegin [30], no inhibitory potential against FXIIIa was found. However, as described earlier, it is expected that the N-terminal fragment enhances the binding between tridegin and FXIIIa. In contrast, the C-terminal fragment exhibits inhibitory activity toward FXIIIa [30].

2.2. Structure Elucidation by NMR Spectroscopy

The structure of isomer B[C19S,C25S] [29] was determined by combining the two structures of the peptide fragments. First, these structures were elucidated by natural abundance 2D NMR spectroscopy, recording TOCSY (total correlation spectroscopy), NOESY (nuclear Overhauser enhancement spectroscopy), COSY (correlated spectroscopy), and 13C-HSQC (heteronuclear single quantum coherence) spectra and calculated by Cyana [35] based on the distance and dihedral angle restraints. Within the signal assignment of the NMR spectra, the spectral resolution of the 2D NMR spectra (TOCSY and NOESY) at a magnetic field of 16.4 T allowed a clear resonance assignment for both peptide fragments (Figure 2, Supplementary Table S3).
The TOCSY, COSY, and HSQC spectra of the C-terminal fragment showed a high resolution and well-separated signals, which could be easily associated with the individual spin systems of the 29 amino acids. A subsequent “sequential walk” within the HN–Hα region of the NOESY and TOCSY spectra (Figure 2a) resulted in a sequence-specific assignment of the signals, despite the difficulty of the occurrence of six prolines. Proline residues lead to an interruption of the sequential walk due to the lack of HN. The assignments within the NOESY spectrum were used to solve the solution structures of the C-terminal fragment via Cyana [35]. The NOESY spectrum already revealed few correlations between amino acids that were separated by more than two amino acids within the sequence. This resulted in a very small number of distance restraints for the structure calculation, which was reflected in a flexible structure of the C-terminal fragment (Figure 3, Supplementary Table S4) with a backbone root mean square deviation (RMSD) of 8.54 ± 2.83 Å.
The NMR spectra of the N-terminal fragment contained more signals due to the higher number of amino acids (37 amino acids), which were successfully attributed to the amino acid positions by assigning the amino acid spin system and subsequent sequential walk within the HN–Hα range (Figure 2b), as it was already performed for the C-terminal fragment. The two disulfide bonds as well as possible secondary structure elements within the N-terminal fragment result in a more rigid structure compared to the C-terminal fragment, enabling NOESY correlations between distant amino acids within the peptide sequence. In this context, NOESY correlations between the amino acids Trp8–Gln36, Pro15–Phe33, Cys17–Phe33, Arg16–Tyr30, Trp18–Cys31, and Ser19–Tyr30 confirmed the disulfide bond connectivity Cys5–Cys37 and Cys17–Cys31 of the N-terminal fragment, which was already verified by the aforementioned MS/MS analysis of the chymotryptic digest. The numerous NOESY correlations led to a higher set of distance restraints, which, in addition to the dihedral angle restraints, were included in the calculation of the structure via Cyana [35]. In contrast to the C-terminal fragment, three secondary structure elements stabilized by the disulfide bonds Cys5–Cys37 and Cys17–Cys31 were identified within the N-terminal fragment (RMSD of residues 5–37: 1.20 ± 0.39 Å), e.g., one helix (Glu24–Tyr30) and two β-turns (His9–Ile12, Trp18–Ala21) (Figure 3, Supplementary Table S4). Especially the fact that the two secondary structure elements in the region Trp18–Tyr30 result in a close proximity of the serine residues at position 19 and 25 is remarkable. This may indicate that the disulfide bond between Cys19–Cys25, which occurs in the three-disulfide-bonded tridegin isomers [30], is formed as a consequence of the secondary structure elements. Regarding this, it may be concluded that during the folding process, this bond is only established after the formation of the secondary structure elements.
To determine the isomer B[C19S,C25S] structure [29], the distance and dihedral angle restraints of the N- and the C-terminal fragments were combined and the complete structure ensemble (100 structures, Figure 4) of the tridegin variant was calculated using Cyana [35]. Thereby, the determined secondary structure elements of the N-terminal fragment as well as the flexibility of the C-terminal fragment were preserved. The possibility of the C-terminal fragment (residue 38–66) to undergo large conformational changes was characterized by an RMSD of 14.14 ± 2.76 Å for the combined isomer B[C19S,C25S] structure [29] with reference to the lowest energy structure (Figure 4a). Additionally, within the structure ensemble of isomer B[C19S,C25S] [29] a massive proportion of this observed variation (12.15 ± 3.27 Å) resulted from the flexible C-terminal fragment (residues 38–66) while the N-terminal fragment alone (residues 1–37) had a backbone RMSD of 3.77 ± 1.42 Å. However, since the calculated structure ensemble of isomer B[C19S,C25S] [29] is only a combination of the independently resolved NMR structures of the N- and C-terminal fragment, the structures were further investigated via molecular dynamics simulations.

2.3. MD-Based Analysis of Isomer B[C19S,C25S]

For the investigation of the aforementioned isomer B[C19S,C25S] structure [29], the model with the lowest energy from the 100-member NMR ensemble of the calculated isomer B[C19S,C25S] structures [29] was subjected to a 300-ns-long all-atom molecular dynamics. The prime motivation behind this effort was to produce a solvent-equilibrated ensemble of the isomer B[C19S,C25S] structure [29], which was originally produced by combining the independently resolved fragments (N-terminal and C-terminal fragments, as described above) determined by NMR spectroscopy. Additionally, the final structure from this MD simulation was also used as an input for molecular docking experiments on the crystal structure of FXIIIa (PDB ID: 4KTY) [16].
As already discussed in Section 2.2, analyses conducted of the NMR ensemble of the combined (N-terminal and C-terminal fragments) isomer B[C19S,C25S] structure [29] (Figure 4a) strongly indicated that the elongated C-terminal segment could potentially undergo large conformational changes upon further equilibration to attain its fully folded conformation. It was therefore clear that, prior to molecular docking experiments, MD-based refinement of this structure was necessary. To this end, a 300-ns-long all-atom MD simulation was performed on the NMR structure with the lowest energy (Figure 4b), and the resulting trajectory was analyzed for structural changes over the course of the simulation. A large conformational rearrangement marked by an 18.50 Å change in backbone RMSD was observed within the first 15 ns of the 300 ns simulation (Figure 4c). This was characterized by a rapid collapse of the flexible C-terminal fragment (residues 38–66) around the disulfide-bonded N-terminal segment (residues 1–37). Remarkably, after this initial dramatic change, the remainder of the simulation saw the structure steadily equilibrate deviating only by an average of 3.33 ± 0.47 Å with respect to the structure at 15 ns. The large deviation between the NMR starting structure and the final simulation structure can clearly be explained due to the fact that, in the NMR structure, the N- and C-terminal segments were resolved as separate entities and were combined as a final step to create the full-length isomer B[C19S,C25S] structure [29]. We believe that our simulation captures the full extent of conformational changes the structure would have undergone if the N-terminal and C-terminal segments could possibly have been resolved as a single entity via NMR spectroscopy and that the structures obtained from any part of the trajectory after the initial 15 ns would surely fall within the low-energy conformations of isomer B[C19S,C25S] [29] in its physiological environment. To further confirm this, we show that the conformational ensemble for the entire structure produced via the MD simulation (containing 2850 snapshots) has an RMSD of 3.33 ± 0.47 Å, which is comparable to the RMSD (3.77 ± 1.42 Å) of the 100-member NMR ensemble of the N-terminal segment of isomer B[C19S,C25S] [29]. We also note that the α-helical motif seen in the NMR starting structure (Glu24–Tyr30) was consistently found through the production phase of the MD simulation. The final snapshot of this 300 ns simulation was used as an input structure for subsequent molecular docking simulations. A video of this MD simulation is provided as supplementary content. Structural alignments of the final snapshot of the simulation with a previously computationally modeled version of the same peptide [29] using the MUSTANG algorithm [36] revealed that although the general fold was preserved, the NMR-based model (current study) varied from the earlier computational model (Figure 5b–d) [29]. This is understandable based on the fact that the earlier model was produced by mutating its three-disulfide-bonded counterpart, which again had its origins in computational modeling. Interestingly though, the parts of the C-terminal segment, seemed to align well against each other (Figure 5b–d). The current produced model has the clear advantage that its origin is from NMR spectroscopy. As a means to confirm that this MD-derived structure from the 300 ns simulation was reproducible, a second, independent 1000-ns-long simulation was conducted. The equilibrated structures from this longer simulation aligned well with the structure produced by the 300 ns simulation in this study (Supplementary Figure S2).

2.4. Molecular Docking of Full-Length Tridegin Variant on FXIIIa

Docking of the full-length model of isomer B[C19S,C25S] [29] (derived from the final 300 ns simulation snapshot) (Figure 5a) onto FXIIIa resulted in one pose amongst the top-ten docking poses that satisfies the inhibitory characteristic of tridegin (i.e., it binds next to the catalytic site) (Supplementary Figure S3). All remaining docks cluster around the β-barrel domain of FXIIIa, which is located significantly far from the FXIIIa catalytic site, and are therefore irrelevant to our analysis. On an atomic/residual level, the important docking pose of isomer B[C19S,C25S] [29] shows the C-terminal region (residues 53–60) participating in several non-covalent interactions in and around the catalytic site (Figure 6, Supplementary Table S5–S7). The original non-proteolytically activated crystal structure from Stieler et al. [16] presents a hydrophobic tunnel at the active site formed by the aromatic interactions between the Trp279 and Trp370 residues of FXIIIa (Figure 6). The two ends of this tunnel serve as the entry points for the FXIIIa substrates, lysine and glutamine, through which they reach the catalytic Cys314 at the bottom of the cavity that the tunnel is formed on top of. The original non-proteolytically activated crystal structure of FXIIIa (PDB ID: 4KTY [16]) was stabilized using the FXIIIa inhibitor ZED1301, which irreversible bound to the catalytic Cys314, thereby occupying the hydrophobic tunnel and preventing access for competitive substrate [16]. The mode of inhibition for the full-length model of isomer B[C19S,C25S] [29] is very similar to that of ZED1301 since residues from the C-terminal region not only participate in hydrophobic interactions with FXIIIa residues that would form the hydrophobic tunnel, but some tridegin residues such as Arg53 and Phe57 also stabilize the C-terminal region by interacting with the critical Trp279 and Trp370 FXIIIa residues (Supplementary Figure S4, Supplementary Table S5–S7). Therefore, the C-terminal region appears to competitively occupy the substrate-entry hydrophobic tunnel, thereby preventing substrate access to the catalytic Cys314 similar to ZED1301. An improvement compared to the inhibitor ZED1301 might be that the model structure of isomer B[C19S,C25S] [29] also potentially blocks the catalytic dyad of FXIIIa (His342, Glu401) which is responsible for incorporating the substrate lysine into the hydrophobic tunnel. This is suspected since the inhibitor ZED1301 does not reach the catalytic dyad, whereas the tridegin analogue rests on the residues His342 and Glu401. In addition to the ZED1301 comparison, the docking results of the NMR-based model structure was aligned to a previously reported docking with the computationally modeled version of the same peptide (Supplementary Figure S5) [29]. In this regard, both dockings showed a similar binding into the hydrophobic tunnel of FXIIIa, despite the overall structural variations between the two (Figure 5b–d). However, there is one significant difference between the two complexes, which is the orientation of the C-terminal part within the hydrophobic pocket. While the NMR-based structure guide through the hydrophobic tunnel showed that the C-terminus of the peptide is on the Q-site (site of the tunnel where glutamine normally binds) of FXIIIa, the C-terminus of the previously reported structure [29] was observed on the K-site (site where lysine binds) close to the catalytic dyad. The differences in the orientation can be attributed to two main things: (a) firstly the differences in the N-terminal structure (residue 1–37) in both isomer B[C19S,C25S] models [29], whose role in this interaction is to orient the C-terminal region within the FXIIIa hydrophobic pocket, and (b) methodological differences in the implementation of docking between the current and previous study. The similarity in binding location though further underlines the importance of the C-terminal fragment in binding and subsequent inhibition of FXIIIa.
Although the Arg16 residue from the tridegin isomer B[C19S,C25S] N-terminal part [29] of the NMR-based model also appears to interact with FXIIIa in this docking pose, the effect of the N-terminal region (residue 1–37) might only be complimentary. The interaction may only play a role in the initial establishment of contact between isomer B[C19S,C25S] [29] and FXIIIa but is unlikely to play any role in the subsequent inhibition. This was confirmed by a docking were only the N-terminal fragment was docked to FXIIIa. Similar to isomer B[C19S,C25S] [29], when we dock only the NMR-based N-terminal fragment structure of the tridegin analogue on FXIIIa, one docking pose amongst top-ten docked close to the catalytic site. However, unlike the full-length isomer B[C19S,C25S] [29] model, the N-terminal fragment structure binds only superficially on top of the catalytic site (Supplementary Figure S6). In this docking pose, neither the hydrophobic tunnel nor any hydrophobic access is disrupted, which is also a hint for the reduced influence of the N-terminal part on the inhibition to FXIIIa.

3. Materials and Methods

3.1. Materials

Fmoc-amino acids were purchased from Orpegen Peptide Chemicals (Heidelberg, Germany) and Novabiochem (Schwalbach, Germany). HBTU, resins, and further chemicals for solid-phase peptide synthesis, including reagent-grade N-methylmorpholine, piperidine, trifluoroacetic acid (TFA), and N,N-dimethylformamide (DMF) were purchased from IRIS Biotech (Marktredwitz, Germany), Sigma-Aldrich Chemie GmbH (Munich, Germany), Alfa Aesar (Karlsruhe, Germany), Abcr GmbH (Karlsruhe, Germany), VWR International (Darmstadt, Germany), and FLUKA Chemika (Seelze, Germany). Solvents (acetonitrile, water, methanol, ethyl acetate, diethyl ether) and chemicals (iodine, acetic acid) used for peptide purification and oxidation were obtained from VWR International (Darmstadt, Germany) and Fisher Scientific (Schwerte, Germany).

3.2. Peptide Synthesis and Purification

The synthesis of the linear peptide precursors according to an automated standard Fmoc SPPS protocol using a ResPep SL peptide synthesizer from Intavis Bioanalytical Instruments GmbH (Cologne, Germany) as well as the subsequent peptide cleavage from the resin and the purification step via RP-HPLC were carried out as previously described [29].
The purity of the linear peptides was >95%, which was assessed by analytical RP-HPLC on a Shimadzu LC-20AD system equipped with a Vydac 218TP column (C18, 250 × 4.6 mm, 5 µm particle size, 300 Å pore size). The gradient elution system contained 0.1% TFA in water (eluent A) and 0.1% TFA in acetonitrile (eluent B). Peptide elution was achieved with a gradient of 20% to 50% eluent B in 30 min and a flow rate of 1 mL/min. The peptides were detected at λ = 220 nm.

3.3. Selective Oxidation of the N-Terminal Fragment

The selective oxidation of the two-disulfide-bonded N-terminal peptide Lys1–Cys37 proceeded in a one-pot reaction slightly modified to the already described protocol [29]. In brief, the linear precursors (1 eq.) were dissolved in 100% acetic acid to yield a final concentration of 0.05 mM. The oxidation of the first disulfide bond with 1.1 eq. iodine (0.1 M in methanol) was stirred at room temperature and under argon atmosphere. After 15 min, 13.9 eq. iodine (0.1 M in methanol) and water was added to the reaction mixture in order to decrease the acetic acid concentration to 70%. The reaction was stirred for a further 43 h (second oxidation). The oxidation reaction was stopped by the removal of iodine by extraction with the same volume of ethyl acetate (three times). The peptide-containing aqueous solutions were combined, freeze-dried, and purified via RP-HPLC.

3.4. Peptide Characterization and Determination of Disulfide Connectivity

Peptide characterization was achieved by means of analytical RP-HPLC (see above for the linear precursor), mass spectrometry, and amino acid analysis.
The peptide concentration as well as the amino acid composition were analyzed using an LC 3000 system from Eppendorf-Biotronik (Hamburg, Germany). Peptide hydrolysis was carried out in 6 N HCl at 110 °C in sealed tubes for 24 h. Subsequently, the hydrolyzed peptides were dried in a vacuum concentrator and redissolved. The respective amino acid concentrations were determined by comparison with an amino acid standard solution (Laborservice Onken, Gründau, Germany).
Peptide masses were determined using matrix-assisted laser desorption/ionization (MALDI) mass spectrometry. As matrix 2,5-dihydroxyacetophenone (2,5-DHAP) was used according to the manufacture’s guideline (Bruker Daltonics, Bremen, Germany). MALDI mass spectra were produced on an UltrafleXtreme instrument (Bruker Daltonics, Bremen, Germany).
The disulfide connectivity of the N-terminal fragment was determined as described recently [29,30]. Therefore, a micrOTOF-Q III device (Bruker Daltonics, Bremen, Germany) was used to measure electrospray ionization (ESI) mass and tandem-mass spectra.

3.5. Enzyme Activity Assay

The inhibitory potential of the N-terminal tridegin fragment toward FXIIIa was performed in an FXIIIa isopeptidase activity assay using the fluorogenic substrate H-Tyr(3-NO2)-Glu(NH-(CH2)4-NH-Abz)-Val-Lys-Val-Ile-NH2 as described previously [30].

3.6. NMR Spectroscopy and Structure Prediction

The natural abundance NMR experiments for the 1H, 15N, and 13C chemical shift assignments were performed at 293 K in 50 mM sodium phosphate pH 6.3 (95% H2O/5% D2O) on a Bruker Avance III HD 700 MHz Cryo spectrometer. The peptides were measured at a concentration of 3 mM, and the backbone as well as the side chain atoms were assigned via a combination of 2D [1H, 1H]-TOCSY, [1H, 1H]-NOESY, [1H, 1H]-COSY, and [1H, 13C]-HSQC spectra using water suppression. The spectra were processed with TOP SPIN 4.0.6 (Bruker) and analyzed using CcpNmr Analysis (Collaborative Computing Project for NMR). Distance constraints were extracted form [1H,1H]-NOESY spectra acquired with a mixing time of 150 ms and a recycle delay of 1.5 s. Upper-limit-distance constraints were calibrated according to their intensities in the NOESY spectra. To calculate the structure of the fragments based on the chemical shift data the program Cyana [35] was used. Cyana [35] was also used to calculate a hypothetical structure of the full-length tridegin analogue out of the combined distance and dihedral restraints from the NMR experiments of the N- and C-terminal fragments. The 100 structures with the lowest energies were selected to represent the NMR-solution structures.

3.7. Molecular Dynamics (MD) Simulations

Molecular dynamics (MD) simulations in this study were carried out using the Gromacs 2018 package [37,38]. The first model of the 100-member NMR ensemble of the full-length model of isomer B[C19S,C25S] [29] was used as the starting structure of the 300-ns-long simulation. The peptide was placed in the center of a cubic simulation cell, with the edges of the cube separated by at least 12 Å from all atoms of the peptide. The peptide was solvated using the TIP3P water model [39] with additional Na+ and Cl− counter ions added to achieve a physiological salt concentration of 0.9%, while maintaining a zero net charge on the entire system. The solvated system was subjected to 5000 steps of steepest-descents energy minimization with the Amber-ff14sb [40] force field, which was used to describe atomic motions for all parts of the simulation. The energy-minimized system was first subjected to a temperature equilibration in the NVT (constant number of atoms, volume, temperature) ensemble for 2 ns, with temperature maintained by the velocity-rescaled variant of the Berendsen thermostat at 300 K [41]. Subsequently a 2 ns simulation to equilibrate the pressure of the system at 1 atm was conducted in the NPT (constant number of atoms, pressure, temperature) ensemble aided by the Parrinello–Rahman barostat [42,43]. During both the temperature and pressure equilibration runs, all heavy atoms were position-restrained by the LINCS algorithm [44]. The production run was conducted for 300 ns using a 2 fs timestep. Long-range interactions were cutoff at 10 Å, and the electrostatics were described by the particle-mesh Ewald method [45,46]. Periodic boundary conditions were employed in the simulations, the effects of which were adjusted prior to conducting analyses on the trajectory. Trajectory analysis, namely the backbone root mean square deviation (RMSD) computation was conducted in VMD 1.9.3. [47], which was also used in the creation of molecular graphics. Snapshots were written at a 100 ps interval to disc, resulting in a total of 3000 snapshots collected for analysis, from this simulation.

3.8. Blind Docking Studies on the FXIIIa Crystall Structure (PDB ID: 4KTY)

The full-length model of isomer B[C19S,C25S] [29] was used for docking onto the modified and simulation-equilibrated crystal structure of FXIIIa originally derived from the source PDB ID: 4KTY [16]. The PDB structure used as a receptor for docking was a simulation-averaged structure of the non-proteolytically activated FXIIIa, which had been modified, i.e., missing regions/loop filled, all heteroatoms/water molecules removed, and subjected to classical molecular dynamic simulation. The details of this simulation have been reported earlier [15]. The full-length model obtained as the final snapshot of the 300 ns MD simulation (described in Section 3.7) was used as the ligand for the docking simulation. Docking was performed on the Hdock server (http://hdock.phys.hust.edu.cn/) in a blind fashion and under the default conditions of the server [48]. The Hdock server is based on a hybrid algorithm of template-based modeling and ab initio free docking and is currently ranked favorably amongst the topmost automated docking servers in the recent CASP (critical assessment of protein structure prediction) competitions [49]. Only the top-ten docking poses generated by the server were closely inspected. The docking pose that strongly agreed with the experimental data was subjected to a docking refinement protocol on the Haddock 2.2 webserver (https://milou.science.uu.nl/services/HADDOCK2.2/haddockserver-refinement.html) that allows for flexible refinement in explicit solvent of the docked complex [50]. The model with the best Haddock score from the output clusters generated out of the refinement protocol was finally inspected in terms of interatomic interactions on the protein interaction calculator webserver (http://pic.mbu.iisc.ernet.in/) [51]. Similarly docking was also performed with the energetically minimum NMR structure of the well-defined N-terminal fragment as the ligand and the above-described FXIIIa structure as the receptor.

4. Conclusions

In summary, we have experimentally determined for the first time, the structure of a two-disulfide-bonded tridegin variant. This structure contained a flexible C-terminal segment and a more rigid N-terminal part which showed distinct secondary structure elements besides being stabilized by the two disulfide bonds Cys5–Cys37 and Cys17–Cys31. The cysteine-rich N-terminal fragment also shows that the two serine residues Ser19 and Ser25, which are substituted by the cysteines Cys19 and Cys25 in the three-disulfide-bonded tridegin isomers, were in close proximity to a helix and a β-turn secondary structure element, which suggests that the secondary structure elements are responsible for the formation of this specific disulfide bond. Therefore, it can be assumed for the disulfide bond Cys19–Cys25 in the three-disulfide-bonded tridegin isomer that it is generated after the formation of the secondary structure elements within the folding process. With regard to the structure–activity relationships, it can be assumed that the low structural influence of the disulfide bond results in an already experimentally described comparable binding behavior and inhibitory potential to FXIIIa for the three-disulfide-bonded isomer B [28] and the two-disulfide-bonded isomer B[C19S,C25S] [29].
The full-length structure of isomer B[C19S,C25S] [29] generated from the NMR-based results of the N-terminal and C-terminal fragment was subsequently refined by MD simulations and used for molecular docking to FXIIIa. Thereby, it could be demonstrated that the isomer B[C19S,C25S] [29] binds in close proximity to the active site resulting in a blockage of the catalytic Cys314 for substrates as described earlier [29]. The isomer B[C19S,C25S] [29], especially the C-terminal region, penetrates into the hydrophobic tunnel where the active site of FXIIIa is located and competitively inhibits the transglutaminase functionality of FXIIIa by blocking the catalytic triad and catalytic dyad. The interaction with the hydrophobic tunnel of FXIIIa demonstrates how important the flexible C-terminal region is for the inhibition mode of isomer B[C19S,C25S] [29] and that the rigid N-terminal segment only plays a minor role for the inhibition by enhancing the binding to FXIIIa.
In conclusion, the deeper insight into the structure–activity relationships of the FXIIIa inhibitor tridegin from the docking studies will prospectively support the future development of anticoagulant-therapy-relevant compounds targeting FXIIIa.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/22/2/880/s1.

Author Contributions

The study has been designed by D.I., T.S. (Thomas Schmitz) performed synthesis as well as functional and structural analysis of the peptides. B.N. performed the MS measurements and T.S. (Torsten Steinmetzer) provided his expertise regarding the FXIIIa activity assay. T.S. (Thomas Schmitz) and C.A.B. performed the FXIIIa activity assay. T.S. (Thomas Schmitz) and O.O. analyzed the NMR spectra and calculated the peptide structures. A.A.P.G. and A.B. carried out the computational studies. Data analysis and interpretation was carried out by all authors. The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

The Bruker ultrafleXtreme TOF/TOF instrument (to M. Engeser) and the Bruker Avance III HD 700 MHz Cryo spectrometer (to S. Nozinovic) were funded by the University of Bonn, the Ministry of Innovation, Science, and Research of North-Rhine Westfalia and the DFG. This study was financially supported by the German Foundation of Heart Research and by the University of Bonn within university science program HSPII (to D.I.).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We are grateful to D. Menche and S. Nozinovic (University of Bonn) as well as R. Hartmann and D. Willbold (FZ Jülich) for measuring the 2D NMR spectra of the peptides on their instruments. This work was supported by the German Foundation of Heart Research to D.I., which is gratefully acknowledged. The FLI (Fritz-Lipmann-Institute) is member of the Leibniz Association (WGL) and financially supported by the Federal Government of Germany and the State of Thuringia (to O.O.).

Conflicts of Interest

The authors declare no competing financial interest.

Abbreviations

AcmAcetamidomethyl
AcOHAcetic acid
ArgArginine
AspAspartic acid
BPTIBovine pancreatic trypsin inhibitor
CASPcritical assessment of protein structure prediction
COSYCorrelated spectroscopy
CysCysteine
2,5-DHAP2,5-Dihydroxyacetophenone
DMFDimethylformamide
ESIElectrospray ionization
FmocFluorenylmethoxycarbonyl
FXIIIBlood coagulation factor XIII
FXIIIaActivated blood coagulation factor XIII
GluGlutamic acid
HBTU2-(1H-Benzotriazol-1-yl)-1,1,3,3-tetramethyluronium hexafluorophosphate
HSQCHeteronuclear single quantum coherence
HisHistidine
Ile Isoleucine
LysLysine
MALDIMatrix-assisted laser desorption/ionization
MDMolecular dynamics
MSMass spectrometry
MS/MSTandem mass spectrometry
NMRNuclear magnetic resonance
NOESYNuclear Overhauser enhancement spectroscopy
NPTconstant number of atoms, pressure, temperature
NSGMNon-saccharide glucosaminoglycan mimetics
NVTconstant number of atoms, volume, temperature
PhePhenylalanine
ProProline
RMSDRoot mean square deviation
RP-HPLCReversed-phase high-performance liquid chromatography
SerSerine
SPPSSolid phase peptide synthesis
TFATrifluoroacetic acid
TOCSYTotal correlation spectroscopy
TrtTrityl
TrpTryptophane
TyrTyrosine
VMDVisual molecular dynamics

References

  1. Whitaker, I.S.; Rao, J.; Izadi, D.; Butler, P.E. Historical Article: Hirudo medicinalis: Ancient origins of, and trends in the use of medicinal leeches throughout history. Br. J. Oral Maxillofac. Surg. 2004, 42, 133–137. [Google Scholar] [CrossRef] [Green Version]
  2. Kvist, S.; Manzano-Marín, A.; de Carle, D.; Trontelj, P.; Siddall, M.E. Draft genome of the European medicinal leech Hirudo medicinalis (Annelida, Clitellata, Hirudiniformes) with emphasis on anticoagulants. Sci. Rep. 2020, 10, 9885. [Google Scholar] [CrossRef] [PubMed]
  3. Salzet, M. Anticoagulants and inhibitors of platelet aggregation derived from leeches. FEBS Lett. 2001, 492, 187–192. [Google Scholar] [CrossRef] [Green Version]
  4. Greinacher, A.; Warkentin, T.E. The direct thrombin inhibitor hirudin. Thromb. Haemost. 2008, 99, 819–829. [Google Scholar] [CrossRef]
  5. Markwardt, F. Hirudin as alternative anticoagulant—A historical review. Semin. Thromb. Hemost. 2002, 28, 405–414. [Google Scholar] [CrossRef]
  6. Callas, D.D.; Hoppensteadt, D.; Fareed, J. Comparative studies on the anticoagulant and protease generation inhibitory actions of newly developed site-directed thrombin inhibitory drugs. Efegatran, argatroban, hirulog, and hirudin. Semin. Thromb. Hemost. 1995, 21, 177–183. [Google Scholar] [CrossRef]
  7. Arsenault, K.A.; Hirsh, J.; Whitlock, R.P.; Eikelboom, J.W. Direct thrombin inhibitors in cardiovascular disease. Nat. Rev. Cardiol. 2012, 9, 402–414. [Google Scholar] [CrossRef]
  8. Lee, C.J.; Ansell, J.E. Direct thrombin inhibitors. Br. J. Clin. Pharmacol. 2011, 72, 581–592. [Google Scholar] [CrossRef]
  9. Griffin, J.H. Blood coagulation. The thrombin paradox. Nature 1995, 378, 337–338. [Google Scholar] [CrossRef]
  10. Shaw, C.M.; O’Hanlon, D.M.; McEntee, G.P. Venous thrombosis. Am. J. Surg. 2003, 186, 167–168. [Google Scholar] [CrossRef]
  11. Lane, D.A.; Philippou, H.; Huntington, J.A. Directing thrombin. Blood 2005, 106, 2605–2612. [Google Scholar] [CrossRef] [PubMed]
  12. Muszbek, L.; Bereczky, Z.; Bagoly, Z.; Komáromi, I.; Katona, É. Factor XIII: A coagulation factor with multiple plasmatic and cellular functions. Physiol. Rev. 2011, 91, 931–972. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Weisel, J.W.; Litvinov, R.I. Fibrin Formation, Structure and Properties. Subcell. Biochem. 2017, 82, 405–456. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gupta, S.; Biswas, A.; Akhter, M.S.; Krettler, C.; Reinhart, C.; Dodt, J.; Reuter, A.; Philippou, H.; Ivaskevicius, V.; Oldenburg, J. Revisiting the mechanism of coagulation factor XIII activation and regulation from a structure/functional perspective. Sci. Rep. 2016, 6, 30105. [Google Scholar] [CrossRef] [PubMed]
  15. Singh, S.; Nazabal, A.; Kaniyappan, S.; Pellequer, J.-L.; Wolberg, A.S.; Imhof, D.; Oldenburg, J.; Biswas, A. The Plasma Factor XIII Heterotetrameric Complex Structure: Unexpected Unequal Pairing within a Symmetric Complex. Biomolecules 2019, 9, 765. [Google Scholar] [CrossRef] [Green Version]
  16. Stieler, M.; Weber, J.; Hils, M.; Kolb, P.; Heine, A.; Büchold, C.; Pasternack, R.; Klebe, G. Structure of active coagulation factor XIII triggered by calcium binding: Basis for the design of next-generation anticoagulants. Angew. Chem. Int. Ed. Engl. 2013, 52, 11930–11934. [Google Scholar] [CrossRef]
  17. Aleman, M.M.; Byrnes, J.R.; Wang, J.-G.; Tran, R.; Lam, W.A.; Di Paola, J.; Mackman, N.; Degen, J.L.; Flick, M.J.; Wolberg, A.S. Factor XIII activity mediates red blood cell retention in venous thrombi. J. Clin. Investig. 2014, 124, 3590–3600. [Google Scholar] [CrossRef] [Green Version]
  18. Byrnes, J.R.; Wolberg, A.S. Newly-Recognized Roles of Factor XIII in Thrombosis. Semin. Thromb. Hemost. 2016, 42, 445–454. [Google Scholar] [CrossRef] [Green Version]
  19. Duval, C.; Allan, P.; Connell, S.D.A.; Ridger, V.C.; Philippou, H.; Ariëns, R.A.S. Roles of fibrin α- and γ-chain specific cross-linking by FXIIIa in fibrin structure and function. Thromb. Haemost. 2014, 111, 842–850. [Google Scholar] [CrossRef]
  20. Wolberg, A.S. Fibrinogen and factor XIII: Newly recognized roles in venous thrombus formation and composition. Curr. Opin. Hematol. 2018, 25, 358–364. [Google Scholar] [CrossRef]
  21. Avery, C.A.; Pease, R.J.; Smith, K.; Boothby, M.; Buckley, H.M.; Grant, P.J.; Fishwick, C.W.G. (±) cis-Bisamido epoxides: A novel series of potent FXIII-A inhibitors. Eur. J. Med. Chem. 2015, 98, 49–53. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Iwata, Y.; Tago, K.; Kiho, T.; Kogen, H.; Fujioka, T.; Otsuka, N.; Suzuki-Konagai, K.; Ogita, T.; Miyamoto, S. Conformational analysis and docking study of potent factor XIIIa inhibitors having a cyclopropenone ring. J. Mol. Graph. Model. 2000, 18, 591–599. [Google Scholar] [CrossRef]
  23. Schmitz, T.; Bäuml, C.A.; Imhof, D. Inhibitors of blood coagulation factor XIII. Anal. Biochem. 2020, 605, 113708. [Google Scholar] [CrossRef] [PubMed]
  24. Al-Horani, R.A.; Kar, S. Factor XIIIa inhibitors as potential novel drugs for venous thromboembolism. Eur. J. Med. Chem. 2020, 200, 112442. [Google Scholar] [CrossRef] [PubMed]
  25. Pasternack, R.; Büchold, C.; Jähnig, R.; Pelzer, C.; Sommer, M.; Heil, A.; Florian, P.; Nowak, G.; Gerlach, U.; Hils, M. Novel inhibitor ZED3197 as potential drug candidate in anticoagulation targeting coagulation FXIIIa (F13a). J. Thromb. Haemost. 2020, 18, 191–200. [Google Scholar] [CrossRef] [Green Version]
  26. Al-Horani, R.A.; Karuturi, R.; Lee, M.; Afosah, D.K.; Desai, U.R. Allosteric Inhibition of Factor XIIIa. Non-Saccharide Glycosaminoglycan Mimetics, but Not Glycosaminoglycans, Exhibit Promising Inhibition Profile. PLoS ONE 2016, 11, e0160189. [Google Scholar] [CrossRef] [Green Version]
  27. Finney, S.; Seale, L.; Sawyer, R.T.; Wallis, R.B. Tridegin, a new peptidic inhibitor of factor XIIIa, from the blood-sucking leech Haementeria ghilianii. Biochem. J. 1997, 324 Pt 3, 797–805. [Google Scholar] [CrossRef] [Green Version]
  28. Bäuml, C.A.; Paul George, A.A.; Schmitz, T.; Sommerfeld, P.; Pietsch, M.; Podsiadlowski, L.; Steinmetzer, T.; Biswas, A.; Imhof, D. Distinct 3-disulfide-bonded isomers of tridegin differentially inhibit coagulation factor XIIIa: The influence of structural stability on bioactivity. Eur. J. Med. Chem. 2020, 201, 112474. [Google Scholar] [CrossRef]
  29. Bäuml, C.A.; Schmitz, T.; Paul George, A.A.; Sudarsanam, M.; Hardes, K.; Steinmetzer, T.; Holle, L.A.; Wolberg, A.S.; Pötzsch, B.; Oldenburg, J.; et al. Coagulation Factor XIIIa Inhibitor Tridegin: On the Role of Disulfide Bonds for Folding, Stability, and Function. J. Med. Chem. 2019, 62, 3513–3523. [Google Scholar] [CrossRef]
  30. Böhm, M.; Bäuml, C.A.; Hardes, K.; Steinmetzer, T.; Roeser, D.; Schaub, Y.; Than, M.E.; Biswas, A.; Imhof, D. Novel insights into structure and function of factor XIIIa-inhibitor tridegin. J. Med. Chem. 2014, 57, 10355–10365. [Google Scholar] [CrossRef]
  31. Böhm, M.; Kühl, T.; Hardes, K.; Coch, R.; Arkona, C.; Schlott, B.; Steinmetzer, T.; Imhof, D. Synthesis and functional characterization of tridegin and its analogues: Inhibitors and substrates of factor XIIIa. ChemMedChem 2012, 7, 326–333. [Google Scholar] [CrossRef] [PubMed]
  32. Chang, J.-Y. Diverse pathways of oxidative folding of disulfide proteins: Underlying causes and folding models. Biochemistry 2011, 50, 3414–3431. [Google Scholar] [CrossRef] [PubMed]
  33. Paul George, A.A.; Heimer, P.; Maaß, A.; Hamaekers, J.; Hofmann-Apitius, M.; Biswas, A.; Imhof, D. Insights into the Folding of Disulfide-Rich μ-Conotoxins. ACS Omega 2018, 3, 12330–12340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Liu, H.; Boudreau, M.A.; Zheng, J.; Whittal, R.M.; Austin, P.; Roskelley, C.D.; Roberge, M.; Andersen, R.J.; Vederas, J.C. Chemical synthesis and biological activity of the neopetrosiamides and their analogues: Revision of disulfide bond connectivity. J. Am. Chem. Soc. 2010, 132, 1486–1487. [Google Scholar] [CrossRef]
  35. Günther, P.; Buchner, L. Combined automated NOE assignment and structure calculation with CYANA. J. Biomol. NMR 2015, 62, 453–471. [Google Scholar] [CrossRef]
  36. Konagurthu, A.S.; Whisstock, J.C.; Stuckey, P.J.; Lesk, A.M. MUSTANG: A multiple structural alignment algorithm. Proteins 2006, 64, 559–574. [Google Scholar] [CrossRef]
  37. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1, 19–25. [Google Scholar] [CrossRef] [Green Version]
  38. Van der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J.C. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef]
  39. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  40. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [Green Version]
  41. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 14101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  43. Nosé, S.; Klein, M.L. Constant pressure molecular dynamics for molecular systems. Mol. Phys. 1983, 50, 1055–1076. [Google Scholar] [CrossRef]
  44. Hess, B.; Bekker, H.; Berendsen, H.J.C.; Fraaije, J.G.E.M. LINCS: A linear constraint solver for molecular simulations. J. Comput. Chem. 1997, 18, 1463–1472. [Google Scholar] [CrossRef]
  45. Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089–10092. [Google Scholar] [CrossRef] [Green Version]
  46. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577–8593. [Google Scholar] [CrossRef] [Green Version]
  47. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  48. Yan, Y.; Tao, H.; He, J.; Huang, S.-Y. The HDOCK server for integrated protein-protein docking. Nat. Protoc. 2020, 15, 1829–1852. [Google Scholar] [CrossRef]
  49. Moult, J.; Fidelis, K.; Kryshtafovych, A.; Schwede, T.; Tramontano, A. Critical assessment of methods of protein structure prediction (CASP)--round x. Proteins 2014, 82 (Suppl. 2), 1–6. [Google Scholar] [CrossRef] [Green Version]
  50. Van Zundert, G.C.P.; Rodrigues, J.P.G.L.M.; Trellet, M.; Schmitz, C.; Kastritis, P.L.; Karaca, E.; Melquiond, A.S.J.; van Dijk, M.; de Vries, S.J.; Bonvin, A.M.J.J. The HADDOCK2.2 Web Server: User-Friendly Integrative Modeling of Biomolecular Complexes. J. Mol. Biol. 2016, 428, 720–725. [Google Scholar] [CrossRef] [Green Version]
  51. Tina, K.G.; Bhadra, R.; Srinivasan, N. PIC: Protein Interactions Calculator. Nucleic Acids Res. 2007, 35, W473–W476. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. (a) Derivation of the two fragments from the two-disulfide-bonded tridegin isomer B[C19S,C25S] [29]. The N-terminal fragment consists of amino acid 1–37 and contains the disulfide bonds Cys5–Cys37 and Cys17–Cys31. The C-terminal fragment is composed of amino acids 38–66 without a cysteine residue; (b) HPLC elution profiles of the stepwise oxidation strategy applied for the synthesis of the N-terminal fragment. During the first oxidation (1.1 eq. iodine, 100% AcOH) the first disulfide bond between Cys5 and Cys37 was formed (1). Subsequently the second disulfide bond between Cys17 and Cys31 was built (2) by increasing the amount of iodine (15 eq.) and adjusting the concentration of AcOH (70%); (c) MS/MS analysis of the N-terminal fragment after digestion with chymotrypsin. The disulfide bond connectivity was confirmed by the characteristic MS/MS fragments C(CW) AF and KLLP C(IPQC) KEW.
Figure 1. (a) Derivation of the two fragments from the two-disulfide-bonded tridegin isomer B[C19S,C25S] [29]. The N-terminal fragment consists of amino acid 1–37 and contains the disulfide bonds Cys5–Cys37 and Cys17–Cys31. The C-terminal fragment is composed of amino acids 38–66 without a cysteine residue; (b) HPLC elution profiles of the stepwise oxidation strategy applied for the synthesis of the N-terminal fragment. During the first oxidation (1.1 eq. iodine, 100% AcOH) the first disulfide bond between Cys5 and Cys37 was formed (1). Subsequently the second disulfide bond between Cys17 and Cys31 was built (2) by increasing the amount of iodine (15 eq.) and adjusting the concentration of AcOH (70%); (c) MS/MS analysis of the N-terminal fragment after digestion with chymotrypsin. The disulfide bond connectivity was confirmed by the characteristic MS/MS fragments C(CW) AF and KLLP C(IPQC) KEW.
Ijms 22 00880 g001
Figure 2. Backbone HN–Hα region of the 150 ms NOESY (nuclear Overhauser enhancement spectroscopy) spectrum of the C-terminal fragment (a) and the N-terminal fragment (b) in 95% H2O/5% D2O with backbone connectivities labeled.
Figure 2. Backbone HN–Hα region of the 150 ms NOESY (nuclear Overhauser enhancement spectroscopy) spectrum of the C-terminal fragment (a) and the N-terminal fragment (b) in 95% H2O/5% D2O with backbone connectivities labeled.
Ijms 22 00880 g002
Figure 3. NMR structure of the rigid N-terminal fragment (left) and the flexible C-terminal fragment (right) presented as ensembles of the 10 structures with the lowest energies. grey: flexible C-terminal fragment; yellow: disulfide bonds Cys5–Cys37 and Cys17–Cys31; blue: helix; green: β-turns.
Figure 3. NMR structure of the rigid N-terminal fragment (left) and the flexible C-terminal fragment (right) presented as ensembles of the 10 structures with the lowest energies. grey: flexible C-terminal fragment; yellow: disulfide bonds Cys5–Cys37 and Cys17–Cys31; blue: helix; green: β-turns.
Ijms 22 00880 g003
Figure 4. (a) Figure depicts 20 structures (green cartoons) from the 100-member NMR ensemble of full-length isomer B[C19S,C25S] [29] aligned on the N-terminal residues (1–37). The lowest energy model 1 is shown as red cartoons. (b) The lowest energy model from the NMR ensemble of isomer B[C19S,C25S] [29] used as input in the 300 ns molecular dynamics (MD) simulation. (c) Backbone root mean square deviation (RMSD) trace from the 300 ns MD simulation of isomer B[C19S,C25S] [29] with respect to its starting structure. The initial 15 ns of this simulation is considered as the equilibration phase (shaded red) and the remainder of the simulation (shaded green) is considered the production phase of the simulation. (d) Figure presents 20 equidistant structures from between 150 and 300 ns of the simulation as an ensemble (white to blue gradient cartoons) with the starting structure (red cartoon) shown as an indicator of the extent of structural change especially at the C-terminal (residues 38–66). (e) The final snapshot of the 300 ns MD simulation (blue cartoon) shown along with the starting structure (red transparent cartoon).
Figure 4. (a) Figure depicts 20 structures (green cartoons) from the 100-member NMR ensemble of full-length isomer B[C19S,C25S] [29] aligned on the N-terminal residues (1–37). The lowest energy model 1 is shown as red cartoons. (b) The lowest energy model from the NMR ensemble of isomer B[C19S,C25S] [29] used as input in the 300 ns molecular dynamics (MD) simulation. (c) Backbone root mean square deviation (RMSD) trace from the 300 ns MD simulation of isomer B[C19S,C25S] [29] with respect to its starting structure. The initial 15 ns of this simulation is considered as the equilibration phase (shaded red) and the remainder of the simulation (shaded green) is considered the production phase of the simulation. (d) Figure presents 20 equidistant structures from between 150 and 300 ns of the simulation as an ensemble (white to blue gradient cartoons) with the starting structure (red cartoon) shown as an indicator of the extent of structural change especially at the C-terminal (residues 38–66). (e) The final snapshot of the 300 ns MD simulation (blue cartoon) shown along with the starting structure (red transparent cartoon).
Ijms 22 00880 g004
Figure 5. (a) The final 20 snapshots of the 300 ns MD simulation of the full-length model of isomer B[C19S,C25S] [29]. (b) Computationally derived structure of the same peptide analogue from an earlier study [29]. (c) Final snapshot from the 300 ns simulation from the current study. (d) Structural alignment of the final structure from the 300 ns MD simulation from the current study (blue) with the computationally derived structure of the same peptide (orange) from an earlier study [29]. The alignment indicates that the two models essentially have varying structures while the underlying fold is preserved. The computationally derived structure from earlier work has a more compact fold owing to the fact that it was created based on its three-disulfide-bonded parent, where residues Cys19 and Cys25 were linked by a disulfide bond. A fair alignment between the trailing end of the C-terminal segment is observed. The alignment was carried out using the MUSTANG algorithm [36].
Figure 5. (a) The final 20 snapshots of the 300 ns MD simulation of the full-length model of isomer B[C19S,C25S] [29]. (b) Computationally derived structure of the same peptide analogue from an earlier study [29]. (c) Final snapshot from the 300 ns simulation from the current study. (d) Structural alignment of the final structure from the 300 ns MD simulation from the current study (blue) with the computationally derived structure of the same peptide (orange) from an earlier study [29]. The alignment indicates that the two models essentially have varying structures while the underlying fold is preserved. The computationally derived structure from earlier work has a more compact fold owing to the fact that it was created based on its three-disulfide-bonded parent, where residues Cys19 and Cys25 were linked by a disulfide bond. A fair alignment between the trailing end of the C-terminal segment is observed. The alignment was carried out using the MUSTANG algorithm [36].
Ijms 22 00880 g005
Figure 6. Molecular docking analysis of the interaction of isomer B[C19S,C25S] [29] with the catalytic site of FXIIIa. Left: Crystal structure of the calcium-activated FXIIIa-monomer (PDB ID: 4KTY [16], grey surface) docked with one of the top-ten-ranked conformations of isomer B[C19S,C25S] [29] (cyan ribbon), which interacts with the catalytic site of FXIIIa. Right: Zoom-in from different directions (lysine access site (top, K-site) and glutamine access site (bottom, Q-site)) of the FXIIIa active site interacting with isomer B[C19S,C25S] [29]. The active site residues of FXIIIa (Cys314, His373, Asp396, Trp279, and Trp370) are shown as blue sticks enclosed by their molecular surfaces (violet translucent surface).
Figure 6. Molecular docking analysis of the interaction of isomer B[C19S,C25S] [29] with the catalytic site of FXIIIa. Left: Crystal structure of the calcium-activated FXIIIa-monomer (PDB ID: 4KTY [16], grey surface) docked with one of the top-ten-ranked conformations of isomer B[C19S,C25S] [29] (cyan ribbon), which interacts with the catalytic site of FXIIIa. Right: Zoom-in from different directions (lysine access site (top, K-site) and glutamine access site (bottom, Q-site)) of the FXIIIa active site interacting with isomer B[C19S,C25S] [29]. The active site residues of FXIIIa (Cys314, His373, Asp396, Trp279, and Trp370) are shown as blue sticks enclosed by their molecular surfaces (violet translucent surface).
Ijms 22 00880 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Schmitz, T.; Paul George, A.A.; Nubbemeyer, B.; Bäuml, C.A.; Steinmetzer, T.; Ohlenschläger, O.; Biswas, A.; Imhof, D. NMR-Based Structural Characterization of a Two-Disulfide-Bonded Analogue of the FXIIIa Inhibitor Tridegin: New Insights into Structure–Activity Relationships. Int. J. Mol. Sci. 2021, 22, 880. https://doi.org/10.3390/ijms22020880

AMA Style

Schmitz T, Paul George AA, Nubbemeyer B, Bäuml CA, Steinmetzer T, Ohlenschläger O, Biswas A, Imhof D. NMR-Based Structural Characterization of a Two-Disulfide-Bonded Analogue of the FXIIIa Inhibitor Tridegin: New Insights into Structure–Activity Relationships. International Journal of Molecular Sciences. 2021; 22(2):880. https://doi.org/10.3390/ijms22020880

Chicago/Turabian Style

Schmitz, Thomas, Ajay Abisheck Paul George, Britta Nubbemeyer, Charlotte A. Bäuml, Torsten Steinmetzer, Oliver Ohlenschläger, Arijit Biswas, and Diana Imhof. 2021. "NMR-Based Structural Characterization of a Two-Disulfide-Bonded Analogue of the FXIIIa Inhibitor Tridegin: New Insights into Structure–Activity Relationships" International Journal of Molecular Sciences 22, no. 2: 880. https://doi.org/10.3390/ijms22020880

APA Style

Schmitz, T., Paul George, A. A., Nubbemeyer, B., Bäuml, C. A., Steinmetzer, T., Ohlenschläger, O., Biswas, A., & Imhof, D. (2021). NMR-Based Structural Characterization of a Two-Disulfide-Bonded Analogue of the FXIIIa Inhibitor Tridegin: New Insights into Structure–Activity Relationships. International Journal of Molecular Sciences, 22(2), 880. https://doi.org/10.3390/ijms22020880

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