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.

Counterions play a significant role in DNA structure and function, and molecular dynamics (MD) simulations offer the prospect of detailed description of the dynamical structure of ions at the molecular level. However, the motions of mobile counterions are notably slow to converge in MD on DNA. Obtaining accurate and reliable MD simulations requires knowing just how much sampling is required for convergence of each of the properties of interest. To address this issue, MD on a d(CGCGAATTCGCG) duplex in a dilute aqueous solution of water and 22 Na ؉ counterions was performed until convergence was achieved. The calculated first shell ion occupancies and DNA-Na ؉ radial distribution functions were computed as a function of time to assess convergence, and compared with relaxation times of the DNA internal parameters shift, slide, rise, tilt, roll, and twist. The sequence dependence of fractional occupancies of ions in the major and minor grooves of the DNA is examined, and the possibility of correlation between ion proximity and DNA minor groove widths is investigated. C ounterion structures and motions have been implicated in recent ideas about sequence effects on DNA structure axis curvature and ligand-induced bending (1-7). In each of these theories, the local interactions of ions with polyionic DNA complement the effects described by counterion condensation (CC) theory (8) and the Poisson Boltzmann (PB) equation (9,10). Although the effect of ions and ionic strength on DNA structural and thermodynamics properties are implied from diverse experiments, obtaining a fully detailed molecular model of the ions interacting with DNA based on these results is usually not possible. Oligonucleotide crystal structures reveal only the few ions that are ordered and can be unequivocally assigned. The positions of ions around DNA are generally underdetermined in experiments based on biophysical methods. Recent studies requiring this information have turned to large-scale molecular dynamics (MD) simulations to obtain computational models (11). However, in MD modeling, the complex aggregate of oligonucleotide, water, counterions, and coions is slow to fully stabilize, and ion motions are a rate-determining step in total convergence (12). Current reviews of MD on DNA (12)(13)(14) indicate that all simulations to date are based on considerably shorter trajectories. Therefore, we initiated a project aimed at obtaining demonstrably converged results on ion structures and motions. A related question is the sensitivity of the fast internal motions of the DNA to ion convergence. Our analysis is based on an MD simulation on the prototype B-form duplex d(CGC-GAATTCGCG) in a dilute aqueous solution of water and Na ϩ counterions carried out on a supercomputer until the point of convergence could be reliably determined. The trajectories are used to study the sequence dependence of ion distributions, the DNA-Na ϩ radial distribution functions, and the sensitivity of groove widths to ion proximity. The simulations are used as a basis for a compa...

The amino acid sequence of a protein governs its function. We used bulk competition and focused deep sequencing to investigate the effects of all ubiquitin point mutants on yeast growth rate. Many aspects of ubiquitin function have been carefully studied, which enabled interpretation of our growth analyses in light of a rich structural, biophysical and biochemical knowledge base. In one highly sensitive cluster on the surface of ubiquitin almost every amino acid substitution caused growth defects. In contrast, the opposite face tolerated virtually all possible substitutions. Surface locations between these two faces exhibited intermediate mutational tolerance. The sensitive face corresponds to the known interface for many binding partners. Across all surface positions, we observe a strong correlation between burial at structurally characterized interfaces and the number of amino acid substitutions compatible with robust growth. This result indicates that binding is a dominant determinant of ubiquitin function. In the solvent inaccessible core of ubiquitin all positions tolerated a limited number of substitutions, with hydrophobic amino acids especially interchangeable. Some mutations null for yeast growth were previously shown to populate folded conformations indicating that for these mutants subtle changes to conformation caused functional defects. The most sensitive region to mutation within the core was located near the C-terminus that is a focal binding site for many critical binding partners. These results indicate that core mutations may frequently cause functional defects through subtle disturbances to structure or dynamics.

A common strategy of metabolic engineering is to increase the endogenous supply of precursor metabolites to improve pathway productivity. The ability to further enhance heterologous production of a desired compound may be limited by the inherent capacity of the imported pathway to accommodate high precursor supply. Here, we present engineered diterpenoid biosynthesis as a case where insufficient downstream pathway capacity limits highlevel levopimaradiene production in Escherichia coli. To increase levopimaradiene synthesis, we amplified the flux toward isopentenyl diphosphate and dimethylallyl diphosphate precursors and reprogrammed the rate-limiting downstream pathway by generating combinatorial mutations in geranylgeranyl diphosphate synthase and levopimaradiene synthase. The mutant library contained pathway variants that not only increased diterpenoid production but also tuned the selectivity toward levopimaradiene. The most productive pathway, combining precursor flux amplification and mutant synthases, conferred approximately 2,600-fold increase in levopimaradiene levels. A maximum titer of approximately 700 mg∕L was subsequently obtained by cultivation in a benchscale bioreactor. The present study highlights the importance of engineering proteins along with pathways as a key strategy in achieving microbial biosynthesis and overproduction of pharmaceutical and chemical products.GGPP synthase | levopimaradiene synthase | metabolic engineering | Escherichia coli | molecular reprogramming M etabolic engineering is the enabling technology for the manipulation of organisms to synthesize high-value compounds of both natural and heterologous origin (1-4). In the case of heterologous production, well-characterized microorganisms are used as production hosts because targeted optimization can be performed using widely available genetic tools and synthetic biology frameworks (5, 6). One important application of engineered microbial systems is geared toward the synthesis of terpenoid natural products (7-9). Terpenoids represent one of the largest classes of secondary metabolites that includes pharmaceuticals, cosmetics, and potential biofuels candidates (8,(10)(11)(12). Metabolic engineering approaches to produce terpenoids in microbial systems such as Escherichia coli and yeast have commonly focused on increasing the precursor flux into the heterologous terpenoid pathway by rerouting endogenous isoprenoid metabolism (13-15). These engineering strategies have relied heavily on changing the enzyme concentrations in the product pathway.Many properties of a metabolic pathway, however, are not limited solely by the enzyme concentration, as is particularly true for the terpenoid pathway. In nature, terpenoid biosynthesis is regulated at multiple metabolic branch points to create large structural and functional diversity (16-18). In the major metabolic branch point in terpenoid biosynthesis, the prenyl transferases and terpenoid synthases catalyze the formation of a wide range of structurally diverse acyclic and cyclic terp...

Molecular dynamics (MD) simulations including water and counterions on B-DNA oligomers containing all 136 unique tetranucleotide basepair steps are reported. The objective is to obtain the calculated dynamical structure for at least two copies of each case, use the results to examine issues with regard to convergence and dynamical stability of MD on DNA, and determine the significance of sequence context effects on all unique dinucleotide steps. This information is essential to understand sequence effects on DNA structure and has implications on diverse problems in the structural biology of DNA. Calculations were carried out on the 136 cases embedded 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 simulations were carried out using a well-defined state-of-the-art MD protocol, the AMBER suite of programs, and the parm94 force field. In a previous article (Beveridge et al. 2004. Biophysical Journal. 87:3799-3813), the research design, details of the simulation protocol, and informatics issues were described. Preliminary results from 15 ns MD trajectories were presented for the d(CpG) step in all 10 unique sequence contexts. The results indicated the sequence context effects to be small for this step, but revealed that MD on DNA at this length of trajectory is subject to surprisingly persistent cooperative transitions of the sugar-phosphate backbone torsion angles alpha and gamma. In this article, we report detailed analysis of the entire trajectory database and occurrence of various conformational substates and its impact on studies of context effects. The analysis reveals a possible direct correspondence between the sequence-dependent dynamical tendencies of DNA structure and the tendency to undergo transitions that "trap" them in nonstandard conformational substates. The difference in mean of the observed basepair step helicoidal parameter distribution with different flanking sequence sometimes differs by as much as one standard deviation, indicating that the extent of sequence effects could be significant. The observations reveal that the impact of a flexible dinucleotide such as CpG could extend beyond the immediate basepair neighbors. The results in general provide new insight into MD on DNA and the sequence-dependent dynamical structural characteristics of DNA.

HIV-1 protease is one of the major antiviral targets in the treatment of patients infected with HIV-1. The nine FDA approved HIV-1 protease inhibitors were developed with extensive use of structure-based drug design, thus the atomic details of how the inhibitors bind are well characterized. From this structural understanding the molecular basis for drug resistance in HIV-1 protease can be elucidated. Selected mutations in response to therapy and diversity between clades in HIV-1 protease have altered the shape of the active site, potentially altered the dynamics and even altered the sequence of the cleavage sites in the Gag polyprotein. All of these interdependent changes act in synergy to confer drug resistance while simultaneously maintaining the fitness of the virus. New strategies, such as incorporation of the substrate envelope constraint to design robust inhibitors that incorporate details of HIV-1 protease’s function and decrease the probability of drug resistance, are necessary to continue to effectively target this key protein in HIV-1 life cycle.

Conformational selection and induced fit are well-known contributors to ligand binding and allosteric effects in proteins. Molecular dynamics (MD) simulations now enable the theoretical study of protein-ligand binding in terms of ensembles of interconverting microstates and the population shifts characteristic of "dynamical allostery." Here we investigate protein-ligand binding and allostery based on a Markov state model (MSM) with states and rates obtained from all-atom MD simulations. As an exemplary case, we consider the single domain protein par-6 PDZ with and without ligand and allosteric effector. This is one of the smallest proteins in which allostery has been experimentally observed. In spite of the increased complexity intrinsic to a statistical ensemble perspective, we find that conformational selection and induced fit mechanisms can be readily identified in the analysis. In the nonallosteric pathway, MD-MSM shows that PDZ binds ligand via conformational selection. However, the allosteric pathway requires an activation step that involves a conformational change induced by the allosteric effector Cdc42. Once in the allosterically activated state, we find that ligand binding can proceed by conformational selection. Our MD-MSM model predicts that allostery in this and possibly other systems involves both induced fit and conformational selection, not just one or the other.

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.