Access to

Topoisomerase 1 (Top1) enzymes regulate DNA superhelicity by forming covalent cleavage complexes that undergo controlled rotation. Substitution of nucleoside analogs at the +1 position of the DNA duplex relative to the Top1 cleavage site inhibits DNA religation. The reduced efficiency for Top1-mediated religation contributes to the anticancer activity of widely used anticancer drugs including fluoropyrimidines and gemcitabine. In the present study, we report that mismatched base pairs at the +1 position destabilize the duplex DNA components for a model Top1 cleavage complex formation even though one duplex component does not directly include a mismatched base pair. Molecular dynamics simulations reveal G-dU and G-FdU mismatched base pairs, but not a G-T mismatched base pair, increase flexibility at the Top1 cleavage site, and affect coupling between the regions required for the religation reaction to occur. These results demonstrate that substitution of dT analogs into the +1 position of the non-scissile strand alters the stability and flexibility of DNA contributing to the reduced efficiency for Top1-mediated DNA religation. These effects are inherent in the DNA duplex and do not require formation of the Top1:DNA complex. These results provide a biophysical rationale for the inhibition of Top1-mediated DNA religation by nucleotide analog substitution

Magnetic fields and density functional theory

A major focus of this dissertation is the development of functionals for the magnetic susceptibility and the chemical shielding within the context of magnetic field density functional theory (BDFT). These functionals depend on the electron density in the absence of the field, which is unlike any other treatment of these responses. There have been several advances made within this theory. The first of which is the development of local density functionals for chemical shieldings and magnetic susceptibilities. There are the first such functionals ever proposed. These parameters have been studied by constructing functionals for the current density and then using the Biot-Savart equations to obtain the responses. In order to examine the advantages and disadvantages of the local functionals, they were tested numerically on some small molecules

Autoinhibitory mechanisms of ERG studied by molecular dynamics simulations

ERG, an ETS-family transcription factor, acts as a regulator of differentiation of early hematopoietic cells. It contains an autoinhibitory domain, which negatively regulates DNA-binding. The mechanism of autoinhibitory is still illusive. To understand the mechanism, we study the dynamical properties of ERG protein by molecular dynamics simulations. These simulations suggest that DNA binding autoinhibition associates with the internal dynamics of ERG. Specifically, we find that (1), The N-C terminal correlation in the inhibited ERG is larger than that in uninhibited ERG that contributes to the autoinhibition of DNA-binding. (2), DNA-binding changes the property of the N-C terminal correlation from being anti-correlated to correlated, that is, changing the relative direction of the correlated motions and (3), For the Ets-domain specifically, the inhibited and uninhibited forms exhibit essentially the same dynamics, but the binding of the DNA decreases the fluctuation of the Ets-domain. We also find from PCA analysis that the three systems, even with quite different dynamics, do have highly similar free energy surfaces, indicating that they share similar conformations

The molecular origin of the MMR-dependent apoptosis pathway from dynamics analysis of MutSα-DNA complexes

<div><p>The cellular response to DNA damage signaling by mismatch-repair (MMR) proteins is incompletely understood. It is generally accepted that MMR-dependent apoptosis pathway in response to DNA damage detection is independent of MMR’s DNA repair function. In this study, we investigate correlated motions in response to the binding of mismatched and platinum cross-linked DNA fragments by MutSα, as derived from 50 ns molecular dynamics simulations. The protein dynamics in response to the mismatched and damaged DNA recognition suggests that MutSα signals their recognition through independent pathways providing evidence for the molecular origin of the MMR-dependent apoptosis. MSH2 subunit is indicated to play a key role in signaling both mismatched and damaged DNA recognition; localized and collective motions within the protein allow identifying sites on the MSH2 surface possible involved in recruiting proteins responsible for downstream events. Unlike in the mismatch complex, predicted key communication sites specific for the damage recognition are on the list of known cancer-causing mutations or deletions. This confirms MSH2’s role in signaling DNA damage-induced apoptosis and suggests that defects in MMR alone is sufficient to trigger tumorigenesis, supporting the experimental evidence that MMR-damage response function could protect from the early occurrence of tumors. Identifying these particular communication sites may have implications for the treatment of cancers that are not defective for MMR, but are unable to function optimally for MMR-dependent responses following DNA damage such as the case of resistance to cisplatin.</p> <p>An animated Interactive 3D Complement (I3DC) is available in Proteopedia at <a href="http://proteopedia.org/w/Journal:JBSD:15" target="_blank">http://proteopedia.org/w/Journal:JBSD:15</a></p> </div

Constant pH molecular dynamics using continuous titration coordinates. Proteins 2004;56: 738–752

ABSTRACT In this work, we explore the question of whether pK a calculations based on a microscopic description of the protein and a macroscopic description of the solvent can be implemented to examine conformationally dependent proton shifts in proteins. To this end, we introduce a new method for performing constant-pH molecular dynamics (PHMD) simulations utilizing the generalized Born implicit solvent model. This approach employs an extended Hamiltonian in which continuous titration coordinates propagate simultaneously with the atomic motions of the system. The values adopted by these coordinates are modulated by potentials of mean force of isolated titratable model groups and the pH to control the proton occupation at particular sites in the polypeptide. Our results for four different proteins yield an absolute average error of ϳ1.6 pK units, and point to the role that thermally driven relaxation of the protein environment in the vicinity of titrating groups plays in modulating the local pK a , thereby influencing the observed pK 1/2 values. While the accuracy of our method is not yet equivalent to methods that obtain pK 1/2 values through the ad hoc scaling of electrostatics, the present approach and constant pH methods in general provide a useful framework for studying pHdependent phenomena. Further work to improve our model to approach quantitative agreement with experiment is outlined. Proteins 2004;56:738 -752

Probing light chain mutation effects on thrombin via molecular dynamics simulations and machine learning

<p>Thrombin is a key component for chemotherapeutic and antithrombotic therapy development. As the physiologic and pathologic roles of the light chain still remain vague, here, we continue previous efforts to understand the impacts of the disease-associated single deletion of LYS9 in the light chain. By combining supervised and unsupervised machine learning methodologies and more traditional structural analyses on data from 10 μs molecular dynamics simulations, we show that the conformational ensemble of the ΔK9 mutant is significantly perturbed. Our analyses consistently indicate that LYS9 deletion destabilizes both the catalytic cleft and regulatory functional regions and result in some conformational changes that occur in tens to hundreds of nanosecond scaled motions. We also reveal that the two forms of thrombin each prefer a distinct binding mode of a Na⁺ ion. We expand our understanding of previous experimental observations and shed light on the mechanisms of the LYS9 deletion associated bleeding disorder by providing consistent but more quantitative and detailed structural analyses than early studies in literature. With a novel application of supervised learning, i.e. the decision tree learning on the hydrogen bonding features in the wild-type and ΔK9 mutant forms of thrombin, we predict that seven pairs of critical hydrogen bonding interactions are significant for establishing distinct behaviors of wild-type thrombin and its ΔK9 mutant form. Our calculations indicate the LYS9 in the light chain has both localized and long-range allosteric effects on thrombin, supporting the opinion that light chain has an important role as an allosteric effector.</p

All-atom molecular dynamics comparison of disease-associated zinc fingers

<p>An important regulatory domain of NF-B Essential Modulator (NEMO) is a ubiquitin-binding zinc finger, with a tetrahedral CYS3HIS1 zinc-coordinating binding site. Two variations of NEMO’s zinc finger are implicated in various disease states including ectodermal dysplasia and adult-onset glaucoma. To discern structural and dynamical differences between these disease states, we present results of 48-s of molecular dynamics simulations for three zinc finger systems each in two states, with and without zinc-bound and correspondingly appropriate cysteine thiol/thiolate configurations. The wild-type protein, often studied for its role in cancer, maintains the most rigid and conformationally stable zinc-bound configuration compared with the diseased counterparts. The glaucoma-related protein has persistent loss of secondary structure except within the dominant conformation. Conformational overlap between wild-type and glaucoma isoforms indicate a competitive binding mechanism may be substantial in the malfunctioning configuration, while the alpha-helical disruption of the ectodermal dysplasia suggests a loss of binding selectivity is responsible for aberrant function.</p