Can we predict the conformational preference of amides?

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

A new series of acylhydrazones derived from metribuzin with modulated herbicidal activity

This paper reports the preparation and herbicidal evaluation of a small library of acylhydrazones based on the synthetic herbicide metribuzin. The hydrazone linkage easily obtained by reaction of metribuzin with aliphatic and aromatic aldehydes, masks efficiently the exocyclic amino group, thereby altering significantly H-bonding with the receptor and increasing the lipophilicity relative to the parent herbicide. The structures of all compounds, including key stereochemical issues on conformation and E/Z configuration around the C[bond, double bond]N bond were thoroughly elucidated by spectroscopic methods, and unambiguously corroborated by X-ray diffraction analysis. The herbicidal assays using an aliphatic and an aromatic acylhydrazone were performed on tomato and rapeseed plants grown in greenhouse. Our results demonstrate, regardless of rate application, that such acylhydrazone formulations do not alter the selectivity of metribuzin. Moreover, the herbicide activity was even higher in the alkyl derivative than that achieved by commercial metribuzin, thus suggesting that this substance can be applied with no need of combination with chemical coadjuvants, unlike most formulations of commercially available herbicides. Therefore, the study shows the promising effect of chemical derivatization of a common herbicide as metribuzin, to improve the herbicide activity without compromising selectivity, and allowing the farmers its use in crop protection safely and effectively.

The Role of Allylic Strain for Conformational Control in Medicinal Chemistry

It is axiomatic in medicinal chemistry that optimization of the potency of a small molecule at a macromolecular target requires complementarity between the ligand and target. In order to minimize the conformational penalty on binding, both enthalpically and entropically, it is therefore preferred to have the ligand preorganized in the bound conformation. In this Perspective, we highlight the role of allylic strain in controlling conformational preferences. Allylic strain was originally described for carbon-based allylic systems, but the same principles apply to other types of structure with sp² or pseudo-sp² arrangements. These systems include benzylic (including heteroaryl methyl) positions, amides, N-aryl groups, aryl ethers, and nucleotides. We have derived torsion profiles from small molecule X-ray structures for these systems. Through multiple examples, we show how these effects have been applied in drug discovery and how they can be used prospectively to influence conformation in the design process.

Triarylmethanes and their Medium-Ring Analogues by Unactivated Truce-Smiles Rearrangement of Benzanilides

Intramolecular nucleophilic aromatic substitution (Truce–Smiles rearrangement) of the anions of 2‐benzyl benzanilides leads to triarylmethanes in an operationally simple manner. The reaction succeeds even without electronic activation of the ring that plays the role of electrophile in the S_NAr reaction, being accelerated instead by the preferred conformation imposed by the tertiary amide tether. The amide substituent of the product may be removed or transformed into alternative functional groups. A ring‐expanding variant (n to n+4) of the reaction provided a route to doubly benzo‐fused medium ring lactams of 10 or 11 members. Hammett analysis returned a ρ value consistent with the operation of a partially concerted reaction mechanism.

Structure-Based Design of a Macrocyclic PROTAC

Constraining a molecule in its bioactive conformation via macrocyclization represents an attractive strategy to rationally design functional chemical probes. While this approach has been applied to enzyme inhibitors or receptor antagonists, to date it remains unprecedented for bifunctional molecules that bring proteins together, such as PROTAC degraders. Herein, we report the design and synthesis of a macrocyclic PROTAC by adding a cyclizing linker to the BET degrader MZ1. A co-crystal structure of macroPROTAC-1 bound in a ternary complex with VHL and the second bromodomain of Brd4 validated the rational design. Biophysical studies revealed enhanced discrimination between the second and the first bromodomains of BET proteins. Despite a 12-fold loss of binary binding affinity for Brd4, macroPROTAC-1 exhibited cellular activity comparable to MZ1. Our findings support macrocyclization as an advantageous strategy to enhance PROTAC degradation potency and selectivity between homologous targets.

Correcting the record: the dimers and trimers of trans-N-methylacetamide

Raman jet spectroscopy reveals three N-methylacetamide molecules organizing into a ring structure, previously overlooked in computations.

Conformational features of secondary N-cyclopropyl amides

NMR studies in conjunction with ab initio calculations revealed unexpected conformational behavior of N-cyclopropylacetamide (1). This secondary amide displays 16-19% E-rotamer (cis) around the carbonyl-nitrogen bond in apolar solvents, in contrast to other aliphatic secondary acetamides in which significant E-rotamer populations are rare due to steric contacts between the substituents on the amide bond. In addition, 1 adopts an ortho conformation around the N-cPr bond instead of the anti conformation generally preferred by secondary acetamides. This distinct conformational behavior was also observed for other secondary N-cyclopropyl amides.

New quinazoline derivatives for telomeric G-quadruplex DNA: effects of an added phenyl group on quadruplex binding ability

To improve the selectivity of indoloquinoline or benzofuroquinoline derivatives, we previously reported several quinazoline derivatives [17]. These compounds could mimic a tetracyclic aromatic system through intramolecular hydrogen bond. Studies showed that these quinazoline derivatives were effective and selective telomeric G-quadruplex ligands. With this encouragement, here we synthesized a series of N-(2-(quinazolin-2-yl)phenyl)benzamide (QPB) compounds as modified quinazoline derivatives. In this modification, a phenyl group was introduced to the aromatic core. The evaluation results showed that part of QPB derivatives had stronger binding ability and better selectivity for telomeric G-quadruplex DNA than LZ-11, the most potential compound of reported quinazoline derivatives. Furthermore, telomerase inhibition of QPB derivatives and their cellular effects were studied.

Disubstituted quinazoline derivatives as a new type of highly selective ligands for telomeric G-quadruplex DNA

A series of 2,4-disubstituted quinazoline derivatives found to be a new type of highly selective ligand to bind with telomeric G-quadruplex DNA, and their biological properties were reported for the first time.Their interactions with telomeric G-quadruplex DNA were evaluated by using fluorescence resonance energy transfer (FRET) melting assay, circular dichroism (CD) spectroscopy, surface plasmon resonance (SPR), nuclear magnetic resonance (NMR), and molecular modeling. Our results showed that these derivatives could well recognize G-quadruplex and have high selectivity toward G-quadruplex over duplex DNA. The structure-activity relationships (SARs) study revealed that the disubstitution of quinazoline and the length of the amide side chain were important for its interaction with the G-quadruplex. Furthermore, telomerase inhibition of the quinazoline derivatives and their cellular effects were studied.

Biological and biophysical properties of the histone deacetylase inhibitor suberoylanilide hydroxamic acid are affected by the presence of short alkyl groups on the phenyl ring

Inhibition of histone deacetylases (HDACs) leads to growth arrest, differentiation, or apoptosis of tumor cell lines, suggesting HDACs as promising targets for cancer therapy. At present, only one HDAC inhibitor (HDACi) is used in therapy: suberoylanilide hydroxamic acid (SAHA). Here, we describe the synthesis and biological evaluation of a new series of compounds derived from SAHA by substituting short alkyl chains at various positions of the phenyl ring. Such modifications induced variable effects ranging from partial loss of activity to increased potency. Through molecular modeling, we describe a possible interaction between HDAC7 proline 809, a residue that is strictly conserved within class 2 enzymes only, and the amide group of HDACi, while nuclear magnetic resonance experiments indicated that dimethyl m-substitution may stabilize the inhibitor in the active site. Our data provide novel information on the structure-activity relationship of HDACi and suggest new ways for developing second generation SAHA-like molecules.

A new model for mapping the peptide backbone: predicting proton chemical shifts in proteins

This paper describes a methodology that correlates experimental chemical shifts (at the alpha proton) of proteins with their geometrical data (both dihedral angles and distances) obtained from 13 representative proteins, which are taken from the Protein Data Bank (PDB) and the BioMagRes Data Bank (BMRB). To this end, the experimentally measured proton chemical shifts of simple amides have been correlated with DFT-based calculated structures, at the B3PW91/6-31G* level. This results in a series of mathematical relationships, which are extrapolated to the above-mentioned proteins giving rise to a modified equation for such skeleta. It is relevant to note that the equation is also supported by a clear comparison with NMR data of a protein beyond the chosen set, such as insulin, even with lower errors. The model also relates the dependence of chemical shifts on hydrophobic and anisotropic effects at the amino acid residues.

Schiff bases from D-glucosamine and aliphatic ketones

Despite the comprehensive literature and enormous versatility of chiral imines derived from aminosugars and aldehydes, the corresponding counterparts generated from ketones remain an underestimated research subject. Filling in the gap, this manuscript sheds light on the synthetic and structural aspects of such substances and updates the few antecedents reported so far.

Estimating the "steric clash" at cis peptide bonds

To account for the scarcity of cis peptide bonds in proteins, especially in nonproline (or secondary amide) cases, a steric-clash argument is often put forward, in a scheme where the R lateral chains are facing parallel one another, and the backbone is kept in an "all- trans"-like arrangement. Although such a steric conflict can be partly relieved through proper adjustment of the backbone dihedral angles, one can try to estimate its associated energy cost. To this end, quantum-chemistry approaches using a differential-torsion protocol and bond-separation-energy analyses are applied to N-ethyl propionamide CH3-CH2-CO-NH-CH2-CH3, regarded as a model capable of exhibiting C beta...C beta interaction as in alanine succession. The calculations provide an increment of 9 kcal/mol, quite close to that obtained in the nearly isostere (gsg) rotamer of n-hexane (10 kcal/mol), suggesting the local effects induced by methyl-methyl contact are similar in both cases. Analogous treatments on larger radicals as encountered in leucine or phenylalanine dimers do not change this increment much, which therefore defines the basic reference per-plaque quota to be overcome along all- cis chains. Explicit modeling indicated it can be reduced by up to a factor of 4.

Anisotropic and hydrogen bonding effects in phenylglyoxamides and mandelamides: theoretical and NMR conformational evaluation

Interesting anisotropic effects were observed for phenylglyoxamides and their respective mandelamides. Such effects were observed in experimental (1)H and (13)C NMR (in CDCl(3), CD(3)OD, and DMSO-d(6) solvents) and in some cases with good correlation to theoretical (1)H and (13)C NMR DFT-GIAO (B3LYP/6-311++G**//B3LYP/6-31G*) calculations. A systematic conformational analysis of these compounds was performed in a two-step methodology, using PM3 and DFT (B3LYP/6-31G*) calculations; with good accomplishment and computational time economy. It was observed that intramolecular hydrogen bonding plays a significant role in the conformation of such compounds. Finally, a geminal nonequivalence of an N-CH(2) moiety, in one of the alkyl side chain (R1 = R2), was found for the tertiary mandelamides studied.

Fluoroolefins as peptide mimetics: a computational study of structure, charge distribution, hydration, and hydrogen bonding

The design of peptide mimetic compounds is greatly facilitated by the identification of functionalities that can act as peptide replacements. The fluoroalkene moiety has recently been employed for that purpose. The purpose of this work is to characterize prototypical fluoroalkenes (fluoroethylene and 2-fluoro-2-butene) with respect to key properties of peptides (amides) including structure, charge distribution, hydration, and hydrogen bonding. The results are compared to those obtained for model peptides (formamide, N-methylacetamide). Calculations have been carried out at the MP2 and B3LYP levels of theory with the 6-311++G(2d,p) and 6-311++G(2d,2p) basis sets. The results suggest that the fluoroalkene is similar in steric requirements to a peptide bond but that there is less charge separation. Calculations of the hydration free energies with the PCM bulk continuum solvent model indicate that the fluoroalkene has much smaller hydration free energies than an amide but that the difference in solvation free energy for cis and trans isomers is comparable. In studies of complexes with water molecules, the fluoroalkene is found to engage in interactions that are analogous to backbone hydrogen-bonding interactions that govern many properties of natural peptides and proteins but with smaller interaction energies. In addition, key structural differences are noted when the fluoroalkene is playing the role of hydrogen-bond acceptor which may have implications in binding, aggregation, and conformational preferences in fluoroalkene peptidomimetics. The issue of cooperativity in hydrogen-bonding interactions in complexes with multiple waters has also been investigated. The fluoroalkene is found to exhibit cooperative effects that mirror those of the peptide but are smaller in magnitude. Thus, pairwise addivitity of interactions appears to more adequately describe the fluoroalkenes than the peptides they are intended to mimic.

Probing helix propensity of monomers within a helical oligomer

A simple strategy is proposed to assess the propensity of a given monomer to follow or not follow a particular helical scheme and to study helix reversal phenomena within helical oligomers. It consists of placing a monomer having a low helix propensity between two conformationally stable helical segments. Helix reversion then occurs preferentially at the site of this monomer, leading to the formation of isomers having P (right-handed) or M (left-handed) helicities at each of the two helical segments. The proportion between the P-P/M-M and P-M isomers is indicative of the stereochemical relations between the inserted monomer and the helical frame. Thus, xylylene or carboxylic acid anhydride spacers have been introduced between two helical oligoamides of 8-amino-2-quinolinecarboxylic acid. Both these spacers presumably lack some of the structural features that confer quinoline units with a high helix propensity. Only one species is observed in solution in the case of an anhydride spacer. This species was shown by x-ray crystallography to be a racemic mixture of P-P and M-M helices. Unexpectedly, the anhydride is consistently incorporated within helical oligoamides. For the xylylene spacer, the P-P/M-M racemate and P-M meso compound are in equal proportions in chloroform, showing that this spacer does not have a propensity to adopt any helical conformation in this solvent. However, the equilibria between the various isomers are shifted in toluene, where one species largely prevails. This species was shown by x-ray crystallography to be the P-P/M-M racemate. Molecular dynamics simulations are consistent with these solution data.

All-cis Cyclic Peptides

Amide bonds -NH-CO- preferentially exist in trans conformations, the cis conformation being thermodynamically unfavored with respect to the trans by about 2 kcal/mol. Yet, the main reason most proteins or peptides cannot be made from cis-peptide plaques only lies in that connecting them into open chains appears to be sterically impracticable. It is possible, however, to build all-cis cyclic peptides in which all cis-plaques are efficiently locked. The present work examines, through quantum calculations, the structural and energetic issues associated with these peculiar arrangements. Systematic exploration at DFT-B3LYP level of the potential-energy surfaces for all-cis cyclopolyglycines cG(n)(c) (n = 2-10,15), and to a lesser extent, for all-cis cyclopolyalanines and all-cis cyclopolyphenylalanines confirms that all these structures are true minima. Optimal ring size occurs around eight peptide units, resulting in planar cG7(c), cG8(c), and cG9(c). In smaller systems, the ring strain is relieved through nonplanar cup-like distortions, particularly in cG6(c). From 10 peptide units and beyond, the ring framework distorts into a saddle-edge shape. These molecules disclose some molecular flexibility, as combinatorial tilting of the plaques may give sets of minima close in energy. Indexes based on isodesmic reactions are used to estimate the energy for joining all-cis or all-trans plaques into cyclic peptides. One of them, the mean plaque-junction energy (MPJE) suggests that within sensible sizes from six peptide units and beyond, all-cis plaque association is almost equally favorable as all-trans one. The frame of radiating cis-amide bonds can be considered as defining a new kind of peptidic material, endowed with specific self-assembling properties.

Computational modelling of the open-state Kv1.5 ion channel block by bupivacaine

Binding of R(+)-bupivacaine to open-state homology models of the mammalian K(v)1.5 membrane ion channel is studied using automated docking and molecular dynamics (MD) methods. Homology models of K(v)1.5 are built using the 3D structures of the KcsA and MthK channels as a template. The packing of transmembrane (TM) alpha-helices in the KcsA structure corresponds to a closed channel state. Opening of the channel may be reached by a conformational transition yielding a bent structure of the internal S6 helices. Our first model of the K(v) open state involves a PVP-type of bending hinge in the internal helices, while the second model corresponds to a Gly-type of bending hinge as found in the MthK channel. Ligand binding to these models is probed using the common local anaesthetic bupivacaine, where blocker binding from the intracellular side of the channel is considered. Conformational properties and partial atomic charges of bupivacaine are determined from quantum mechanical HF/6-31G* calculations with inclusion of solvent effects. The automated docking and MD calculations for the PVP-bend model predict that bupivacaine could bind either in the central cavity or in the PVP region of the channel pore. Linear interaction energy (LIE) estimates of the binding free energies for bupivacaine predict strongest binding to the PVP region. Surprisingly, no binding is predicted for the Gly-bend model. These results are discussed in light of electrophysiological data which show that the K(v)1.5 channel is unable to close when bupivacaine is bound.

A novel thiourea-based dual fluorescent anion receptor with a rigid hydrazine spacer

[structure: see text] A neutral receptor with a rigid hydrazine spacer, N-p-(dimethylamino)benzamido-N'-phenylthiourea, was prepared, and its dual fluorescence in acetonitrile was found to show response toward the presence of anions such as AcO(-), H(2)PO(4)(-), HSO(4)(-), Br(-), Cl(-), F(-), and ClO(4)(-) with high sensitivity and selectivity toward AcO(-).