The key role of atom types, reference states, and interaction cutoff radii in the knowledge-based method: new variational approach

We present a variational method to derive knowledge-based potentials. The method is based on an optimization procedure of objective variables: atom types, reference states, and interaction cutoff radii. We suggest and apply new unsymmetrical reference states. The cutoff radii and atom types are optimized to improve docking accuracy of the corresponding potentials. The atom types are varied along an atom type tree, with 6 root and 49 top atom types, and the set of 18 optimal atom types is obtained. We demonstrate strong dependence between the choice of atom types and the docking accuracy of the potentials derived with these atom types. The averaged root-mean square deviations (RMSDs) of the ligand docked positions relative to the experimentally determined positions decrease when the elements C, N, O are split into the optimal types.     

Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparison

A method for the prediction of hydrogen positions in proteins is presented. The method is based on the knowledge of the heavy atom positions obtained, for instance, from X-ray crystallography. It employs an energy minimization limited to the environment of the hydrogen atoms bound to a common heavy atom or to a single water molecule. The method is not restricted to proteins and can be applied without modification to nonpolar hydrogens and to nucleic acids. The method has been applied to the neutron diffraction structures of trypsin ribonuclease A, and bovine pancreatic trypsin inhibitor. A comparison of the constructed and the observed hydrogen positions shows few deviations except in situations in which several energetically similar conformations are possible. Analysis of the potential energy of rotation of Lys amino and Ser, Thr, Tyr hydroxyl groups reveals that the conformations of lowest intrinsic torsion energies are statistically favored in both the crystal and the constructed structures.     

On the use of one-step perturbation to investigate the dependence of NOE-derived atom–atom distance bound violations of peptides upon a variation of force-field parameters

The method of one-step perturbation can be used to predict from a single molecular dynamics simulation the values of observable quantities as functions of variations in the parameters of the Hamiltonian or biomolecular force field used in the simulation. The method is used to predict violations of nuclear overhauser effect (NOE) distance bounds measured in nuclear magnetic resonance (NMR) experiments by atom–atom distances of the NOE atom pairs when varying force-field parameters. Predictions of NOE distance bound violations between different versions of the GROMOS force field for a hexa-β-peptide in solution show that the technique works for rather large force-field parameter changes as well as for very different NOE bound violation patterns. The effect of changing individual force-field parameters on the NOE distance bound violations of the β-peptide and an α-peptide was investigated too. One-step perturbation, which in this case is equivalent to reweighting configurations, constitutes an efficient technique to predict many values of different quantities from a single conformational ensemble for a particular system, which makes it a powerful force-field development technique that easily reduces the number of required separate simulations by an order of magnitude.     

Ribozyme catalysis via orbital steering

Orbital steering is invoked to explain how the three-dimensional structure of a small self-cleaving RNA, the hammerhead ribozyme, both prevents and enhances RNA autocatalysis. Within the conserved catalytic core of the ribozyme, the position of the 2' oxygen atom of the G8 ribose is observed to be aligned almost perfectly with the phosphorus atom and the 5' oxygen atom of the adjacent A9 phosphate group for self-cleavage via an in-line attack mechanism. Despite this apparent near-perfect atomic positioning, no cleavage takes place. The explanation proposed is that a network of hydrogen bonds in the ribozyme core orients or steers the orbitals containing the electron lone pairs of the attacking nucleophile (the 2' oxygen atom) away from the A9 phosphorus atom, eliminating overlap with the vacant phosphorus d-orbitals despite the near-perfect in-line positioning of the oxygen atom, thus preventing catalysis. Because of the near-perfect atomic positioning of the 2' oxygen atom relative to the phosphate group, orbital steering effects in this case are fortuitously uncoupled from conformational, distance and orientation effects, allowing an assessment of the catalytic power due purely to orbital steering. In contrast, a conformational change at the cleavage site is required to bring the 2' oxygen atom and the scissile phosphate group into atomic positions amenable to an in-line attack mechanism. In addition, the conformationally changed structure must then steer the lone-pair orbitals of the correctly positioned 2' oxygen atom toward the scissile phosphorus atom in order for cleavage to take place. We estimate that fulfillment of each of these two required changes may contribute separately an approximately 1000-fold rate enhancement, potentially accounting for a significant fraction of the catalytic power of this ribozyme. Orbital steering therefore appears to be a general phenomenon that may help to explain catalysis in both ribozymes and protein enzymes in a unified manner.     

The differentiation of isomeric biological compounds using collision-induced dissociation of ions generated by fast atom bombardment

Though fast atom bombardment ionization makes possible the ionization and molecular weight determination of polar or thermally labile biological compounds, the resulting mass spectra commonly give few or no fragment ions which would allow detailed structural analysis. In particular, isomeric compounds often give identical spectra. Collision-induced dissociation of ions resulting from fast atom bombardment ionization is shown to be a powerful combination which can differentiate isomeric substances. The technique is applied to isomeric bile acid salts and steroid conjugates and is capable of differentiating structural isomers which have similar fast atom bombardment mass spectra. A range of isomeric cyclic nucleotides is also shown to be amenable to the method. Sensitivity limits are examined and the unequivocal identification of two 3',5'-cyclic nucleotides isolated from living systems is demonstrated.     

Is a nitrogen atom an important pharmacophoric element in sigma ligand binding?

A lingering question in sigma receptor ligand development is whether a nitrogen atom serves as an important pharmacophoric element in binding affinity. To address this question, we have synthesized several phenylalkylpiperidines and phenylalkylpiperazines and demonstrated that removal of the N atom from a typical phenylalkylpiperidine led to little or no binding to sigma receptors. In addition, where two N atoms occur in a compound, such as with phenylalkylpiperazines, the N atom on the longer alkyl chain appears to be more important. Thus, based on this study, the N atom is an important pharmacophoric element in the binding of phenylalkylpiperidines and phenylalkylpiperazines to sigma receptors.     

Membrane structure characterization using variable-period x-ray standing waves.

The variable-period x-ray standing wave (XSW) technique is emerging as a powerful tool for studying membrane structure. However, two significant problems arise when the method is used to characterize membranes of thickness dL < 100 A. First, the surface roughness, sigma(r), of the supporting reflecting mirror convolutes with the intrinsic half-width of the marker atom distribution in the membrane, sigma(in), and contributes to an apparent half-width, sigma, which is measured in the XSW experiment. Here we show how the latter terms are related quantitatively [sigma(in) = (sigma2 - sigma(r)2)(1/2)], such that rough mirrors give rise to larger marker atom distribution widths, sigma, and how the required quantity sigma(in) can be determined in the XSW measurement. Second, when the mean position of the marker atom layer, (z), is close to one or both boundaries of the membrane, its distribution function is truncated at the boundary. In such cases, we show why marker atom distribution should be expressed in terms of its first and second moments. We also demonstrate by numerical simulations of realistic samples how the physical parameters, sigma(r), sigma, (z), and dL, affect x-ray reflectivity and fluorescence yield profiles as an aid in their interpretation.     

Predicting immunoglobulin-like hypervariable loops

A two-stage method is developed to search the conformational space of small protein segments for low energy structures. Central features of the method are efficient procedures for generating small, eight-backbone atom, local moves in Cartesian coordinates and for introducing geometric constraints in adaptable Monte Carlo procedures. This allows natural implementation of an adaptive simulated annealing algorithm, which achieves an effective trade-off between speed and acceptance ratio. The method is applied to the calculation of various immunoglobulin loops. We also develop data base derived rules for identifying constraint conditions, and show that the incorporation of an identified side-chain constraint allows a 1.2 Å all-backbone atom rms deviation prediction of a 9 residue long L1 loop. © 1994 John Wiley & Sons, Inc.     

Question 9: Quantum Self-Assembly and Photoinduced Electron Tunneling in Photosynthetic Systems of Artificial Minimal Living Cells

Natural and artificial living cells and their substructures are self-assembling, due to electron correlation interactions among biological and water molecules, which lead to attractive dispersion forces and hydrogen bonds. Dispersion forces are weak intermolecular forces that arise from the attractive force between quantum multipoles. A hydrogen bond is a special type of quantum attractive interaction that exists between an electronegative atom and a hydrogen atom bonded to another electronegative atom; and this hydrogen atom exist in two quantum states. The best method to simulate these dispersion forces and hydrogen bonds is to perform quantum mechanical non-local density functional potential calculations of artificial minimal living cells consisting of around 1,000 atoms. The cell systems studied are based on peptide nucleic acid and are 3.0–4.2 nm in diameter. The electron tunneling and associated light absorption of the most intense transitions, as calculated by the time dependent density functional theory method, differs from spectroscopic experiments by only 0.2–0.3 nm, which is within the value of experiment errors. This agreement implies that the quantum mechanically self-assembled structures of artificial minimal living cells very closely approximate realistic ones.     

Protein-protein recognition: method for finding complementary surfaces of interacting proteins

A new method is described for searching for complementary surfaces of protein molecules from known coordinates of their non-hydrogen atoms. Each atom is assigned an arbitrary feature which is specific of its interactions with other atoms. The plots representing surfaces are generated. For each pair of surfaces a number of coincidences is calculated which increases as the number of contacts is increased between these atom pairs whose contribution to the energy of interaction is considered to be essential. The results obtained for the well-known autoassociation of insulin show the applicability of the method for the prediction of possible reaction modes between macromolecules.     

A new electronic-topological investigation of the relationship between chemical structure and ambergris odour.

An electronic-topological approach has been used to define an active ambergris fragment (AAF) which correctly describes the presence (or absence) of the ambergris odour of all 181 compounds investigated. The AAF consists of one oxygen atom and three carbon atoms (alpha, beta, gamma) which are separated by certain key distances and which possess certain atomic charges. The C(alpha) atom must bear at least one hydrogen atom (H(alpha)) which is located at a certain distance from one of the unshared electronic pairs of the oxygen atom.     

A Molecular Electrostatic Potential Mapping Study of Some Fluoroquinolone Anti-Bacterial Agents

Molecular electrostatic potential (MEP) maps of some fluoroquinolones having varying degrees of activity against the bacterium Staphylococcus Aureus have been studied using the optimized hybridization displacement charges (HDC) combined with Löwdin charges obtained by the AM1 method. The roles of different substitutions at the N₁-position in the parent quinolone ring have been studied. The conformation of the carboxylic group attached to the quinolone ring was shown to be such that there is an intramolecular hydrogen bonding between the hydrogen atom of this group and the oxygen atom of the carbonyl group of the quinolone moiety. The carbonyl oxygen atom of the quinolone moiety, hydroxyl oxygen atom of the carboxylic group and the terminal nitrogen atom of the piperazin ring attached to the quinolone ring appear to be involved in the action of the drugs through electrostatic interactions while the N₁-alkyl substituents seem to be involved in the same through hydrophobic interactions.     

An algorithm to filter out packing arrangements based on steric clashes

This document outlines the use of an algorithm to filter out impossible crystal-packing arrangements based on steric considerations. Within an exhaustive grid search frame, the space sample is reduced by analysis of spherical areas where atom pairs from different rigid units might clash.This technique finds areas in the state space where the global energy minimum might lie. The minimum can then be found by the usual methods of molecular modeling restricted to these particular areas.Only a tiny fraction of atom pair distances need to be tested; usually a single quantity on average per one state of model space! For example, a crystal of three rigid molecules, each containing 12 atoms, has 3×12×12=432 atom pairs just in one unit cell but our method needs to test on average 1 to 4 atom pairs per state.Using modern computers, about 10^12–15 models can be tested within several hours or days. For example, a crystal model with six rotational degrees of freedom (two rigid molecules in the unit cell) each with step 3° can be tested in a few hours on a 1-GHz x86 processor-based machine.The method presented here has been implemented in the SUPRAMOL program.     

Design,Synthesis and Evaluation of Fe-S Targeted Adenosine 5′-Phosphosulfate Reductase Inhibitors

Adenosine 5′-phosphosulfate reductase (APR) is an iron-sulfur enzyme that is vital for survival of Mycobacterium tuberculosis during dormancy and is an attractive target for the treatment of latent tuberculosis (TB) infection. The 4Fe-4S cluster is coordinated to APR by sulfur atoms of four cysteine residues, is proximal to substrate, adenosine 5′-phopsphosulfate (APS), and is essential for catalytic activity. Herein, we present an approach for the development of a new class of APR inhibitors. As an initial step, we have employed an improved solid-phase chemistry method to prepare a series of N⁶-substituted adenosine analogues and their 5′-phosphates as well as adenosine 5′-phosphate diesters bearing different Fe and S binding groups, such as thiols or carboxylic and hydroxamic acid moieties. Evaluation of the resulting compounds indicates a clearly defined spacing requirement between the Fe-S targeting group and adenosine scaffold and that smaller Fe-S targeting groups are better tolerated. Molecular docking analysis suggests that the S atom of the most potent inhibitor may establish a favorable interaction with an S atom in the cluster. In summary, this study showcases an improved solid-phase method that expedites the preparation of adenosine and related 5′-phosphate derivatives and presents a unique Fe-S targeting strategy for the development of APR inhibitors.     

Complexation of uranium by cells and S-layer sheets of Bacillus sphaericus JG-A12

Bacillus sphaericus JG-A12 is a natural isolate recovered from a uranium mining waste pile near the town of Johanngeorgenstadt in Saxony, Germany. The cells of this strain are enveloped by a highly ordered crystalline proteinaceous surface layer (S-layer) possessing an ability to bind uranium and other heavy metals. Purified and recrystallized S-layer proteins were shown to be phosphorylated by phosphoprotein-specific staining, inductive coupled plasma mass spectrometry analysis, and a colorimetric method. We used extended X-ray absorption fine-structure (EXAFS) spectroscopy to determine the structural parameters of the uranium complexes formed by purified and recrystallized S-layer sheets of B. sphaericus JG-A12. In addition, we investigated the complexation of uranium by the vegetative bacterial cells. The EXAFS analysis demonstrated that in all samples studied, the U(VI) is coordinated to carboxyl groups in a bidentate fashion with an average distance between the U atom and the C atom of 2.88 +/- 0.02 A and to phosphate groups in a monodentate fashion with an average distance between the U atom and the P atom of 3.62 +/- 0.02 A. Transmission electron microscopy showed that the uranium accumulated by the cells of this strain is located in dense deposits at the cell surface.     

Thiaisoleucine and protein synthesis.

Thiaisoleucine is an isoleucine analogue having the gamma-methylene group of the valerianic carbon chain substituted by a sulphur atom. It has been demonstrated that thiaisoleucine is activated and transferred to tRNAIle by rat liver aminoacyl-tRNA synthetase and inhibits isoleucine incorporation into polypeptides in protein synthesizing systems from rat liver or rabbit reticulocytes, whereas it does not affect either leucine incorporation or ribosome run-off or polypeptide chain elongation rate. All tests were performed in comparison with O-methyl-threonine, an isoleucine analogue with the gamma-methylene group substituted by an oxygen atom. In all the reactions studied, both thiaisoleucine and O-methyl-threonine act as competitive inhibitors of isoleucine. With respect to O-methyl-threonine, thiaisoleucine shows higher activity as an isoleucine inhibitor.     

Ab-initio refinement of an orbital-centred force field for biomolecules. test cases including peptides,a sulphonamide and modelling of DNA helices

The application of the Oribital-centred Force Field (OFF) to the calculation of the conformational behaviour of biomolecules is further considered. This method includes classical type potential functions between atom centres, but also specifies the location of non-core orbital lobes and their strength of electrostatic interaction with other lobes or atom centres. It is shown that the OFF method is not particularly advantageous for peptide systems, since although non-core orbital lobes are present with typical electrostatic strength, they do not dominate the potential surface. For other systems, however, of which sulphonamide is used as an example, the interactions between non-core orbitals are the major contribution, as was earlier shown to be the case for the nucleic acid backbone. An interesting test case is the dioxy system OXO, for which a priori the OFF approach may be inadequate because of an anomeric effect discussed by other authors interested in saccharide systems. The possibility that an intrinsic rotational potential may be required is considered but remains by no means obvious (the need for such potentials is normally absent when the OFF approach is used). Finally, the applications of the OFF approach are illustrated and tested by reference to the conformational analysis of 5′-AMP and to the computer-aided model building of DNA helices.     

Use of flavin photochemistry to probe intraprotein and interprotein electron transfer mechanisms

Photoexcitation of flavin analogs generates the lowest triplet state (via intersystem crossing from the first excited singlet state) in the nanosecond time domain and with high quantum efficiency. The triplet, being a strong oxidant, can abstract a hydrogen atom (or an electron) from a reduced donor in a diffusion-controlled reaction. If the donor is a redox protein, the oxidation process can be used to initiate an electron transfer sequence involving either intramolecular or intermolecular reactions. If the donor is an organic compound such as EDTA, the neutral flavin semiquinone will be produced by H atom abstraction; this is a strong reductant and can subsequently transfer a hydrogen atom (or an electron) to an oxidized redox protein, thereby again initiating a sequence of intramolecular or intermolecular processes. If flavin photoexcitation is accomplished using a pulsed laser light source, the initiation of these protein electron transfer reactions can be made to occur in the nanosecond to microsecond time domain, and the sequence of events can be followed by time-resolved spectrophotometry to obtain rate constants and thus mechanistic information. The present paper describes this technology, and selected examples of its use in the investigation of redox protein mechanisms are given.     

Incorporation of fluorotryptophan into triostin antibiotics by Streptomyces triostinicus

The quinoxaline chromophores of the antibiotics produced by Streptomyces triostinicus are derived from tryptophan. Protoplasts of this organism made novel products when they were incubated with DL-5-fluorotryptophan or DL-6-fluorotryptophan. When added to batch cultures of the organism, DL-5-fluorotryptophan, at concentrations as low as 10 microM, inhibited both mycelial growth and triostin production, but gave rise to novel products. These have been characterized, using fast atom bombardment mass spectrometry, as novel triostins in which one or both of the quinoxaline rings contain an atom of fluorine. The chromatographic properties of the triostins arising from the incorporation of DL-5-fluorotryptophan are very similar to those of triostins containing chlorine or bromine at position 6 of the quinoxaline ring; they are clearly different from those having a chlorine atom at position 7. Accordingly, it is suggested that the carbon atom at position 5 of the indole ring of tryptophan ends up at position 6 of the quinoxaline ring system in triostins A and C.     

Spontaneous resolution of binary copper(II) complexes with racemic dipeptides: crystal structures of glycyl-L-alpha-amino-n-butyrato copper(II) monohydrate, glycyl-D-valinato copper(II) hemihydrate, and glycyl-L-valinato copper(II) hemihydrate

Copper(II) complexes with glycyl-DL-alpha-amino-n-butyric acid (H2gly-DL-but), glycyl-DL-valine (H2gly-DL-val), glycyl-DL-norleucine (H2gly-DL-norleu), glycyl-DL-threonine (H2gly-DL-thr), glycyl-DL-serine (H2gly-DL-ser), glycyl-DL-phenylalanine (H2gly-DL-phe), and glycyl-L-valine (H2gly-L-val), have been prepared and characterized by IR, powder diffuse reflection, CD and ORD spectra, and magnetic susceptibility measurements, and by single-crystal X-ray diffraction. The crystal structures of the copper complex with H2gly-DL-but, the copper complex with H2gly-DL-val, and [Cu(gly-L-val)]n.0.5nH2O have been determined by a single-crystal X-ray diffraction method. As for the structure of the copper complex with H2gly-DL-but, the configuration around the asymmetric carbon atom is similar to that of [Cu(gly-L-val)]n.0.5nH2O. Therefore it is concluded that the copper complex with H2gly-DL-but is [Cu(gly-L-but)]n.nH2O. On the contrary, as for the structure of the copper complex with H2gly-DL-val, the configuration around the asymmetric carbon atom is different from that of [Cu(gly-L-val)]n.0.5nH2O. Therefore it is concluded that the copper complex with H2gly-dl-val is [Cu(gly-D-val)]n.0.5nH2O. So during the crystallization of the copper(II) complexes with H2gly-DL-but and H2gly-DL-val, spontaneous resolution has been observed; the four complexes have separated as [Cu(gly-D-but)]n.nH2O, [Cu(gly-L-but)]n.nH2O, [Cu(gly-D-val)]n.0.5nH2O, and [Cu(gly-L-val)]n.0.5nH2O, respectively. [Cu(gly-L-but)]n.nH2O is orthorhombic with the space group P2(1)2(1)2(1). [Cu(gly-D-val)]n.0.5nH2O and [Cu(gly-L-val)]n.0.5nH2O are monoclinic with the space group C2. In these complexes, the copper atom is in a square-pyramidal geometry, ligated by a peptide nitrogen atom, an amino nitrogen atom, a carboxyl oxygen atom, and a carboxyl oxygen atom and a peptide oxygen atom from neighboring molecules. So these complexes consist of a two-dimensional polymer chain bridged by a carboxyl oxygen atom and a peptide oxygen atom from neighboring molecules. The axial oxygen atom is located above the basal plane and the side chain of an amino acid is located below it. These polymer chains consist of only one or the other type of optical isomers; no racemic dipeptides are found. Therefore, spontaneous resolution has been observed in the crystallization of copper(II) complexes with H2gly-DL-but and H2gly-DL-val. The crystal structure of [Cu(gly-D-val)]n.0.5nH2O agrees almost completely with that of [Cu(gly-L-val)]n.0.5nH2O, except for the configuration around the asymmetric carbon atom.