In molecular simulations with fixed-charge force fields, the choice of partial atomic charges influences numerous computed physical properties, including binding free energies. Many molecular mechanics force fields specify how nonbonded parameters should be determined, but various choices are often available for how these charges are to be determined for arbitrary small molecules. Here, we compute hydration free energies for a set of 44 small, neutral molecules in two different explicit water models (TIP3P and TIP4P-Ew) to examine the influence of charge model on agreement with experiment. Using the AMBER GAFF force field for nonbonded parameters, we test several different methods for obtaining partial atomic charges, including two fast methods exploiting semiempirical quantum calculations and methods deriving charges from the electrostatic potentials computed with several different levels of ab initio quantum calculations with and without a continuum reaction field treatment of solvent. We find that the best charge sets give a root-mean-square error from experiment of roughly 1 kcal/mol. Surprisingly, agreement with experimental hydration free energies does not increase substantially with increasing level of quantum theory, even when the quantum calculations are performed with a reaction field treatment to better model the aqueous phase. We also find that the semiempirical AM1-BCC method for computing charges works almost as well as any of the more computationally expensive ab initio methods and that the root-mean-square error reported here is similar to that for implicit solvent models reported in the literature. Further, we find that the discrepancy with experimental hydration free energies grows substantially with the polarity of the compound, as does its variation across theory levels.

We relate the dynamical behavior of molecular liquids confined in mesoscopic cylindrical pores to the thermodynamic properties, heat capacity and density and to the static structure by combining different experimental methods (H-NMR, calorimetry, elastic and inelastic neutron scattering, numerical simulations). The crystallization process is greatly reduced or avoided by confinement under standard cooling conditions, instead a glass transition temperature T(g) at the 1000s time scale can be observed. The pore averaged local structure of the confined liquid is not noticeably affected when "excluded-volume" corrections are carefully applied, but follows the density changes reflected by the Bragg peak intensities of the porous matrices. The pore size dependence of T(g) is dominated by two factors, surface interaction and finite-size effect. For the smallest pores ([Formula: see text], [Formula: see text] being the van der Waals radius of a molecule), one observes an increase of T(g) and a broadening of the transition region, related to the interaction with the surface that induces a slowing-down of the molecules close to the wall. This is confirmed by neutron scattering experiments and molecular-dynamics simulations at shorter time scales and higher temperatures, which indicate a remaining fraction of frozen molecules. For larger pore sizes, taking the decrease of density under confinement conditions into account, a decrease of T(g) is observed. This could be related to finite-size effects onto the putative cooperativity length that is often invoked to explain glass formation. However, no quantitative determination of this length (not to mention its T-dependence) can be extracted, since the interaction with the wall itself introduces an additional length that adds to the complexity of the problem.

TABLE 1:: Experimental a and Computed b Hydration Free Energies for Molecules in this Study (in kcal/mol) -4.4

DNA damages induced by oxidative intrastrand cross-links have been the subject of intense research during the past decade. Yet, the currently available experimental protocols used to isolate such lesions only allow to get structural information about linked dinucleotides. The detailed structure of the damaged DNA macromolecule has remained elusive. In this study we generated in silico the most frequent oxidative intrastrand cross-link adduct, G[8,5-Me]T, embedded in a solvated DNA dodecamer by means of quantum mechanics/molecular mechanics (QM/MM) Car-Parrinello simulations. The free energy of activation required to bring the reactant close together and to form the C-C covalent-bond is estimated to be ~10 kcal/mol. We observe that the G[8,5-Me]T tandem lesion is accommodated with almost no perturbation of the Watson-Crick hydrogen-bond network and induces bend and unwinding angles of ~20° and 8°, respectively. This rather small structural distortion of the DNA macromolecule compared to other well characterized intrastrand cross-links, such as cyclobutane pyrimidines dimers or cisplatin-DNA complex adduct, is a probable rationale for the known lack of efficient repair of oxidative damages.

The reaction of singlet molecular oxygen with purine DNA bases is investigated by computational means. We support the formation of a transient endoperoxide for guanine and by classical molecular dynamics simulations we demonstrate that the formation of this adduct does not affect the B-helicity. We thus identify the guanine endoperoxide as a key intermediate, confirming a low-temperature nuclear magnetic resonance proof of its existence, and we delineate its degradation pathway, tracing back the preferential formation of 8-oxoguanine versus spiro-derivates in B-DNA. Finally, the latter oxidized 8-oxodGuo product exhibits an almost barrierless reaction profile, and hence is found, coherently with experience, to be much more reactive than guanine itself. On the contrary, in agreement with experimental observations, singlet-oxygen reactivity onto adenine is kinetically blocked by a higher energy transition state.

Protein crystals with their precise, periodic array of functional building blocks have potential applications in biomaterials, sensing and catalysis. This paper describes how a highly-porous crystalline framework of a cationic redox protein and an anionic macrocycle can be modulated by a small cationic effector. Ternary composites of protein (~13 kDa), calix[8]arene (~1.5 kDa) and effector (~0.2 kDa) formed distinct crystalline architectures, dependent on the effector concentration and the crystallization technique. A combination of X-ray crystallography and DFT calculations was used to decipher the framework variations, which appear to depend on a calixarene conformation change mediated by the effector. This "switch" calixarene was observed in three states, each of which is associated with a different interaction network. Two structures obtained by co-crystallization with the effector contained an additional protein "pillar", resulting in framework duplication and decreased porosity. These results suggest how protein assembly can be engineered by supramolecular host-guest interactions.

Benzophenone, the parent of the diarylketone family, is a versatile compound commonly used as a UV blocker. It may also trigger triplet-based DNA photosensitization. Therefore, benzophenone is involved in DNA photodamage induction. In the absence of experimentally resolved structure, the mechanism of DNA damage production remains elusive. Employing a hybrid quantum mechanics/molecular mechanics approach, here we address the spin transfer mechanism between this drug and proximal thymine, that is, the DNA nucleobase most prone to suffer triplet damages.

Clustered apurinic/apyrimidinic (AP; abasic) DNA lesions produced by ionizing radiation are by far more cytotoxic than isolated AP lesion entities. The structure and dynamics of a series of seven 23-bp oligonucleotides featuring simple bistranded clustered damage sites, comprising of two AP sites, zero, one, three or five bases 3′ or 5′ apart from each other, were investigated through 400 ns explicit solvent molecular dynamics simulations. They provide representative structures of synthetically engineered multiply damage sites-containing oligonucleotides whose repair was investigated experimentally (Nucl. Acids Res. 2004, 32:5609-5620; Nucl. Acids Res. 2002, 30: 2800–2808). The inspection of extrahelical positioning of the AP sites, bulge and non Watson–Crick hydrogen bonding corroborates the experimental measurements of repair efficiencies by bacterial or human AP endonucleases Nfo and APE1, respectively. This study provides unprecedented knowledge into the structure and dynamics of clustered abasic DNA lesions, notably rationalizing the non-symmetry with respect to 3′ to 5′ position. In addition, it provides strong mechanistic insights and basis for future studies on the effects of clustered DNA damage on the recognition and processing of these lesions by bacterial or human DNA repair enzymes specialized in the processing of such lesions.