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 hepatitis C virus (HCV) envelope glycoproteins E1 and E2 are critical in viral attachment and cell fusion, and studies of these proteins may provide valuable insights into their potential uses in vaccines and antiviral strategies. Progress has included elucidating the crystal structures of portions of their ectodomains, as well as many other studies of hypervariable regions, stem regions, glycosylation sites, and the participation of E1/E2 in viral fusion with the endosomal membrane. The available structural data have shed light on the binding sites of cross-neutralizing antibodies. A large amount of information has been discovered concerning heterodimerization, including the roles of transmembrane domains, disulfide bonding, and heptad repeat regions. The possible organization of higher order oligomers within the HCV virion has also been evaluated on the basis of experimental data. In this review, E1/E2 structure and function is discussed, and some important issues requiring further study are highlighted.

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.

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.