A simplified homology‐model builder toward highly protein‐like structures: An inspection of restraining potentials

A homology model builder using simple restraining potentials based on spline‐interpolated quadratic functions is developed and interfaced with CHARMM package. The continuity and stability of the potential function were validated, and the parameters were optimized using the CASP7 targets. The performance of the model builder was benchmarked to the Modeller program using the template‐based modeling targets in CASP9. The benchmark results show that, while our builder yields the structures with slightly lower packing, backbone, and template modeling scores, our models show much better protein‐like scores in terms of normalized discrete optimized protein energy, dipolar distance‐scaled finite‐ideal gas reference, Molprobity clash, Ramachandran appearance Z‐score, and rotamer Z‐score. As our model builder is interfaced with CHARMM, it is advantageous to directly use other CHARMM functionality and energy functions to refine the model structures or to use the models for other computational studies using CHARMM. © 2012 Wiley Periodicals, Inc.     

CHARMM‐GUI ligand reader and modeler for CHARMM force field generation of small molecules

Reading ligand structures into any simulation program is often nontrivial and time consuming, especially when the force field parameters and/or structure files of the corresponding molecules are not available. To address this problem, we have developed Ligand Reader & Modeler in CHARMM‐GUI. Users can upload ligand structure information in various forms (using PDB ID, ligand ID, SMILES, MOL/MOL2/SDF file, or PDB/mmCIF file), and the uploaded structure is displayed on a sketchpad for verification and further modification. Based on the displayed structure, Ligand Reader & Modeler generates the ligand force field parameters and necessary structure files by searching for the ligand in the CHARMM force field library or using the CHARMM general force field (CGenFF). In addition, users can define chemical substitution sites and draw substituents in each site on the sketchpad to generate a set of combinatorial structure files and corresponding force field parameters for throughput or alchemical free energy simulations. Finally, the output from Ligand Reader & Modeler can be used in other CHARMM‐GUI modules to build a protein‐ligand simulation system for all supported simulation programs, such as CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Ligand Reader & Modeler is available as a functional module of CHARMM‐GUI at http://www.charmm-gui.org/input/ligandrm . © 2017 Wiley Periodicals, Inc.     

Producing DFT/MM enzyme reaction trajectories from SCC‐DFTB/MM driving forces to probe the underlying electronics of a glycosyltransferase reaction

The SCC‐DFTB/MIO/CHARMM free energy surface for a glycosyltransferase, TcTS, is benchmarked against a DFT/MM reaction trajectory using the same CHARMM MM force field ported to the NWChem package. The popular B3LYP functional, against which the MIO parameter set was parameterized is used to optimize TS structures and run DFT reaction dynamics. A novel approach was used to generate reaction forces from a SCC‐DFTB/MIO/CHARMM reaction surface to drive B3LYP/6‐31G/MM and B3LYP/6‐31G(d)/MM reaction trajectories. Although TS structures compare favorably, differences stemming primarily from a minimal basis set approximation prevented a successful 6‐31G(d) FEARCF reaction dynamics trajectory. None the less, the dynamic evolution of the B3LYP/6‐31G/MM‐computed electron density provided an opportunity to perform NBO analysis along the reaction trajectory. Here, we illustrate that a successful ab initio reaction trajectory is computationally accessible when the underlying potential energy function of the semi‐empirical method used to produce driving forces is sufficiently close to the ab initio potential. © 2017 Wiley Periodicals, Inc.     

Performances of Five Representative Force Fields on Gaseous Amino Acids with Di erent Termini

There is a growing interest in the study of structures and properties of biomolecules in gas phase. Applications of force fields are highly desirable for the computational efficiency of the gas phase study. To help the selection of force fields, the performances of five repre-sentative force fields for gaseous neutral, protonated, deprotonated and capped amino acids are systematically examined and compared. The tested properties include relative conforma-tional energies, energy differences between cis and trans structures, the number and strength of predicted hydrogen bonds, and the quality of the optimized structures. The results of BHandHLYP/6-311++G(d,p) are used as the references. GROMOS53A6 and ENCADS are found to perform poorly for gaseous biomolecules, while the performance of AMBER99SB, CHARMM27 and OPLSAA/L are comparable when applicable. Considering the general availability of the force field parameters, CHARMM27 is the most recommended, followed by OPLSAA/L, for the study of biomolecules in gas phase.     

Computational Insights into the Catalytic Mechanism of Is-PETase: An Enzyme Capable of Degrading Poly(ethylene) Terephthalate

Is-PETase has become an enzyme of significant interest due to its ability to catalyse the degradation of polyethylene terephthalate (PET) at mesophilic temperatures. We performed hybrid quantum mechanics and molecular mechanics (QM/MM) at the DSD-PBEP86-D3/ma-def2-TZVP/CHARMM27//rev-PBE-D3/dev2-SVP/CHARMM level to calculate the energy profile for the degradation of a suitable PET model by this enzyme. Very low overall barriers are computed for serine protease-type hydrolysis steps (as low as 34.1 kJ mol⁻¹). Spontaneous deprotonation of the final product, terephthalic acid, with a high computed driving force indicates that product release could be rate limiting.     

Optimization of CHARMM force field parameters for the chalcone fragment

In this work,we developed the CHARMM all-atom force field parameters for the nonstandard biological residue chalcone,followed by the standard protocol for the CHARMM27 force field development.Target data were generated via ab initio calculations at the MP2/6-31G* and HF/6-31G* levels.The reference data included interaction energies between water and the model compound F(a fragment of chalcone).Bond,angle,and torsion parameters were derived from the ab initio calculations and renormalized to maintain compatibility with the existing CHARMM27 parameters of standard residues.The optimized CHARMM parameters perform well in reproducing the target data.We expect that the extension of the CHARMM27 force field parameters for chalcone will facilitate the molecular simulation studies of the reaction mechanism of intramolecular cyclization of chalcone catalyzed by chalcone isomerase.     

Thermoelasticity of stretched elastomers under small torsions: Effects of crosslinking,extension, and swelling

By using a simple version of Treloar's torsion method the internal energy contribution M_e to the restoring couple M has been determined for stretched testpieces subjected to small torsions. Values obtained for M_e/M decreased significantly with increasing extension. Corresponding decreases were observed on swelling, but no specific effect of swelling liquids could be detected. The degree of crosslinking had no effect on M_e/M. Thus, M_e/M at small strains appears to reflect a characteristic molecular property, largely independent of local environment, network density, or type of deformation, as Flory has maintained. The present values are generally consistent with previously published data, except in the case of polyethylene, for which a much larger negative value is obtained, i.e., ?1.6 compared to ?0.5. The reason for this discrepancy is not known.     

The compressive electrical field electrostrictive coefficient M33 of electroactive polymer composites and its saturation versus electrical field,polymer thickness,frequency, and fillers

Electroactive polymers are widely studied because of their large electrical‐field‐induced strain. Their flexibility and their ability to be deposited on large surfaces make them promising candidates as electroactive materials for actuators or energy‐harvesting devices. For actuation purposes, the material efficiency is directly related to the electrical‐field‐related electrostrictive coefficient M₃₃ through S₃₃ = M₃₃E², where S₃₃ is the electrical‐field‐induced strain and E is the applied electrical field. Numerous studies concern the increase of M₃₃, but very few have been devoted to its saturations versus electrical field. To this end, the present paper describes the variation of M₃₃ versus thickness, composition, frequency, and electrical field for polyurethane‐based composites. Based on the saturation of the electrical‐field‐induced polarization within the studied polymer composites, a model of the M₃₃ behavior was also proposed, and it was found to show a good agreement with the experimental data. In addition, this model predicts the dielectric constant and the saturation electrical field to be the key parameters ruling the M₃₃ saturation. Copyright © 2011 John Wiley & Sons, Ltd.     

The CHARMM–TURBOMOLE interface for efficient and accurate QM/MM molecular dynamics,free energies,and excited state properties

The quantum mechanical (QM)/molecular mechanical (MM) interface between Chemistry at HARvard Molecular Mechanics (CHARMM) and TURBOMOLE is described. CHARMM provides an extensive set of simulation algorithms, like molecular dynamics (MD) and free energy perturbation, and support for mature nonpolarizable and Drude polarizable force fields. TURBOMOLE provides fast QM calculations using density functional theory or wave function methods and excited state properties. CHARMM–TURBOMOLE is well‐suited for extended QM/MM MD simulations using first principles methods with large (triple‐ζ) basis sets. We demonstrate these capabilities with a QM/MM simulation of Mg²⁺(aq), where the MM outer sphere water molecules are represented using the SWM4‐NDP Drude polarizable force field and the ion and inner coordination sphere are represented using QM PBE, PBE0, and MP2 methods. The relative solvation free energies of Mg²⁺ and Zn²⁺ were calculated using thermodynamic integration. We also demonstrate the features for excited state properties. We calculate the time‐averaged solution absorption spectrum of indole, the emission spectrum of the indole excited state, and the electronic circular dichroism spectrum of an oxacepham. © 2014 Wiley Periodicals, Inc.     

QuanPol: A full spectrum and seamless QM/MM program

The quantum chemistry polarizable force field program (QuanPol) is implemented to perform combined quantum mechanical and molecular mechanical (QM/MM) calculations with induced dipole polarizable force fields and induced surface charge continuum solvation models. The QM methods include Hartree–Fock method, density functional theory method (DFT), generalized valence bond theory method, multiconfiguration self‐consistent field method, Møller–Plesset perturbation theory method, and time‐dependent DFT method. The induced dipoles of the MM atoms and the induced surface charges of the continuum solvation model are self‐consistently and variationally determined together with the QM wavefunction. The MM force field methods can be user specified, or a standard force field such as MMFF94, Chemistry at Harvard Molecular Mechanics (CHARMM), Assisted Model Building with Energy Refinement (AMBER), and Optimized Potentials for Liquid Simulations‐All Atom (OPLS‐AA). Analytic gradients for all of these methods are implemented so geometry optimization and molecular dynamics (MD) simulation can be performed. MD free energy perturbation and umbrella sampling methods are also implemented. © 2013 Wiley Periodicals, Inc.     

A truncated Newton minimizer adapted for CHARMM and biomolecular applications

We report the adaptation of the truncated Newton minimization package TNPACK for CHARMM and biomolecular energy minimization. TNPACK is based on the preconditioned linear conjugate–gradient technique for solving the Newton equations. The structure of the problem—sparsity of the Hessian—is exploited for preconditioning. Experience with the new version of TNPACK is presented on a series of molecular systems of biological and numerical interest: alanine dipeptide (N-methyl-alanyl-acetamide), a dimer of N-methyl-acetamide, deca-alanine, mellitin (26 residues), avian pancreatic polypeptide (36 residues), rubredoxin (52 residues), bovine pancreatic trypsin inhibitor (58 residues), a dimer of insulin (99 residues), and lysozyme (130 residues). Detailed comparisons among the minimization algorithms available in CHARMM, particularly those used for large-scale problems, are presented along with new mathematical developments in TNPACK. The new TNPACK version performs significantly better than ABNR, the most competitive minimizer in CHARMM, for all systems tested in terms of CPU time when curvature information (Hessian/vector product) is calculated by a finite-difference of gradients (the numeric option of TNPACK). The remaining derivative quantities are, however, evaluated analytically in TNPACK. The CPU gain is 50% or more (speedup factors of 1.5 to 2.5) for the largest molecular systems tested and even greater for smaller systems (CPU factors of 1 to 4 for small systems and 1 to 5 for medium systems). TNPACK uses curvature information to escape from undesired configurational regions and to ensure the identification of true local minima. It converges rapidly once a convex region is reached and achieves very low final gradient norms, such as of order 10^?8, with little additional work. Even greater overall CPU gains are expected for large-scale minimization problems by making the architectures of CHARMM and TNPACK more compatible with respect to the second-derivative calculations. © 1994 by John Wiley & Sons, Inc.     

On the accuracy of the many-points local procedures

As until now proposed in the literature, in the local energy calculations we can distinguish “few-points” procedures, in which the number M of configurational points is strictly related to the number N of trial functions used, and statistical “many-points” procedures, in which the number M of points can be arbitrarily increased. In this paper we demonstrate that the energy errors resulting from a “many-points” calculation M points/N functions (M > N) can be connected in a simple way with the errors of the (^M_N) partial calculations N points/N functions. This suggests a possible approach for the problem of the choice of the configurational points to be introduced in the calculation, and leads to a simple interpretation of the numerical meaning of the error associated with the ordinary Ritz energy. Numerical examples on the hydrogen atom are reported.     

On the use of conformationally dependent geometry trends from ab initio dipeptide studies to refine potentials for the empirical force field CHARMM

Previous 4-21G ab initio geometry optimizations of various conformations of the model dipeptides (N-acetyl N'methyl amides) of glycine (GLY) and the alanine (ALA) have been used to help refine the empirical force constants and equilibrium geometry in the CHARMM force field for peptides. Conformationally dependent geometry trends from ab initio calculations and positions of energy minima on the ab initio energy surfaces have been used as guides in the parameter refinement, leading to modifications in the bond stretch, angle bending, and some torsional parameters. Preliminary results obtained with these refined empirical parameters are presented for the protein Crambin. Results for the cyclic (Ala-Pro-DPhe)₂ are compared with those from other calculations. It seems that the dihedral angle fit achieved by the new parameters is significantly improved compared with results from force fields whose derivation does not include ab initio geometry trends.     

Thermoplastic resins,polyquats, and membranes based on epoxies

High-molecular-weight branched-addition polymers made from oligomeric solid epoxy resins and amines carrying two active hydrogen atoms are described. The reaction with aliphatic amines like ethanolamine, 2-ethylhexylamine, or piperazine can easily be controlled carrying out the polymerization in about 15% solutions of boiling 2-methoxyethanol. Less basic derivatives of aniline and epoxy resins are processed at 200°C in the molten state. Polymers with M _n up to 33,000 and M _w to 720,000 are obtained. Properties of films and injection moldings are described. The tertiary nitrogen atoms present in the polymer chain can be quaternized with alkylating agents, like dimethylsulfate, methyliodide, or trimethyl phosphate. Thin membranes cast from the polyquats thus obtained exhibit a good flow rate for water and a high salt rejection under reverse osmosis conditions.     

A Computational Study on Closed-Shell Molecular Hexafluorides MF6 (M=S,Se, Te,Po, Xe,Rn, Cr,Mo, W,U) – Molecular Structure,Anharmonic Frequency Calculations,and Prediction of the NdF6 Molecule

Quantum chemical methods were used to study the molecular structure and anharmonic IR spectra of the experimentally known closed-shell molecular hexafluorides MF₆ (M=S, Se, Te, Xe, Mo, W, U). First, the molecular structures and harmonic frequencies were investigated using Density Functional Theory (DFT) with all-electron basis sets and explicitly considering the influence of spin-orbit coupling. Second, anharmonic frequencies and IR intensities were calculated with the CCSD(T) coupled cluster method and compared, where available, with IR spectra recorded by us. These comparisons showed satisfactory results. The anharmonic IR spectra provide means for identifying experimentally too little studied or unknown MF₆ molecules with M=Cr, Po, Rn. To the best of our knowledge, we predict the NdF₆ molecule for the first time and show it to be a true local minimum on the potential energy surface. We used intrinsic bond orbital (IBO) analyses to characterize the bonding situation in comparison with the UF₆ molecule.     

CHARMM‐GUI Membrane Builder toward realistic biological membrane simulations

CHARMM‐GUI Membrane Builder, http://www.charmm‐gui.org/input/membrane , is a web‐based user interface designed to interactively build all‐atom protein/membrane or membrane‐only systems for molecular dynamics simulations through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builder that allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types, such as phosphoinositides, cardiolipin (CL), sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance‐based algorithm for faster initial ion displacement, (5) CHARMM inputs for P₂₁ image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: CL lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM‐GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments. © 2014 Wiley Periodicals, Inc.     

A QM-MM interface between CHARMM and TURBOMOLE: implementation and application to systems in bulk phase and biologically active systems

The implementation of a hybrid QM-MM approach combining ab initio and density functional methods of TURBOMOLE with the molecular mechanics program package CHARMM is described. An interface has been created to allow data exchange between the two applications. With this method the efficient multiprocessor capabilities of TURBOMOLE can be utilized with CHARMM running as a single processor application. Therefore, features of nonparallel running code in CHARMM like the TRAVEL module for locating saddle points or VIBRAN for the calculation of second derivatives can be exploited by running the CPU intensive QM calculations in parallel. To test the methodology, several small systems are studied with both Hartree-Fock and density functional methods and varying QM-MM boundaries. Also, the computationally efficient RI-J method has been examined for use in QM-MM applications. A B(12) cofactor containing cobalt has been studied, to examine systems with a large QM region and transition metals. All tested methods perform satisfactory in comparison with pure quantum calculations. Additionally, algorithms for the characterization of saddle points have been tested for their potential use in QM-MM problems. The TRAVEL module of CHARMM has been applied to the Menshutkin reaction in the condensed phase, and a saddle point was located. This saddle point was verified by calculation of a steepest descent path connecting educt, transition state, and product, and by calculation of vibrational modes.     

Density-functional calculations of MnC(M = Fe, Co, Ni, Cu, n = 1–6) clusters

The structural, energetic and magnetic properties of M_nC(M = Fe, Co, Ni, Cu, n = 1–6) clusters are systematically investigated by density-functional calculations. We found that the ground-state geometrical structures of M_nC clusters are different from those of pure M_n+1 clusters. Fe₄C, Ni₂C and Cu₄C possess relatively higher stabilities. Doping of a C atom enhances the binding energy of M_n clusters, and the binding energies of Fe_n-C, Co_n-C and Ni_n-C are stronger than that of Cu_n-C.     

CHARMMing: a new, flexible web portal for CHARMM

A new web portal for the CHARMM macromolecular modeling package, CHARMMing (CHARMM interface and graphics, http://www.charmming.org), is presented. This tool provides a user-friendly interface for the preparation, submission, monitoring, and visualization of molecular simulations (i.e., energy minimization, solvation, and dynamics). The infrastructure used to implement the web application is described. Two additional programs have been developed and integrated with CHARMMing: GENRTF, which is employed to define structural features not supported by the standard CHARMM force field, and a job broker, which is used to provide a portable method for using grid and cluster computing with CHARMMing. The use of the program is described with three proteins: 1YJP , 1O1O , and 1UFY . Source code is provided allowing CHARMMing to be downloaded, installed, and used by supercomputing centers and research groups that have a CHARMM license. Although no software can replace a scientist's own judgment and experience, CHARMMing eases the introduction of newcomers to the molecular modeling discipline by providing a graphical method for running simulations.     

Fortran IV programs for the penultimate copolymerization model

Three programs have been written for calculations involving use of the penultimate copolymerization model. The first computes the penultimate reactivity ratios from composition-conversion data, without constraints, at any conversion. A nonlinear leastsquares technique using Marquardt's algorithm is employed. The second program computes the four optimum starting monomer feed ratios, M₁⁰/M₂⁰ which should be used by the experimenter from the penultimate reactivity ratios. These optimum feed ratios are obtained by choosing the conditions necessary to minimize the determinant of the variance-covariance matrix. The input for the first program includes estimates of known values of the penultimate reactivity ratios. By using these two programs sequentially the experimenter has an optimized experimental approach toward evaluating penultimate reactivity ratios at any conversion. Finally, a program has been provided to calculate composition–conversion data, given penultimate reactivity ratios.