Dynamic behavior of DNA base pairs containing 8-oxoguanine

The process by which DNA repair enzymes recognize and selectively excise damaged bases in duplex DNA is fundamental to our mechanistic understanding of these critical biological reactions. 8-Oxoguanine (8-oxoG) is the most common form of oxidative DNA damage; unrepaired, this lesion generates a G:C-->T:A mutation. Central to the recognition and repair of DNA damage is base extrusion, a process in which the damaged base lesion or, in some cases, its partner disengages from the helix and is bound to the enzyme's active site where base excision takes place. The conformation adopted by 8-oxoG in duplex DNA is affected by the base positioned opposite this lesion; conformational changes may also take place when the damaged base binds to its cognate repair enzyme. We performed unrestrained molecular dynamics simulations for several 13-mer DNA duplexes. Oligomers containing G:C and 8oxoG:C pairs adopted Watson-Crick geometries in stable B-form duplexes; 8oxoG showed increased local and global flexibility and a reduced barrier to base extrusion. Duplexes containing the G:A mismatch showed much larger structural fluctuations and failed to adopt a well-defined structure. For the 8oxoG:A mismatch that is recognized by the DNA glycosylase MutY, the damaged nucleoside underwent spontaneous and reproducible anti-->syn transitions. The syn conformation is thermodynamically preferred. Steric hindrance and unfavorable electrostatics associated with the 8oxoG O8 atom in the anti conformation were the major driving forces for this transition. Transition events follow two qualitatively different pathways. The overall anti-->syn transition rate and relative probability of the two transition paths were dependent on local sequence context. These simulations indicate that both the dynamic and equilibrium behavior of the duplex change as a result of oxidation; these differences may provide valuable new insight into the selective action of enzymes on damaged DNA.     

Use of the tubulin bound paclitaxel conformation for structure-based rational drug design

A new computational docking protocol has been developed and used in combination with conformational information inferred from REDOR-NMR experiments on microtubule bound 2-(p-fluorobenzoyl)paclitaxel to delineate a unique tubulin binding structure of paclitaxel. A conformationally constrained macrocyclic taxoid bearing a linker between the C-14 and C-3'N positions has been designed and synthesized to enforce this "REDOR-taxol" conformation. The novel taxoid SB-T-2053 inhibits the growth of MCF-7 and LCC-6 human breast cancer cells (wild-type and drug resistant) on the same order of magnitude as paclitaxel. Moreover, SB-T-2053 induces in vitro tubulin polymerization at least as well as paclitaxel, which directly validates our drug design process. These results open a new avenue for drug design of next generation taxoids and other microtubule-stabilizing agents based on the refined structural information of drug-tubulin complexes, in accordance with typical enzyme-inhibitor medicinal chemistry precepts.     

Cumulant functions in optical coherence theory

Cumulant functions are introduced to describe the statistical state of a radiation field. These functions are simply related to the optical coherence functions but have some interesting features. It is shown that if the cumulant functions of all orders greater than some numberN ₀ vanish then they also vanish for all orders greater than 2. Thermal field is the only field having this property. This property holds whether the field is described by a classical stochastic process or by a quantum density operator. Further the particular operator ordering used in defining these cumulant functions for the quantized field affects only the second order cumulant function. To describe the statistical state of a vector field such as partially polarized or unpolarized radiation, one would need to introduce cumulant tensors.     

Gas phase activation energy for unimolecular dissociation of biomolecular ions determined by focused RAdiation for gaseous multiphoton ENergy transfer (FRAGMENT)

We present a novel approach for the determination of activation energy for the unimolecular dissociation of a large (>50 atoms) ion, based on measurement of the unimolecular dissociation rate constant as a function of continuous-wave CO(2) laser intensity. Following a short ( approximately 1 s) induction period, CO(2) laser irradiation produces an essentially blackbody internal energy distribution, whose 'temperature' varies inversely with laser intensity. The only currently available method for measuring such activation energies is blackbody infrared radiative dissociation (BIRD). Compared with BIRD, FRAGMENT: (a) eliminates the need to heat the surrounding ion trap and vacuum chamber to each of several temperatures (each requiring hours for temperature equilibration); (b) offers a three-fold wider range of effective blackbody temperature; and (c) extends the range of applications to include initially cold ions (e.g., gas-phase H/D exchange). Our FRAGMENT-determined activation energy for dissociation of protonated bradykinin, 1.2 +/- 0.1 eV, agrees within experimental error to the value, 1.3 +/- 0.1 eV, previously reported by Williams et al. from BIRD experiments. Copyright 1999 John Wiley & Sons, Ltd.     

The open structure of a multi-drug-resistant HIV-1 protease is stabilized by crystal packing contacts

The introduction of HIV-1 protease (HIV-PR) inhibitors has led to a dramatic increase in patient survival; however, these gains are threatened by the emergence of multi-drug-resistant strains. Design of inhibitors that overcome resistance would be greatly facilitated by deeper insight into the mechanistic events associated with binding of substrates and inhibitors, as well as an understanding of the effects of resistance mutations on the structure and dynamic behavior of HIV-PR. We previously reported a series of simulations that provide a model for HIV-PR dynamics, with spontaneous conversions between the bound and unbound crystal forms upon addition or removal of an inhibitor. Importantly, the unbound protease transiently sampled a third fully open state that permits entry to the active site, unlike both crystallographic forms. Recently, a crystal structure of unbound HIV-PR was reported for the MDR 769 isolate (PDB: 1TW7); unlike all previous experimental structures, the binding pocket is open. It is suggested that drug resistance in this strain arises at least in part from the inability of inhibitors to induce closing. We carried out simulations of the MDR 769 HIV-PR mutant and observed that the reported structure is unstable in solution and rapidly adopts the semi-open conformation observed for the unbound wild-type protease in solution. Further analysis suggests that the wide-open structure observed for MDR 769 arises not from sequence variation, but instead is an artifact from crystal packing. Thus, despite being the first experimental structure to reveal flap opening sufficient for substrate access to the active site, this structure may not be directly relevant to studies of inhibitor entry or to the cause of HIV-PR drug resistance.     

Energetic preference of 8-oxoG eversion pathways in a DNA glycosylase

Base eversion is a fundamental process in the biochemistry of nucleic acids, allowing proteins engaged in DNA repair and epigenetic modifications to access target bases in DNA. Crystal structures reveal end points of these processes, but not the pathways involved in the dynamic process of base recognition. To elucidate the pathway taken by 8-oxoguanine during base excision repair by Fpg, we calculated free energy surfaces during eversion of the damaged base through the major and minor grooves. The minor groove pathway and free energy barrier (6-7 kcal/mol) are consistent with previously reported results (Qi, Y.; Spong, M. C.; Nam, K.; Banerjee, A.; Jiralerspong, S.; Karplus, M.; Verdine, G. L. Nature 2009, 462, 762.) However, eversion of 8-oxoG through the major groove encounters a significantly lower barrier (3-4 kcal/mol) more consistent with experimentally determined rates of enzymatic sliding during lesion search (Blainey, P. C.; van Oijent, A. M.; Banerjee, A.; Verdine, G. L.; Xie, X. S. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 5752.). Major groove eversion has been suggested for other glycosylases, suggesting that in addition to function, dynamics of base eversion may also be conserved.     

Secondary structure bias in generalized Born solvent models: comparison of conformational ensembles and free energy of solvent polarization from explicit and implicit solvation

The effects of the use of three generalized Born (GB) implicit solvent models on the thermodynamics of a simple polyalanine peptide are studied via comparing several hundred nanoseconds of well-converged replica exchange molecular dynamics (REMD) simulations using explicit TIP3P solvent to REMD simulations with the GB solvent models. It is found that when compared to REMD simulations using TIP3P the GB REMD simulations contain significant differences in secondary structure populations, most notably an overabundance of alpha-helical secondary structure. This discrepancy is explored via comparison of the differences in the electrostatic component of the free energy of solvation (DeltaDeltaG(pol)) between TIP3P (via thermodynamic Integration calculations), the GB models, and an implicit solvent model based on the Poisson equation (PE). The electrostatic components of the solvation free energies are calculated using each solvent model for four representative conformations of Ala10. Since the PE model is found to have the best performance with respect to reproducing TIP3P DeltaDeltaG(pol) values, effective Born radii from the GB models are compared to effective Born radii calculated with PE (so-called perfect radii), and significant and numerous deviations in GB radii from perfect radii are found in all GB models. The effect of these deviations on the solvation free energy is discussed, and it is shown that even when perfect radii are used the agreement of GB with TIP3P DeltaDeltaG(pol) values does not improve. This suggests a limit to the optimization of the effective Born radius calculation and that future efforts to improve the accuracy of GB models must extend beyond such optimizations.     

The Amber biomolecular simulation programs

We describe the development, current features, and some directions for future development of the Amber package of computer programs. This package evolved from a program that was constructed in the late 1970s to do Assisted Model Building with Energy Refinement, and now contains a group of programs embodying a number of powerful tools of modern computational chemistry, focused on molecular dynamics and free energy calculations of proteins, nucleic acids, and carbohydrates.     

Solution structure of HIV-1 protease flaps probed by comparison of molecular dynamics simulation ensembles and EPR experiments

The introduction of multidrug treatment regimens has dramatically prolonged the progression and survival of AIDS patients. However, the success of the long-term treatment has been hindered by strains of HIV that are increasingly resistant to inhibitors of targets such as HIV protease (HIV PR). Therefore, the need for a thorough understanding of the structure and dynamics of HIV PR and how these are altered in resistant mutants is crucial for the design of more effective treatments. Crystal structures of unbound HIV PR show significant heterogeneity and often have extensive crystal packing interactions. Recent site-directed spin labeling (SDSL) and double electron-electron resonance (DEER) spectroscopy studies characterized flap conformations in HIV-1 protease in an inhibited and uninhibited form and distinguished the extent of flap opening in an unbound form. However, the correlation between EPR-measured interspin distances and structural/dynamic features of the flaps has not been established. In this report, we link EPR-based data and 900 ns of MD simulation in explicit water to gain insight into the ensemble of conformations sampled by HIV PR flaps in solution, both in the presence and in the absence of an FDA-approved HIV PR inhibitor.     

Coupling of replica exchange simulations to a non-Boltzmann structure reservoir

Computing converged ensemble properties remains challenging for large biomolecules. Replica exchange molecular dynamics (REMD) can significantly increase the efficiency of conformational sampling by using high temperatures to escape kinetic traps. Several groups, including ours, introduced the idea of coupling replica exchange to a pre-converged, Boltzmann-populated reservoir, usually at a temperature higher than that of the highest temperature replica. This procedure reduces computational cost because the long simulation times needed for extensive sampling are only carried out for a single temperature. However, a weakness of the approach is that the Boltzmann-weighted reservoir can still be difficult to generate. We now present the idea of employing a non-Boltzmann reservoir, whose structures can be generated through more efficient conformational sampling methods. We demonstrate that the approach is rigorous and derive a correct statistical mechanical exchange criterion between the reservoir and the replicas that drives Boltzmann-weighted probabilities for the replicas. We test this approach on the trpzip2 peptide and demonstrate that the resulting thermal stability profile is essentially indistinguishable from that obtained using very long (>100 ns) standard REMD simulations. The convergence of this reservoir-aided REMD is significantly faster than for regular REMD. Furthermore, we demonstrate that modification of the exchange criterion is essential; REMD simulations using a standard exchange function with the non-Boltzmann reservoir produced incorrect results.