Conformational Investigations in Flexible Molecules Using Orientational NMR Constraints in Combination with 3J-Couplings and NOE Distances

Design of rigid protein-protein interaction inhibitors enables targeting of undruggable Mcl-1

The structure-based design of small-molecule inhibitors targeting protein-protein interactions (PPIs) remains a huge challenge as the drug must bind typically wide and shallow protein sites. A PPI target of high interest for hematological cancer therapy is myeloid cell leukemia 1 (Mcl-1), a prosurvival guardian protein from the Bcl-2 family. Despite being previously considered undruggable, seven small-molecule Mcl-1 inhibitors have recently entered clinical trials. Here, we report the crystal structure of the clinical-stage inhibitor AMG-176 bound to Mcl-1 and analyze its interaction along with clinical inhibitors AZD5991 and S64315. Our X-ray data reveal high plasticity of Mcl-1 and a remarkable ligand-induced pocket deepening. Nuclear Magnetic Resonance (NMR)-based free ligand conformer analysis demonstrates that such unprecedented induced fit is uniquely achieved by designing highly rigid inhibitors, preorganized in their bioactive conformation. By elucidating key chemistry design principles, this work provides a roadmap for targeting the largely untapped PPI class more successfully.

Structures Controlled by Entropy: The Flexibility of Strychnine as Example

To study the flexibility of strychnine, we performed molecular dynamics simulations with orientational tensorial constraints (MDOC). Tensorial constraints are derived from nuclear magnetic resonance (NMR) interaction tensors, for instance, from residual dipolar couplings (RDCs). Used as orientational constraints, they rotate the whole molecule and molecular parts with low rotational barriers. Since the NMR parameters are measured at ambient temperatures, orientational constraints generate conformers that populate the whole landscape of Gibbs free energy. In MDOC, structures are populated that are not only controlled by energy but by the entropy term TΔS of the Gibbs free energy. In the case of strychnine, it is shown that ring conformers are populated, which has not been discussed in former investigations. These conformer populations are not only in accordance with RDCs but fulfill nuclear Overhauser effect (NOE)-derived distance constraints and ³J_HH couplings as well.

Weizhouochrones: Gorgonian-Derived Symmetric Dimers and Their Structure Elucidation Using Anisotropic NMR Combined with DP4+ Probability and CASE-3D

Natural product dimers have intriguing structural features and often have remarkable pharmacological activities. We report here two uncommon marine gorgonian-derived symmetric dimers, weizhouochrones A (1) and B (2), with indenone-derived monomers, that were isolated from the coral Anthogorgia ochracea collected from the South China Sea. These dimers are difficult targets for structure elucidation that solely relies upon conventional NMR data such as NOEs and J-couplings. Here, to explore the application of emerging methods on the structure elucidation of challenging molecules, we explored a number of different anisotropic and computational NMR approaches. The measurements of anisotropic NMR parameters of weizhouochrone A, including residual dipolar couplings (RDCs) and residual chemical shift anisotropy (RCSA), allowed us to successfully determine the planar structure and its relative configuration. This result was corroborated by a computational NMR analysis based on DP4+ probability and computer-assisted 3D structure elucidation (CASE-3D).

Statistical evaluation of simulated NMR data of flexible molecules

A new probability score-named χ-probability-is introduced for evaluating the fit of mixed NMR datasets to calculate molecular model ensembles, in order to answer challenging structural questions such as the determination of stereochemical configurations. Similar to the DP4 parameter, the χ-probability is based on Bayes theorem and expresses the probability that an experimental NMR dataset fits to a given individual within a finite set of candidate structures or configurations. Here, the χ-probability is applied to single out the correct configuration in four example cases, with increasing complexity and conformational mobility. The NMR data (which include RDCs, NOE distances and ³J couplings) are calculated from MDOC (Molecular Dynamics with Orientational Constraints) trajectories and are investigated against experimentally measured data. It is demonstrated that this approach singles out the correct stereochemical configuration with probabilities more than 98%, even for highly mobile molecules. In more demanding cases, a decisive χ-probability test requires that the datasets include high-quality NOE distances in addition to RDC values.

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.

TITANIA: Model-Free Interpretation of Residual Dipolar Couplings in the Context of Organic Compounds

Residual dipolar couplings (RDCs) become increasingly important as additional NMR parameters in the structure elucidation of organic compounds but are usually used in fitting procedures to discriminate between (computed) structures that are in accordance with RDCs and others that can be ruled out. Thus, the determination of configurations requires prior structural information. The direct use of RDCs as restraints to construct structures based on RDCs has only recently begun also in organic compounds. No protocol has been published though that uses the vector and dynamics information available in multialignment data sets directly for the joint determination of conformation and configuration of organic compounds. This is proposed in the current study. We show that by employing these data, even a flat or random start structure converges into the correctly configured structure when employing multiple alignment data sets in our iterative procedure. The requirements in terms of the number of RDCs and alignment media necessary are discussed in detail.

Nuclear magnetic resonance free ligand conformations and atomic resolution dynamics

Knowledge of free ligand conformational preferences (energy minima) and conformational dynamics (rotational energy barriers) of small molecules in solution can guide drug design hypotheses and help rank ideas to bias syntheses towards more active compounds. Visualization of conformational exchange dynamics around torsion angles, by replica exchange with solute tempering molecular dynamics (REST-MD), gives results in agreement with high-resolution 1 H nuclear magnetic resonance (NMR) spectra and complements free ligand conformational analyses. Rotational energy barriers around individual bonds are comparable between calculated and experimental values, making the in-silico method relevant to ranking prospective design ideas in drug discovery programs, particularly across a series of analogs. Prioritizing design ideas, based on calculations and analysis of measurements across a series, efficiently guides rational discovery towards the "right molecules" for effective medicines.

Determination of Configuration and Conformation of a Reserpine Derivative with Seven Stereogenic Centers Using Molecular Dynamics with RDC-Derived Tensorial Constraints*

NMR-based determination of the configuration of complex molecules containing many stereocenters is often not possible using traditional NOE data and coupling patterns. Making use of residual dipolar couplings (RDCs), we were able to determine the relative configuration of a natural product containing seven stereocenters, including a chiral amine lacking direct RDC data. To identify the correct relative configuration out of 32 possible ones, experimental RDCs were used in three different approaches for data interpretation: by fitting experimental data based singular value decomposition (SVD) using a single alignment tensor and either (i) a single conformer or (ii) multiple conformers, or alternatively (iii) using molecular dynamics simulations with tensorial orientational constraints (MDOC). Even though in all three approaches one and the same configuration could be selected and clear discrimination between possible configurations was achieved, the experimental data was not fully satisfied by the methods based on single tensor approaches. While these two approaches are faster, only MDOC is able to fully reproduce experimental results, as the obtained conformational ensemble adequately covers the conformational space necessary to describe the molecule with inherent flexibility.

The simulation of NMR data of flexible molecules: sagittamide A as an example for MD simulations with orientational constraints

The recently developed MDOC (Molecular Dynamics with Orientational Constraints) simulation is applied for the first time to a fully flexible molecule. MDOC simulations aim to single out the naturally existing configuration of molecules and to elucidate conformer populations. The performance of the method was first demonstrated on a well-studied test case, the five-membered ring lactone (α-methylene-γ-butyrolactone). In the case of sagittamide A, one-bond ¹H-¹³C residual dipolar couplings (RDC) are used as orientational constraints that reorient the molecule or parts of it. In addition, NOE distances and ³J scalar couplings are used as constraints. Five possible configurations of sagittamide A (labelled a to e) are considered. One experimental RDC value per flexible unit was available and this was not sufficient to single out one valid configuration. The problem could be solved by including NOE distances as well as ³J couplings as complementary constraints into the MDOC simulations. In accordance with former investigations, we confirmed the configuration a for the natural product. A detailed analysis of conformers of the central chain of 6 chiral carbon atoms could be given by inspecting the MDOC trajectory. The relative abundance of these conformers is crucial in fulfilling all three sets of constraints.