Assessing the Ability of Spectroscopic Methods to Determine the Difference in the Folding Propensities of Highly Similar β-Hairpins

Going Viral: An Investigation into the Chameleonic Behaviour of Antiviral Compounds

The ability to adjust conformations in response to the polarity of the environment, i.e. molecular chameleonicity, is considered to be important for conferring both high aqueous solubility and high cell permeability to drugs in chemical space beyond Lipinski's rule of 5. We determined the conformational ensembles populated by the antiviral drugs asunaprevir, simeprevir, atazanavir and daclatasvir in polar (DMSO-d₆ ) and non-polar (chloroform) environments with NMR spectroscopy. Daclatasvir was fairly rigid, whereas the first three showed large flexibility in both environments, that translated into major differences in solvent accessible 3D polar surface area within each conformational ensemble. No significant differences in size and polar surface area were observed between the DMSO-d₆ and chloroform ensembles of these three drugs. We propose that such flexible compounds are characterized as "partial molecular chameleons" and hypothesize that their ability to adopt conformations with low polar surface area contributes to their membrane permeability and oral absorption.

Hydrocarbon Macrocycle Conformer Ensembles and 13 C-NMR Spectra

NMR as a routine analytical method provides important three-dimensional structure information of compounds in solution. Here we apply the recently released CRENSO computational workflow for the automated generation of conformer ensembles to the quantum mechanical calculation of ¹³ C-NMR spectra of a series of flexible cycloalkanes up to C₂₀ H₄₀ . We evaluate the computed chemical shifts in comparison with corresponding experimental data in chloroform. It is shown that accurate and properly averaged theoretical NMR data can be obtained in about a day of computation time on a standard workstation computer. The excellent agreement between theory and experiment enables one to deduce preferred conformations of large, non-rigid macrocycles under ambient conditions from our automated procedure.

Pushing the Limits of Characterising a Weak Halogen Bond in Solution

Detection and characterisation of very weak, non-covalent interactions in solution is inherently challenging. Low affinity, short complex lifetime and a constant battle against entropy brings even the most sensitive spectroscopic methods to their knees. Herein we introduce a strategy for the accurate experimental description of weak chemical forces in solution. Its scope is demonstrated by the detailed geometric and thermodynamic characterisation of the weak halogen bond of a non-fluorinated aryl iodide and an ether oxygen (0.6 kJ mol^-1 ). Our approach makes use of the entropic advantage of studying a weak force intramolecularly, embedded into a cooperatively folding system, and of the combined use of NOE- and RDC-based ensemble analyses to accurately describe the orientation of the donor and acceptor sites. Thermodynamic constants (ΔG, ΔH and ΔS), describing the specific interaction, were derived from variable temperature chemical shift analysis. We present a methodology for the experimental investigation of remarkably weak halogen bonds and other related weak forces in solution, paving the way for their improved understanding and strategic use in chemistry and biology.

Employing complementary spectroscopies to study the conformations of an epimeric pair of side-chain stapled peptides in aqueous solution

Understanding the conformational preferences of free ligands in solution is often necessary to rationalize structure-activity relationships in drug discovery. Herein, we examine the conformational behavior of an epimeric pair of side-chain stapled peptides that inhibit the FAD dependent amine oxidase lysine specific demethylase 1 (LSD1). The peptides differ only at a single stereocenter, but display a major difference in binding affinity. Their Raman optical activity (ROA) spectra are most likely dominated by the C-terminus, obscuring the analysis of the epimeric macrocycle. By employing NMR spectroscopy, we show a difference in conformational behavior between the two compounds and that the LSD1 bound conformation of the most potent compound is present to a measurable extent in aqueous solution. In addition, we illustrate that Molecular Dynamics (MD) simulations produce ensembles that include the most important solution conformations, but that it remains problematic to identify relevant conformations with no a priori knowledge from the large conformational pool. Furthermore, this work highlights the importance of understanding the scope and limitations of the available techniques for conducting conformational analyses. It also emphasizes the importance of conformational selection of a flexible ligand in molecular recognition.

A Chemical Biology Approach to Understanding Molecular Recognition of Lipid II by Nisin(1-12): Synthesis and NMR Ensemble Analysis of Nisin(1-12) and Analogues

Natural products that target lipid II, such as the lantibiotic nisin, are strategically important in the development of new antibacterial agents to combat the rise of antimicrobial resistance. Understanding the structural factors that govern the highly selective molecular recognition of lipid II by the N-terminal region of nisin, nisin(1-12), is a crucial step in exploiting the potential of such compounds. In order to elucidate the relationships between amino acid sequence and conformation of this bicyclic peptide fragment, we have used solid-phase peptide synthesis to prepare two novel analogues of nisin(1-12) in which the dehydro residues have been replaced. We have carried out an NMR ensemble analysis of one of these analogues and of the wild-type nisin(1-12) peptide in order to compare the conformations of these two bicyclic peptides. Our analysis has shown the effects of residue mutation on ring conformation. We have also demonstrated that the individual rings of nisin(1-12) are pre-organised to an extent for binding to the pyrophosphate group of lipid II, with a high degree of flexibility exhibited in the central amide bond joining the two rings.

Conformational Sampling of Macrocyclic Drugs in Different Environments: Can We Find the Relevant Conformations?

Conformational flexibility is a major determinant of the properties of macrocycles and other drugs in beyond rule of 5 (bRo5) space. Prediction of conformations is essential for design of drugs in this space, and we have evaluated three tools for conformational sampling of a set of 10 bRo5 drugs and clinical candidates in polar and apolar environments. The distance-geometry based OMEGA was found to yield ensembles spanning larger structure and property spaces than the ensembles obtained by MOE-LowModeMD (MOE) and MacroModel (MC). Both MC and OMEGA but not MOE generated different ensembles for polar and apolar environments. All three conformational search methods generated conformers similar to the crystal structure conformers for 9 of the 10 compounds, with OMEGA performing somewhat better than MOE and MC. MOE and OMEGA found all six conformers of roxithromycin that were identified by NMR in aqueous solutions, whereas only OMEGA sampled the three conformers observed in chloroform. We suggest that characterization of conformers using molecular descriptors, e.g., the radius of gyration and polar surface area, is preferred to energy- or root-mean-square deviation-based methods for selection of biologically relevant conformers in drug discovery in bRo5 space.

Halogen Bonding: A Powerful Tool for Modulation of Peptide Conformation

Halogen bonding is a weak chemical force that has so far mostly found applications in crystal engineering. Despite its potential for use in drug discovery, as a new molecular tool in the direction of molecular recognition events, it has rarely been assessed in biopolymers. Motivated by this fact, we have developed a peptide model system that permits the quantitative evaluation of weak forces in a biologically relevant proteinlike environment and have applied it for the assessment of a halogen bond formed between two amino acid side chains. The influence of a single weak force is measured by detection of the extent to which it modulates the conformation of a cooperatively folding system. We have optimized the amino acid sequence of the model peptide on analogues with a hydrogen bond-forming site as a model for the intramolecular halogen bond to be studied, demonstrating the ability of the technique to provide information about any type of weak secondary interaction. A combined solution nuclear magnetic resonance spectroscopic and computational investigation demonstrates that an interstrand halogen bond is capable of conformational stabilization of a β-hairpin foldamer comparable to an analogous hydrogen bond. This is the first report of incorporation of a conformation-stabilizing halogen bond into a peptide/protein system, and the first quantification of a chlorine-centered halogen bond in a biologically relevant system in solution.