Atropisomeric xanthines: Synthesis, stereodynamics and absolute configuration

During the thesis period a new class of atropisomeric xanthine derivatives has been studied. We decided to focus our attention on these purine bases because of their various biological activities, that could play an important role in the discovery of new bioactive atropisomers. The synthesized compounds bear an Aryl-N chiral axis in position 1 of the xanthine scaffold, around which the rotation is prevented by the presence of bulky ortho substituents. Through a retro synthetic analysis we synthesized three atropisomeric structures bearing in position 1 of the purine scaffold respectively an o-tolyl, o-nitrophenyl and a 1-naphthyl group. The conformational studies by DFT simulations showed that the interconversion energy barrier between the two available skewed conformations is higher enough to obtain thermally stable atropisomers. After the separation of the atropisomers, the experimental energy of interconversion was investigated by means of kinetic studies following the thermal racemization process using an enantioselective HPLC column. The absolute configuration of each atropisomer was assigned by experimental ECD analysis and TD-DFT simulations of the ECD spectra.

Synthesis and dynamic study of atropisomeric compounds containing boron-carbon bond

In this thesis we studied the stereodynamic behavior of 1,2-azaborines variously substituted on boron (7a, 7b, 13). Depending on the hindrance of the asymmetric aryl substituent the resulting conformations could be stereolabile or configurationally stable. Through dynamic NMR and lineshape simulation, the energy rotational barriers of the different conformers are obtained. When the barrier is higher than 22-23 kcal/mol stable atropisomers that are fisically separable could be obtained (case of compound 13) and the free activation energy barrier is determinable by kinetic analysis. Absolute configuration of two atropisomers were assigned by comparison between computational calculations and experimental ECD. Isosteric compound 21 is then synthesized in order to compare the rotational barrier around B-Caryl with the one around Cnaphth-Caryl bond.

Novel octreotide dicarba-analogues with high affinity and different selectivity for somatostatin receptors

A limited set of novel octreotide dicarba-analogues with non-native aromatic side chains in positions 7 and/or 10 were synthesized. Their affinity toward the ssts1-5 was determined. Derivative 4 exhibited a pan-somatostatin activity, except sst4, and derivative 8 exhibited high affinity and selectivity toward sst5. Actually, compound 8 has similar sst5 affinity (IC50 4.9 nM) to SRIF-28 and octreotide. Structure-activity relationships suggest that the Z geometry of the double-bond bridge is that preferred by the receptors. The NMR study on the conformations of these compounds in SDS(-d25) micelles solution shows that all these analogues have the pharmacophore beta-turn spanning Xaa7-D-Trp8-Lys9-Yaa10 residues. Notably, the correlation between conformation families and affinity data strongly indicates that the sst5 selectivity is favored by a helical conformation involving the C-terminus triad, while a pan-SRIF mimic activity is based mainly on a conformational equilibrium between extended and folded conformational states.

The contribution of surface residues to membrane binding and ligand transfer by the -tocopherol transfer protein ( -TTP)

Previous work has shown that the -tocopherol transfer protein ( -TTP) can bind to vesicular or immobilized phospholipid membranes. Revealing the molecular mechanisms by which -TTP associates with membranes is thought to be critical to understanding its function and role in the secretion of tocopherol from hepatocytes into the circulation. Calculations presented in the Orientations of Proteins in Membranes database have provided a testable model for the spatial arrangement of -TTP and other CRAL-TRIO family proteins with respect to the lipid bilayer. These calculations predicted that a hydrophobic surface mediates the interaction of -TTP with lipid membranes. To test the validity of these predictions, we used site-directed mutagenesis and examined the substituted mutants with regard to intermembrane ligand transfer, association with lipid layers and biological activity in cultured hepatocytes. Substitution of residues in helices A8 (F165A and F169A) and A10 (I202A, V206A and M209A) decreased the rate of intermembrane ligand transfer as well as protein adsorption to phospholipid bilayers. The largest impairment was observed upon mutation of residues that are predicted to be fully immersed in the lipid bilayer in both apo (open) and holo (closed) conformations such as Phe165 and Phe169. Mutation F169A, and especially F169D, significantly impaired -TTP-assisted secretion of -tocopherol outside cultured hepatocytes. Mutation of selected basic residues (R192H, K211A, and K217A) had little effect on transfer rates, indicating no significant involvement of nonspecific electrostatic interactions with membranes.

The binding mode of side chain- and C3-modified epothilones to tubulin

The tubulin-binding mode of C3- and C15-modified analogues of epothilone A (Epo A) was determined by NMR spectroscopy and computational methods and compared with the existing structural models of tubulin-bound natural Epo A. Only minor differences were observed in the conformation of the macrocycle between Epo A and the C3-modified analogues investigated. In particular, 3-deoxy- (compound 2) and 3-deoxy-2,3-didehydro-Epo A (3) were found to adopt similar conformations in the tubulin-binding cleft as Epo A, thus indicating that the 3-OH group is not essential for epothilones to assume their bioactive conformation. None of the available models of the tubulin-epothilone complex is able to fully recapitulate the differences in tubulin-polymerizing activity and microtubule-binding affinity between C20-modified epothilones 6 (C20-propyl), 7 (C20-butyl), and 8 (C20-hydroxypropyl). Based on the results of transferred NOE experiments in the presence of tubulin, the isomeric C15 quinoline-based Epo B analogues 4 and 5 show very similar orientations of the side chain, irrespective of the position of the nitrogen atom in the quinoline ring. The quinoline side chain stacks on the imidazole moiety of beta-His227 with equal efficiency in both cases, thus suggesting that the aromatic side chain moiety in epothilones contributes to tubulin binding through strong van der Waals interactions with the protein rather than hydrogen bonding involving the heteroaromatic nitrogen atom. These conclusions are in line with existing tubulin polymerization and microtubule-binding data for 4, 5, and Epo B.

Molecular Dynamics Simulation of a PNA·DNA·PNA Triple Helix in Aqueous Solution

Molecular dynamics simulations have been used to explore the conformational flexibility of a PNA·DNA·PNA triple helix in aqueous solution. Three 1.05 ns trajectories starting from different but reasonable conformations have been generated and analyzed in detail. All three trajectories converge within about 300 ps to produce stable and very similar conformational ensembles, which resemble the crystal structure conformation in many details. However, in contrast to the crystal structure, there is a tendency for the direct hydrogen-bonds observed between the amide hydrogens of the Hoogsteen-binding PNA strand and the phosphate oxygens of the DNA strand to be replaced by water-mediated hydrogen bonds, which also involve pyrimidine O2 atoms. This structural transition does not appear to weaken the triplex structure but alters groove widths and so may relate to the potential for recognition of such structures by other ligands (small molecules or proteins). Energetic analysis leads us to conclude that the reason that the hybrid PNA/DNA triplex has quite different helical characteristics from the all-DNA triplex is not because the additional flexibility imparted by the replacement of sugar−phosphate by PNA backbones allows motions to improve base-stacking but rather that base-stacking interactions are very similar in both types of triplex and the driving force comes from weak but definate conformational preferences of the PNA strands.

Quantum Mechnical Investion of Cyclic 3',5'- Adenosine Monophosphate, the Second Hormonal Messenger

Full geometry optimizations using the PM3, AM1, 3-21G∗/HF and 6-31G∗/HF levels of theory were conducted on the syn and anti conformations of cyclic3′,5′-adenosine monophosphate (cAMP). Comparison of the anti crystal structures with the semiempirical and ab initio results revealed that the ab initio results agree well with the experimental results. The results of semiempirical calculations are in qualitative agreement with experimental and ab initio values, with the exception of the glycosyl torsion angle for the anti conformer. Sugar puckering, which is not handled properly by semiempirical methods for unconstrained sugars, nucleosides, nucleotides and nucleotide base pairs, is modeled reasonably well by the semiempirical methods for cAMP. This improvement results from the constraints introduced by the cyclization of AMP to form the phosphodiester.

Molecular Dynamics Simulations of the d(TAT) Triple Helix

Molecular dynamics (MD) simulations have been used to study the dynamical and time-averaged characteristics of the DNA triple helix d(T)10âd(A)10âd(T)10. The structures sampled during the trajectory resemble closely the B-type model for the DNA triplex proposed on the basis of NMR data, although there are some subtle differences. Alternative P- and A-type conformations for the triplex, suggested from X-ray experiments, are not predicted to contribute significantly to the structure of the DNA triplex in solution. Comparison with the best available experimental data supports the correctnes of the MD-generated structures. The analysis of the collected data gives a detailed picture of the characteristics of triple-helix DNA. A new and interesting pattern of hydration, specific for triplex DNA, is an important observation. The results suggest that molecular dynamics can be useful for the study of novel nucleic acid structures.

Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.

First-principle studies of intermolecular and intramolecular catalysis of protonated cocaine

We have performed a series of first-principles electronic structure calculations to examine the reaction pathways and the corresponding free energy barriers for the ester hydrolysis of protonated cocaine in its chair and boat conformations. The calculated free energy barriers for the benzoyl ester hydrolysis of protonated chair cocaine are close to the corresponding barriers calculated for the benzoyl ester hydrolysis of neutral cocaine. However, the free energy barrier calculated for the methyl ester hydrolysis of protonated cocaine in its chair conformation is significantly lower than for the methyl ester hydrolysis of neutral cocaine and for the dominant pathway of the benzoyl ester hydrolysis of protonated cocaine. The significant decrease of the free energy barrier, ∼4 kcal/mol, is attributed to the intramolecular acid catalysis of the methyl ester hydrolysis of protonated cocaine, because the transition state structure is stabilized by the strong hydrogen bond between the carbonyl oxygen of the methyl ester moiety and the protonated tropane N. The relative magnitudes of the free energy barriers calculated for different pathways of the ester hydrolysis of protonated chair cocaine are consistent with the experimental kinetic data for cocaine hydrolysis under physiologic conditions. Similar intramolecular acid catalysis also occurs for the benzoyl ester hydrolysis of (protonated) boat cocaine in the physiologic condition, although the contribution of the intramolecular hydrogen bonding to transition state stabilization is negligible. Nonetheless, the predictability of the intramolecular hydrogen bonding could be useful in generating antibody-based catalysts that recruit cocaine to the boat conformation and an analog that elicited antibodies to approximate the protonated tropane N and the benzoyl O more closely than the natural boat conformer might increase the contribution from hydrogen bonding. Such a stable analog of the transition state for intramolecular catalysis of cocaine benzoyl-ester hydrolysis was synthesized and used to successfully elicit a number of anticocaine catalytic antibodies.

Thermodynamics of the Hydroxyl Radical Addition to Isoprene

Oxidation of isoprene by the hydroxyl radical leads to tropospheric ozone formation. Consequently, a more complete understanding of this reaction could lead to better models of regional air quality, a better understanding of aerosol formation, and a better understanding of reaction kinetics and dynamics. The most common first step in the oxidation of isoprene is the formation of an adduct, with the hydroxyl radical adding to one of four unsaturated carbon atoms in isoprene. In this paper, we discuss how the initial conformations of isoprene, s-trans and s-gauche, influences the pathways to adduct formation. We explore the formation of pre-reactive complexes at low and high temperatures, which are often invoked to explain the negative temperature dependence of this reaction’s kinetics. We show that at higher temperatures the free energy surface indicates that a pre-reactive complex is unlikely, while at low temperatures the complex exists on two reaction pathways. The theoretical results show that at low temperatures all eight pathways possess negative reaction barriers, and reaction energies that range from −36.7 to −23.0 kcal·mol−1. At temperatures in the lower atmosphere, all eight pathways possess positive reaction barriers that range from 3.8 to 6.0 kcal·mol−1 and reaction energies that range from −28.8 to −14.4 kcal·mol−1.

A closer look at the high affinity benzodiazepine binding site on GABAA receptors

Ligands of the benzodiazepine binding site of the GABA(A) receptor come in three flavors: positive allosteric modulators, negative allosteric modulators and antagonists all of which can bind with high affinity. The GABA(A) receptor is a pentameric protein which forms a chloride selective ion channel and ligands of the benzodiazepine binding site stabilize three different conformations of this protein. Classical benzodiazepines exert a positive allosteric effect by increasing the apparent affinity of channel opening by the agonist γ-aminobutyric acid (GABA). We concentrate here on the major adult isoform, the α(1)β(2)γ(2) GABA(A) receptor. The classical binding pocket for benzodiazepines is located in a subunit cleft between α(1) and γ(2) subunits in a position homologous to the agonist binding site for GABA that is located between β(2) and α(1) subunits. We review here approaches to this picture. In particular, point mutations were performed in combination with subsequent analysis of the expressed mutant proteins using either electrophysiological techniques or radioactive ligand binding assays. The predictive power of these methods is assessed by comparing the results with the predictions that can be made on the basis of the recently published crystal structure of the acetylcholine binding protein that shows homology to the N-terminal, extracellular domain of the GABA(A) receptor. In addition, we review an approach to the question of how the benzodiazepine ligands are positioned in their binding pocket. We also discuss a newly postulated modulatory site for benzodiazepines at the α(1)/β(2) subunit interface, homologous to the classical benzodiazepine binding pocket.

Conformationally Controlled Electron Delocalization in n-Type Rods: Synthesis, Structure, and Optical, Electrochemical, and Spectroelectrochemical Properties of Dicyanocyclophanes

A series of dicyanobiphenyl-cyclophanes 1-6 with various pi-backbone conformations and characteristic n-type semiconductor properties is presented. Their synthesis, optical, structural, electrochemical, spectroelectrochemical, and packing properties are investigated. The X-ray crystal structures of all n-type rods allow the systematic correlation of structural features with physical properties. In addition, the results are supported by quantum mechanical calculations based on density functional theory. A two-step reduction process is observed for all n-type rods, in which the first step is reversible. The potential gap between the reduction processes depends linearly on the cos(2) value of the torsion angle phi between the pi-systems. Similarly, optical absorption spectroscopy shows that the vertical excitation energy of the conjugation band correlates with the cos(2) value of the torsion angle phi. These correlations demonstrate that the fixed intramolecular torsion angle phi is the dominant factor determining the extent of electron delocalization in these model compounds, and that the angle phi measured in the solid-state structure is a good proxy for the molecular conformation in solution. Spectroelectrochemical investigations demonstrate that conformational rigidity is maintained even in the radical anion form. In particular, the absorption bands corresponding to the SOMO-LUMO+i transitions are shifted bathochromically, whereas the absorption bands corresponding to the HOMO-SOMO transition are shifted hypsochromically with increasing torsion angle phi.

Position-dependent effects on stability in tricyclo-DNA modified oligonucleotide duplexes

A series of oligodeoxyribonucleotides and oligoribonucleotides containing single and multiple tricyclo(tc)-nucleosides in various arrangements were prepared and the thermal and thermodynamic transition profiles of duplexes with complementary DNA and RNA evaluated. Tc-residues aligned in a non-continuous fashion in an RNA strand significantly decrease affinity to complementary RNA and DNA, mostly as a consequence of a loss of pairing enthalpy DeltaH. Arranging the tc-residues in a continuous fashion rescues T(m) and leads to higher DNA and RNA affinity. Substitution of oligodeoxyribonucleotides in the same way causes much less differences in T(m) when paired to complementary DNA and leads to substantial increases in T(m) when paired to complementary RNA. CD-spectroscopic investigations in combination with molecular dynamics simulations of duplexes with single modifications show that tc-residues in the RNA backbone distinctly influence the conformation of the neighboring nucleotides forcing them into higher energy conformations, while tc-residues in the DNA backbone seem to have negligible influence on the nearest neighbor conformations. These results rationalize the observed affinity differences and are of relevance for the design of tc-DNA containing oligonucleotides for applications in antisense or RNAi therapy.

An engineered disulfide bond reversibly traps the IgE-Fc3-4 in a closed, nonreceptor binding conformation

IgE antibodies interact with the high affinity IgE Fc receptor, FcεRI, and activate inflammatory pathways associated with the allergic response. The IgE-Fc region, comprising the C-terminal domains of the IgE heavy chain, binds FcεRI and can adopt different conformations ranging from a closed form incompatible with receptor binding to an open, receptor-bound state. A number of intermediate states are also observed in different IgE-Fc crystal forms. To further explore this apparent IgE-Fc conformational flexibility and to potentially trap a closed, inactive state, we generated a series of disulfide bond mutants. Here we describe the structure and biochemical properties of an IgE-Fc mutant that is trapped in the closed, non-receptor binding state via an engineered disulfide at residue 335 (Cys-335). Reduction of the disulfide at Cys-335 restores the ability of IgE-Fc to bind to its high affinity receptor, FcεRIα. The structure of the Cys-335 mutant shows that its conformation is within the range of previously observed, closed form IgE-Fc structures and that it retains the hydrophobic pocket found in the hinge region of the closed conformation. Locking the IgE-Fc into the closed state with the Cys-335 mutation does not affect binding of two other IgE-Fc ligands, omalizumab and DARPin E2_79, demonstrating selective blocking of the high affinity receptor binding.