Lincomycin and clindamycin conformations. A fragment shared by macrolides, ketolides and lincosamides determined from TRNOE ribosome-bound conformations

Two important lincosamide antibiotics, lincomycin and clindamycin were studied in the complex state with the bacterial ribosome after a conformational analysis by 1H and 13C NMR spectroscopy and molecular modelling of the unbound molecules. Lincosamide-ribosome interactions were investigated using two-dimensional transferred nuclear Overhauser effect spectroscopy (TRNOESY), resulting in a bound structure compatible with the experimental NMR data. The results compared with the conformational analysis of the substrates in solution indicate that specific conformations are preferred in the bound state. Clindamycin, the more bioactive antibiotic studied, displayed a stronger NMR response than lincomycin showing that in lincosamide-ribosome interactions, a low affinity binding level is associated to the tight binding one and is related to biological activity. This study shows that conformation plays an essential role for the low affinity binding site. Superimposition of lincosamide, macrolide and ketolide bound structures exhibited conformational similarities in a particular fragment which is in agreement with a hypothesis of partial overlapping lincosamide and macrolide binding sites.     

Solution structure of a common substrate mimetic of cyclooxygenase-downstream synthases bound to an engineered thromboxane A2 synthase using a high-resolution NMR technique

Understanding the docking mechanism of the common substrate, prostaglandin H(2) (PGH(2)), into the active sites of different cyclooxygenase(COX)-downstream synthases is a key step toward uncovering the molecular basis of the isomerization of PGH(2) to different prostanoids. A high-resolution NMR spectroscopy was used to determine the conformational changes and solution 3D structure of U44069, a PGH(2) analogue, bound to one of the COX-downstream synthases-an engineered thromboxane A(2) synthase (TXAS). The dynamic binding was clearly observed by (1)D NMR titration. The detailed conformational change and 3D structure of U44069 bound to the TXAS were demonstrated by 2D (1)H NMR experiments using transferred NOEs. Through the assignments for the 2D (1)H NMR spectra, TOCSY, DQF-COSY, NOESY, and the structural calculations based on the NOE constraints, they demonstrated that the widely open conformation with a triangle shape of the free U44069 changed to a compact structure with an oval shape when bound to the TXAS. The putative substrate-binding pocket of the TXAS model fits the conformation of the TXAS-bound U44069 appropriately, but could not fit the free form of U44069. It was the first to provide structural information for the dynamic docking of the PGH(2) mimic of the TXAS in solution, and to imply that PGH(2) undergoes conformational changes when bound to different COX-downstream synthases, which may play important roles in the isomerization of PGH(2) to different prostanoids. The NMR technique can be used as a powerful tool to determine the conformations of PGH(2) bound to other COX-downstream synthases.     

Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin.

Multiple molecular dynamics (MD) simulations of crambin with different initial atomic velocities are used to sample conformations in the vicinity of the native structure. Individual trajectories of length up to 5 ns sample only a fraction of the conformational distribution generated by ten independent 120 ps trajectories at 300 K. The backbone atom conformational space distribution is analyzed using principal components analysis (PCA). Four different major conformational regions are found. In general, a trajectory samples only one region and few transitions between the regions are observed. Consequently, the averages of structural and dynamic properties over the ten trajectories differ significantly from those obtained from individual trajectories. The nature of the conformational sampling has important consequences for the utilization of MD simulations for a wide range of problems, such as comparisons with X-ray or NMR data. The overall average structure is significantly closer to the X-ray structure than any of the individual trajectory average structures. The high frequency (less than 10 ps) atomic fluctuations from the ten trajectories tend to be similar, but the lower frequency (100 ps) motions are different. To improve conformational sampling in molecular dynamics simulations of proteins, as in nucleic acids, multiple trajectories with different initial conditions should be used rather than a single long trajectory.     

Conformational analysis and molecular dynamics simulations of maltose

Constrained conformational energy minimizations have been used to calculate an adiabatic (Φ, ψ) potential energy surface for the disaccharide β-maltose. The inclusion of molecular flexibility in the conformational energy analysis of the disaccharide was found to significantly lower the barriers to conformational transitions, as has been observed previously for other systems. Several low energy wells were identified on the adiabatic surface which differ in energy by small amounts and with low absolute barriers separating them, indicating the possibility of a non-negligible equilibrium population distribution in each well. If such a distribution of conformations existed in the physical system, the conformation observed by NMR NOE measurements would thus be a “virtual” conformation. Molecular dynamics simulations of the motions of this molecule in vacuum were also conducted and indicate that the rate of relaxation of the molecule to the adiabatic surface may be slower than the typical timescale of conformational fluctuations. This effect is apparently due to an unphysical persistence of hydrogen bond patterns in vacuum which does not occur in aqueous solution. Trajectories undergoing transitions between wells were calculated and the effects of such conformational transitions upon the ensemble mean structure, such as might be observed in an NMR experiment, were demonstrated.     

Conformational analysis of Amphotericin B

Within a theoretical approach to the problem of antifungal action of Amphotericin B (AmB), a conformational analysis of the neutral and zwitterionic form of this antibiotic in vacuo was performed by the MM2P and AM1 methods. The analysis was carried out with regard to the mutual orientation of the macrolidic and glycosidic fragments of the molecule, which is defined by the phi and psi steric angles. This orientation defines the overall shape of the molecule and is postulated to be important for the antifungal action of the drug. As a result of the MM2P calculations, phi, psi steric energy and population maps were prepared. Several conformers were found on these maps but only two of them (one each for the zwitterionic and the neutral forms of the antibiotic) were previously observed experimentally for isolated molecules. Our other calculated conformers were not observed experimentally but we propose that they may also appear in the AmB channel structure. The results of our conformational analysis were compared with experimental NMR data (nuclear Overhauser effects between selected hydrogen atoms) obtained previously. New structural information obtained for AmB in the present work will be useful for building a molecular model of AmB-target interactions as well as for designing new derivatives of AmB.     

Convergence in peptide folding simulation: multiple trajectories of a potential AIDS pharmacophore

To examine the conformational properties in aqueous solution of a 15-residue peptide that is a potential pharmacophore for AIDS vaccine development, molecular dynamics simulations were performed in water starting from structures determined experimentally in three different organic solvents. Convergence characteristics of the simulation are examined in Cartesian and conformational spaces. In addition, novel analysis tools are employed including a multidimensional scaling method to represent the distance between trajectory frames. As these methods are based on a variety of physical parameters, they provide a useful cross-check on the structural convergence. Theoretical two-dimensional (2D) 1H-NMR spectra are also generated. These are superficially quite different in appearance, demonstrating that backbone similarities difficult to identify by visual inspection of 2D NMR data can be revealed using the methods described here.     

NMR and molecular modeling studies of the interaction between wheat germ agglutinin and the beta-D-GlcpNAc-(1-->6)-alpha-D-Manp epitope present in glycoproteins of tumor cells

The beta-D-GlcpNAc-(1-->6)-alpha-D-Manp disaccharide is a constituent of highly branched cell-surface glycoconjugates that are malignancy markers. The conformational preference of the disaccharide beta-D-GlcpNAc-(1-->6)-alpha-D-Manp-OMe in solution has been studied by molecular modeling and NMR spectroscopy including 1D (1)H,(1)H T-ROESY experiments and analysis of (3)J(H,H) of the hydroxymethyl group being part of the glycosidic linkage of the disaccharide, which revealed the relative populations of the omega torsion angle as gt = 0.60, gg = 0.35, and tg = 0.05. Good agreement was obtained between the effective proton-proton distances from the experiment and those obtained by molecular modeling when the flexibility at the omega torsion angle was taken into account. Molecular modeling of the disaccharide in the binding sites of the lectin wheat germ agglutinin indicates that several conformations could be adopted in the bound state. (1)H NMR and transfer NOESY experiments confirmed that binding took place, and trans-glycosidic proton-proton interactions indicated that a conformational preference was present in the bound state, as observed by the relative change of the NOEs from H1' to H6(pro-R) and H6(pro-S). STD NMR experiments showed that binding occurred in the region of the N-acetyl group of the terminal sugar residue. In addition, the O-methyl group received saturation transfer because of the proximity to the protein. (1)H,(1)H NOEs indicated that the two methyl groups were close in space, as observed in only one of the predicted bound conformations. Experimental and theoretical data therefore agree that one conformation with a gt conformation of the hydroxymethyl group and a negative sign for the psi torsion angle is indeed selected by the lectin upon binding.     

BI-BII transitions in B-DNA.

Molecular modelling is used to study the conformational and energetic aspects of BI-BII transitions within the backbone of a B-DNA dodecamer d(CATGACGTCATG) whose fine structure has previously been determined by molecular modelling combined with NMR spectroscopy. It is shown that while the dodecamer under investigation does not contain any BII junctions, the central CpG step can most easily undergo the transition. More generally, it is also found that the base sequence and hence the backbone geometry of a DNA segment, strongly influences both the conformational impact of the transition, the associated energy barrier and the stability of the resulting BII state.     

Interaction of a <Emphasis Type="Italic">Salmonella enteritidis</Emphasis> O-antigen octasaccharide with the phage P22 tailspike protein by NMR spectroscopy and docking studies

The tailspike protein P22 recognizes an octasaccharide derived from the O-antigen polysaccharide of Salmonella enteritidis in a shallow groove and molecular docking successfully identifies this binding region on the protein surface. Analysis by 2D ¹H,¹H-T-ROESY and transferred NOESY NMR experiments indicate that the bound octasaccharide ligand has a conformation similar to that observed in solution. The results from a saturation transfer difference NMR experiment show that a large number of protons in the octasaccharide are in close contact with the protein as a result of binding. A comparison of the crystal structure of the complex and a molecular dynamics simulation of the octasaccharide with explicit water molecules suggest that only minor conformational changes are needed upon binding to the tailspike protein.     

Synthesis, absolute configuration, conformational analysis and binding affinity properties of enantiomeric forms of DAU 5750, a novel M1–M3 muscarinic receptor antagonist

Both the enantiomeric forms of DAU 5750, a novel muscarinic receptor antagonist, have been synthesized in order to assess the relevance of configurational/conformational features for high affinity binding to muscarinic receptor subtypes. The attribution of absolute stereochemistry and conformational analysis by means of molecular modelling and NMR techniques are also reported.     

Protein flexibility: multiple molecular dynamics simulations of insulin chain B

Multiple molecular dynamics simulations totaling more than 100 ns were performed on chain B of insulin in explicit solvent at 300 K and 400 K. Despite some individual variations, a comparison of the protein dynamics of each simulation showed similar trends and most structures were consistent with NMR experimental values, even at the elevated temperature. The importance of packing interactions in determining the conformational transitions of the protein was observed, sometimes resulting in conformations induced by localized hydrophobic interactions. The high temperature simulation generated a more diverse range of structures with similar elements of secondary structure and populated conformations to the simulations at room temperature. A broad sampling of the conformational space of insulin chain B illustrated a wide range of conformational states with many transitions at room temperature in addition to the conformational states observed experimentally. The T-state conformation associated with insulin activity was consistently present and a possible mechanism of behavior was suggested.     

Simulation study on the disordered state of an Alzheimer's beta amyloid peptide Abeta(12 36) in water consisting of random-structural, beta-structural, and helical clusters

The monomeric Alzheimer's beta amyloid peptide, Abeta, is known to adopt a disordered state in water at room temperature, and a circular dichroism (CD) spectroscopy experiment has provided the secondary-structure contents for the disordered state: 70% random, 25% beta-structural, and 5% helical. We performed an enhanced conformational sampling (multicanonical molecular dynamics simulation) of a 25-residue segment (residues 12-36) of Abeta in explicit water and obtained the conformational ensemble over a wide temperature range. The secondary-structure contents calculated from the conformational ensemble at 300 degrees K reproduced the experimental secondary-structure contents. The constructed free-energy landscape at 300 degrees K was not plain but rugged with five clearly distinguishable clusters, and each cluster had its own characteristic tertiary structure: a helix-structural cluster, two beta-structural clusters, and two random-structural clusters. This indicates that the contribution from the five individual clusters determines the secondary-structure contents experimentally measured. The helical cluster had a similarity with a stable helical structure for monomeric Abeta in 2,2,2-trifluoroethanol (TFE)/water determined by an NMR experiment: The positions of helices in the helical cluster were the same as those in the NMR structure, and the residue-residue contact patterns were also similar with those of the NMR structure. The cluster-cluster separation in the conformational space indicates that free-energy barriers separate the clusters at 300 degrees K. The two beta-structural clusters were characterized by different strand-strand hydrogen-bond (H-bond) patterns, suggesting that the free-energy barrier between the two clusters is due to the H-bond rearrangements.     

Solution structure by nuclear magnetic resonance of the two lantibiotics 97518 and NAI‐107

Lantibiotics 97518 and NAI‐107, produced by the related genera Planomonospora and Microbispora respectively, are members of a family of nisin‐related compounds. They represent promising compounds to treat infections caused by multiresistant Gram‐positive pathogens. Despite their similar structure and a similar antibacterial spectrum, the two lantibiotics exhibit significant differences in their potency. To gain an insight into the structure–activity relationships, their conformational properties in solution are determined by NMR. After carrying out an NOE analysis of 2D ¹H NMR spectra, high‐resolution 3D structures are determined using molecular dynamics simulations. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

An improved ultrafast 2D NMR experiment: Towards atom-resolved real-time studies of protein kinetics at multi-Hz rates

Multidimensional NMR spectroscopy is a well-established technique for the characterization of structure and fast-time-scale dynamics of highly populated ground states of biological macromolecules. The investigation of short-lived excited states that are important for molecular folding, misfolding and function, however, remains a challenge for modern biomolecular NMR techniques. Off-equilibrium real-time kinetic NMR methods allow direct observation of conformational or chemical changes by following peak positions and intensities in a series of spectra recorded during a kinetic event. Because standard multidimensional NMR methods required to yield sufficient atom-resolution are intrinsically time-consuming, many interesting phenomena are excluded from real-time NMR analysis. Recently, spatially encoded ultrafast 2D NMR techniques have been proposed that allow one to acquire a 2D NMR experiment within a single transient. In addition, when combined with the SOFAST technique, such ultrafast experiments can be repeated at high rates. One of the problems detected for such ultrafast protein NMR experiments is related to the heteronuclear decoupling during detection with interferences between the pulses and the oscillatory magnetic field gradients arising in this scheme. Here we present a method for improved ultrafast data acquisition yielding higher signal to noise and sharper lines in single-scan 2D NMR spectra. In combination with a fast-mixing device, the recording of ¹H–¹⁵N correlation spectra with repetition rates of up to a few Hertz becomes feasible, enabling real-time studies of protein kinetics occurring on time scales down to a few seconds.     

Analysis of protein-carbohydrate interaction at the lower size limit of the protein part (15-mer peptide) by NMR spectroscopy,electrospray ionization mass spectrometry,and molecular modeling

Structural analysis of minimally sized lectins will offer insights into fundamentals of intermolecular recognition and potential for biomedical applications. We thus moved significantly beyond the natural limit of lectin size to determine the structure of synthetic mini-lectins in solution, their carbohydrate selectivity and the impact of ligand binding on their conformational behavior. Using three disaccharide (Thomsen-Friedenreich antigen; Gal beta 1,3GalNAc alpha 1,R)-binding pentadecapeptides without internal disulfide bridges as role models, we successfully tested a combined strategy with different techniques of NMR spectroscopy, electrospray ionization mass spectrometry, and molecular modeling. In solution, the peptides invariably displayed flexibility with rather limited restrictions, shown by NMR experiments including nearly complete resonance assignments and molecular dynamics simulations. The occurrence of aromatic/nonpolar amino acids in the sequence did not lead to formation of a hydrophobic core known from microbial chitinase modules. Selectivity of disaccharide binding was independently observed by mass spectrometry and NMR analysis. Specific ligand interaction yielded characteristic NMR signal alterations but failed to reduce conformational flexibility significantly. We have thereby proven effectiveness of our approach to analyze even low-affinity interactions (not restricted to carbohydrates as ligands). It will be useful to evaluate the impact of rational manipulation of lead peptide sequences.     

Sequence-specific resonance assignment and conformational analysis of subtilin by 2D NMR.

Subtilin, a 32-amino acid peptide with potent antimicrobial activity, has been isolated from Bacillus subtilis ATCC6633. The chemical structure has been confirmed by the unambiguous sequence-specific assignment of its 1H NMR spectrum. Detailed NMR analysis revealed that subtilin is a rather flexible molecule; the only observed conformational contraints were those imposed by the cyclic structures created by the lanthionine and 3-methyllanthionine residues. These results suggest that in aqueous solution subtilin and the homologous peptide nisin have similar conformations.     

Design,synthesis and conformational analysis of hGM-CSF[13-31]-Gly-Pro-Gly-[103-116]

On the basis of the X-ray structure and results from structure--activity relationship studies, the following GM–CSF analogue was designed and synthesized by solid-phase methodology: hGM–CSF[13-31]-Gly-Pro-Gly-[103–116]-NH₂. This analogue was constructed to comprise helices A and D of the native hGM–CSF, covalently linked in an antiparallel orientation by the tripeptide spacer Gly-Pro-Gly, which is known as a turn-inducing sequence. The conformational analysis of the analogue by CD spectroscopy revealed an essentially random structure in water, while α-helix formation was observed upon addition of TFE. In 40% TFE the helix content was ∼45%. By two-dimensional NMR experiments in 1:1 water/trifluoroethanol mixture two helical sequences were identified comprising the segments corresponding to helix A and helix D. In addition to medium-range NOESY connectivities, a long-range cross-peak was found involving the leucine residues at positions 13 and 35. Based on the experimentally derived data (54 NOEs), the structure was refined by restrained molecular dynamics simulations over 120 ps at various temperatures. A representative conformation derived from the computer simulation is mainly characterized by two helical segments connected by a loop region. The overall three-dimensional structure of the analogue is comparable to the X-ray structure of hGM–CSF in that helices A and D are oriented in an antiparallel fashion, forming a two α-helix bundle. Nevertheless, there are small differences in the topology of the helices between the solution structure of the designed analogue and the X-ray structure of hGM–CSF. The possible implications of these conformational features at the effects of biological activity are discussed. © 1997 European Peptide Society and John Wiley & Sons, Ltd. J. Pep. Sci.3: 323–335 No. of Figures: 10. No. of Tables: 5. No. of References: 46