Algorithm for normal mode analysis with general internal coordinates

A technique for performing normal vibrational analysis for biological macromolecules using general internal coordinates is proposed. The technique is based on the conventional algorithm for calculating the second derivatives of potential and kinetic energies using intramolecular dihedral angles, intermolecular translation, and rotation as variables [Braun, W. et al., J Phys Soc Jpn 1984, 53, 3269]. We extend the algorithm to include more general internal coordinates, bond stretching, angle bending, and so forth, without assuming two-body interactions. The essential point is the separation of the variables for potential functions and vibrational analysis. With our technique, we can arbitrarily choose any combination of internal coordinates as variables, free from the functional form of potential energy. We can analyze complex systems such as a multiple molecular system including solvents or a transition state of chemical reactions. In addition, mixed use of the potentials of molecular mechanics and quantum chemistry is possible.     

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.     

Internal‐to‐Cartesian back transformation of molecular geometry steps using high‐order geometric derivatives

In geometry optimizations and molecular dynamics calculations, it is often necessary to transform a geometry step that has been determined in internal coordinates to Cartesian coordinates. A new method for performing such transformations, the high‐order path‐expansion (HOPE) method, is here presented. The new method treats the nonlinear relation between internal and Cartesian coordinates by means of automatic differentiation. The method is reliable, applicable to any system of internal coordinates, and computationally more efficient than the traditional method of iterative back transformations. As a bonus, the HOPE method determines not just the Cartesian step vector but also a continuous step path expressed in the form of a polynomial, which is useful for determining reaction coordinates, for integrating trajectories, and for visualization. © 2013 Wiley Periodicals, Inc.     

Expressing kinetic‐energy for vibration–rotation motions in general rotating systems of axes and introducing quasirectilinear vibrational coordinates to simplify Hamiltonian forms

To understand how the internal and rotational motions of a polyatomic system depend on which rotating system of axes is selected, we derived the explicit form of the atomic velocities determined by an observer stationed on the general rotating system of axes. Using the derived velocities, we formulated the kinetic energy expression for vibration–rotation motions with respect to the rotating system of axes. From this expression, we clarified covariant metric tensors under zero angular momentum, which have been confused with an erroneous expression even in the professional literature, and the relationship between the kinetic energy expression and the rotating system of axes. Furthermore, to simplify the Hamiltonian form, we introduced quasirectilinear vibrational coordinates to describe the Hamiltonian. The resulting Hamiltonian form is superior to those of the previous studies in that the kinetic and potential energy expressions are simple and the vibrational frequencies are independent of the original internal coordinates used. In fact, we show that its application for three examples is useful. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 83: 22–29, 2001     

Exploring the effect of anharmonicity of molecular vibrations on thermodynamic properties

Thermodynamic properties of selected small and medium size molecules were calculated using harmonic and anharmonic vibrational frequencies. Harmonic vibrational frequencies were obtained by normal mode analysis, whereas anharmonic ones were calculated using the vibrational self-consistent field (VSCF) method. The calculated and available experimental thermodynamic data for zero point energy, enthalpy, entropy, and heat capacity are compared. It is found that the anharmonicity and coupling of molecular vibrations can play a significant role in predicting accurate thermodynamic quantities. Limitations of the current VSCF method for low frequency modes have been partially removed by following normal mode displacements in internal, rather than Cartesian, coordinates.     

The use of internal cartesian coordinates for describing molecular vibrations

An analysis of the influence of isotope substitution on the system of electronic-nuclear equations for an arbitrary molecular system was used as a basis for formulating invariance conditions with respect to isotope substitution of the potential energy surface written in the Cartesian coordinates rigidly bound with the center of mass of the molecule (internal Cartesian coordinates). This property of the potential function obviates the necessity of using curvilinear natural coordinates, which can be replaced by Cartesian coordinates, in theoretical studies of the vibrational spectra of molecules and their isotopomers and in solving the direct and inverse anharmonic problems. An equation for the quantum-mechanical Hamiltonian of a normal molecule in internal Cartesian coordinates was obtained.     

Scaled in Cartesian Coordinates Ab Initio Molecular Force Fields of DNA Bases: Application to Canonical Pairs

The model of Regularized Quantum Mechanical Force Field (RQMFF) was applied to the joint treatment of ab initio and experimental vibrational data of the four primary nucleobases using a new algorithm based on the scaling procedure in Cartesian coordinates. The matrix of scaling factors in Cartesian coordinates for the considered molecules includes diagonal elements for all atoms of the molecule and off-diagonal elements for bonded atoms and for some non-bonded atoms (1–3 and some 1–4 interactions). The choice of the model is based on the results of the second-order perturbation analysis of the Fock matrix for uncoupled interactions using the Natural Bond Orbital (NBO) analysis. The scaling factors obtained within this model as a result of solving the inverse problem (regularized Cartesian scale factors) of adenine, cytosine, guanine, and thymine molecules were used to correct the Hessians of the canonical base pairs: adenine–thymine and cytosine–guanine. The proposed procedure is based on the block structure of the scaling matrix for molecular entities with non-covalent interactions, as in the case of DNA base pairs. It allows avoiding introducing internal coordinates (or coordinates of symmetry, local symmetry, etc.) when scaling the force field of a compound of a complex structure with non-covalent H-bonds.     

Ferreting out correlations from trajectory data

Thermally driven materials characterized by complex energy landscapes, such as proteins, exhibit motions on a broad range of space and time scales. Principal component analysis (PCA) is often used to extract modes of motion from protein trajectory data that correspond to coherent, functional motions. In this work, two other methods, maximum covariance analysis (MCA) and canonical correlation analysis (CCA) are formulated in a way appropriate to analyze protein trajectory data. Both methods partition the coordinates used to describe the system into two sets (two measurement domains) and inquire as to the correlations that may exist between them. MCA and CCA provide rotations of the original coordinate system that successively maximize the covariance (MCA) or correlation (CCA) between modes of each measurement domain under suitable constraint conditions. We provide a common framework based on the singular value decomposition of appropriate matrices to derive MCA and CCA. The differences between and strengths and weaknesses of MCA and CCA are discussed and illustrated. The application presented here examines the correlation between the backbone and side chain of the peptide met-enkephalin as it fluctuates between open conformations, found in solution, to closed conformations appropriate to when it is bound to its receptor. Difficulties with PCA carried out in Cartesian coordinates are found and motivate a formulation in terms of dihedral angles for the backbone atoms and selected atom distances for the side chains. These internal coordinates are a more reliable basis for all the methods explored here. MCA uncovers a correlation between combinations of several backbone dihedral angles and selected side chain atom distances of met-enkephalin. It could be used to suggest residues and dihedral angles to focus on to favor specific side chain conformers. These methods could be applied to proteins with domains that, when they rearrange upon ligand binding, may have correlated functional motions or, for multi-subunit proteins, may exhibit correlated subunit motions.     

The impact of approximate VSCF schemes and curvilinear coordinates on the anharmonic vibrational frequencies of formamide and thioformamide

The accuracy of the vibrational self-consistent field (VSCF) method for the computation of anharmonic vibrational frequencies in the infrared (IR) spectrum of formamide and thioformamide is investigated. The importance of triple potentials in the commonly used hierarchical expansion of the potential energy surface (PES) is studied in detail. The PES is expanded in terms of Cartesian as well as internal coordinate normal mode displacements. It is found that triples play an important role when using rectilinear coordinates. A VSCF computation based on rectilinear displacements exhibits serious shortcomings which are only remedied by a large vibrational configuration interaction (VCI) treatment including triple potentials. These limitations are partially removed when using curvilinear coordinates. The merits and disadvantages of either type of displacements for the generation of the PES are discussed.     

First-principle calculation of reduced masses in vibrational analysis using generalized internal coordinates: some crucial aspects and examples

In this paper we present and analyze the most essential aspects of reduced masses along generalized internal coordinates. The definition of reduced masses in the internal coordinate formalism is established through the Wilson G-matrix concept and includes sophisticated relations between internal and Cartesian coordinates. Moreover, reduced masses in internal coordinates are, in general, no longer constant but coordinate-dependent. Based on the approach presented earlier [Stare, J.; Balint-Kurti, G. G. J. Phys. Chem. A 2003, 107, 7204-7214] and on our experience with reduced masses discussed in this paper, we have developed a robust program for the calculation of Wilson G-matrix elements and their functional coordinate dependence. The approach is based on the first principles and can be used in virtually any (internal) coordinate set. Since the program allows for projection of any kind of nuclear motion on the selected internal coordinates, the method is particularly suitable for ab initio or DFT potential energy functions calculated by partial geometry optimization. Moreover, reduced masses obtained by this program can be used as a decision tool for selecting the most appropriate internal coordinates for the considered vibrational problem and for the inclusion or omission of the kinetic coupling terms in the vibrational Hamiltonian.     

The choice of internal coordinates in complex chemical systems

This article presents several considerations for the appropriate choice of internal coordinates in various complex chemical systems. The appropriate and black box recognition of internal coordinates is of fundamental importance for the extension of internal coordinate algorithms to all fields where previously Cartesian coordinates were the preferred means of geometry manipulations. Such fields range from local and global geometry optimizations to molecular dynamics as applied to a wide variety of chemical systems. We present a robust algorithm that is capable to quickly determine the appropriate choice of internal coordinates in a wide range of atomic arrangements. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Representing and selecting vibrational angular momentum states for quasiclassical trajectory chemical dynamics simulations

Linear molecules with degenerate bending modes have states, which may be represented by the quantum numbers N and L. The former gives the total energy for these modes and the latter identifies their vibrational angular momentum jz. In this work, the classical mechanical analog of the N,L-quantum states is reviewed, and an algorithm is presented for selecting initial conditions for these states in quasiclassical trajectory chemical dynamics simulations. The algorithm is illustrated by choosing initial conditions for the N = 3 and L = 3 and 1 states of CO2. Applications of this algorithm are considered for initial conditions without and with zero-point energy (zpe) included in the vibrational angular momentum states and the C-O stretching modes. The O-atom motions in the x,y-plane are determined for these states from classical trajectories in Cartesian coordinates and are compared with the motion predicted by the normal-mode model. They are only in agreement for the N = L = 3 state without vibrational angular momentum zpe. For the remaining states, the Cartesian O-atom motions are considerably different from the elliptical motion predicted by the normal-mode model. This arises from bend-stretch coupling, including centrifugal distortion, in the Cartesian trajectories, which results in tubular instead of elliptical motion. Including zpe in the C-O stretch modes introduces considerable complexity into the O-atom motions for the vibrational angular momentum states. The short-time O-atom motions for these trajectories are highly irregular and do not appear to have any identifiable characteristics. However, the O-atom motions for trajectories integrated for substantially longer period of times acquire unique properties. With C-O stretch zpe included, the long-time O-atom motion becomes tubular for trajectories integrated to approximately 14 ps for the L = 3 states and to approximately 44 ps for the L = 1 states.     

Rotational symmetry of the molecular potential energy in the Cartesian coordinates

We consider the molecular Born-Oppenheimer potential energy as a function of atomic Cartesian coordinates and discuss the non-stationary Hessian properties arising due to rotational symmetry. A connection with the extended Hessian theory is included. New applications of Cartesian representation for examining and correcting raw numerical Hessian data and a simple formulation of harmonic vibrational analysis of partially optimized systems are proposed. Exemplary calculations for the porphyrin molecule with an internal proton transfer are reported. We also develop the normal transformation method to incorporate the rotational symmetry into the approximate analytical potentials, which are parametrized in the Cartesian coordinates. The transformation converts the coordinates from the space fixed frame to the frame which translates and rotates with the molecule and is determined by the Eckart conditions. New simple analytical formulas for the first and second derivatives of the transformed potential are derived. This fast method can be used to calculate the potential and its derivatives in the simulations of chemical reaction dynamics in the space fixed Cartesian frame without the need to constrain the molecular rotation or to define the local non-redundant internal coordinates.     

Calculating energy levels of isomerizing tetra-atomic molecules. II. The vibrational states of acetylene and vinylidene

A general, full-dimensional computational method for the accurate calculation of rotationally and vibrationally excited states of tetra-atomic molecules is further developed. The resulting computer program may be run in serial and parallel modes and is particularly appropriate for molecules executing wide-amplitude motions and isomerizations. An application to the isomerizing acetylene/vinylidene system is presented. Large-scale calculations using a coordinate system based on orthogonal satellite vectors have been performed in six dimensions and vibrational term values and wave functions for acetylene and vinylidene states up to approximately 23 000 cm(-1) above the potential minimum have been determined. This has permitted the characterization of acetylene and vinylidene states at and above the isomerization barrier. These calculations employ more extensive vibrational basis sets and hence consider a much higher density of states than in any variational calculations reported hitherto for this system. Comparison of the calculated density of states with that determined empirically suggests that our calculations are the most realistic achieved for this system to date. Indeed more states have been converged than in any previous study of this system. Calculations on lower lying excited states of acetylene based on HC-CH diatom-diatom coordinates give nearly identical results to those based on orthogonal satellite vectors. Comparisons are also made with calculations based on HH-CC diatom-diatom coordinates.     

Using redundant internal coordinates to optimize equilibrium geometries and transition states

A redundant internal coordinate system for optimizing molecular geometries is constructed from all bonds, all valence angles between bonded atoms, and all dihedral angles between bonded atoms. Redundancies are removed by using the generalized inverse of the G matrix; constraints can be added by using an appropriate projector. For minimizations, redundant internal coordinates provide substantial improvements in optimization efficiency over Cartesian and nonredundant internal coordinates, especially for flexible and polycyclic systems. Transition structure searches are also improved when redundant coordinates are used and when the initial steps are guided by the quadratic synchronous transit approach. © 1996 by John Wiley & Sons, Inc.     

Construction of the free energy landscape of biomolecules via dihedral angle principal component analysis

A systematic approach to construct a low-dimensional free energy landscape from a classical molecular dynamics (MD) simulation is presented. The approach is based on the recently proposed dihedral angle principal component analysis (dPCA), which avoids artifacts due to the mixing of internal and overall motions in Cartesian coordinates and circumvents problems associated with the circularity of angular variables. Requiring that the energy landscape reproduces the correct number, energy, and location of the system's metastable states and barriers, the dimensionality of the free energy landscape (i.e., the number of essential components) is obtained. This dimensionality can be determined from the distribution and autocorrelation of the principal components. By performing an 800 ns MD simulation of the folding of hepta-alanine in explicit water and using geometric and kinetic clustering techniques, it is shown that a five-dimensional dPCA energy landscape is a suitable and accurate representation of the full-dimensional landscape. In the second step, the dPCA energy landscape can be employed (e.g., in a Langevin simulation) to facilitate a detailed investigation of biomolecular dynamics in low dimensions. Finally, several ways to visualize the multidimensional energy landscape are discussed.     

The PyPES library of high quality semi‐global potential energy surfaces

In this article, we present a Python‐based library of high quality semi‐global potential energy surfaces for 50 polyatomic molecules with up to six atoms. We anticipate that these surfaces will find widespread application in the testing of new potential energy surface construction algorithms and nuclear ro‐vibrational structure theories. To this end, we provide the ability to generate the energy derivatives required for Taylor series expansions to sixth order about any point on the potential energy surface in a range of common coordinate systems, including curvilinear internal, Cartesian, and normal mode coordinates. The PyPES package, along with FORTRAN, C, MATLAB and Mathematica wrappers, is available at http://sourceforge.net/projects/pypes-lib . © 2015 Wiley Periodicals, Inc.     

CAMERRA: An analysis tool for the computation of conformational dynamics by evaluating residue–residue associations

A computational method which extracts the dominant motions from an ensemble of biomolecular conformations via a correlation analysis of residue–residue contacts is presented. The algorithm first renders the structural information into contact matrices, then constructs the collective modes based on the correlated dynamics of a selected set of dynamic contacts. Associated programs can bridge the results for further visualization using graphics software. The aim of this method is to provide an analysis of conformations of biopolymers from the contact viewpoint. It may assist a systematical uncovering of conformational switching mechanisms existing in proteins and biopolymer systems in general by statistical analysis of simulation snapshots. In contrast to conventional correlation analyses of Cartesian coordinates (such as distance covariance analysis and Cartesian principal component analysis), this program also provides an alternative way to locate essential collective motions in general. Herein, we detail the algorithm in a stepwise manner and comment on the importance of the method as applied to decoding allosteric mechanisms. © 2018 Wiley Periodicals, Inc.     

Transfer of molecular property tensors in cartesian coordinates: A new algorithm for simulation of vibrational spectra

A direct transfer of Cartesian molecular force fields (FF) and electric property tensors is tested on model systems and compared to transfer in internal coordinates with an aim to improve simulation of vibrational spectra for larger molecules. This Cartesian transformation can be implemented easily and offers greater flexibility in practical computations. It can be also applied for transfer of anharmonic derivatives. The results for model calculations of the force field and vibrational frequencies for N-methylacetamide show that our method removes errors associated with numerical artifacts caused by nonlinearity of the otherwise required Cartesian to internal coordinate transformation. For determination of IR absorption and vibrational circular dichroism intensities, atomic polar and axial tensors were also transferred in the Cartesian representation. For the latter, which are dependent upon the magnetic dipole operator, a distributed origin gauge is used to avoid an origin dependence. Comparison of the results of transferring ab initio FF and intensity parameters from an amide dimer fragment onto a tripeptide with those from a conventionally determined tripeptide FF document some limitations of the transfer method and its possible applications in the vibrational spectroscopy. Finally, application to determination of the FF and spectra for helical heptapeptide are presented and compared to experimental results. © 1997 by John Wiley & Sons, Inc.