Rapid parameterization of small molecules using the force field toolkit

The inability to rapidly generate accurate and robust parameters for novel chemical matter continues to severely limit the application of molecular dynamics simulations to many biological systems of interest, especially in fields such as drug discovery. Although the release of generalized versions of common classical force fields, for example, General Amber Force Field and CHARMM General Force Field, have posited guidelines for parameterization of small molecules, many technical challenges remain that have hampered their wide‐scale extension. The Force Field Toolkit (ffTK), described herein, minimizes common barriers to ligand parameterization through algorithm and method development, automation of tedious and error‐prone tasks, and graphical user interface design. Distributed as a VMD plugin, ffTK facilitates the traversal of a clear and organized workflow resulting in a complete set of CHARMM‐compatible parameters. A variety of tools are provided to generate quantum mechanical target data, setup multidimensional optimization routines, and analyze parameter performance. Parameters developed for a small test set of molecules using ffTK were comparable to existing CGenFF parameters in their ability to reproduce experimentally measured values for pure‐solvent properties (<15% error from experiment) and free energy of solvation (±0.5 kcal/mol from experiment). © 2013 Wiley Periodicals, Inc.     

Non-isothermal pyrolysis of used lubricating oil and the catalytic effect of carbon-based nanomaterials on the process performance

Pyrolysis is a commonly used method for the recovery of used lubricating oil (ULO), which should be kinetically improved by a catalyst, due to its high level of energy consumption. In this research, the catalytic effects of carbon nanotube (CNT) and graphene nanoplatelets on the pyrolysis of ULO were studied through thermogravimetric analysis. First, the kinetic parameters of ULO pyrolysis including activation energy were calculated to be 170.12 and 167.01 kJ mol^?1 by FWO and KAS methods, respectively. Then, the catalytic effects of CNT and graphene nanoplatelets on pyrolysis kinetics were studied. While CNT had a negligible effect on the pyrolysis process, graphene nanoplatelets significantly reduced the temperature of maximum conversion during pyrolysis from 400 to 350 °C, due to high thermal conductivity and homogenous heat transfer in the pyrolysis process. On the other hand, graphene nanoplatelets maximized the rate of conversion of highly volatile components at lower temperatures (<?100 °C), which was mainly due to the high affinity of these components toward graphene nanoplatelets and also the effect of nanoplatelets’ edges which have free tails and can bond with other molecules. Moreover, graphene nanoplatelets decreased the activation energy of the conversion to 154.48 and 152.13 kJ mol^?1 by FWO and KAS methods, respectively.

     

Classical force field parameters for the heme prosthetic group of cytochrome c

Accurate force fields are essential for describing biological systems in a molecular dynamics simulation. To analyze the docking of the small redox protein cytochrome c (cyt c) requires simulation parameters for the heme in both the reduced and oxidized states. This work presents parameters for the partial charges and geometries for the heme in both redox states with ligands appropriate to cyt c. The parameters are based on both protein X-ray structures and ab initio density functional theory (DFT) geometry optimizations at the B3LYP/6-31G* level. The simulations with the new parameter set reproduce the geometries of the X-ray structures and the interaction energies between water and heme prosthetic group obtained from B3LYP/6-31G* calculations. The parameter set developed here will provide new insights into docking processes of heme containing redox proteins.     

Rationalizing the generation of broad spectrum antibiotics with the addition of a positive charge

Antibiotic resistance of Gram-negative bacteria is largely attributed to the low permeability of their outer membrane (OM). Recently, we disclosed the eNTRy rules, a key lesson of which is that the introduction of a primary amine enhances OM permeation in certain contexts. To understand the molecular basis for this finding, we perform an extensive set of molecular dynamics (MD) simulations and free energy calculations comparing the permeation of aminated and amine-free antibiotic derivatives through the most abundant OM porin of E. coli, OmpF. To improve sampling of conformationally flexible drugs in MD simulations, we developed a novel, Monte Carlo and graph theory based algorithm to probe more efficiently the rotational and translational degrees of freedom visited during the permeation of the antibiotic molecule through OmpF. The resulting pathways were then used for free-energy calculations, revealing a lower barrier against the permeation of the aminated compound, substantiating its greater OM permeability. Further analysis revealed that the amine facilitates permeation by enabling the antibiotic to align its dipole to the luminal electric field of the porin and form favorable electrostatic interactions with specific, highly-conserved charged residues. The importance of these interactions in permeation was further validated with experimental mutagenesis and whole cell accumulation assays. Overall, this study provides insights on the importance of the primary amine for antibiotic permeation into Gram-negative pathogens that could help the design of future antibiotics. We also offer a new computational approach for calculating free-energy of processes where relevant molecular conformations cannot be efficiently captured.

A rapid pathway sampling method combining Monte Carlo and graph theory, developed to describe permeation pathways through outer membrane porins, can distinguish between structurally similar analogs with different permeabilities.     

Interplay of mechanical and binding properties of Fibronectin type I

Fibronectins (FNs) are a major component of the extracellular matrix (ECM), and provide important binding sites for a variety of ligands outside and on the surface of the cell. Similar to other ECM proteins, FNs are consistently subject to mechanical stress in the ECM. Therefore, it is important to study their structure and binding properties under mechanical stress and understand how their binding and mechanical properties might affect each other. Although certain FN modules have been extensively investigated, no simulation studies have been reported for the FN type I (Fn1) domains, despite their prominent role in binding of various protein modules to FN polymers in the ECM. Using equilibrium and steered molecular dynamics simulations, we have studied mechanical properties of Fn1 modules in the presence or the absence of a specific FN-binding peptide (FnBP). We have also investigated how the binding of the FnBP peptide to Fn1 might be affected by tensile force. Despite the presence of disulfide bonds within individual Fn1 modules that are presumed to prevent their extension, it is found that significant internal structural changes within individual modules are induced by the forces applied in our simulations. These internal structural changes result in significant variations in the accessibility of different residues of the Fn1 modules, which affect their exposure, and, thus, the binding properties of the Fn1 modules. Binding of the FnBP appears to reduce the flexibility of the linker region connecting individual Fn1 modules (exhibited in the form of reduced fluctuation and motion of the linker region), both with regard to bending and stretching motions, and hence stabilizes the inter-domain configuration under force. Under large tensile forces, the FnBP peptide unbinds from Fn1. The results suggest that Fn1 modules in FN polymers do contribute to the overall extension caused by force-induced stretching of the polymer in the ECM, and that binding properties of Fn1 modules can be affected by mechanically induced internal protein conformational changes in spite of the presence of disulfide bonds which were presumed to completely abolish the capacity of Fn1 modules to undergo extension in response to external forces.     

The effect of the protein environment on the structure and charge distribution of the retinal Schiff base in bacteriorhodopsin

The effects of the amino acid side chains of the binding pocket of bacteriorhodopsin (bR) and of a water molecule on the structure of the retinal Schiff base have been studied using Becke3LYP/6-31G* level of density functional theory. A model protonated Schiff base structure including six conjugated double bonds and methyl substituents was optimized in the presence of several amino acid side chains and of a water molecule, separately. The Schiff base structure was also calculated in the form of a neutral species. At each optimized complex geometry the atomic charges of the model Schiff base were calculated using Mulliken population analysis. In agreement with previously proposed counterion(s) of the protonated retinal Schiff base in bR, the results show that Asp₈₅ and Asp₂₁₂, which are present in the form of negatively charged groups, have significantly large effects on the structure and electronic configuration of both unprotonated and protonated model Schiff bases. The presence of a water molecule in the vicinity of the Schiff base demonstrates significant effects which are comparable to those of aspartate groups. Other side chains studied did not show any significant effect in this direction. Apart from the aspartate groups and the water molecule, in none of the other complexes studied are the atomic charges and the bond alternation of the model Schiff base significantly influenced by the presence of the neighboring amino acids. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 10 December 1998     

Resolution enhancement in solid-state NMR of oriented membrane proteins by anisotropic differential linebroadening

We demonstrate that a significant improvement in the spectral resolution may be achieved in solid-state NMR experiments of proteins in inhomogeneously disordered oriented lipid bilayers. Using 1H homonuclear decoupling instead of standard 1H heteronuclear decoupling, the 15N line widths may be reduced by up to seven times for such samples. For large oriented membrane proteins, such resolution enhancements may be crucial for assignment and structural interpretation.     

Toward theoretical analysis of long-range proton transfer kinetics in biomolecular pumps

Motivated by the long-term goal of theoretically analyzing long-range proton transfer (PT) kinetics in biomolecular pumps, researchers made a number of technical developments in the framework of quantum mechanics-molecular mechanics (QM/MM) simulations. A set of collective reaction coordinates is proposed for characterizing the progress of long-range proton transfers; unlike previous suggestions, the new coordinates can describe PT along highly nonlinear three-dimensional pathways. Calculations using a realistic model of carbonic anhydrase demonstrated that adiabatic mapping using these collective coordinates gives reliable energetics and critical geometrical parameters as compared to minimum energy path calculations, which suggests that the new coordinates can be effectively used as reaction coordinate in potential of mean force calculations for long-range PT in complex systems. In addition, the generalized solvent boundary potential was implemented in the QM/MM framework for rectangular geometries, which is useful for studying reactions in membrane systems. The resulting protocol was found to produce water structure in the interior of aquaporin consistent with previous studies including a much larger number of explicit solvent and lipid molecules. The effect of electrostatics for PT through a membrane protein was also illustrated with a simple model channel embedded in different dielectric continuum environments. The encouraging results observed so far suggest that robust theoretical analysis of long-range PT kinetics in biomolecular pumps can soon be realized in a QM/MM framework.