CHARMM Force Field Parameters for Nitroalkanes and Nitroarenes

Interfacial Behavior of Cs+, K+, Na+, and Rb+ Extraction in the Presence of Dibenzo-18-Crown-6 from the Nitrobenzene-Water Biphasic System: Experimental, Quantum Chemical, and Molecular Dynamic Studies

Extraction of metal ions (i.e., Cs⁺, K⁺, Na⁺, and Rb⁺) in the presence of ionophore such as dibenzo-18-crown-6 (DB18C6) from the nitrobenzene-water biphasic system is reported by COSMO-RS (conductor-like screening model for real solvents) predictions, molecular dynamics simulation, along with experimental validation. The predicted values of selectivity as obtained for the Na⁺-DB18C6 complex were 4.571, 4.877, and 4.947 at 298.15, 308.15, and 318.15 K, respectively. This was then confirmed by the experimental distribution coefficient (D) as obtained in the diluent systems along with by varying the metal ion to crown ether ligand (M-L) mole ratios: 10:1 (0.1 M M⁺ and 0.01 M DB18C6), 1:1 (0.01 M M⁺ and 0.01 M DB18C6), and 1:10 (0.001 M M⁺ and 0.01 M DB18C6). The experimentally determined values of D _Na (i.e., 0.059, 0.060, and 0.056) were found to be very large as compared to the values of D _Cs (i.e., 0.001, 0.010, and 0.024) in the nitrobenzene phase. It indicates an excellent extraction ability of DB18C6 for Na⁺. The rate of phase separation for the Cs⁺NO₃ ^- system was slow as compared to other metal ion systems. The binding energies, free energies, and nonbonded interaction energies of the complexed metal ion in solution were calculated with both explicit and implicit solvent models. A higher interaction energy between Na⁺-DB18C6 complex and nitrobenzene was observed (i.e., -289.92 in the explicit model and -143.12 kcal/mol in the implicit model) when compared with other metal ions (i.e., Cs⁺, K⁺, and Rb⁺).

Modelling water diffusion in plasticizers: development and optimization of a force field for 2,4-dinitroethylbenzene and 2,4,6-trinitroethylbenzene

A classical all-atom force field has been developed for 2,4,6-trinitroethylbenzene and 2,4-dinitroethylbenzene and applied in molecular dynamics simulations of the two pure and two mixed plasticizer systems. Bonding parameters and partial charges were derived through electronic and geometry optimization of the single molecules. The other required parameters were derived from values already available in the literature for generic nitro aromatic compounds, which were adjusted to reproduce to a high level of accuracy the densities of 2,4-dinitroethylbenzene, 2,4,6-trinitroethylbenzene and the energetic plasticizers K10 and R8002. This force field has been applied to both K10 and R8002, which when used as plasticizers form an energetic binder with nitrocellulose. Nitrocellulose decomposes in storage, under varying conditions, but in particular where it may become increasingly dry. Following the derivation of the force field, we have therefore applied it to calculate water diffusion coefficients for each of the different materials at 298 K and 338 K, thereby providing a starting point for understanding water behaviour in a nitrocellulose binder.

A classical all-atom force field has been developed for the plasticizer molecules 2,4,6-trinitroethylbenzene and 2,4-dinitroethylbenzene which can be used to investigate properties and energetic output of nitrocellulose-based propellants and binders.

Biochemical characterization of recombinant Enterovirus 71 3C protease with fluorogenic model peptide substrates and development of a biochemical assay

Enterovirus 71 (EV71), a primary pathogen of hand, foot, and mouth disease (HFMD), affects primarily infants and children. Currently, there are no effective drugs against HFMD. EV71 3C protease performs multiple tasks in the viral replication, which makes it an ideal antiviral target. We synthesized a small set of fluorogenic model peptides derived from cleavage sites of EV71 polyprotein and examined their efficiencies of cleavage by EV71 3C protease. The novel peptide P08 [(2-(N-methylamino)benzoyl) (NMA)-IEALFQGPPK(DNP)FR] was determined to be the most efficiently cleaved by EV71 3C protease, with a kinetic constant kcat/Km of 11.8 ± 0.82 mM(-1) min(-1). Compared with literature reports, P08 gave significant improvement in the signal/background ratio, which makes it an attractive substrate for assay development. A Molecular dynamics simulation study elaborated the interactions between substrate P08 and EV71 3C protease. Arg39, which is located at the bottom of the S2 pocket of EV71 3C protease, may participate in the proteolysis process of substrates. With an aim to evaluate EV71 3C protease inhibitors, a reliable and robust biochemical assay with a Z' factor of 0.87 ± 0.05 was developed. A novel compound (compound 3) (50% inhibitory concentration [IC50] = 1.89 ± 0.25 μM) was discovered using this assay, which effectively suppressed the proliferation of EV 71 (strain Fuyang) in rhabdomyosarcoma (RD) cells with a highly selective index (50% effective concentration [EC50] = 4.54 ± 0.51 μM; 50% cytotoxic concentration [CC50] > 100 μM). This fast and efficient assay for lead discovery and optimization provides an ideal platform for anti-EV71 drug development targeting 3C protease.

The zero-multipole summation method for estimating electrostatic interactions in molecular dynamics: analysis of the accuracy and application to liquid systems

In the preceding paper [I. Fukuda, J. Chem. Phys. 139, 174107 (2013)], the zero-multipole (ZM) summation method was proposed for efficiently evaluating the electrostatic Coulombic interactions of a classical point charge system. The summation takes a simple pairwise form, but prevents the electrically non-neutral multipole states that may artificially be generated by a simple cutoff truncation, which often causes large energetic noises and significant artifacts. The purpose of this paper is to judge the ability of the ZM method by investigating the accuracy, parameter dependencies, and stability in applications to liquid systems. To conduct this, first, the energy-functional error was divided into three terms and each term was analyzed by a theoretical error-bound estimation. This estimation gave us a clear basis of the discussions on the numerical investigations. It also gave a new viewpoint between the excess energy error and the damping effect by the damping parameter. Second, with the aid of these analyses, the ZM method was evaluated based on molecular dynamics (MD) simulations of two fundamental liquid systems, a molten sodium-chlorine ion system and a pure water molecule system. In the ion system, the energy accuracy, compared with the Ewald summation, was better for a larger value of multipole moment l currently induced until l ≲ 3 on average. This accuracy improvement with increasing l is due to the enhancement of the excess-energy accuracy. However, this improvement is wholly effective in the total accuracy if the theoretical moment l is smaller than or equal to a system intrinsic moment L. The simulation results thus indicate L ∼ 3 in this system, and we observed less accuracy in l = 4. We demonstrated the origins of parameter dependencies appearing in the crossing behavior and the oscillations of the energy error curves. With raising the moment l we observed, smaller values of the damping parameter provided more accurate results and smoother behaviors with respect to cutoff length were obtained. These features can be explained, on the basis of the theoretical error analyses, such that the excess energy accuracy is improved with increasing l and that the total accuracy improvement within l ⩽ L is facilitated by a small damping parameter. Although the accuracy was fundamentally similar to the ion system, the bulk water system exhibited distinguishable quantitative behaviors. A smaller damping parameter was effective in all the practical cutoff distance, and this fact can be interpreted by the reduction of the excess subset. A lower moment was advantageous in the energy accuracy, where l = 1 was slightly superior to l = 2 in this system. However, the method with l = 2 (viz., the zero-quadrupole sum) gave accurate results for the radial distribution function. We confirmed the stability in the numerical integration for MD simulations employing the ZM scheme. This result is supported by the sufficient smoothness of the energy function. Along with the smoothness, the pairwise feature and the allowance of the atom-based cutoff mode on the energy formula lead to the exact zero total-force, ensuring the total-momentum conservations for typical MD equations of motion.

Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types

A significant modification to the additive all-atom CHARMM lipid force field (FF) is developed and applied to phospholipid bilayers with both choline and ethanolamine containing head groups and with both saturated and unsaturated aliphatic chains. Motivated by the current CHARMM lipid FF (C27 and C27r) systematically yielding values of the surface area per lipid that are smaller than experimental estimates and gel-like structures of bilayers well above the gel transition temperature, selected torsional, Lennard-Jones and partial atomic charge parameters were modified by targeting both quantum mechanical (QM) and experimental data. QM calculations ranging from high-level ab initio calculations on small molecules to semiempirical QM studies on a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer in combination with experimental thermodynamic data were used as target data for parameter optimization. These changes were tested with simulations of pure bilayers at high hydration of the following six lipids: DPPC, 1,2-dimyristoyl-sn-phosphatidylcholine (DMPC), 1,2-dilauroyl-sn-phosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-phosphatidylcholine (POPC), 1,2-dioleoyl-sn-phosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-sn-phosphatidylethanolamine (POPE); simulations of a low hydration DOPC bilayer were also performed. Agreement with experimental surface area is on average within 2%, and the density profiles agree well with neutron and X-ray diffraction experiments. NMR deuterium order parameters (S(CD)) are well predicted with the new FF, including proper splitting of the S(CD) for the aliphatic carbon adjacent to the carbonyl for DPPC, POPE, and POPC bilayers. The area compressibility modulus and frequency dependence of (13)C NMR relaxation rates of DPPC and the water distribution of low hydration DOPC bilayers also agree well with experiment. Accordingly, the presented lipid FF, referred to as C36, allows for molecular dynamics simulations to be run in the tensionless ensemble (NPT), and is anticipated to be of utility for simulations of pure lipid systems as well as heterogeneous systems including membrane proteins.