Study of flexible joints and permanent bends in DNA fragments by brownian dynamics simulations

The use of Brownian dynamics simulations to investigate the presence of structural (kinks) and dynamic (bulges) anomalies in short DNA stretches is analyzed in connection with a string-of-beads model. A scaling method to choose the hydrodynamic translational and rotational parameters of the beads is proposed and tested on straight, kinked and bulged DNA fragments 17 nm long. The model reproduces the rigid-body rotational diffusion for the straight DNA and for the fluorescence polarization anisotropy decay of the kinked and bulged DNAs the model predicts a different behavior which is found experimentally. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 23 November 1998     

Peptides in membranes: assessment of environmental effects via simulations using an implicit solvation model

A recently developed implicit solvation model is applied to Monte Carlo simulations of peptides in bilayer-mimetic and polar environments. The model employs the formalism of atomic solvation parameters and reproduces experimental data. Solvent effects on the␣structure of the following peptides were studied: 20-residue poly-Leu and poly-Val, transmembrane helix A of bacteriorhodopsin, magainin2. It was shown that a␣membrane-like environment considerably promotes α-helix formation (all the peptides were found to be α-helical), while simulations in water reveal helix distortion. Consistency of the results with experimental data and further implications of the model are discussed. Received: 24 April 1998 / Accepted: 3 September 1998 / Published online: 10 December 1998     

The spin-forbidden reaction CH(2∏) + N2→ HCN + N(4S) revisited I. Ab initio study of the potential energy surfaces

High-level ab initio electronic structure theories have been applied to investigate the detailed reaction mechanism of the spin-forbidden reaction CH(²∏) + N₂ → HCN + N(⁴S). The G2M(RCC) calculations provide accurate energies for the intermediates and transition states involved in the reaction, whereas the B3LYP/6-311G(d,p) method overestimates the stability of some intermediates by as much as about 10 kcal/mol. A few new structures have been found for both the doublet and quartet electronic states, which are mainly involved in the dative pathways. However, due to the higher energies of these structures, the dominant mechanism remains the one involving the C ₂ intersystem-crossing step. The C ₂ minima on the seam of crossing (MSX) structures and the spin-orbit coupling between the doublet and quartet electronic states are rather close to those found in previous studies. Vibrational frequencies orthogonal to the normal of the seam which have been applied in a separate publication to calculate the rate of the CH(²∏) + N₂ → HCN + N(⁴S) reaction with a newly proposed nonadiabatic transition-state theory for spin-forbidden reactions have been calculated at the MSX from first principles. Received: 23 June 1998 / Accepted: 21 September 1998 / Published online: 8 February 1999     

The structure and energetics of cryolite melts

A method of solvation energy computation is proposed for ions and molecules in the environment of an ionic melt, based on the approximation of the ionic melt as an ideal conductor. The method is used to compute equilibrium constants of some equilibria in cryolite melt. Theoretically obtained results predict that aluminium is bound in tetrafluorocomplexes AlF₄ ⁻. Received: 16 March 1998 / Accepted: 19 June 1998 / Published online: 7 October 1998     

Computation of stationary points via a homotopy method

A homotopy method is presented that locates both minimizers and saddle points of energy functions in an efficient manner. In contrast to other methods, it makes possible the exploration of large parts of potential energy surfaces. Along a homotopy path stationary points of odd and even order occur alternately. A path tracing procedure requiring only gradients and at most one evaluation of the Hessian matrix is given. Test results on a model potential and three MINDO/3 potentials are reported. Received: 6 May 1996 / Accepted: 2 April 1998 / Published online: 23 June 1998     

Subduction of Q-conjugacy representations and characteristic monomials for combinatorial enumeration

A new method of combinatorial enumeration is presented. The subduction of Q-conjugacy representations gives a characteristic subduction table and a characteristic monomial table. A cycle index is defined on the basis of such monomials and used for combinatorial enumeration of isomers. Received: 10 October 1997 / Accepted: 13 February 1998 / Published online: 17 June 1998     

An implementation of the conductor-like screening model of solvation within the Amsterdam density functional package

The conductor-like screening model (COSMO) of solvation has been implemented in the Amsterdam density functional program with maximum flexibility in mind. Four cavity definitions have been incorporated. Several iterative schemes have been tested for solving the COSMO equations. The biconjugate gradient method proves to be both robust and memory-conserving. The interaction between the surface charges and the electron density may be calculated by integrating over either the fitted or exact density, or by calculating the molecular potential. A disk-smearing algorithm is applied in the former case to avoid singularities. Several self-consistent field/COSMO coupling schemes were examined in an attempt to reduce computational effort. A gradient-preserving algorithm for removing outlying charge has been implemented. Preliminary optimized radii are given. Applications to the benzene oxide-oxepin valence tautomerization and to glycine conformation are presented. Received: 13 November 1998 / Accepted: 16 December 1998 / Published online: 16 March 1999     

Electron affinities, excitation energies and ionization potentials of the transition metals (I) Sc and Ti

The electron affinities of the Sc and Ti atoms have been obtained by configuration interaction calculations. Energy convergence with respect to the systematic expansion of both the one-electron basis set and the configuration space was investigated for valence electrons, and the inclusion of correlation contributions from core electrons and relativistic effects gave the electron affinities of 0.181 eV and 0.163 eV for Sc and Ti, respectively. These are in excellent agreement with the observed values of 0.189 ± 0.020 eV and 0.080 eV. The same approach was applied for the first excited states and positive ions of both atoms. Excellent agreement with the experimental results was also obtained for these states. Received: 16 February 1998 / Accepted: 2 April 1998 / Published online: 23 June 1998     

Sequence classification of water channels and related proteins in view of functional predictions

We have worked with a classification method based upon a notion of probabilistic similarity or “likelihood of similarity” between aligned sequences. One important parameter, among others, affecting the sequence similarities and hence the classification results is the amino acid similarity matrix. We present a method for choosing the most adapted matrix to classify protein sequences. This method has been applied to the transmembrane channels of the major intrinsic protein (MIP) family. At present, two functional subgroups have been well characterized in this family: (1) specific water transport by the aquaporins and (2) small neutral solutes transport. The aim of the present study is to show the usefulness of the classification method in the prediction of sequence segments important for substrate selectivity. Moreover, we show that this method can also be used to predict the function of undetermined MIP proteins. The method could be applied to other protein families as well. Received: 24 April 1998 / Accepted: 4 August 1998 / Published online: 2 November 1998     

Extracting parameters for base-pair level models of DNA from molecular dynamics simulations

 A method is described to extract a complete set of sequence-dependent energy parameters for a rigid base-pair model of DNA from molecular dynamics (MD) simulations. The method is properly consistent with equilibrium statistical mechanics and leads to effective inertia parameters for the base-pair units as well as stacking and stiffness parameters for the base-pair junctions. We give explicit formulas that yield a complete set of base-pair model parameters in terms of equilibrium averages that can be estimated from a time series generated in an MD simulation. The expressions to be averaged depend strongly both on the choice of coordinates used to describe rigid-body orientations and on the choice of strain measures at each junction. Received: 12 July 2000 / Accepted: 5 January 2001 / Published online: 3 May 2001     

Adaptive umbrella sampling of the potential energy: modified updating procedure of the umbrella potential and application to peptide folding

Adaptive umbrella sampling of the potential energy is used as a search method to determine the structures and thermodynamics of peptides in solution. It leads to uniform sampling of the potential energy, so as to combine sampling of low-energy conformations that dominate the properties of the system at room temperature with sampling of high-energy conformations that are important for transitions between different minima. A modification of the procedure for updating the umbrella potential is introduced to increase the number of transitions between folded and unfolded conformations. The method does not depend on assumptions about the geometry of the native state. Two peptides with 12 and 13 residues, respectively, are studied using the CHARMM polar-hydrogen energy function and the analytical continuum solvent potential for treatment of solvation. In the original adaptive umbrella sampling simulations of the two peptides, two and six transitions occur between folded and unfolded conformations, respectively, over a simulation time of 10 ns. The modification increases the number of transitions to 6 and 12, respectively, in the same simulation time. The precision of estimates of the average effective energy of the system as a function of temperature and of the contributions to the average effective energy of folded conformations obtained with the adaptive methods is discussed. Received: 11 July 1998 / Accepted: 22 September 1998 / Published online: 17 December 1998     

Using simplified protein representation as a reference potential for all-atom calculations of folding free energy

An effective approach for evaluating folding free-energy surfaces of explicit all-atom models is developed and examined. This approach is based on using the potential of a simplified protein model as a reference potential for calculating the free energy of the corresponding explicit model. Preliminary results are presented for the folding free energy of a 12-residue helix. The potential of the method for studies of protein-folding processes is discussed, emphasizing the ability to determine the difference between the results of simplified and explicit models. This can help in establishing the validity of simplified folding models. Received: 30 December 1998 / Accepted: 27 January 1999 / Published online: 5 May 1999     

Implementation of the IMOMM methodology for performing combined QM/MM molecular dynamics simulations and frequency calculations

A technique for implementing the integrated molecular orbital and molecular mechanics (IMOMM) methodology developed by Maseras and Morokuma that is used to perform combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, frequency calculations and simulations of macromolecules including explicit solvent is presented. Although the IMOMM methodology is generalized to any coordinate system, the implementation first described by Maseras and Morokuma requires that the QM and MM gradients be transformed into internal coordinates before they are added together. This coordinate transformation can be cumbersome for macromolecular systems and can become ill-defined during the course of a molecular dynamics simulation. We describe an implementation of the IMOMM method in which the QM and MM gradients are combined in the cartesian coordinate system, thereby avoiding potential problems associated with using the internal coordinate system. The implementation can be used to perform combined QM/MM molecular dynamics simulations and frequency calculations within the IMOMM framework. Finally, we have examined the applicability of thermochemical data derived from IMOMM framework. Finally, we have examined the applicability of thermochemical data derived from IMOMM frequency calculations. Received: 11 May 1998 / Accepted: 14 August 1998 / Published online: 16 November 1998     

Chemical potential inequality principle

Regional density functional theory has been extended to treat irreversible thermodynamic electronic processes for application to adiabatic electron-transfer processes of chemical reactions. Onsager's local equilibrium hypothesis is slightly modified to take into account the quantum mechanical nature of the electron. The quantum mechanical interference effect has been demonstrated to be included in the entropy production rate formula associated with electron transfer through an interface. A new formula for the determination of the transition state of a chemical reaction has been postulated that corresponds to the maximum of the regional electron transferability. A quantum mechanical law of mass action has been established and applied to prove the regional electrochemical potential inequality principle. Received: 1 July 1998 / Accepted: 2 September 1998 / Published online: 8 February 1999     

Structure effects on isomerization pathways and vibrational spectra of epoxysaccharides: a numerical study

Results of a comparative analysis of conformational possibilities of the hexopyranose ring of six epoxysaccharides differing from each other by the position of the oxirane ring within the limits of the hexapyranose ring and having different orientations of substituents and different positions of the oxirane ring with respect to the skeleton plane of the molecules are presented. Numerical simulations based on the Wiberg and Boyd method made it possible to determine all the stationary forms in which anhydropyranose rings can exist. The effect of various structural factors on the character of conformational transformations, heights of transition barriers, and the energy of stationary forms has been investigated. Normal vibrational modes of the stationary forms of the compounds were calculated using molecular mechanics. Based on results of our simulations, we predict a strong effect of steric factors on the vibrational spectra of sugar epoxides. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 17 December 1998     

Implicit two-phase solvation model as a tool to assess conformation and energetics of proteins in membrane-mimetic media

 A heterogeneous implicit membrane-mimetic model is applied to simulations of membrane proteins. The model employs atomic solvation parameters for gas–water and gas–cyclohexane transfer. It is used to analyze structure, energetics, and orientation with respect to the bilayer of two polypeptides with different modes of membrane binding – hydrophobic segment of human glycophorin A (GpA) and cytotoxin II from Naja naja oxiana snake venom (CTX). The native state of GpA represents a transmembrane (TM) α helix, while CTX is a water-soluble protein, which is able to interact with the cell membrane. The conformational space of the polypeptides was explored in Monte Carlo simulations. The results show that the most stable conformers of GpA represent a TM α helix. They are additionally stabilized by an applied TM voltage. The results also show that CTX inserts with its three loops, does not cross the hydrophobic layer, and stays partially immersed in the membrane. This agrees well with the experimental data, thus confirming the validity of the solvation model. Received: 13 June 2000 / Accepted: 15 September 2000 / Published online: 19 January 2001     

Transition-state structures for describing the enzyme-catalyzed mechanisms of rubisco

The carboxylation and oxygenation processes of a model substrate, 3,4-dihydroxy-2-pentanone, have been theoreticaly characterized as a set of steps, mimicking the corresponding reactions of D-ribulose-1,5-bisphosphate catalyzed by rubisco. A theoretical characterization is carried out of transition-state structures and possible molecular intermediates represented as saddle points of index 1 and minimum energy structures, respectively. The quantum chemical characterization, at the HF/3-21G calculation level, of these stationary points is used to rationalize and to discuss both catalyzed sequences. The reported set of these stationary points maps out most experimental aspects of the reaction pathways for the real system. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 10 December 1998     

Visual valence bond rules for chemical reactions

A symmetry-adaptation rule of the valence bond structure for concerted reactions was developed within the bonded tableau valence bond formalism. According to a symmetry analysis of the valence bond structure segments accounting for the reaction, one can predict whether a chemical process is favored or unfavored. This method is based on conceptual resonance theory and the visual valence bond approach, without carrying out any explicitly theoretical calculations to know orbital details. Furthermore, by imposing a phase factor on each bonding pair, namely, the phase alternation postulate, the mechanisms of the concerted reactions can be easily outlined. These rules have been applied to organic and inorganic reactions including the participation of biradicals and species with multi-reference character. Received: 7 September 1998 / Accepted: 9 November 1998 / Published online: 1 February 1999     

Electric and magnetic multipole integrals in STO basis sets: an analytical approach

A new method for the evaluation of one- and two-centre magnetic and electric multipole integrals for Slater-type functions is presented. The method is strictly analytical in that no approximations of any kind are involved. Two simple functions, ℐ₁ ^aug and ℐ₂ ^aug, are introduced, which employ only functions that are well known in electronic structure theory. With the use of augmentation exponents these functions apply to multipole integrals as well as other one-electron integrals, such as nuclear attraction integrals. The proposed method includes the analytic determination of derivatives of the integrals with respect to atomic displacements. Some illustrative test calculations are presented and compared to results from the literature. Received: 20 April 1998 / Accepted: 13 October 1998 / Published online: 1 February 1999     

On the automatic restricted-step rational-function-optimization method

The rational function optimization algorithm is one of the widely used methods to search stationary points on surfaces. However, one of the drawbacks of this method is the step reduction procedure to deal with the overstepping problem. We present and comment on a method such that the step obtained from the solution of the rational function equations possesses the desired correct length. The analysis and discussion of the method is mainly centered on the location and optimization of transition states. Received: 18 June 1998 / Accepted: 17 September 1998 / Published online: 23 November 1998