Next Article in Journal
Application of Alginate Hydrogels for Next-Generation Articular Cartilage Regeneration
Previous Article in Journal
In Vivo Rapid Investigation of CRISPR-Based Base Editing Components in Escherichia coli (IRI-CCE): A Platform for Evaluating Base Editing Tools and Their Components
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

NMR and MD Analysis of the Bonding Interaction of Vancomycin with Muramyl Pentapeptide

by
Rafał Ślusarz
*,
Barbara Dmochowska
,
Justyna Samaszko-Fiertek
,
Krzysztof Brzozowski
and
Janusz Madaj
Faculty of Chemistry, University of Gdańsk, Wita Stwosza 63, 80-308 Gdańsk, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(3), 1146; https://doi.org/10.3390/ijms23031146
Submission received: 14 December 2021 / Revised: 17 January 2022 / Accepted: 18 January 2022 / Published: 20 January 2022
(This article belongs to the Section Macromolecules)

Abstract

:
The article describes an NMR spectroscopy study of interactions between vancomycin and a muramyl pentapeptide in two complexes: vancomycin and a native muramyl pentapeptide ended with D-alanine (MPP-D-Ala), and vancomycin and a modified muramyl pentapeptide ended with D-serine (MPP-D-Ser). The measurements were made in a 9:1 mixture of H2O and D2O. The obtained results confirmed the presence of hydrogen bonds previously described in the literature. At the same time, thanks to the pentapeptide model used, we were able to prove the presence of two more hydrogen bonds formed by the side chain amino group of L-lysine and oxygen atoms from the vancomycin carboxyl and amide groups. This type of interaction has not been described before. The existence of these hydrogen bonds was confirmed by the 1H NMR and molecular modeling. The formation of these bonds incurs additional through-space interactions, visible in the NOESY spectrum, between the protons of the L-lysine amino group and a vancomycin-facing hydrogen atom in the benzylic position. The presence of such interactions was also confirmed by molecular dynamics trajectory analysis.

1. Introduction

Vancomycin is one of the most important glycoprotein antibiotics and, thanks to its effectiveness, it is often called a last resort antibiotic. For over 60 years it has been used to treat infections with bacteria such as: Staphylococci, Streptococci, Enterococci, and Pneumococci as well as Corynebacterium, Listeria, Bacillus spp., Clostridia. It was first isolated by Dr. E.C. Kornfield, an organic chemist at Eli Lilly, from a fungus called Streptomyces orientalis. However, since then, many strains of bacteria resistant to vancomycin treatment have been isolated and described in the literature. In the mid-1980s, vancomycin-resistant S. aureus (VRSA) [1,2,3] and vancomycin-resistant S. epidermidis (VRSE) [4,5] strains appeared. The emergence of these vancomycin-resistant strains necessitates the search for new antibiotics that would be more effective in their treatment. The discovery of these vancomycin-resistant strains necessitates the search for new antibiotics that would be more effective in treating them. In the search for new, more effective antibiotics, the knowledge of their mechanism of action is extremely helpful. In the case of vancomycin, the most likely mechanism is the formation of a complex of vancomycin with the last two D-alanine residues in muramyl pentapeptide sequence [6,7]. This blocks the possibility of peptidoglycan cross-linking and, consequently, the expansion of the cell wall of gram-positive bacteria. The vancomycin-muramyl complex is stabilized by the hydrophobic van der Waals interactions and five hydrogen bonds formed between the L-lysine and subsequent two D-alanine residues, and the vancomycin binding pocket formed by the cyclic heptapeptide (Figure 1). To counteract this, bacteria have developed a defense mechanism involving the exchange of a terminal amino acid from D-alanine to D-serine or lactate. Especially in the latter case, this substitution causes significant decrease in effectiveness, which is explained by a weakening of the network of previously formed hydrogen bonds [8,9,10].
On the other hand, studies of the activity of vancomycin analogues modified in the sugar part show increased effectiveness, as pointed out by Min Ge and co-authors [11]. As they suggest, the mechanism of action of such antibiotics could be different from the interaction with the terminal peptidoglycan fragment. This possibility is also confirmed by our research on the activity of aglycone against some reference strains, which, despite the presence of full structure of the cyclic heptapeptide, shows a significantly reduced effectiveness [12]. An additional mechanism of vancomycin’s action has been described in the literature that involves the inhibition of ribonucleic acid synthesis by changing the permeability of the bacterial cell membrane [13]. This all shows that more research is still needed to further elucidate the effects of vancomycin as an antibiotic.

2. Results and Discussion

Research aimed at elucidating the mechanism of the interaction of vancomycin with peptidoglycan forming the gram-positive bacteria cell wall has been conducted for many years. Initially, attempts were made to determine these interactions based on X-ray analysis [14]. Soon after the structure of the vancomycin molecule was fully elucidated, the groundbreaking work of Williams’ team appeared, in which they used NMR spectroscopy to study the interaction of D-Ala-D-Ala dipeptide with vancomycin [15]. Further work that tried to explain this mechanism used increasingly complex research models with an increasingly longer peptidoglycan chain, such as the Whitesides study [16], in which they used a tripeptide including lysine. In our study on this mechanism, we decided to use the muramyl pentapeptide derivative, the synthesis, and characteristics of which we described in our earlier work [17].

2.1. NMR Study

We used an NMR spectroscopy to analyze the interactions. As an environment, we decided to use H2O and D2O mixture in a volume ratio of 9 to 1. Although most of the vancomycin NMR analyzes, including the first full assignment of signals [18,19,20] and probably the newest [21], were performed with DMSO, this system differs from what we observe in nature. Although the spectrum of vancomycin in water does not have such well separated signals, the use of high-resolution 700 MHz 1H NMR spectra, in combination with the solvent mixture used, turned out to have great advantages. The group of H-N amide signals and peptidoglycan Cα-H proton signals were completely separated, which made it possible to accurately track changes in their chemical shift on the spectrum and changes in the structure of the signals. We recorded the NMR spectra of two muramyl pentapeptides. The first is native pentapeptide with C-terminal D-alanine (MPP-D-Ala) and the second with D-serine (MPP-D-Ser) suggested in the literature [22]. For both pentapeptides, the spectra were also recorded after the addition of vancomycin in a molar ratio of 1:1. Moreover, for the native pentapeptide MPP-D-Ala, the spectrum was recorded after the addition of vancomycin in a molar ratio of 1:0.5. It allowed us to observe interesting effects confirming the influence of vancomycin concentration on the change of the structure of peptidoglycan signals.
At the highest values of chemical shifts, signals from the 3 OH groups located on the aromatic ring of vancomycin are observed. Although these signals are broad, we can observe clear changes in their position as well. This is clearly seen in Figure 2, showing this fragment of the spectrum.
Much more characteristic changes can be observed in the region where the muramyl pentapeptide chain amide protons signals occur, which is shown in Figure 3.
Analysis of these spectra shows that both the L-Ala amide proton 37 (for atom numbering see Figure 1) and the muramic acid NHAc groups proton 39, undergo relatively slight changes in their positions. Considering the suggested mechanism of the interaction of vancomycin and muramyl pentapeptide, in which two D-alanine residues and oxygen atom from the L-lysine carbonyl group are involved, it is not surprising.
The notable change in the value of the chemical shift was observed for the amide proton 15 signal of the D-alanine amino acid residue penultimate in the pentapeptide chain. This is also seen in the MPP-D-Ser with vancomycin spectra. In the pure MPP-D-Ala spectrum, the signal of the amide proton of this D-alanine appears as the clearly separated two doublets (Δδ 0.065 ppm). When vancomycin was added, the two signals overlapped. It is even more spectacularly visible on the example of MMP-D-Ser. The presence of a hydroxyl group in D-serine causes the amide proton 15 signal to appear in the spectrum at higher δ values. When vancomycin is added, it clearly moves up the field. This behavior can be explained by Figure 4.
If the oxygen atom 17 in the amide group begins to form a hydrogen bond, the electron density on it decreases. As a result, the order of the C-N partial double bond decreases, and the amide proton moves up the field. Such a situation is observed in the case of the carboxylic oxygen atom 17 of L-lysine and the amide hydrogen atom 15 of penultimate D-alanine, which leads to a shift of the signal of this proton observed in the spectrum. The analysis of the shift in the HSQC spectrum of the proton alpha 19 signal (Figure 5) confirms the participation of this oxygen atom in the formation of hydrogen bond. The acidity of the proton in the alpha position to the carbonyl carbon is well known in organic chemistry. The change in the electron density on this oxygen atom because of hydrogen bond formation reduces the acidity of this proton, which in the NMR spectrum results in its shift towards lower δ values. At the same time there is a shift of the carbon atom 18 towards higher δ values. We can observe it both in the spectra of the mixture of vancomycin with MPP-D-Ala and MPP-D-Ser, as shown in Figure 5.
Analogously to oxygen atom 17 of amide group of L-lysine, terminal D-alanine, or D-serine carboxyl oxygen atoms 1 and 2 are involved in the hydrogen bond formation with vancomycin. As in the case of L-lysine, the spectra show analogous shifts of alpha proton 5. The presence of a hydroxyl group in D-serine makes this shift more spectacular in this case. Fragments of 1H NMR spectra illustrating this phenomenon are shown in Figure 6.
In addition to the two carboxyl oxygen atoms of the muramyl pentapeptide, the amide hydrogen 8 of the terminal D-alanine or D-serine is also involved in binding to vancomycin. According to the literature [23] for these atoms, this should result in a signal shift down the field. The analysis of the spectra in Figure 3 confirms the expected changes, but they have a smaller range and, similarly to the previous one, due to the specifics of the spectra, they are better visible for D-serine.
All these observations presented above confirm the contribution of carboxy oxygen atoms and amide hydrogen in the formation of a complex with vancomycin presented in the literature. However, further analysis of the spectra indicates the formation of one more hydrogen bond, previously not described in the literature. Observing the position of the proton signal of the L-lysine amino group, it can be clearly seen that after adding vancomycin it moves down the field (Figure 7). The analysis of spectra with different vancomycin concentration (b and c) turns out to be useful in this case. It is clearly seen that the signal resulting from the interaction increases with increasing vancomycin concentration.
An analogous shift of the signal of protons of this group in the 1H NMR spectrum is observed in the case of MPP-D-Ser (Figure 7f).
The question of which vancomycin group is hydrogen bonded by the L-lysine side chain amino group could be answered by analysis of the NOESY spectrum (Figure 8). It shows a cross-peak indicating the through-space interaction of the protons of the amino group of L-lysine and the benzyl proton 42 of vancomycin (in Figure 9 it is marked in green) resulting from the hydrogen bond formed between the amino group of L-lysine and the oxygen atoms 40 and 41 of the carboxyl and amide groups of vancomycin, as shown in Figure 9.
While reviewing the literature on the interaction of vancomycin with muramyl pentapeptide, we did not find information about this type of interaction. This may result from the previously used research models, in which this group was ignored or even protected with an acetyl group [24].

2.2. MD Study

The possibility of this binding did not come as a surprise to us. During the work on the material that we had previously published [25], we observed the formation of such a complex in molecular dynamics (MD) simulations. In MD, after initial geometry optimization, the distance between L-lysine side chain amine group and vancomycin C-terminal part was favorable for hydrogen bond formation.
During the MD simulations we monitored existence of the expected hydrogen bonds between vancomycin and muramyl pentapeptide MPP-D-Ala (marked in Figure 1 and Figure 9) and, additionally, the distances between L-lysine side chain amine group and any of the potential proton acceptors: carboxyl oxygen atom from vancomycin C-terminus or carboxyl oxygen atom from its neighboring amide group (atoms 40 and 41, respectively; for atom numbering see Figure 9). From the MD run, a multiplots of those distances were prepared (Figure 10).
Lengths of newly described interactions were measured as a virtual distance between center of mass of L-lysine side chain amine group (involving four atoms) and center of mass of two oxygen atoms from the vancomycin C-terminal carboxyl group (distance 1) or previously described amine group and the oxygen atom 41 included in the amide group next to the vancomycin C-terminus (distance 2).
Initial 100 ps of the MD run (TDMD; see the Section 3) in Figure 10 is ended with red vertical, dashed line and is not considered as a part of actual MD therefore the distances between interacting parties in this region of the plots are not discussed here.
From the beginning of free, unrestrained MD, both distances between newly discovered centers of interactions are kept just below 4 Å, which is suitable for hydrogen bond formation. After 0.6 ns of MD these identified hydrogen bonds are disturbed and broken—the values of monitored distances are growing above 5 Å. Moreover, from the plots in Figure 10 one can see, that the distance 1 is even more often preserved on constant level than distance 2 which confirms higher adaptative conformational flexibility of vancomycin C-terminal carboxyl group then amide group involved in vancomycin macrocyclic peptide backbone.
The favorable distances between interacting parties are not preserved during whole MD which suggests that the interacting parties are observed only in some local energetic minimum induced by the initial TDMD. One can observe the presence of the hydrogen bonds discussed but the MD run in our investigation does not confirm their presence during the whole MD. It is caused by the full conformational freedom of all interacting parties and is expected. The confirmation of this can be found in Figure 11.
In Figure 11 one can see control distances measured between known and expected points of interactions between vancomycin and pentapeptide fragment complexed during vancomycin active bonding to the peptidoglycan cell wall [14,15,26] resulting in formation of five hydrogen bonds. If the line representing any distance is being plot below the contractual 5 Å limit, we assume (and can visually confirm) that the hydrogen bond is present in this moment of our MD simulation. Presence of these hydrogen bonds during MD simulation confirms proper model building and stable vancomycin-muramyl pentapeptide interaction.
Before 0.6 ns of MD simulation all monitored distances are kept at the expected level—below 5 Å—excluding distance 3. The distance 3 values are erratic and only for two thirds of the 0.6 ns of MD discussed is kept below 5 Å. This can be explained by the fact that this hydrogen bond is formed between hydrogen atom from the vancomycin C-terminal neighboring amide group and L-lysine amide group oxygen atom. The hydrogen atom discussed belongs to the amide group involved in hydrogen bond formation described as distance 2 in Figure 10. Wherever muramyl pentapeptide L-lysine residue is being pulled apart from vancomycin (from its C-terminus) both distances are increased—the hydrogen bonds represented by distances 2 and 3 are disrupted and possibly, temporarily terminated.
Where the control distances are increased above our contractual 5 Å threshold the hydrogen bonds are being broken. In the simulated complex this is caused by freeing the muramyl pentapeptide from the vancomycin binding cavity partially (between 0.6 and 0.8 ns of MD) or as a whole (after 0.8 ns of MD). We did not continue the MD after this point and did not impose any restraints to keep the ligand in the position meeting crystallographic criteria in active vancomycin-ligand complex.
From the analysis of the trajectory of complex simulated one more critical distance got our attention. This is the distance between L-lysine side chain amine group protons and aliphatic hydrogen atom 42 (HB). This specific hydrogen atom appeared to induce proton-proton interaction observed in NOESY spectrum (Figure 8).
We prepared a plot illustrating how close those protons were during whole MD and presented it in Figure 12.
The distance analysis shows that at almost any point in the MD run at least one, most of the time two hydrogen atoms forming L-lysine amine group, is in the near vicinity to the aliphatic hydrogen atom 42 discussed. This makes it possible to influence the proton chemical shifts as observed in the NOESY spectra marked in green in Figure 8.
The interactions are possible, and we confirmed their existence in manual inspection of the MD trajectory snapshots excluding 190–230, 500–550 and above 690 ps MD regions. In those specific intervals the distances were greater, and the possible proton–proton interactions could be weaker or impossible due to too great distance. For most of the MD time discussed in this work (in the 0.1–0.6 ps interval of the MD it is almost 72% of the MD time) and even after the time when we no longer observe new, here described hydrogen bonds (after 0.6 ns of the MD time), the possible proton–proton interactions were present.
From the MD snapshots in which we could confirm presence of all discussed hydrogen bonds formed between muramyl pentapeptide and vancomycin we selected the complex with lowest total energy and presented it in Figure 13.

3. Materials and Methods

3.1. NMR Measurements

Amounts of vancomycin (12.1 mg, 8.1 μmol); vancomycin (6.1 mg, 4.1 μmol) and MPP-D-Ala (6.4 mg, 8.1 μmol); vancomycin (12.1 mg, 8.1 μmol) and MPP-D-Ala (6.4 mg, 8.1 μmol); vancomycin (13.1 mg, 8.9 μmol) and MPP-D-Ser (7.2 mg, 8.9 μmol) were dissolved in 0.5 mL of a mixture of H2O and D2O (vol. 9:1).
For the samples studied the following spectra were recorded: 1H and 2D TOCSY [27], NOESY [28,29] and ROESY [30] with over 4000 data points in f2 at mixing times of 80, 200 and 200 ms, respectively. Additionally for shift confirmation HSQC correlation spectrum was performed using standard sequence [31,32,33,34,35,36]. All measurements were carried out at Bruker 700 MHz spectrometer. All spectra were recorded at controlled temperature of 298 K using TXI inverse probe. To reflect natural conditions mixture of H2O and D2O at the ratio of 9:1 was used. Water signal suppression was achieved by the excitation sculpting protocol [37]. Obtained spectra were processed and analyzed with the usage of TopSpin 3.2 (Bruker BioSpin GmbH, Germany) software.

3.2. MD Simulations

For the MD simulations we have parameterized all non-standard groups and fragments for the AMBER [38] ff03.r1 [39,40] and glycam06g [41,42] force fields as described in detail in [25]. In the crystallographic complex of vancomycin with di-acetyl-Lys-D-Ala-D-Ala ligand (PDB ID 1FVM, chain O) the ligand structure has been replaced with muramyl pentapeptide derivative described in [17,25] of the sequence: Mur-Ala-D-iGlu-Lys-D-Ala-D-Ala. Vancomycin vancosamine and the Lys side chain amine groups were protonated.
The pre-equilibrated water boxes, TIP3P [43] model, were added to the computational system. The rectangular water box containing interacting complex and water molecules size was 10 Å in each direction, which resulted the computational system containing about 1800 water molecules.
Geometry of the complex was optimized, and energy minimized. After initial geometry optimization, the hydrogen bond restraints defined to meet the crystallographic criteria were imposed to provide mutual adjustment of vancomycin and the muramyl pentapeptide in water solution.
Initial equilibration with adjustment was driven in the 300 K with the force constant 1 kcal/mol to achieve the crystallographic criteria in this Target Driven Molecular Dynamics (TDMD). In TDMD the interacting parties are guided towards a final (target) structure by means of steering forces in selected interacting points, which in this case were defined as expected hydrogen bonds. The initial TDMD was carried out for 100 ps and the force constant was applied for first 50 ps and then gradually reduced to 0 kcal/mol in the last 50 ps of TDMD.
Subsequently the simulation complex was submitted to the isothermal-isobaric (300 K) molecular dynamics (MD) in the AMBER suite of programs. For the MD we used SHAKE algorithm [44] with default time step 0.001 ps, Berendsen [45] bath coupling (1 ps) and 9.0 Å cut-off for non-bonded interactions, with PME [46,47,48].
Snapshots from the trajectory were taken every 1000 steps (1 ps) and analyzed over whole 1 ns of the run in the CPPTRAJ [49] module.
Total vancomycin-muramyl pentapeptide complex energy (AMBER’s EAMBER energy) was calculated by SANDER module of AMBER suite of programs.

4. Conclusions

It can be stated that the use of high-resolution NMR allowed to obtain spectra of vancomycin complexes with MPP-D-Ala and MPP-D-Ser in the mixture of solvents: H2O and D2O. We were able to confirm the presence of hydrogen bonds suggested in the literature. The analysis of the chemical shifts of the alpha and amide protons of the muramyl pentapeptides used was extremely helpful. In addition, the model we used revealed two more, previously not described, hydrogen bonds formed with vancomycin by the amino group of L-lysine from muramyl pentapeptide and an additional proton–proton through-space interaction (NOESY spectrum) affecting vancomycin-muramyl pentapeptide complex formation. We confirmed this model in in silico modeling using molecular dynamics and trajectory analysis.

Author Contributions

Conceptualization, J.M.; Investigation, R.Ś. and K.B.; Data curation, R.Ś., B.D., J.S.-F. and K.B.; Writing—original draft preparation, R.Ś. and J.M.; Writing—review & editing, R.Ś. and J.M.; Supervision, J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Polish Ministry of Science and Higher Education DS 531-T100-D 501-21.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The computational time in the Academic Computer Center in Gdansk CI TASK, Poland is acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fridkin, S.K.; Hageman, J.; McDougal, L.K.; Mohammed, J.; Jarvis, W.R.; Perl, T.M.; Tenover, F.C.; Vancomycin-Intermediate Staphylococcus aureus Epidemiology Study Group. Epidemiological and Microbiological Characterization of Infections Caused byStaphylococcus aureuswith Reduced Susceptibility to Vancomycin, United States, 1997–2001. Clin. Infect. Dis. 2003, 36, 429–439. [Google Scholar] [CrossRef] [Green Version]
  2. Hiramatsu, K. Vancomycin-resistant Staphylococcus aureus: A new model of antibiotic resistance. Lancet Infect. Dis. 2001, 1, 147–155. [Google Scholar] [CrossRef]
  3. Howden, B.P.; Davies, J.K.; Johnson, P.D.R.; Stinear, T.P.; Grayson, M.L. Reduced Vancomycin Susceptibility in Staphylococcus aureus, Including Vancomycin-Intermediate and Heterogeneous Vancomycin-Intermediate Strains: Resistance Mechanisms, Laboratory Detection, and Clinical Implications. Clin. Microbiol. Rev. 2010, 23, 99–139. [Google Scholar] [CrossRef] [Green Version]
  4. Biavasco, F.; Vignaroli, C.; Varaldo, P.E. Glycopeptide Resistance in Coagulase-Negative Staphylococci. Eur. J. Clin. Microbiol. Infect. Dis. 2000, 19, 403–417. [Google Scholar] [CrossRef]
  5. Strahilevitz, J.; Rubinstein, E. Novel agents for resistant Gram-positive infections—a review. Int. J. Infect. Dis. 2002, 6, S38–S46. [Google Scholar] [CrossRef] [Green Version]
  6. Perkins, H.R. Specificity of combination between mucopeptide precursors and vancomycin or ristocetin. Biochem. J. 1969, 111, 195–205. [Google Scholar] [CrossRef] [Green Version]
  7. Barna, J.C.J.; Williams, D.H. The structure and mode of action of glycopeptide antibiotics of the vancomycin group. Annu. Rev. Microbiol. 1984, 38, 339–357. [Google Scholar] [CrossRef]
  8. Bugg, T.D.H.; Wright, G.; Dutka-Malen, S.; Arthur, M.; Courvalin, P.; Walsh, C.T. Molecular basis for vancomycin resistance in Enterococcus faecium BM4147: Biosynthesis of a depsipeptide peptidoglycan precursor by vancomycin resistance proteins VanH and VanA. Biochemistry 1991, 30, 10408–10415. [Google Scholar] [CrossRef]
  9. Arthur, M.; Courvalin, P. Genetics and mechanisms of glycopeptide resistance in enterococci. Antimicrob. Agents Chemother. 1993, 37, 1563–1571. [Google Scholar] [CrossRef] [Green Version]
  10. Walsh, C.T. Vancomycin Resistance: Decoding the Molecular Logic. Science 1993, 261, 308–309. [Google Scholar] [CrossRef]
  11. Ge, M.; Chen, Z.; Russell, H.; Onishi, H.R.; Kohler, J.; Silver, L.L.; Kerns, R.; Fukuzawa, S.; Thompson, C.; Kahne, D. Vancomycin Derivatives That Inhibit Peptidoglycan Biosynthesis Without Binding d -Ala- d -Ala. Science 1999, 284, 507–511. [Google Scholar] [CrossRef] [PubMed]
  12. Samaszko-Fiertek, J.; Szulc, M.; Dmochowska, B.; Jaśkiewicz, M.; Kamysz, W.; Ślusarz, R.; Madaj, J. Influence of Carbohydrate Residues on Antibacterial Activity of Vancomycin. Lett. Org. Chem. 2020, 17, 287–293. [Google Scholar] [CrossRef]
  13. Jordan, D.C.; Inniss, W.E. Selective Inhibition of Ribonucleic Acid Synthesis in Staphylococcus aureus by Vancomycin. Nature 1959, 184, 1894–1895. [Google Scholar] [CrossRef]
  14. Sheldrick, G.M.; Jones, P.G.; Kennard, O.; Williams, D.H.; Smith, G.A. Structure of vancomycin and its complex with acetyl-D-alanyl-D-alanine. Nature 1978, 271, 223–225. [Google Scholar] [CrossRef]
  15. Williams, D.H.; Williamson, M.P.; Butcher, D.W.; Hammond, S.J. Detailed binding sites of the antibiotics vancomycin and ristocetin A: Determination of intermolecular distances in antibiotic/substrate complexes by use of the time-dependent NOE. J. Am. Chem. Soc. 1983, 105, 1332–1339. [Google Scholar] [CrossRef]
  16. Rao, J.; Whitesides, G.M. Tight Binding of a Dimeric Derivative of Vancomycin with Dimeric l-Lys-d-Ala-d-Ala. J. Am. Chem. Soc. 1997, 119, 10286–10290. [Google Scholar] [CrossRef]
  17. Szulc, M.; Samaszko-Fiertek, J.; Slusarz, R.; Kowalska, K.; Sikorski, A.; Madaj, J. Solid phase synthesis of two muramyl pentapeptide derivatives. J. Carbohydr. Chem. 2016, 35, 161–171. [Google Scholar] [CrossRef]
  18. Williams, D.H.; Kalman, J.R. Structural and mode of action studies on the antibiotic vancomycin. Evidence from 270-MHz proton magnetic resonance. J. Am. Chem. Soc. 1977, 99, 2768–2774. [Google Scholar] [CrossRef]
  19. Williamson, M.P.; Williams, D.H. Structure revision of the antibiotic vancomycin. Use of nuclear Overhauser effect difference spectroscopy. J. Am. Chem. Soc. 1981, 103, 6580–6585. [Google Scholar] [CrossRef]
  20. Pearce, C.M.; Williams, D.H. Complete assignment of the 13C NMR spectrum of vancomycin. J. Chem. Soc. Perkin Trans. 2 1995, 153–157. [Google Scholar] [CrossRef]
  21. Tian, Y.; Chong, X.; Yao, S.; Xu, M. Spectral Data Analysis and Identification of Vancomycin Hydrochloride. Front. Chem. 2021, 9, 782. [Google Scholar] [CrossRef] [PubMed]
  22. Sujatha, S.; Praharaj, I. Glycopeptide Resistance in Gram-Positive Cocci: A Review. Interdiscip. Perspect. Infect. Dis. 2012, 2012, 1–10. [Google Scholar] [CrossRef] [Green Version]
  23. Hibbert, F.; Emsley, J. Hydrogen Bonding and Chemical Reactivity. In Advances in Physical Organic Chemistry; Bethell, D., Ed.; Academic Press: San Diego, CA, USA, 1990; Volume 26, pp. 255–379. [Google Scholar]
  24. Yao, N.-H.; He, W.-Y.; Lam, K.S.; Liu, G. Conformational Studies of Resin-Bound Vancomycin and the Complex of Vancomycin and Ac2-l-Lys-d-Ala-d-Ala. J. Comb. Chem. 2004, 7, 123–129. [Google Scholar] [CrossRef]
  25. Ślusarz, R.; Szulc, M.; Madaj, J. Molecular modeling of Gram-positive bacteria peptidoglycan layer, selected glycopeptide antibiotics and vancomycin derivatives modified with sugar moieties. Carbohydr. Res. 2014, 389, 154–164. [Google Scholar] [CrossRef] [PubMed]
  26. Kannan, R.; Harris, C.M.; Harris, T.M.; Waltho, J.P.; Skelton, N.J.; Williams, D.H. Function of the amino sugar and N-terminal amino acid of the antibiotic vancomycin in its complexation with cell wall peptides. J. Am. Chem. Soc. 1988, 110, 2946–2953. [Google Scholar] [CrossRef]
  27. Bax, A.; Davis, D.G. MLEV-17-based two-dimensional homonuclear magnetization transfer spectroscopy. J. Magn. Reson. (1969) 1985, 65, 355–360. [Google Scholar] [CrossRef]
  28. Jeener, J.; Meier, B.; Bachmann, P.; Ernst, R.R. Investigation of exchange processes by two-dimensional NMR spectroscopy. J. Chem. Phys. 1979, 71, 4546–4553. [Google Scholar] [CrossRef]
  29. Macura, S.; Ernst, R.R. Elucidation of cross relaxation in liquids by two-dimensional N.M.R. spectroscopy. Mol. Phys. 1980, 41, 95–117. [Google Scholar] [CrossRef]
  30. Bax, A.; Davis, D.G. Practical aspects of two-dimensional transverse NOE spectroscopy. J. Magn. Reson. (1969) 1985, 63, 207–213. [Google Scholar] [CrossRef]
  31. Palmer, A.G.; Cavanagh, J.; Wright, P.; Rance, M. Sensitivity improvement in proton-detected two-dimensional heteronuclear correlation NMR spectroscopy. J. Magn. Reson. (1969) 1991, 93, 151–170. [Google Scholar] [CrossRef]
  32. Kay, L.; Keifer, P.; Saarinen, T. Pure absorption gradient enhanced heteronuclear single quantum correlation spectroscopy with improved sensitivity. J. Am. Chem. Soc. 1992, 114, 10663–10665. [Google Scholar] [CrossRef]
  33. Schleucher, J.; Schwendinger, M.; Sattler, M.; Schmidt, P.; Schedletzky, O.; Glaser, S.; Griesinger, C. A general enhancement scheme in heteronuclear multidimensional NMR employing pulsed field gradients. J. Biomol. NMR 1994, 4, 301–306. [Google Scholar] [CrossRef]
  34. Fukushi, E. Advanced NMR Approaches for a Detailed Structure Analysis of Natural Products. Biosci. Biotechnol. Biochem. 2006, 70, 1803–1812. [Google Scholar] [CrossRef] [Green Version]
  35. Zwahlen, C.; Legault, P.; Vincent, S.J.F.; Greenblatt, J.; Konrat, R.; Kay, L.E. Methods for Measurement of Intermolecular NOEs by Multinuclear NMR Spectroscopy: Application to a Bacteriophage λ N-Peptide/boxB RNA Complex. J. Am. Chem. Soc. 1997, 119, 6711–6721. [Google Scholar] [CrossRef]
  36. Boyer, R.D.; Johnson, R.; Krishnamurthy, K. Compensation of refocusing inefficiency with synchronized inversion sweep (CRISIS) in multiplicity-edited HSQC. J. Magn. Reson. 2003, 165, 253–259. [Google Scholar] [CrossRef] [PubMed]
  37. Hwang, T.-L.; Shaka, A. Water Suppression That Works. Excitation Sculpting Using Arbitrary Wave-Forms and Pulsed-Field Gradients. J. Magn. Reson. Ser. A 1995, 112, 275–279. [Google Scholar] [CrossRef]
  38. Case, D.A.; Betz, R.M.; Cerutti, D.S.; Cheatham, I.T.E.; Darden, T.A.; Duke, R.E.; Giese, T.J.; Gohlke, H.; Goetz, A.W.; Homeyer, N.; et al. AMBER 16. University of California, San Francisco, CA, USA. 2016. Available online: https://ambermd.org/doc12/Amber16.pdf (accessed on 13 December 2021).
  39. Duan, Y.; Wu, C.; Chowdhury, S.; Lee, M.C.; Xiong, G.; Zhang, W.; Yang, R.; Cieplak, P.; Luo, R.; Lee, T.; et al. A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations. J. Comput. Chem. 2003, 24, 1999–2012. [Google Scholar] [CrossRef]
  40. Lee, M.C.; Duan, Y. Distinguish protein decoys by Using a scoring function based on a new AMBER force field, short molecular dynamics simulations, and the generalized born solvent model. Proteins Struct. Funct. Bioinform. 2004, 55, 620–634. [Google Scholar] [CrossRef] [PubMed]
  41. Kirschner, K.N.; Yongye, A.B.; Tschampel, S.M.; González-Outeiriño, J.; Daniels, C.R.; Foley, B.L.; Woods, R.J. GLYCAM06: A generalizable biomolecular force field. Carbohydrates. J. Comput. Chem. 2008, 29, 622–655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Foley, B.L.; Tessier, M.B.; Woods, R.J. Carbohydrate force fields. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2011, 2, 652–697. [Google Scholar] [CrossRef]
  43. Price, D.J.; Brooks, C.L. A modified TIP3P water potential for simulation with Ewald summation. J. Chem. Phys. 2004, 121, 10096–10103. [Google Scholar] [CrossRef] [PubMed]
  44. Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H.J.C. Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comput. Phys. 1977, 23, 327–341. [Google Scholar] [CrossRef] [Green Version]
  45. Berendsen, H.J.C.; Postma, J.P.M.; Van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684. [Google Scholar] [CrossRef] [Green Version]
  46. Ewald, P.P. Die Berechnung optischer und elektrostatischer Gitterpotentiale. Ann. der Phys. 1921, 369, 253–287. [Google Scholar] [CrossRef] [Green Version]
  47. Darden, T.; Perera, L.; Li, L.; Pedersen, L. New tricks for modelers from the crystallography toolkit: The particle mesh Ewald algorithm and its use in nucleic acid simulations. Structure 1999, 7, R55–R60. [Google Scholar] [CrossRef] [Green Version]
  48. Salomon-Ferrer, R.; Götz, A.W.; Poole, D.; Le Grand, S.; Walker, R.C. Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. J. Chem. Theory Comput. 2013, 9, 3878–3888. [Google Scholar] [CrossRef]
  49. Roe, D.R.; Cheatham, T.E., III. PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theory Comput. 2013, 9, 3084–3095. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The structure of the complex formed by the peptidoglycan with the vancomycin binding pocket.
Figure 1. The structure of the complex formed by the peptidoglycan with the vancomycin binding pocket.
Ijms 23 01146 g001
Figure 2. Fragment of overlapping 1H NMR spectra of: (a) MPP-D-Ser, (b) vancomycin, (c) MPP-D-Ala.
Figure 2. Fragment of overlapping 1H NMR spectra of: (a) MPP-D-Ser, (b) vancomycin, (c) MPP-D-Ala.
Ijms 23 01146 g002
Figure 3. Fragments of the 1H NMR spectra of: (a) MPP-D-Ala, (b) mixture of MPP-D-Ala and vancomycin in a ratio of 1:0.5, (c) mixture of MPP-D-Ala and vancomycin in a ratio 1:1, (d) MPP-D-Ser, (e) mixture of MPP-D-Ser and vancomycin in a ratio of 1:1.
Figure 3. Fragments of the 1H NMR spectra of: (a) MPP-D-Ala, (b) mixture of MPP-D-Ala and vancomycin in a ratio of 1:0.5, (c) mixture of MPP-D-Ala and vancomycin in a ratio 1:1, (d) MPP-D-Ser, (e) mixture of MPP-D-Ser and vancomycin in a ratio of 1:1.
Ijms 23 01146 g003
Figure 4. Diagram of hydrogen bond formation by the L-lysine amide group oxygen atom.
Figure 4. Diagram of hydrogen bond formation by the L-lysine amide group oxygen atom.
Ijms 23 01146 g004
Figure 5. Fragments of HSQC spectra of 1:1 molar ratio of vancomycin mixtures with: (a) MPP-D-Ala, (b) MPP-D-Ser. Arrows show the shift of the signal after adding vancomycin.
Figure 5. Fragments of HSQC spectra of 1:1 molar ratio of vancomycin mixtures with: (a) MPP-D-Ala, (b) MPP-D-Ser. Arrows show the shift of the signal after adding vancomycin.
Ijms 23 01146 g005
Figure 6. Fragments of the 1H NMR spectra of: (a) MPP-D-Ala, (b) mixture of MPP-D-Ala and vancomycin in a ratio of 1:0.5, (c) mixture of MPP-D-Ala and vancomycin in a ratio 1:1, (d) MPP-D-Ser, (e) mixture of MPP-D-Ser and vancomycin in a ratio of 1:1.
Figure 6. Fragments of the 1H NMR spectra of: (a) MPP-D-Ala, (b) mixture of MPP-D-Ala and vancomycin in a ratio of 1:0.5, (c) mixture of MPP-D-Ala and vancomycin in a ratio 1:1, (d) MPP-D-Ser, (e) mixture of MPP-D-Ser and vancomycin in a ratio of 1:1.
Ijms 23 01146 g006
Figure 7. Fragments of the 1H NMR spectra of: (a) vancomycin, (b) MPP-D-Ala and vancomycin in molar ratio 1:1, (c) MPP-D-Ala and vancomycin in molar ratio 1:0.5, (d) MPP-D-Ala, (e) vancomycin, (f) MPP-D-Ser and vancomycin in molar ratio 1:1, (g) MPP-D-Ser.
Figure 7. Fragments of the 1H NMR spectra of: (a) vancomycin, (b) MPP-D-Ala and vancomycin in molar ratio 1:1, (c) MPP-D-Ala and vancomycin in molar ratio 1:0.5, (d) MPP-D-Ala, (e) vancomycin, (f) MPP-D-Ser and vancomycin in molar ratio 1:1, (g) MPP-D-Ser.
Ijms 23 01146 g007
Figure 8. NOESY spectrum of the mixture of MPP-D-Ala and vancomycin in molar ratio 1:1 (the green circle marks the cross-peak indicating the through-space interaction of the protons of the amino group of L-lysine and the benzyl proton 42 of vancomycin).
Figure 8. NOESY spectrum of the mixture of MPP-D-Ala and vancomycin in molar ratio 1:1 (the green circle marks the cross-peak indicating the through-space interaction of the protons of the amino group of L-lysine and the benzyl proton 42 of vancomycin).
Ijms 23 01146 g008
Figure 9. The proposed structure of the muramyl pentapeptide-vancomycin complex (MPP-D-Ala), considering the additional hydrogen bond (marked in red) between the amino group of lysine and the oxygen of the carboxyl and amide groups of vancomycin.
Figure 9. The proposed structure of the muramyl pentapeptide-vancomycin complex (MPP-D-Ala), considering the additional hydrogen bond (marked in red) between the amino group of lysine and the oxygen of the carboxyl and amide groups of vancomycin.
Ijms 23 01146 g009
Figure 10. Newly described interactions’ measured lengths. The interacting parts are L-lysine side chain amine group and vancomycin oxygen atom 40 (distance 1) or L-lysine side chain amine group and vancomycin oxygen atom 41 (distance 2).
Figure 10. Newly described interactions’ measured lengths. The interacting parts are L-lysine side chain amine group and vancomycin oxygen atom 40 (distance 1) or L-lysine side chain amine group and vancomycin oxygen atom 41 (distance 2).
Ijms 23 01146 g010
Figure 11. Control distances confirming hydrogen bonds presence in vancomycin-pentapeptide complex [14,15,26]. The interacting parts are hydrogen atom from the vancomycin C-terminal neighboring amide group and L-lysine backbone (amide) oxygen atom (distance 3), pentapeptide C-terminal D-alanine amide hydrogen and mid-vancomycin amide oxygen atom (distance 4), pentapeptide C-terminal D-alanine carboxyl group and three subsequent amide hydrogen atoms from the vancomycin macrocyclic peptide backbone (distances 5, 6 and 7)—all marked as a black dashed lines in Figure 1 and Figure 9.
Figure 11. Control distances confirming hydrogen bonds presence in vancomycin-pentapeptide complex [14,15,26]. The interacting parts are hydrogen atom from the vancomycin C-terminal neighboring amide group and L-lysine backbone (amide) oxygen atom (distance 3), pentapeptide C-terminal D-alanine amide hydrogen and mid-vancomycin amide oxygen atom (distance 4), pentapeptide C-terminal D-alanine carboxyl group and three subsequent amide hydrogen atoms from the vancomycin macrocyclic peptide backbone (distances 5, 6 and 7)—all marked as a black dashed lines in Figure 1 and Figure 9.
Ijms 23 01146 g011
Figure 12. Distances between vancomycin hydrogen atom 42 (HB) and every L-lysine side chain amine group (HZ1, HZ2 or HZ3).
Figure 12. Distances between vancomycin hydrogen atom 42 (HB) and every L-lysine side chain amine group (HZ1, HZ2 or HZ3).
Ijms 23 01146 g012
Figure 13. Stereo view (crossed eyes representation) of selected vancomycin-muramyl pentapeptide (MPP-D-Ala) MD snapshot with hydrogen bonds marked in yellow.
Figure 13. Stereo view (crossed eyes representation) of selected vancomycin-muramyl pentapeptide (MPP-D-Ala) MD snapshot with hydrogen bonds marked in yellow.
Ijms 23 01146 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ślusarz, R.; Dmochowska, B.; Samaszko-Fiertek, J.; Brzozowski, K.; Madaj, J. NMR and MD Analysis of the Bonding Interaction of Vancomycin with Muramyl Pentapeptide. Int. J. Mol. Sci. 2022, 23, 1146. https://doi.org/10.3390/ijms23031146

AMA Style

Ślusarz R, Dmochowska B, Samaszko-Fiertek J, Brzozowski K, Madaj J. NMR and MD Analysis of the Bonding Interaction of Vancomycin with Muramyl Pentapeptide. International Journal of Molecular Sciences. 2022; 23(3):1146. https://doi.org/10.3390/ijms23031146

Chicago/Turabian Style

Ślusarz, Rafał, Barbara Dmochowska, Justyna Samaszko-Fiertek, Krzysztof Brzozowski, and Janusz Madaj. 2022. "NMR and MD Analysis of the Bonding Interaction of Vancomycin with Muramyl Pentapeptide" International Journal of Molecular Sciences 23, no. 3: 1146. https://doi.org/10.3390/ijms23031146

APA Style

Ślusarz, R., Dmochowska, B., Samaszko-Fiertek, J., Brzozowski, K., & Madaj, J. (2022). NMR and MD Analysis of the Bonding Interaction of Vancomycin with Muramyl Pentapeptide. International Journal of Molecular Sciences, 23(3), 1146. https://doi.org/10.3390/ijms23031146

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