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].
2. NMR Study
Authors used an NMR spectroscopy to analyze the interactions. As an environment, they decided to use H
2O and D
2O 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 authors 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. Authors 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, they can observe clear changes in their position as well. This is clearly seen in Figure 2, showing this fragment of the spectrum.
Figure 2. Fragment of overlapping 1H NMR spectra of: (a) MPP-D-Ser, (b) vancomycin, (c) MPP-D-Ala.
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.
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.
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.
Figure 4. Diagram of hydrogen bond formation by the L-lysine amide group oxygen atom.
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. Authors 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.
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.
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.
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.
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.
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.
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.
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 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.
While reviewing the literature on the interaction of vancomycin with muramyl pentapeptide, authors 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].
3. MD Study
During the work on the material that authors had previously published
[25], they 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 authors 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).
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).
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 Materials and Methods section) 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.
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.
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). Authors 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).
Authors prepared a plot illustrating how close those protons were during whole MD and presented it in Figure 12.
Figure 12. Distances between vancomycin hydrogen atom 42 (HB) and every L-lysine side chain amine group (HZ1, HZ2 or HZ3).
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 authors 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 they 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 authors could confirm presence of all discussed hydrogen bonds formed between muramyl pentapeptide and vancomycin they selected the complex with lowest total energy and presented it in Figure 13.
Figure 13. Stereo view (crossed eyes representation) of selected vancomycin-muramyl pentapeptide (MPP-D-Ala) MD snapshot with hydrogen bonds marked in yellow.
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. Authors 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 they 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. They confirmed this model in in silico modeling using molecular dynamics and trajectory analysis.