The solution structure of a copper(II) compound of a new cyclic octapeptide by EPR spectroscopy and force field calculations

A new cyclic octapeptide, cyclo(Ile-Ser-(Gly)Thz-Ile-Thr-(Gly)Thz) (PatN), related to patellamide A, has been synthesized and reacted with copper(II) and base to form mono- and dinuclear complexes. The coordination environments around copper(TI) have been characterized by EPR spectroscopy. The solution structure of the thermodynamically most stable product, a purple dicopper(TI) compound, has been examined by simulating weakly dipole-dipole coupled EPR spectra based upon structural parameters obtained from force field (MM and MD) calculations. The MM-EPR method produces a saddle-shaped structure for [Cu-2(PatN)(OH2)(6)] that is similar to the known solution structure of patellamide A and the known solid-state structure of [Cu-2(AscidH(2))CO3(OH2)(2)]. Compared with the latter, [Cu-2(PatN)] has no carbonate bridge and a significantly flatter topology. The MM-EPR approach to solution-structure determination for paramagnetic metallopeptides may find wide applications to other metallopeptides and metalloproteins.

Fast docking using the CHARMM force field with EADock DSS.

The prediction of binding modes (BMs) occurring between a small molecule and a target protein of biological interest has become of great importance for drug development. The overwhelming diversity of needs leaves room for docking approaches addressing specific problems. Nowadays, the universe of docking software ranges from fast and user friendly programs to algorithmically flexible and accurate approaches. EADock2 is an example of the latter. Its multiobjective scoring function was designed around the CHARMM22 force field and the FACTS solvation model. However, the major drawback of such a software design lies in its computational cost. EADock dihedral space sampling (DSS) is built on the most efficient features of EADock2, namely its hybrid sampling engine and multiobjective scoring function. Its performance is equivalent to that of EADock2 for drug-like ligands, while the CPU time required has been reduced by several orders of magnitude. This huge improvement was achieved through a combination of several innovative features including an automatic bias of the sampling toward putative binding sites, and a very efficient tree-based DSS algorithm. When the top-scoring prediction is considered, 57% of BMs of a test set of 251 complexes were reproduced within 2 Å RMSD to the crystal structure. Up to 70% were reproduced when considering the five top scoring predictions. The success rate is lower in cross-docking assays but remains comparable with that of the latest version of AutoDock that accounts for the protein flexibility. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011.

SwissParam: a fast force field generation tool for small organic molecules.

The drug discovery process has been deeply transformed recently by the use of computational ligand-based or structure-based methods, helping the lead compounds identification and optimization, and finally the delivery of new drug candidates more quickly and at lower cost. Structure-based computational methods for drug discovery mainly involve ligand-protein docking and rapid binding free energy estimation, both of which require force field parameterization for many drug candidates. Here, we present a fast force field generation tool, called SwissParam, able to generate, for arbitrary small organic molecule, topologies, and parameters based on the Merck molecular force field, but in a functional form that is compatible with the CHARMM force field. Output files can be used with CHARMM or GROMACS. The topologies and parameters generated by SwissParam are used by the docking software EADock2 and EADock DSS to describe the small molecules to be docked, whereas the protein is described by the CHARMM force field, and allow them to reach success rates ranging from 56 to 78%. We have also developed a rapid binding free energy estimation approach, using SwissParam for ligands and CHARMM22/27 for proteins, which requires only a short minimization to reproduce the experimental binding free energy of 214 ligand-protein complexes involving 62 different proteins, with a standard error of 2.0 kcal mol(-1), and a correlation coefficient of 0.74. Together, these results demonstrate the relevance of using SwissParam topologies and parameters to describe small organic molecules in computer-aided drug design applications, together with a CHARMM22/27 description of the target protein. SwissParam is available free of charge for academic users at www.swissparam.ch.

Investigating the spontaneous formation of SDS micelle in aqueous solution using a coarse-grained force field

A 1µs Molecular Dynamic simulation was performed with a realistic model system of Sodium Dodecyl Sulfate (SDS) micelles in aqueous solution, comprising of 360 DS-, 360 Na+ and 90000 water particles. After 300 ns three different micellar shapes and sizes 41, 68 and 95 monomers, were observed. The process led to stabilization in the total number of SDS clusters and an increase in the micellar radius to 2.23 nm, in agreement with experimental results. An important conclusion, is be aware that simulations employed in one aggregate, should be considered as a constraint. Size and shape distribution must be analyzed.

Studies in molecular mechanics using spectroscopic data exact force field and bond in homogeneity

This thesis deals with some studies in molecular mechanic using spectroscopic data. It includes an improvement in the parameter technique for the evaluation of exact force fields, the introduction of a new and simple algebraic method for the force field calculation and a study of asymmetric variation of bonding forces along a bond.

Anharmonic force field of acetylene

The harmonic and anharmonic force field of acetylene has been determined in a least-squares calculation from recently determined data on the spectroscopic constants of various isotopic species (including the vibrational l-doubling constant). A general quadratic and cubic force field was used, but a constrained quartic force field containing only 8 of the 23 possible quartic constants. The results are discussed and compared with earlier work.

Microwave spectrum of CHD2CN and CHD2NC, and the harmonic force field and rz structure of methyl cyanide and isocyanide

The microwave spectra of CHD2CN and CHD2NC have been measured from 18 to 40 GHz; about 20 type A and 30 type C transitions have been observed for each molecule. These have been fitted to a Hamiltonian using 3 rotational constants, and 5 quartic and 4 sextic distortion constants, in the IrS reduction of Watson [in “Vibrational spectra and structure” Vol. 6 (1977)]; the standard error of the fit is 26 kHz. For methyl cyanide the 5 quartic distortion constants have been used to further refine the recent harmonic force field of Duncan et al. [J. Mol. Spectrosc. 69, 123 (1978)], but the changes are small. Finally, for both molecules, the harmonic force field has been used to determine zero point average moments of inertia Iz from the ground state rotational constants for many isotopic species, and these have been used to determine an rz structure. The results are compared with rs structure calculations.

Excited vibrational state microwave spectra, structure, and force-field of trifluorosilane

The J = 2−1 microwave spectrum of six isotopic species of HSiF3 has been observed and assigned in excited states of five of the six fundamental vibrations. The assignment is based on relative intensities, double resonance experiments, and trial anharmonic force constant calculations. Analysis of the spectra leads to experimental values for five of the constants, all three l-doubling constants qt, one Fermi resonance constant φ233, and one zeta constant. The harmonic force field has been refined to all the available data on vibration wavenumbers, centrifugal distortion constants, and zeta constants. The cubic anharmonic force field has been refined to the data on and qt constants, using two models: a valence force model with two cubic force constants for SiH and SiF stretching, and a more sophisticated model. With the help of these calculations, the following equilibrium structure has been determined: re(SiH) = 1.4468(±5) Å, re(SiF) = 1.5624(±1) Å, HSiF = 110.64(±3)°,

Redundant coordinates in harmonic force-field calculations

Redundancy relations between vibrational coordinates may be linear (as for rectilinear coordinates used in deriving a G matrix), or non-linear (as for curvilinear coordinates used in formulating model force fields). It is shown that geometrically defined internal coordinates are necessarily curvilinear. Hence it is shown that linear force constants can occur in model force field calculations involving redundant coordinates, in disagreement with the recent proposal of Gussoni and Zerbi.

The anharmonic force field of HCN and HCP

The quadratic, cubic, and quartic force field of HCN has been calculated by a least squares refinement to fit the most recent observed data on the vibration-rotation constants of HCN, DCN and H13CN. All of the observed parameters are fitted within their standard errors of observation. The corresponding parameters for other isotopic species are calculated. For HCP and DCP the more limited data available have been fitted to an anharmonic force field using constraints based on comparison with HCN. Using this force field the zero-point rotational constants B0 have been corrected to obtain the equilibrium constants Be, and hence the equilibrium structure has been determined to be re(CH) = 1•0692(7)A, and re(CP) = 1•5398(2)A.

Vibration-rotation spectra and the harmonic force field of HOCl

Vibration-rotation spectra of HOCl have been measured at a resolution of 0.05 cm−1 to determine vibration rotation constants, and 35–37 Cl isotope shifts in the vibration frequencies. The spectrum of DOCl has also been recorded, and a preliminary analysis for the band origins has been made. The vibrational frequency data and centrifugal distortion constants have been used to determine the harmonic force field in a least-squares refinement; the force field obtained also gives a good fit to data on the vibrational contributions to the inertial defect. The equilibrium rotational constants of HOCl have been obtained, and an equilibrium structure has been estimated.

An anharmonic force field for SO3

An anharmonic force field for SO3 based on the valence force model has been investigated. The results of extending the model to include some further estimated cubic interaction potential constants have also been investigated. The phenomenological parameters calculated from both model force fields agree with those few values which have been experimentally determined. A calculation of the inertia defect has been made, and thus the value of C0 has been determined. The equilibrium structure has been determined to be: re = 1.4184 ± 0.0010 Å.

Nitrous acid: vibrational frequency shifts due to 15N substitution, and the harmonic force field

Vibration rotation spectra of HO15 NO and DO15 NO have been measured at a resolution of 0•04 cm-1 to determine the isotopic shifts in the vibrational band origins. These have been used together with recently determined data on the vibrational band origins, Coriolis constants, and centrifugal distorition constants, to determine the harmonic force field of both cis and trans nitrous acid in least squares refinement calculations. The results are discussed in relation to recent ab initio calculations, the inertia defects, and the torsional potential function.

The harmonic force field and rz structure of HNCO

The presently available microwave, millimeter wave, and far-infrared data of five isotopic species of isocyanic acid, namely, HNCO, H15NCO, HN13CO, HNC18O, and DNCO, have been used to obtain improved values of the ground-state rotational constants, the five quartic distortion constants, and some higher-order distortion constants in the Ir S reduced Hamiltonian of Watson. The appropriate planarity relation among the quartic centrifugal distortion constants has been imposed in the fitting procedure. The general harmonic force field of isocyanic acid has been determined using all existing data, and assuming a trans bent equilibrium geometry of the molecule with an NCO angle of 170°. Finally an rz structure has been obtained using the Az, Bz, and Cz rotational constants of five isotopic species. The bending of the NCO chain is found to be 8° in the trans configuration.