We describe herein a computationally intensive project aimed at carrying out molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide base sequences. This initiative was undertaken by an international collaborative effort involving nine research groups, the "Ascona B-DNA Consortium" (ABC). Calculations were carried out on the 136 cases imbedded in 39 DNA oligomers with repeating tetranucleotide sequences, capped on both ends by GC pairs and each having a total length of 15 nucleotide pairs. All MD simulations were carried out using a well-defined protocol, the AMBER suite of programs, and the parm94 force field. Phase I of the ABC project involves a total of approximately 0.6 mus of simulation for systems containing approximately 24,000 atoms. The resulting trajectories involve 600,000 coordinate sets and represent approximately 400 gigabytes of data. In this article, the research design, details of the simulation protocol, informatics issues, and the organization of the results into a web-accessible database are described. Preliminary results from 15-ns MD trajectories are presented for the d(CpG) step in its 10 unique sequence contexts, and issues of stability and convergence, the extent of quasiergodic problems, and the possibility of long-lived conformational substates are discussed.

Newly generated actin free barbed ends at the front of motile cells provide sites for actin filament assembly driving membrane protrusion. Growth factors induce a rapid biphasic increase in actin free barbed ends, and we found both phases absent in fibroblasts lacking H+ efflux by the Na-H exchanger NHE1. The first phase is restored by expression of mutant cofilin-H133A but not unphosphorylated cofilin-S3A. Constant pH molecular dynamics simulations and nuclear magnetic resonance (NMR) reveal pH-sensitive structural changes in the cofilin C-terminal filamentous actin binding site dependent on His133. However, cofilin-H133A retains pH-sensitive changes in NMR spectra and severing activity in vitro, which suggests that it has a more complex behavior in cells. Cofilin activity is inhibited by phosphoinositide binding, and we found that phosphoinositide binding is pH-dependent for wild-type cofilin, with decreased binding at a higher pH. In contrast, phosphoinositide binding by cofilin-H133A is attenuated and pH insensitive. These data suggest a molecular mechanism whereby cofilin acts as a pH sensor to mediate a pH-dependent actin filament dynamics.

Actin filament binding by the focal adhesion (FA)-associated protein talin stabilizes cell-substrate adhesions and is thought to be rate-limiting in cell migration. Although F-actin binding by talin is known to be pH-sensitive in vitro, with lower affinity at higher pH, the functional significance of this pH dependence is unknown. Because increased intracellular pH (pHi) promotes cell migration and is a hallmark of metastatic carcinomas, we asked whether it increases FA remodeling through lower-affinity talin-actin binding. Talin contains several actin binding sites, but we found that only the COOH-terminal USH-I/LWEQ module showed pH-dependent actin binding, with lower affinity and decreased maximal binding at higher pH. Molecular dynamics simulations and NMR of this module revealed a structural mechanism for pH-dependent actin binding. A cluster of titratable amino acids with upshifted pKa values, including His-2418, was identified at one end of the five-helix bundle distal from the actin binding site. Protonation of His-2418 induces changes in the conformation and dynamics of the remote actin binding site. Structural analyses of a mutant talin-H2418F at pH 6.0 and 8.0 suggested changes different from the WT protein, and we confirmed that actin binding by talin-H2418F was relatively pH-insensitive. In motile fibroblasts, increasing pHi decreased FA lifetime and increased the migratory rate. However, expression of talin-H2418F increased lifetime 2-fold and decreased the migratory rate. These data identify a molecular mechanism for pH-sensitive actin binding by talin and suggest that FA turnover is pH-dependent and in part mediated by pH-dependent affinity of talin for binding actin. intracellular pH ͉ NHE1 ͉ migration F ocal adhesion (FA) remodeling is a rate-limiting determinant in haptokinetic migration of adherent cells. At the leading edge of migrating cells, FAs undergo rapid cycles of assembly and turnover, creating and disrupting, respectively, sites of traction necessary for forward movement of the cell body. Force generation for traction requires linkage among the extracellular matrix, integrin receptors, and actin filaments. Actin filaments do not directly bind to the cytoplasmic domain of integrins but bind to integrin-associated FA proteins such as talin and vinculin. Although several mechanisms contribute to FA remodeling in migrating cells (1), emerging evidence indicates that talin plays a central role in the dynamic linkage between integrins and actin filaments necessary for cell migration (2, 3). Talin functions in distinct albeit complementary mechanisms that promote FA turnover. First is cleavage of talin by the protease calpain, which also modulates adhesion complex composition and likely signaling functions of talin (4). Second is regulated talin binding to actin filaments, which is proposed to act as a clutch to control FA turnover and membrane protrusion dynamics (3, 5-7). How actin binding by talin is dynamically regulated during cell migration, however, remains undetermined.Previous stu...

Dihydrofolate reductase (DHFR) is the enzyme responsible for the NADPH-dependent reduction of 5,6-dihydrofolate to 5,6,7,8-tetrahydrofolate, an essential cofactor in the synthesis of purines, thymidylate, methionine, and other key metabolites. Because of its importance in multiple cellular functions, DHFR has been the subject of much research targeting the enzyme with anticancer, antibacterial, and antimicrobial agents. Clinically used compounds targeting DHFR include methotrexate for the treatment of cancer and diaminopyrimidines (DAPs) such as trimethoprim (TMP) for the treatment of bacterial infections. DAP inhibitors of DHFR have been used clinically for >30 years and resistance to these agents has become widespread. Methicillin-resistant Staphylococcus aureus (MRSA), the causative agent of many serious nosocomial and community acquired infections, and other gram-positive organisms can show resistance to DAPs through mutation of the chromosomal gene or acquisition of an alternative DHFR termed ''S1 DHFR.'' To develop new therapies for health threats such as MRSA, it is important to understand the molecular basis of DAP resistance. Here, we report the crystal structure of the wild-type chromosomal DHFR from S. aureus in complex with NADPH and TMP. We have also solved the structure of the exogenous, TMP resistant S1 DHFR, apo and in complex with TMP. The structural and thermodynamic data point to important molecular differences between the two enzymes that lead to dramatically reduced affinity of DAPs to S1 DHFR. These differences in enzyme binding affinity translate into reduced antibacterial activity against strains of S. aureus that express S1 DHFR. Proteins 2009; 76:706-717. V V C 2009 Wiley-Liss, Inc.

We present a thermodynamical approach to identify changes in macromolecular structure and dynamics in response to perturbations such as mutations or ligand binding, using an expansion of the Kullback-Leibler Divergence that connects local population shifts in torsion angles to changes in the free energy landscape of the protein. While the Kullback-Leibler Divergence is a known formula from information theory, the novelty and power of our implementation lies in its formal developments, connection to thermodynamics, statistical filtering, ease of visualization of results, and extendability by adding higher-order terms. We present a formal derivation of the Kullback-Leibler Divergence expansion and then apply our method at a first-order approximation to molecular dynamics simulations of four protein systems where ligand binding or pH titration is known to cause an effect at a distant site. Our results qualitatively agree with experimental measurements of local changes in structure or dynamics, such as NMR chemical shift perturbations and hydrogen-deuterium exchange mass spectrometry. The approach produces easy-to-analyze results with low background, and as such has the potential to become a routine analysis when molecular dynamics simulations in two or more conditions are available. Our method is implemented in the MutInf code package and is available on the SimTK website at https://simtk.org/home/mutinf.

A virtual screening protocol has been applied to seek non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs) and its K103N mutant. First, a chemical similarity search on the Maybridge library was performed using known NNRTIs as reference structures. The top-ranked molecules obtained from this procedure plus 26 known NNRTIs were then docked into the binding sites of the wild-type reverse transcriptase (HIV-RT) and its K103N variant (K103N-RT) using Glide 3.5. The top-ranked 100 compounds from the docking for both proteins were post-scored with a procedure using molecular mechanics and continuum solvation (MM-GB/SA). The validity of the virtual screening protocol was supported by (i) testing of the MM-GB/SA procedure, (ii) agreement between predicted and crystallographic binding poses, (iii) recovery of known potent NNRTIs at the top of both rankings, and (iv) identification of top-scoring library compounds that are close in structure to recently reported NNRTI HTS hits. However, purchase and assaying of selected top-scoring compounds from the library failed to yield active anti-HIV agents. Nevertheless, the highest-ranked database compound, S10087, was pursued as containing a potentially viable core. Subsequent synthesis and assaying of S10087 analogues proposed by further computational analysis yielded anti-HIV agents with EC50 values as low as 310 nM. Thus, with the aid of computational tools, it was possible to evolve a false positive into a true active.

Virtual screening of the Maybridge library of ca. 70,000 compounds was performed using a similarity filter, docking, and MM-GB/SA post-processing to seek potential non-nucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs). Though known NNRTIs were retrieved well, purchase and assaying of representative, top-scoring compounds from the library failed to yield any active anti-HIV agents. However, the highest-ranked library compound, oxadiazole 1, was pursued as a potential "near-miss" with the BOMB program to seek constructive modifications. Subsequent synthesis and assaying of several polychloro-analogs did yield anti-HIV agents with EC 50 values as low as 310 nM. The study demonstrates that it is possible to learn from a formally unsuccessful virtual-screening exercise and, with the aid of computational analyses, to evolve efficiently a false positive into a true active. In addition, the need for adequate structure validation was confirmed by the apparent misrepresentation of a purchased compound elsewhere as the present oxadiazole core compound, 16.

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.