New formulation for derivatives of torsion angles and improper torsion angles in molecular mechanics: Elimination of singularities

A new set of formulae is developed for the derivatives of torsion angle energy terms and is introduced into the program CHARMM. These formulae, which are based on derivatives of the torsion angle itself, avoid the singularities introduced by use of the derivatives of the torsion angle cosine. The potential energy can include any differentiable function of the torsion angle and there is no need for a special treatment for cases where planar conformations are not extrema. The resulting code is simpler than the original version and yields correct derivatives in all practical situations. Because the minimum of the torsion energy can be at any angle, the functionality of the existing energy routines is generalized. © 1996 by John Wiley & Sons, Inc.     

T-Analyst: a program for efficient analysis of protein conformational changes by torsion angles

T-Analyst is a user-friendly computer program for analyzing trajectories from molecular modeling. Instead of using Cartesian coordinates for protein conformational analysis, T-Analyst is based on internal bond-angle-torsion coordinates in which internal torsion angle movements, such as side-chain rotations, can be easily detected. The program computes entropy and automatically detects and corrects angle periodicity to produce accurate rotameric states of dihedrals. It also clusters multiple conformations and detects dihedral rotations that contribute hinge-like motions. Correlated motions between selected dihedrals can also be observed from the correlation map. T-Analyst focuses on showing changes in protein flexibility between different states and selecting representative protein conformations for molecular docking studies. The program is provided with instructions and full source code in Perl.     

Transition state geometry in radical abstraction reactions: comparison of interatomic distances in the intersecting parabolas and Morse curves models with quantum-chemical calculations

Interatomic distances in the transition state were estimated for the reactions of radical abstraction: H^· + H₂, H^· + HCl, H^· + CH₄, N^·H₂ + NH₃, HO^· + H₂O, HO₂ ^· + HOOH, and C^·H₃ + SiH₄. The calculation was performed by the quantum-chemical density functional method or coupled clusters method (QCH), as well as by the methods of intersecting parabolas (IPM) and Morse curves (IMM), using experimental data (activation energies and reaction enthalpies). The results of the latter two methods are close to the quantum-chemical calculation and differ only by the increment a: r(IPM or IMM) = a + r(QCH), where a = –4.5·10^–12 m for IPM and a = +1.9·10^–12 m for IMM.     

PLUMS: a program for the rapid optimization of focused libraries

PLUMS is a new method to perform rational monomer selection for combinatorial chemistry libraries. The algorithm has been developed to optimize focused libraries with specific two-dimensional and/or three-dimensional properties. A preliminary step is the identification of those molecules in the initial virtual library which satisfy the imposed property constraints; we define these molecules as the virtual hits. From the virtual hits, PLUMS generates a starting library, which is the true combinatorial library that includes all the virtual hits. Monomers are then removed in an iterative fashion, thus reducing the size of the library. At each iteration, the worst monomer is removed. Each sublibrary is selected using a global scoring function, which balances effectiveness and efficiency. The iterative process continues until one is left with a library that consists entirely of virtual hits. The optimal library, which is the best compromise between effectiveness and efficiency, can then be selected according to the score. During the iterative process, equivalent solutions may well occur and are taken into account by the algorithm, according to a user-defined parameter. The number of monomers for each substitution site and the size of the library are parameters that can be either optimized or used to constrain the selection. The results obtained on two test libraries are presented. PLUMS was compared with genetic algorithms (GA) and monomer frequency analysis (MFA), which are widely used for monomer selection. For the two test libraries, PLUMS and GA gave equivalent results. MFA is the fastest method, but it can give misleading solutions. Possible advantages and disadvantages of the different methods are discussed.     

Efficient evaluation of sampling quality of molecular dynamics simulations by clustering of dihedral torsion angles and Sammon mapping

A direct conformational clustering and mapping approach for peptide conformations based on backbone dihedral angles has been developed and applied to compare conformational sampling of Met-enkephalin using two molecular dynamics (MD) methods. Efficient clustering in dihedrals has been achieved by evaluating all combinations resulting from independent clustering of each dihedral angle distribution, thus resolving all conformational substates. In contrast, Cartesian clustering was unable to accurately distinguish between all substates. Projection of clusters on dihedral principal component (PCA) subspaces did not result in efficient separation of highly populated clusters. However, representation in a nonlinear metric by Sammon mapping was able to separate well the 48 highest populated clusters in just two dimensions. In addition, this approach also allowed us to visualize the transition frequencies between clusters efficiently. Significantly, higher transition frequencies between more distinct conformational substates were found for a recently developed biasing-potential replica exchange MD simulation method allowing faster sampling of possible substates compared to conventional MD simulations. Although the number of theoretically possible clusters grows exponentially with peptide length, in practice, the number of clusters is only limited by the sampling size (typically much smaller), and therefore the method is well suited also for large systems. The approach could be useful to rapidly and accurately evaluate conformational sampling during MD simulations, to compare different sampling strategies and eventually to detect kinetic bottlenecks in folding pathways.     

A rapid method for molecular shape comparison of medium-sized molecules. Conformational calculations on sandalwood odor,IV

Summary A fast method of a surface comparison of two or more molecules to be matched is presented. The Van der Waals surfaces of molecules are described by points calculated as the intersection of grid lines with the molecular surface. The mean surface of various molecules with the same biological activity can be constructed. It is used for further comparisons with similar molecules lacking this activity. Deviations of any molecular surface from the mean surface can be mapped onto the surface. The method was tested on a distinct group of sandalwood odor molecules and it was shown that such matching and comparison procedures are useful in the investigation of odor structure-activity relationships proposed as CAFD (computer aided fragrance design).Dedicated to Prof. Dr. W. Fleischhacker, on occasion of his 60th anniversary     

Molecular cyclicity and centricity of polycyclic graphs. I. Cyclicity based on resistance distances or reciprocal distances

Rules for molecular cyclicity based on the global indices resulting from reciprocal distances (Harary number, H) or from resistance distances (Kirchhoff number, Kf) were tested in comparison with those elaborated earlier by means of the Wiener index, W. The Harary number and the Wiener number were found to match molecular cyclicity in an almost identical manner. The Kirchhoff number also generally follows cyclicity trends described previously. H is slightly less degenerate than is W, but Kf has practically no degeneracy in the graphs investigated here. Being much more discriminating than the Wiener number (i.e., practically nondegenerate), Kf allowed the formulation of new rules for systems formed from linearly condensed ribbons of even-membered rings with different sizes as well as for branched ribbons. The topological cyclicity patterns are thus reformulated in an extended basis, proceeding from three different graph metrics. © 1994 John Wiley & Sons, Inc.     

Fuzzy clustering as a means of selecting representative conformers and molecular alignments

This paper describes the first application of fuzzy c-means clustering for the selection of representatives from assemblies of conformations or alignments. In case of alignments, their quality is taken into account using a weighted c-means scheme, developed in this work. The performance of fuzzy cluster validity measures, such as compactness, partition function, and entropy, are studied on several examples, but the visual 3D representation of data points is shown to be most beneficial in determining the optimum number of clusters. Fuzzy clustering is expected to perform better than crisp clustering methods in cases where there are a significant number of "outliers", such as in molecular dynamics simulations and molecular alignments.     

Software news and updates. Carma: a molecular dynamics analysis program

A computer program has been developed to aid the analysis of molecular dynamics trajectories. The program is tuned for macromolecular large-scale problems and supports features such as removal of global translations-rotations of the solute, calculation of average distance maps and their corresponding standard deviations, calculation of the variance-covariance and cross-correlation matrices, and principal component analysis of trajectories with the added ability to create artificial trajectories based on selected eigenvectors. Limited graphics (trajectory viewing) capabilities are also available.     

Ab Initio computed molecular structures and energies of the conformers of glucose

Ab initio computations indicate the existence of several stable and some unstable conformers in isolated α and β glucose molecules. All of the lower-energy conformers exhibit a strikingly regular pattern of internal hydrogen bonding. Five such stable structures have been identified for each of the α and β anomers, differing primarily in the orientation of the CH₂OH group. In each conformer, the α anomer is predicted to be lower in energy than the corresponding conformers of β anomer. The difference is about 2 kcal/mol in the 4-31G basis but only 0.4 kcal/mol in the 6-31G* basis. It is found that the electronic contributions to the free energy difference stabilize the α anomer while the nuclear motion contributions stabilize the β anomer. The implications of these predictions and the future investigations required to understand the relative stabilities of the two anomers are pointed out. © 1992 by John Wiley & Sons, Inc.     

Novel biosensor for the rapid measurement of estrogen based on a ligand-receptor interaction.

A bioaffinity sensor was developed aiming at the detection of estrogen. This biosensor system is based on the specific binding of estrogen to its receptor immobilized on a gold disk electrode. The recombinant DNA encoding human estrogen receptor ligand-binding domain was expressed in bacteria using the histidine-tag fusion system. The expression of the fusion protein was under control of a bacteriophage T7 promoter, and the protein was purified under native conditions by affinity chromatography, which is based on a specific interaction between a histidine-tag, located in the N-terminus of the protein, and the Ni(II) chelate adsorbent. The protein was immobilized on an Au-electrode with Ni(II)-mediated chemisorption using a histidine tag and thiol-modified iminodiacetic acid. Cyclic voltammetric measurements showed that the reversible electrochemical reaction of a ferrocyanide/ferricyanide redox couple was suppressed by the presence of estrogen in a concentration-dependent manner. It seems reasonable to suppose that the electrostatic property of the protein layer on the electrode surface was altered by complexation with estrogen. These data suggest that this biosensor is applicable to the evaluation binding activities of the chemicals toward the human estrogen receptor.     

TGSA: A molecular superposition program based on topo‐geometrical considerations

In this article, a new molecular alignment procedure to provide general‐purpose, fast, automatic, and user‐intuitive three‐dimensional molecular alignments is presented. This procedure, called Topo‐Geometrical Superposition Approach (TGSA), is only based on comparisons of atom types and interatomic distances; hence, the procedure can handle large molecular sets within affordable computational costs. The method is able to accurately align 3D structures using the common molecular substructures, as inferred by the bonding pattern (atom correspondences), where present. The algorithm has been implemented into a program named TGSA99, and it has been tested over eight different molecular sets: flavilium salts, amino acids, indole derivatives, AZT, steroids, anilide derivatives, poly‐aromatic‐hydrocarbons, and inhibitors of thrombine. The TGSA algorithm performance is evaluated by means of computational time, number of superposed atoms, and index of fit between the compared structures. © 2000 John Wiley & Sons, Inc. J Comput Chem 22: 255–263, 2001     

Chemiluminescence sensor for the determination of perphenazine based on a molecular imprinted polymer.

A luminol-potassium ferricyanide-perphenazine CL reaction with high sensitivity for the determination of perphenazine was found. A perphenazine molecular imprinted polymer (MIP) was synthesized. Using the perphenazine MIP as a recognition material, a novel molecular imprinting-chemiluminescence (MI-CL) sensor for the determination of perphenazine was made based on the luminol-potassium ferricyanide-perphenazine CL reaction. The sensor displayed good selectivity and high sensitivity. The linear-response range of the sensor was 5.0 x 10(-8) - 1.0 x 10(-5) g/ml with a linear correlation coefficient of 0.9961. The detection limit was 2 x 10(-8) g/ml perphenazine and the relative standard deviation for 1.0 x 10(-6) g/ml perphenazine solution was 3.7% (n = 11). The sensor was applied to the determination of perphenazine in urine samples with satisfactory results.     

SSC: a tool for constructing libraries for systematic screening of conformers

Even a relatively small molecule with 10-20 atoms might have a few local minima, which correspond to different conformers. The number of local minima quickly increases with molecular size and the most common algorithms, driven by calculated forces, frequently identify a minimum, which is closest to the initial structure, rather than the most stable conformer. Here we discuss how to perform a systematic search of the conformational space for a chain-like molecule. Our approach is fully automated and a user has control which chemical bonds will be probed and with which increments. Moreover, whole fragments of the molecule, which are adjacent to each selected rotational bond, are rotated in a properly selected cylindrical coordinate system and unchemical hybridizations and some "clashes" between neighboring groups, which are common when standard Z-matrices are used, are avoided. A library of potentially relevant conformers is created with a tool, which we call SSC, denoting Systematic Screening of Conformers. Each member of the library is prescreened at a predefined level of theory and the most promising conformers are identified. Finally, they are further evaluated at a higher level of theory to identify the most stable structures and their physicochemical properties. As an example, we demonstrate the results of this approach for 2'-deoxycytidine.