Search

The RNA-dependent RNA polymerase (RdRp) of norovirus is an attractive target of antiviral agents aimed at providing protection against norovirusassociated gastroenteritis. Here, we perform molecular dynamics simulations of the crystal structure of norovirus RdRp in complex with several known binders, as well as free-energy simulations by free-energy perturbation (FEP) to determine binding free energies of these molecules relative to the natural nucleotide substrates. We determine experimental EC(50) values and nucleotide incorporation efficiencies for several of these compounds. Moreover, we investigate the mechanism of inhibition of some of these ligands. Using FEP, we screened a virtual nucleotide library with 121 elements for binding to the polymerase and successfully identified two novel chain terminators.

The coupled reference interaction site model-hyper-netted chain (RISM-HNC)/ simulation methodology determines solvation free energies as a function of the set of all radial distribution functions of solvent atoms about atomic solute sites. These functions are determined from molecular dynamics (MD) or Monte Carlo (MC) simulations rather than from solving the RISM and HNC equations iteratively. Previous applications of the method showed that it can predict relative free energies of solvation for small solutes accurately. However, the errors scale with the system size. In this study, we propose the use of the hard-sphere free energy as the reference and a linear response approximation to improve the performance, i.e., accuracy and robustness, of the method, particularly removing the size dependency of the error. The details of the new formalism are presented. To validate the proposed formalism, solvation free energies of N-methylacetamide and methylamine are computed using the new RISM-HNC-based expressions in addition to a linear response expression, which are compared to previous thermodynamic integration and thermodynamic perturbation results performed with the same force field. Additionally, free energies of solvation for cyclohexane, pyridine, benzene and derivatives, and other small organic molecules are calculated and compared to experimental values.

Energy released by the hydrolysis of the high-energy phosphate bond of nucleoside triphosphate (NTP) cofactors is the driving force behind most biological processes. To understand how this energy is used to induce differences in protein structure and function, we examine the transfer of vibrational energy into the nucleotide-bound actin active site immediately after reaction activation. To this end, we perform Born-Oppenheimer molecular dynamics simulations of the active site at the level of density functional theory (DFT) starting at the calculated transition state (TS) structure. Similarly to the mechanism determined in many nucleotide-bound protein systems, the Os-Pγ bond is first elongated. Then, nucleophilic attack of the lytic water on Pγ occurs. Subsequently, protons are transferred in a cycle formed by water molecules, a protein residue, Asp154, and the γ-phosphate group, resulting in the formation of H2PO4(-). To investigate the possible creation of excited vibrational states in the products, power spectra of bond-length autocorrelation functions for relevant bonds within the active site are compared for simulations that start at the TS, at reactants, and at reaction end products. The hydroxyl bond formed in the final proton transfer to the phosphate molecule is observed to exhibit relatively high kinetic energies and large oscillations during reaction. It is also likely that some of the energy released by the reaction is captured by the low-energy stretching vibrations of the phosphoryl bonds of orthophosphate, which oscillate with large amplitudes in nonequilibrium simulations of end products.

An application of the coupled reference interaction site model (RISM)/simulation methodology to the calculation of the potential of mean force (PMF) curve in aqueous solution for the identity nucleophilic substitution reaction Cl(-) + CH(3)Cl is performed. The free energy of activation is calculated to be 27.1 kcal/mol which compares very well with the experimentally determined barrier height of 26.6 kcal/mol. Furthermore, the calculated PMF is almost superimposed with that previously calculated using the computationally rigorous Monte Carlo with importance sampling method (Chandrasekhar, J.; Smith, S. F.; Jorgensen, W. L. J. Am. Chem. Soc. 1985, 107, 154). Using the calculated PMF, a crude estimate of the solvated kinetic transmission coefficient also compares well with that of previous more accurate simulations. These results indicate that the coupled RISM/simulation method provides a cost-effective methodology for studying reactions in solution.

We present an application of the recently proposed coupled reference interation site model (RISM)/molecular dynamics (MD) solvation free-energy method (Freedman, H.; Truong, T. N. Chem. Phys. Lett. 2003, 381, 362; J. Chem. Phys. 2004, 121, 2187) to study the hydration effect upon the tautomeric equilibria of 2-hydroxypyridine and 2-oxopyridine and of cytosine in water. In this methodology, simulation trajectories of solvated systems are computed and averaged to determine radial distribution functions, which are then inserted into a RISM expression for the solvation free energy. The computed differences in free energy of hydration for the two tautomerizations agree well with the experimental data as well as results from previous free-energy simulations. This is particularly encouraging since the coupled RISM/MD method requires a single MD simulation as compared to the more complicated and computationally demanding free-energy simulation methods.

We present a new methodology for computing solvation free energy, which is based upon the reference interaction site model (RISM)/hypernetted chain (HNC) solvation free energy expression, but which substitutes radial distribution functions taken from simulations for those calculated by simultaneous solution of the RISM and HNC equations. Consequently, solvation free energy can be obtained from a single molecular dynamics or Monte Carlo simulation. Here we describe in detail the coupled RISM/simulation approach, and offer some error analysis. Finally we give the results of its application to a set of small test molecules in aqueous solution. The success shown in some of our results demonstrates that the coupled RISM/simulation approach is worth considering further as a potentially useful tool in studies of solvated systems, such as aqueous molecular biosystems.

We propose a new computational methodology that couples reference interaction site model (RISM) and molecular dynamics (MD) or Monte Carlo (MC) simulation methods for determination of solvation free energies. We employ the RISM formulation of solvation free energy. The correlation functions entering this expression are derived from radial distribution functions supplied by MD or MC simulations, instead of coming from simultaneous solution of the RISM and hyper-netted chain equations. We apply this approach to determining free energies of solvation for several small molecules. Our results are in good agreement with experimental values, demonstrating the potential of this method for applications to larger systems.

The cytoskeleton is essential to cell morphology, cargo trafficking, and cell division. As the neuronal cytoskeleton is extremely complex, it is no wonder that a startling number of neurodegenerative disorders (including but not limited to Alzheimer's disease, Parkinson's disease and Huntington's disease) share the common feature of a dysfunctional neuronal cytoskeleton. Recently, concern has been raised about a possible link between anesthesia, post-operative cognitive dysfunction, and the exacerbation of neurodegenerative disorders. Experimental investigations suggest that anesthetics bind to and affect cytoskeletal microtubules, and that anesthesia-related cognitive dysfunction involves microtubule instability, hyper-phosphorylation of the microtubule-associated protein tau, and tau separation from microtubules. However, exact mechanisms are yet to be identified. In this paper the interaction of anesthetics with the microtubule subunit protein tubulin is investigated using computer-modeling methods. Homology modeling, molecular dynamics simulations and surface geometry techniques were used to determine putative binding sites for volatile anesthetics on tubulin. This was followed by free energy based docking calculations for halothane (2-bromo-2-chloro-1,1,1-trifluoroethane) on the tubulin body, and C-terminal regions for specific tubulin isotypes. Locations of the putative binding sites, halothane binding energies and the relation to cytoskeleton function are reported in this paper. NSERC (Natural Sciences and Engineering Research Council); Alberta Cancer Foundation [23783]; Alberta Advanced Education and Technology [TC-LSI-08-09]; Allard Foundation [23676]; Portuguese Foundation for Science and Technology (FCT) info:eu-repo/semantics/publishedVersion

To provide more accurate computational estimates of binding free energies in solution from molecular dynamics (MD) simulations, a separate solvation contribution for the binding ligand is determined from a linear response treatment. We use explicit water coordinates for this term and combine with MM-PBSA (molecular mechanics, Poisson-Boltzmann, and surface area contributions) in a new approach (MM-PB/LRA-SA). To assess this method, application to the binding between theophylline and its derivatives to an RNA aptamer was performed and compared with experimental binding affinities. Explicitly solvated MD trajectories were generated with the same parameter set used in the previous work by Gouda et al., who compared the relative binding of these molecules by both the MM-PBSA and thermodynamic integration methods. Substituting the linear response term for the ligand in the MM-PB/LRA-SA approach led to an improvement upon MM-PBSA when compared with experimental and thermodynamic integration results at approximately twice the computational cost. The balance between accuracy and computational expense achieved using this method suggests potential advantages in applying it in the virtual drug-screening process.

Adenosine triphosphate (ATP) is known to be the main energy currency of the living cell, and is used as a coenzyme to generate energy for many cellular processes through hydrolysis to adenosine diphosphate (ADP), although the mechanism of energy transfer is not well understood. It has been proposed that following hydrolysis of the ATP cofactor bound to a protein, up to two quanta of amide-I vibrational energy are excited and utilized to bring about important structural changes in the protein. To study whether, and how, amide-I vibrational excitations are capable of leading to protein structural changes, we have added components arising from quantum-mechanical amide-I vibrational excitations to the total energy and force terms within a molecular-dynamics simulation. This model is applied to helical deca-alanine as a test case to investigate how its dynamics differs in the presence or absence of an amide-I excitation. We find that the presence of an amide-I excitation can bias the structure toward a more helical state.