Structure cristalline d'une thiazolorifamycine: rifamycine P

Crystal Structure of a Thiazolorifamycine: Rifamycine P The crystal structure of a natural thiazolorifamycin, rifamycin P , as an adduct with EtOH and H₂O, has been determined by X-ray crystallography by a combination of direct methods and Patterson techniques. Block-diagonal least-squares refinement using 5996 independent reflections and 1475 parameters led to an R factor of 0.092 (two independent molecules of rifamycin P). The middle part of the ansa chain, essential for activity against the bacterial DNA-dependent RNA polymerase (DDRP) resembles that of active rifamycins. The four OH groups lie on the same side of the molecule and are almost in the same plane. Interatomic distances between the O-atoms agree well with a spatial model derived from X-ray studies, performed on semisynthetic active rifamycins.     

The conformation of rifamycin S in solution by 1H NMR spectroscopy

The ¹H NMR spectra of rifamycin S in different solvents and at different temperatures strongly suggest the existence of a dominant conformer. The nine dihedral angles of the ansa chain from C-28 to C-19, considered conformationally flexible, were obtained from the vicinal interproton coupling constants by the Karplus equation and the proper alternative for each of them compares well with the corresponding values given by X-ray analysis in solid state.All the possible conformations derivable from the intrinsic alternatives of the NMR method were calculated for the ansa chain between C-28 and C-19 and, by applying geometrical considerations, such as closure of the polygonal path of the ansa and the steric incompatibility between the various atoms of the ansa and of the chromophore, only two of them appeared real. Thus, NMR spectroscopy can be used for studying the conformation of the ansa chain of rifamycins in solution.     

A fast and efficient method to generate biologically relevant conformations

Summary Mutual binding between a ligand of low molecular weight and its macromolecular receptor demands structural complementarity of both species at the recognition site. To predict binding properties of new molecules before synthesis, information about possible conformations of drug molecules at the active site is required, especially if the 3D structure of the receptor is not known. The statistical analysis of small-molecule crystal data allows one to elucidate conformational preferences of molecular fragments and accordingly to compile libraries of putative ligand conformations. A comparison of geometries adopted by corresponding fragments in ligands bound to proteins shows similar distributions in conformation space. We have developed an automatic procedure that generates different conformers of a given ligand. The entire molecule is decomposed into its individual ring and open-chain torsional fragments, each used in a variety of favorable conformations. The latter ones are produced according to the library information about conformational preferences. During this building process, an extensive energy ranking is applied. Conformers ranked as energetically favorable are subjected to an optimization in torsion angle space. During minimization, unfavorable van der Waals interactions are removed while keeping the open-chain torsion angles as close as possible to the experimentally most frequently observed values. In order to assess how well the generated conformers map conformation space, a comparison with experimental data has been performed. This comparison gives some confidence in the efficiency and completeness of this approach. For some ligands that had been structurally characterized by protein crystallography, the program was used to generate sets of some 10 to 100 conformers. Among these, geometries are found that fall convincingly close to the conformations actually adopted by these ligands at the binding site.     

Indor and 1H-NMR. Studies of the Ansa Conformation of Hexahydrorifamycin S

The study of 16, 17, 18, 19, 28, 29-hexahydrorifamycin S by ¹H-NMR. spectroscopy using INDOR at 100 MHz and homodecoupling at 270 MHz has shown that the conformation of the ansa fragment C(20)-C(27) and the position of this fragment relative to the naphthoquinone chromophore are practically the same as those of rifamycin S. This finding rules out the hypothesis that the lower activity of the hydrogenated rifamycins could be due to a different conformation.     

Influence of an Unnatural Amino Acid Side Chain on the Conformational Dynamics of Peptides

In this work, a non‐natural amino acid, H‐propargylglycine‐OH (Pra), is chosen to examine the side‐chain effect on the backbone conformation of small peptides. The conformations of two synthesized Pra‐containing tripeptides, Ac‐Pra‐Pra‐NH₂ (PPTP) and Ac‐Pra‐Ala‐NH₂ (PATP), are examined by infrared (IR) spectroscopy in combination with molecular dynamics (MD) simulations and quantum chemical computations. By analyzing the joint distributions of backbone torsional angles, several significant conformations can be identified for the two tripeptides solvated in D₂O. At room temperature, 44 % of PPTP exists in the α‐α conformation and 33 % of PATP exists in the α‐polyproline‐II conformation. Larger structural inhomogeneity is seen in both cases by MD simulations at elevated temperatures. Thus even a small side chain, such as the propargyl group can significantly alter the peptide backbone conformations. The results suggest that there is no overwhelming conformational propensity of the Pra residue in short peptides. IR spectra simulated in the amide‐I region using two different methods, reasonably reproduce the experimental IR spectra and their temperature dependence.     

Fine grained sampling of residue characteristics using molecular dynamics simulation

In a fine-grained computational analysis of protein structure, we investigated the relationships between a residue's backbone conformations and its side-chain packing as well as conformations. To produce continuous distributions in high resolution, we ran molecular dynamics simulations over a set of protein folds (dynameome). In effect, the dynameome dataset samples not only the states well represented in the PDB but also the known states that are not well represented in the structural database. In our analysis, we characterized the mutual influence among the backbone ?,ψ angles with the first side-chain torsion angles (χ₁) and the volumes occupied by the side-chains. The dependencies of these relationships on side-chain environment and amino acids are further explored. We found that residue volumes exhibit dependency on backbone 2° structure conformation: side-chains pack more densely in extended β-sheet than in α-helical structures. As expected, residue volumes on the protein surface were larger than those in the interior. The first side-chain torsion angles are found to be dependent on the backbone conformations in agreement with previous studies, but the dynameome dataset provides higher resolution of rotamer preferences based on the backbone conformation. All three gauche^?, gauche⁺, and trans rotamers show different patterns of ?,ψ dependency, and variations in χ₁ value are skewed from their canonical values to relieve the steric strains. By demonstrating the utility of dynameomic modeling on the native state ensemble, this study reveals details of the interplay among backbone conformations, residue volumes and side-chain conformations.     

Chi(1) rotamer populations and angles of mobile surface side chains are accurately predicted by a torsion angle database potential of mean force

The equilibrium angles and distributions of chi(1) rotamers for mobile surface side chains of the small, 63-residue, B1 domain of protein L have been calculated from the static crystal structure by rigid body/torsion angle simulated annealing using a torsion angle database potential of mean force and compared to those deduced by Monte Carlo analysis of side chain residual dipolar couplings measured in solution. Good agreement between theory and experiment is observed, indicating that for side chains undergoing rotamer averaging that is fast on the chemical shift time scale, the equilibrium angles and distribution of chi(1) rotamers are largely determined by the backbone phi/psi torsion angles.     

A kinematic view of loop closure

We consider the problem of loop closure, i.e., of finding the ensemble of possible backbone structures of a chain segment of a protein molecule that is geometrically consistent with preceding and following parts of the chain whose structures are given. We reduce this problem of determining the loop conformations of six torsions to finding the real roots of a 16th degree polynomial in one variable, based on the robotics literature on the kinematics of the equivalent rotator linkage in the most general case of oblique rotators. We provide a simple intuitive view and derivation of the polynomial for the case in which each of the three pair of torsional axes has a common point. Our method generalizes previous work on analytical loop closure in that the torsion angles need not be consecutive, and any rigid intervening segments are allowed between the free torsions. Our approach also allows for a small degree of flexibility in the bond angles and the peptide torsion angles; this substantially enlarges the space of solvable configurations as is demonstrated by an application of the method to the modeling of cyclic pentapeptides. We give further applications to two important problems. First, we show that this analytical loop closure algorithm can be efficiently combined with an existing loop-construction algorithm to sample loops longer than three residues. Second, we show that Monte Carlo minimization is made severalfold more efficient by employing the local moves generated by the loop closure algorithm, when applied to the global minimization of an eight-residue loop. Our loop closure algorithm is freely available at http://dillgroup. ucsf.edu/loop_closure/.     

Vibrational and electronic properties of polyaza chain molecules and their metal complexes

The polyaza chain molecules exhibit a quasi planar backbone with all-trans geometry. The chelation of several metallic ions such as copper (II) and zinc (II) constrains different conformations of the chain molecules. The vibrational and electronic properties are typical of the conformation of the polyaza backbone as well as the spin spin exchange between the metallic ions through the azine bonds.     

Protein structure prediction using a combination of sequence homology and global energy minimization: II. Energy functions

A protein energy surface is constructed. Validation is through applications of global energy minimization to surface loops of protein crystal structures. For 9 of 10 predictions, the native backbone conformation is identified correctly. Electrostatic energy is modeled as a pairwise sum of interactions between anisotropic atomic charge densities. Model repulsion energy has a softness similar to that seen in ab initio data. Intrinsic torsional energy is modeled as a sum over pairs of adjacent torsion angles of 2-dimensional Fourier series. Hydrophobic energy is that of a hydration shell model. The remainder of hydration free energy is obtained as the energetic effect of a continuous dielectric medium. Parameters are adjusted to reproduce the following data: a complete set of ab initio energy surfaces, meaning one for each pair of adjacent torsion angles of each blocked amino acid; experimental crystal structures and sublimation energies for nine model compounds; ab initio energies over 1014 conformations of 15 small-molecule dimers; and experimental hydration free energies for 48 model compounds. All ab initio data is at the Hartree–Fock/6–31G* level. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 548–573, 1998     

A Vibrational Raman Optical Activity Study of 1,1′‐Binaphthyl Derivatives

The Raman polarized and vibrational Raman optical activity (VROA) backward spectra are simulated for a series of 2,2′‐substituted 1,1′‐binaphthyl compounds presenting a variety of torsion angles between the two naphthalene rings. The substitution prevents free rotation along this torsion angle and the chirality of these compounds is thus called atropisomerism. However, the rotation is not completely frozen so that two different conformations, namely cisoid and transoid, are found and their Raman and VROA signatures are studied. As expected, the Raman spectra are not very sensitive whereas the VROA spectra present more complex patterns, which evolve as a function of the torsion angle between the two naphthalene groups. In particular, our analysis shows that some modes can be used as a probe for the determination of the torsion angle of these molecules in solution. The contributions of both invariants to the VROA backward intensity are also assessed.     

Systematic drug receptor mapping: A new approach to the analysis of conformational energy calculations of flexible molecules with application to dopaminergic and adrenergic agonists

The computer program PRODIS is used to find low energy conformations of flexible molecules by searching the potential energy surface(s) of one or more torsion angles via rigid rotation. The n-dimensional grid of energy versus torsion angles is then converted to a Boltzman probability distribution, with the probability being represented not as a function of torsion angle, but rather a distance between two atoms. These atoms are chosen by comparison with a known, active analogue in which certain atoms have previously been determined as requirements for drug activity. PRODIS produces a list of low energy conformations, their corresponding interatomic distances and the Boltzman probability for each distance ±0.125, as well as the total probability for each conformation. The user also specifies a target interatomic distance and range (usually derived from a more rigid analogue) for which PRODIS lists all conformations and their Boltzman probability that meet this distance.     

Torsional flexing: Conformational searching of cyclic molecules in biased internal coordinate space

A new stochastic (Monte Carlo) procedure, termed torsional flexing, has been devised for searching the conformational space of cyclic molecules. Torsional flexing causes a local, torsion angle-biased, distortion of a ring bond in a cyclic molecule. Because torsional flexing does not cause large atomic movements, even when it is applied to several bonds simultaneously, subsequent energy minimization generally proceeds rapidly. Nevertheless, the torsional flexing method is prone to generate structures that cross energy barriers so that the structure resulting after energy minimization is frequently a different conformer of the cyclic molecule. Conformational searches on cycloheptadecane, oxobrefeldin A, cyclopenta-L -alanine, and rifamycin SV based upon torsional flexing indicated that torsional flexing is among the best methods yet devised for searching the conformational space of flexible cyclic molecules. © 1993 John Wiley & Sons, Inc.     

Molecular determinants for drug-receptor interactions Part 14. X-ray molecular structure and theoretical conformational (MNDO, MM2) studies of the narcotic dualist buprenorphine hydrochloride

The X-ray molecular structure of buprenorphine hydrochloride [(C₂₉H₄₂O₄N)⁺Cl⁻] was determined. The molecule forms colourless tetragonal crystals; the space group is P4₁2₁2₁ with eight molecules per unit cell of dimensions a = 11.507(2) Å and c = 41.952(9) Å. The structure was determined by direct methods and refined by the full-matrix least-squares procedure to R = 0.056 for 2251 observed reflections. The protonated buprenorphine (BPNH⁺) in the solid state adopts the T-shape form typical of benzomorphan compounds; the substituent group at the N atom is in the equatorial configuration and the cyclopropylmethyl side-chain is almost perpendicular to the piperidine ring. Theoretical studies of the conformations and the rotational energetics of BPNH⁺ were performed using both the semiquantitative MNDO and the force field (MM2) methods. The latter method predicted four low energy conformations about the N-C(21) and C(21)-C(22) bonds. Two of these were more significantly populated (59% and 30%). They may easily interconvert through only a single torsional process. The less-populated (30%) conformation was consistent with that adopted in the solid state. The MNDO method predicted different low energy conformations of the N-cyclopropylmethyl moiety. Neither method reproduced the favoured torsional angles adopted by the structurally analogous naltrexone molecule which has the opposite pharmacological profile.     

Calculation of structural behavior of indirect NMR spin-spin couplings in the backbone of nucleic acids

Calculated indirect NMR spin-spin coupling constants (J-couplings) between (31)P, (13)C, and (1)H nuclei were related to the backbone torsion angles of nucleic acids (NAs), and it was shown that J-couplings can facilitate accurate and reliable structural interpretation of NMR measurements and help to discriminate between their distinct conformational classes. A proposed stepwise procedure suggests assignment of the J-couplings to torsion angles from the sugar part to the phosphodiester link. Some J-couplings show multidimensional dependence on torsion angles, the most prominent of which is the effect of the sugar pucker. J-couplings were calculated in 16 distinct nucleic acid conformations, two principal double-helical DNAs, B- and A-, the main RNA form, A-RNA, as well as in 13 other RNA conformations. High-level quantum mechanics calculations used a baseless dinucleoside phosphate as a molecular model, and the effect of solvent was included. The predicted J-couplings correlate reliably with available experimental data from the literature.     

Determination of peptide backbone torsion angles using double-quantum dipolar recoupling solid-state NMR spectroscopy

Several approaches for utilizing dipolar recoupling solid-state NMR (ssNMR) techniques to determine local structure at high resolution in peptides and proteins have been developed. However, many of these techniques measure only one torsion angle or are accurate for only certain classes of secondary structure. Additionally, the efficiency with which these dipolar recoupling experiments suppress the deleterious effects of chemical shift anisotropy (CSA) at high magnetic field strengths varies. Dipolar recoupling with a windowless sequence (DRAWS) has proven to be an effective pulse sequence for exciting double-quantum (DQ) coherences between adjacent carbonyl carbons along the peptide backbone. By allowing this DQ coherence to evolve, it is possible to measure the relative orientations of the CSA tensors and subsequently use this information to determine the Ramachandran torsion angles phi and psi. Here, we explore the accuracies of the assumptions made in interpreting DQ-DRAWS data and demonstrate their fidelity in measuring torsion angles corresponding to a variety of secondary structures irrespective of hydrogen-bonding patterns. It is shown how a simple choice of isotopic labels and experimental conditions allows accurate measurement of backbone secondary structures without any prior knowledge. This approach is considerably more sensitive for determining structure in helices and has comparable accuracy for beta-sheet and extended conformations relative to other methods. We also illustrate the ability of DQ-DRAWS to distinguish between structures in heterogeneous samples.     

The effects of biasing torsional mutations in a conformational GA

This paper describes the effects of incorporating torsional bias into a conformational Genetic Algorithm (GA) such as that found in the GASP program. Several major conclusions can be drawn. Biasing torsional angles toward values associated with local energy minima increases the rate of convergence of the fitness function (consisting of energy, steric, and pharmacophoric compatibility terms) for a set of molecules, but a definite tradeoff exists between total model energy and the steric and pharmacophoric compatibility terms in the fitness score. Biasing torsions in favor of sets of angles drawn from low-energy conformations does not guarantee low total energy, but biased torsional sampling does generally produce less strained models than does the uniform torsional sampling in classical GASP. Overall, torsionally biased sampling produces good models comprised of energetically favorable ligand conformations.     

Quantum theory of the structure and bonding in proteins: Part 4. The question of the quality of the basis set and of the stability of the extended forms

The conformational energy map of the molecule NH₂CH₂CHO is constructed using both the long basis set and the short basis set of the GAUSSIAN 70 package. It is shown that the two maps are qualitatively the same but that there are quantitative differences between them. Analysis of these and other results suggests that the stability of the extended planar structure is due in the main to the minimization of torsional barriers about the two single bonds of the backbone of the protein chain.There is broad general agreement between the semi-empirical method and its results as used by Scheraga and his colleagues on the one hand and our results from the ab initio computations on the other. Both methods suggest that the details of the conformations of these small peptides, treated as isolated molecules, are determined by hydrogen bonds and torsional barriers.     

Effect of ring torsion on intramolecular vibrational redistribution dynamics of 1,1'-binaphthyl and 2,2'-binaphthyl

The role of ring torsion in the enhancement of intramolecular vibrational energy redistribution (IVR) in aromatic molecules was investigated by conducting excitation and dispersed fluorescence spectroscopy of 1,1'-binaphthyl (1,1'-BN) and 2,2'-BN. The dispersed fluorescence spectra of 1,1'-BN in the origin region of S(1)-S(0) were well resolved, which presented 25-27 cm(-1) gaps of torsional mode in the ground state. The overall profile of the dispersed spectra of 1,1'-BN is similar to that of naphthalene. In contrast, the spectra of 2,2'-BN were not resolved due to the multitude of the active torsional modes. In both cases, dissipative IVR was observed to take place with a relatively small excess vibrational energy: 237.5 cm(-1) for 1,1'-BN and 658 cm(-1) for 2,2'-BN, which clearly shows that ring torsion efficiently enhances the IVR rate. Ab initio and density functional theory calculations with medium-sized basis sets showed that the torsional potential of 1,1'-BN has a very flat minimum over the range of torsional angles from ca. 60° to 120°, whereas that of 2,2'-BN showed two well-defined potential minima at ca. 40° and 140°, in resemblance to the case of biphenyl. In this work, we propose that aromatic molecules be classified into "strong" and "weak" torsional hindrance cases: molecules with strong hindrance case show shorter torsional progressions and more effective IVR dynamics than do those with weak hindrance.     

Head-to-backbone cyclization of peptides on solid support by nucleophilic aromatic substitution

A new versatile synthetic route is presented for the cyclization of tripeptides on solid support using nucleophilic aromatic substitution in the cyclization step. Identification of all conformers within a limit of 3 kcal/mol from the identified global minimum conformations by Monte Carlo conformational searching reveals that five out of six synthesized compounds have well-defined peptide backbone conformational properties. This was determined by clustering the identified conformers against a filter of seven to nine torsion angles in the peptide backbone. Thus, the results meet our goal to find synthetic routes to peptides that are conformationally sufficiently locked to serve as convenient leads for further development of pharmacophoric models. The strategy is based on Fmoc-peptide chemistry on a N-aminoethyl-substituted glycine bound to the commercially available Rink amide PS-resin. After deprotection of the N-terminus of the tripeptide, it is acylated with a fluoronitrobenzoic acid. Subsequently, a Boc group on the N-bound aminoethyl substituent is selectively deprotected allowing cyclization from the head (N-terminus) to the backbone substituent, thereby leading to the desired cyclized tripeptides. A number of representative examples of peptides cyclized by this method have been synthesized and characterized by NMR. Protecting groups that allow the incorporation of side chain functionalized amino acids have been found. Thus, the route provides access to generic libraries of conformationally restricted peptide sequences expressing a range of proteinogenic pharmacophores.