Energetic coupling between DNA bending and base pair opening

Systematic assessment of the flexibility of uracil damaged DNA

Uracil is a common DNA lesion which is recognized and removed by uracil DNA-glycosylase (UDG) as a part of the base excision repair pathway. Excision proceeds by base flipping, and UDG efficiency is thought to depend on the ease of deformability of the bases neighboring the lesion. We used molecular dynamics simulations to assess the flexibility of a large library of dsDNA strands, containing all tetranucleotide motifs with U:A, U:G, T:A or C:G base pairs. Our study demonstrates that uracil damaged DNA largely follows trends in flexibility of undamaged DNA. Measured bending persistence lengths, groove widths, step parameters and base flipping propensities demonstrate that uracil increases the flexibility of DNA, and that U:G base paired strands are more flexible than U:A strands. Certain sequence contexts are more deformable than others, with a key role for the 3' base next to uracil. Flexibilities are large when this base is an A or G, and repressed for a C or T. A 5' T adjacent to the uracil strongly promotes flexibility, but other 5' bases are less influential. DNA bending is correlated to step deformations and base flipping, and bending aids flipping. Our study implies that the link between substrate flexibility and UDG efficiency is widely valid, helps explain why UDG prefers to bind U:G base paired strands, and suggests that the DNA bending angle of the UDG-substrate complex is optimal for base flipping.Communicated by Ramaswamy H. Sarma.

Nucleosome Array Deformation in Chromatin is Sustained by Bending, Twisting and Kinking of Linker DNA

Chromatin in the nucleus undergoes mechanical stresses from different sources during the various stages of cell life. Here a trinucleosome array is used as the minimal model to study the mechanical response to applied stress at the molecular level. By using large-scale, all-atom steered-molecular dynamics simulations, we show that the largest part of mechanical stress in compression is accommodated by the DNA linkers joining pairs of nucleosomes, which store the elastic energy accumulated by the applied force. Different mechanical instabilities (Euler bending, Brazier kinking, twist-bending) can deform the DNA canonical structure, as a function of the increasing force load. An important role of the histone tails in assisting the DNA deformation is highlighted. The overall response of the smallest chromatin fragment to compressive stress leaves the nucleosome assembly with a substantial plastic deformation and localised defects, which can have a potential impact on DNA transcription, downstream signaling pathways, the regulation of gene expression, and DNA repair.

Register-Shifted Structures in Base-Flipped Uracil-Damaged DNA

We report the occurrence of register-shifted structures in simulations of uracil-containing dsDNA. These occur when the 3' base vicinal to uracil is thymine in U:A base-paired DNA. Upon base flipping of uracil, this 3' thymine hydrogen bonds with the adenine across the uracil instead of its complementary base. The register-shifted structure is persistent and sterically blocks re-entry of uracil into the helix stack. Register shifting might be important for DNA repair since the longer exposure of the lesion in register-shifted structures could facilitate enzymatic recognition and repair.

Uracil-DNA glycosylase efficiency is modulated by substrate rigidity

Uracil DNA-glycosylase (UNG) is a DNA repair enzyme that removes the highly mutagenic uracil lesion from DNA using a base flipping mechanism. Although this enzyme has evolved to remove uracil from diverse sequence contexts, UNG excision efficiency depends on DNA sequence. To provide the molecular basis for rationalizing UNG substrate preferences, we used time-resolved fluorescence spectroscopy, NMR imino proton exchange measurements, and molecular dynamics simulations to measure UNG specificity constants (kcat/KM) and DNA flexibilities for DNA substrates containing central AUT, TUA, AUA, and TUT motifs. Our study shows that UNG efficiency is dictated by the intrinsic deformability around the lesion, establishes a direct relationship between substrate flexibility modes and UNG efficiency, and shows that bases immediately adjacent to the uracil are allosterically coupled and have the greatest impact on substrate flexibility and UNG activity. The finding that substrate flexibility controls UNG efficiency is likely significant for other repair enzymes and has major implications for the understanding of mutation hotspot genesis, molecular evolution, and base editing.

Bicyclic Caged Morpholino Oligonucleotides for Optical Gene Silencing

Caged morpholino oligonucleotides (cMOs) are synthetic tools that allow light-inducible gene silencing in live organisms. Previously reported cMOs have utilized hairpin, duplex, and cyclic structures, as well as caged nucleobases. While these antisense technologies enable efficient optical control of RNA splicing and translation, they can have limited dynamic range. A new caging strategy was developed where the two MO termini are conjugated to an internal position through a self-immolative trifunctional linker, thereby generating a bicyclic cMO that is conformationally resistant to RNA binding. The efficacy of this alternative cMO design has been demonstrated in zebrafish embryos and compared to linear MOs and monocyclic constructs.

Markov state models elucidate the stability of DNA influenced by the chiral 5S-Tg base

The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol-1 through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol-1 of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6HO6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson-Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 KBT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

Strand asymmetry influences mismatch resolution during a single-strand annealing

Background

Biases of DNA repair can shape the nucleotide landscape of genomes at evolutionary timescales. The molecular mechanisms of those biases are still poorly understood because it is difficult to isolate the contributions of DNA repair from those of DNA damage.

Results

Here, we develop a genome-wide assay whereby the same DNA lesion is repaired in different genomic contexts. We insert thousands of barcoded transposons carrying a reporter of DNA mismatch repair in the genome of mouse embryonic stem cells. Upon inducing a double-strand break between tandem repeats, a mismatch is generated if the break is repaired through single-strand annealing. The resolution of the mismatch showed a 60–80% bias in favor of the strand with the longest 3′ flap. The location of the lesion in the genome and the type of mismatch had little influence on the bias. Instead, we observe a complete reversal of the bias when the longest 3′ flap is moved to the opposite strand by changing the position of the double-strand break in the reporter.

Conclusions

These results suggest that the processing of the double-strand break has a major influence on the repair of mismatches during a single-strand annealing.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02665-3.

CRISPR-Cas9 bends and twists DNA to read its sequence

In bacterial defense and genome editing applications, the CRISPR-associated protein Cas9 searches millions of DNA base pairs to locate a 20-nucleotide, guide RNA-complementary target sequence that abuts a protospacer-adjacent motif (PAM). Target capture requires Cas9 to unwind DNA at candidate sequences using an unknown ATP-independent mechanism. Here we show that Cas9 sharply bends and undertwists DNA on PAM binding, thereby flipping DNA nucleotides out of the duplex and toward the guide RNA for sequence interrogation. Cryogenic-electron microscopy (cryo-EM) structures of Cas9-RNA-DNA complexes trapped at different states of the interrogation pathway, together with solution conformational probing, reveal that global protein rearrangement accompanies formation of an unstacked DNA hinge. Bend-induced base flipping explains how Cas9 'reads' snippets of DNA to locate target sites within a vast excess of nontarget DNA, a process crucial to both bacterial antiviral immunity and genome editing. This mechanism establishes a physical solution to the problem of complementarity-guided DNA search and shows how interrogation speed and local DNA geometry may influence genome editing efficiency.

The roles of HMGB1-produced DNA gaps in DNA protection and aging biomarker reversal

The endogenous DNA damage triggering an aging progression in the elderly is prevented in the youth, probably by naturally occurring DNA gaps. Decreased DNA gaps are found during chronological aging in yeast. So we named the gaps "Youth-DNA-GAPs." The gaps are hidden by histone deacetylation to prevent DNA break response and were also reduced in cells lacking either the high-mobility group box (HMGB) or the NAD-dependent histone deacetylase, SIR2. A reduction in DNA gaps results in shearing DNA strands and decreasing cell viability. Here, we show the roles of DNA gaps in genomic stability and aging prevention in mammals. The number of Youth-DNA-GAPs were low in senescent cells, two aging rat models, and the elderly. Box A domain of HMGB1 acts as molecular scissors in producing DNA gaps. Increased gaps consolidated DNA durability, leading to DNA protection and improved aging features in senescent cells and two aging rat models similar to those of young organisms. Like the naturally occurring Youth-DNA-GAPs, Box A-produced DNA gaps avoided DNA double-strand break response by histone deacetylation and SIRT1, a Sir2 homolog. In conclusion, Youth-DNA-GAPs are a biomarker determining the DNA aging stage (young/old). Box A-produced DNA gaps ultimately reverse aging features. Therefore, DNA gap formation is a potential strategy to monitor and treat aging-associated diseases.

Dynamics of 5R-Tg Base Flipping in DNA Duplexes Based on Simulations─Agreement with Experiments and Beyond

Damaged or mismatched DNA bases are normally thought to be able to flip out of the helical stack, providing enzymes with access to the faulty genetic information otherwise hidden inside the helix. Thymine glycol (Tg) is one of the most common products of nucleic acid damage. However, the static and dynamic structures of DNA duplexes affected by 5R-Tg epimers are still not clearly understood, including the ability of these to undergo spontaneous base flipping. Structural effects of the 5R-Tg epimers on the duplex DNA are herein studied using molecular dynamics together with reliable DFT based calculations. In comparison with the corresponding intact DNA, the cis-5R,6S-Tg epimer base causes little perturbation to the duplex DNA, and a barrier of 4.9 kcal mol^–1 is obtained by meta-eABF for cis-5R,6S-Tg base flipping out of the duplex DNA, comparable to the 5.4 kcal mol^–1 obtained for the corresponding thymine flipping in intact DNA. For the trans-5R,6R-Tg epimer, three stable local structures were identified, of which the most stable disrupts the Watson–Crick hydrogen-bonded G5/C20 base pair, leading to conformational distortion of the duplex. Interestingly, the relative barrier height of the 5R-Tg flipping is only 1.0 kcal mol^–1 for one of these trans-5R,6R-Tg epimers. Water bridge interactions were identified to be essential for 5R-Tg flipping. The study clearly demonstrates the occurrence of partial trans-5R,6R-Tg epimer flipping in solution.

Using diverse potentials and scoring functions for the development of improved machine-learned models for protein-ligand affinity and docking pose prediction

The advent of computational drug discovery holds the promise of significantly reducing the effort of experimentalists, along with monetary cost. More generally, predicting the binding of small organic molecules to biological macromolecules has far-reaching implications for a range of problems, including metabolomics. However, problems such as predicting the bound structure of a protein-ligand complex along with its affinity have proven to be an enormous challenge. In recent years, machine learning-based methods have proven to be more accurate than older methods, many based on simple linear regression. Nonetheless, there remains room for improvement, as these methods are often trained on a small set of features, with a single functional form for any given physical effect, and often with little mention of the rationale behind choosing one functional form over another. Moreover, it is not entirely clear why one machine learning method is favored over another. In this work, we endeavor to undertake a comprehensive effort towards developing high-accuracy, machine-learned scoring functions, systematically investigating the effects of machine learning method and choice of features, and, when possible, providing insights into the relevant physics using methods that assess feature importance. Here, we show synergism among disparate features, yielding adjusted R² with experimental binding affinities of up to 0.871 on an independent test set and enrichment for native bound structures of up to 0.913. When purely physical terms that model enthalpic and entropic effects are used in the training, we use feature importance assessments to probe the relevant physics and hopefully guide future investigators working on this and other computational chemistry problems.

New machine learning and physics-based scoring functions for drug discovery

Scoring functions are essential for modern in silico drug discovery. However, the accurate prediction of binding affinity by scoring functions remains a challenging task. The performance of scoring functions is very heterogeneous across different target classes. Scoring functions based on precise physics-based descriptors better representing protein–ligand recognition process are strongly needed. We developed a set of new empirical scoring functions, named DockTScore, by explicitly accounting for physics-based terms combined with machine learning. Target-specific scoring functions were developed for two important drug targets, proteases and protein–protein interactions, representing an original class of molecules for drug discovery. Multiple linear regression (MLR), support vector machine and random forest algorithms were employed to derive general and target-specific scoring functions involving optimized MMFF94S force-field terms, solvation and lipophilic interactions terms, and an improved term accounting for ligand torsional entropy contribution to ligand binding. DockTScore scoring functions demonstrated to be competitive with the current best-evaluated scoring functions in terms of binding energy prediction and ranking on four DUD-E datasets and will be useful for in silico drug design for diverse proteins as well as for specific targets such as proteases and protein–protein interactions. Currently, the MLR DockTScore is available at www.dockthor.lncc.br.

Coarse-grained modelling of DNA plectoneme pinning in the presence of base-pair mismatches

Damaged or mismatched DNA bases result in the formation of physical defects in double-stranded DNA. In vivo, defects in DNA must be rapidly and efficiently repaired to maintain cellular function and integrity. Defects can also alter the mechanical response of DNA to bending and twisting constraints, both of which are important in defining the mechanics of DNA supercoiling. Here, we use coarse-grained molecular dynamics (MD) simulation and supporting statistical-mechanical theory to study the effect of mismatched base pairs on DNA supercoiling. Our simulations show that plectoneme pinning at the mismatch site is deterministic under conditions of relatively high force (>2 pN) and high salt concentration (>0.5 M NaCl). Under physiologically relevant conditions of lower force (0.3 pN) and lower salt concentration (0.2 M NaCl), we find that plectoneme pinning becomes probabilistic and the pinning probability increases with the mismatch size. These findings are in line with experimental observations. The simulation framework, validated with experimental results and supported by the theoretical predictions, provides a way to study the effect of defects on DNA supercoiling and the dynamics of supercoiling in molecular detail.

Interactions of replication initiator RctB with single- and double-stranded DNA in origin opening of Vibrio cholerae chromosome 2

Studies of bacterial chromosomes and plasmids indicate that their replication initiator proteins bind to origins of replication at many double-stranded sites and also at AT-rich regions where single-stranded DNA is exposed during origin opening. Single-strand binding apparently promotes origin opening by stabilizing an open structure, but how the initiator participates in this process and the contributions of the several binding sites remain unclear. Here, we show that the initiator protein of Vibrio cholerae specific to chromosome 2 (Chr2) also has single-strand binding activity in the AT-rich region of its origin. Binding is strand specific, depends on repeats of the sequence 5'ATCA and is greatly stabilized in vitro by specific double-stranded sites of the origin. The stability derives from the formation of ternary complexes of the initiator with the single- and double-stranded sites. An IHF site lies between these two kinds of sites in the Chr2 origin and an IHF-induced looping out of the intervening DNA mediates their interaction. Simultaneous binding to two kinds of sites in the origin appears to be a common mechanism by which bacterial replication initiators stabilize an open origin.

Application of a Distance-Dependent Sigmoidal Dielectric Constant to the REMC/SAAP3D Simulations of Chignolin, Trp-Cage, and the G10q Mutant

The replica-exchange Monte Carlo method based on the single amino acid potential (SAAP) force field, i.e., REMC/SAAP3D, was recently developed by our group for the molecular simulation of short peptides. In this study, the method has been improved by applying a distance-dependent dielectric (DDD) constant and extended to the peptides containing D-amino acid (AA) residues. For chignolin (10 AAs), a sigmoidal DDD model reasonably allocated the native-like β-hairpin structure with all-atom root mean square deviation (RMSD) = 2.0 Å as a global energy minimum. The optimal DDD condition was subsequently applied for Trp-cage (20 AAs) and its G10q mutant. The native-like α-rich folded structures with main-chain RMSD = 3.7 and 3.8 Å were obtained as global energy minima for Trp-cage and G10q, respectively. The results suggested that the REMC/SAAP3D method with the sigmoidal DDD model is useful for structural prediction for the short peptides comprised of up to 20 AAs. In addition, the relative contributions of SAAP to the total energy (%SAAP) were evaluated by energetic component analysis. The ratios of %SAAP were about 40 and 20% for chignolin and Trp-cage (or G10q), respectively. It was proposed that SAAP is more important for the secondary structure formation than for the assembly to a higher-order folded structure, in which the attractive van der Waals interaction may play a more important role.

Softening of DNA near melting as disappearance of an emergent property

Near the melting transition the bending elastic constant κ, an emergent property of double-stranded DNA (dsDNA), is shown not to follow the rodlike scaling for small-length N. The reduction in κ with temperature is determined by the denatured bubbles for a continuous transition, e.g., when the two strands are Gaussian, but by the broken bonds near the open end in a Y-like configuration for a first-order transition as for strands with excluded volume interactions. In the latter case, a lever rule is operational, implying a phase coexistence although dsDNA is known to be a single phase.

Critical Sites of DNA Backbone Integrity for Damaged Base Removal by Formamidopyrimidine-DNA Glycosylase

DNA glycosylases, the enzymes that initiate base excision DNA repair, recognize damaged bases through a series of precisely orchestrated movements. Most glycosylases sharply kink the DNA axis at the lesion site and extrude the target base from the DNA double helix into the enzyme's active site. Little attention has been paid so far to the role of the physical continuity of the DNA backbone in allowing the required conformational distortion. Here, we analyze base excision by formamidopyrimidine-DNA glycosylase (Fpg) from substrates keeping all phosphates but containing a nick within three nucleotides of the lesion in either DNA strand. Four phosphoester linkages at the damaged nucleotide and two nucleotides 3' to it were essential for Fpg activity, while the breakage of the others, even at the same critical phosphates, had no effect or even stimulated the reaction. Reduction of the likelihood of hydrogen bonding at the nicks by using dideoxynucleotides as their 3'-terminal groups was more detrimental for the activity. All phosphoester bonds in the complementary strand were dispensable for base excision, but nicks close to the orphaned nucleotide caused early termination of damaged strand cleavage. Elastic network analysis of Fpg-DNA structures showed that the vibrational motions of the critical phosphates are strongly correlated, in part due to the presence of the protein. Overall, our results suggest that mechanical forces propagating along the DNA backbone play a critical role in the correct conformational distortion of DNA by Fpg and possibly by other target base-everting DNA glycosylases.

Emerging Roles of DNA Glycosylases and the Base Excision Repair Pathway

The base excision repair (BER) pathway historically has been associated with maintaining genome integrity by eliminating nucleobases with small chemical modifications. In the past several years, however, BER was found to play additional roles in genome maintenance and metabolism, including sequence-specific restriction modification and repair of bulky adducts and interstrand crosslinks. Central to this expanded biological utility are specialized DNA glycosylases - enzymes that selectively excise damaged, modified, or mismatched nucleobases. In this review we discuss the newly identified roles of the BER pathway and examine the structural and mechanistic features of the DNA glycosylases that enable these functions.

Combinatorial control of gene function with wavelength-selective caged morpholinos

Caged morpholino oligonucleotides (cMOs) are useful research tools in developmental biology because they allow spatiotemporal control of gene expression in whole organisms. While cMOs are usually triggered by light of a single wavelength, the introduction of spectrally distinct chromophores can enable combinatorial regulation of multiple genes. This chapter describes the general principles and methods of wavelength-selective cMO design and synthesis from commercially available reagents. Synthetic protocols for the linkers and the two-step cMO assembly are described in detail, as well as the microinjection and photoactivation techniques. Following these protocols, spectrally separated cyclic cMOs for multiple genes of interest can be prepared, enabling their inhibition in zebrafish embryos and other animal models.

Toxicity and repair of DNA adducts produced by the natural product yatakemycin

Yatakemycin (YTM) is an extraordinarily toxic DNA alkylating agent with potent antimicrobial and antitumor properties and the most recent addition to the CC-1065 and duocarmycin family of natural products. While bulky DNA lesions the size of those produced by YTM are normally removed from the genome by the nucleotide excision repair (NER) pathway, YTM adducts are also a substrate for the bacterial DNA glycosylases AlkD and YtkR2, unexpectedly implicating base excision repair (BER) in their elimination. The reason for the extreme toxicity of these lesions and the molecular basis for how they are eliminated by BER have been unclear. Here, we describe the structural and biochemical properties of YTM adducts responsible for their toxicity, and define the mechanism by which they are excised by AlkD. These findings delineate an alternative strategy for repair of bulky DNA damage and establish the cellular utility of this pathway relative to that of NER.

Free Energy Coupling between DNA Bending and Base Flipping

Free energy simulations are presented to probe the energetic coupling between DNA bending and the flipping of a central thymine in double stranded DNA 13mers. The energetics are shown to depend on the neighboring base pairs, and upstream C or T or downstream C tended to make flipping more costly. Flipping to the major groove side was generally preferred. Bending aids flipping, by pushing the system up in free energy, but for small and intermediate bending angles the two were uncorrelated. At higher bending angles, bending and flipping became correlated, and bending primed the system for base flipping toward the major groove. Flipping of the 6-4 pyrimidine-pyrimidone and pyrimidine dimer photoproducts is shown to be more facile than for undamaged DNA. For the damages, major groove flipping was preferred, and DNA bending was much facilitated in the 6-4 pyrimidine-pyrimidone damaged system. Aspects of the calculations were verified by structural analyses of protein-DNA complexes with flipped bases.

[KCl] Dependence of B-DNA Groove Bending Anisotropy

The energetics of B-DNA bending toward the major and minor grooves were quantified by free energy simulations at four different KCl concentrations. Increased [KCl] led to more flexible DNA, with persistence lengths that agreed well with experimental values. At all salt concentrations, major groove bending was preferred, although preferences for major and minor groove bending were similar for the A-tract containing sequence. Since the phosphate repulsions and DNA internal energy favored minor groove bending, the preference for major groove bending was thought to originate from differences in solvation. Water in the minor groove was tighter bound than water in the major groove, and harder to displace than major groove water, which favored the compression of the major groove upon bending. Higher [KCl] decreased the persistence length for both major and minor groove bending but did not greatly affect the free energy spacing between the minor and major groove bending curves. For sequences without A-tracts, salt affected major and minor bending to nearly the same degree, and did not change the preference for major groove bending. For the A-tract containing sequence, an increase in salt concentration decreased the already small energetic difference between major and minor groove bending. Since salts did not significantly affect the relative differences in bending energetics and hydration, it is likely that the increased bending flexibilities upon salt increase are simply due to screening.

Targeting HIV-1 Envelope Proteins Using a Fragment Discovery All-Atom Computational Algorithm

INTRODUCTION

HIV viral envelope proteins are targets for small inhibitor molecules aimed at disrupting the cellular entry process. Potential peptide-class inhibitor molecules (rDNA drugs) have been previously identified, with mixed results, through biomimicry and phage display experimental methods. Here we describe a new approach based on computational fragment discovery. The method has the potential to not only optimize peptide binding affinity but also to rapidly produce alternative inhibitors against mutated strains.

METHODS

A comprehensive, all-atom implicit solvent method is used to bombard the C-heptad repeat unit of HIV-1 target envelope protein GP41 with single D-amino acid residues as they exist in their native state. A nascent peptide computational search process then identifies potential favorable sequences of attached ligands based on four peptide bond criteria. Finally, dynamic simulations of nascent peptides attached to host targets help refine potential peptide inhibitors for experimental HIV-1 challenge assays and testing.

RESULTS AND DISCUSSION

Initial testing of the method was done using 64,000 total ligands of D-amino acid residues at a total computational time of 0.05 microseconds per ligand, which resulted in several thousand attached ligands. Peptide bond criteria search employing three of the four bond constraints with a tolerance of 20 percent, resulted in four potential peptide inhibitors of 5 to 6 residues in length. Only one of the four peptides demonstrated IC50 values and partial viral inhibition based on cell challenge assays using CEM-SS host cells. That peptide inhibitor also computationally demonstrated long-time attachment and stability to a helical groove in its C-heptad target. This initial testing of peptide fragment discovery against HIV-1 has helped us refine the protocols and identify key areas of improvement.

CONCLUSION

Our methods demonstrate the potential to design efficient peptide inhibitors to viral target proteins based on an all-atom dynamic simulation and using a ligand library as fragments of potential nascent peptides. Our methods can be greatly improved through the use of higher numbers of ligands, increased time of bombardment, and tighter constraints on the peptide bond search step. Our method may be important in the need to rapidly respond to target mutations and to advance multiple targeting methods based multiple peptide inhibitors.

Both DNA global deformation and repair enzyme contacts mediate flipping of thymine dimer damage

The photo-induced cis-syn-cyclobutane pyrimidine (CPD) dimer is a frequent DNA lesion. In bacteria photolyases efficiently repair dimers employing a light-driven reaction after flipping out the CPD damage to the active site. How the repair enzyme identifies a damaged site and how the damage is flipped out without external energy is still unclear. Employing molecular dynamics free energy calculations, the CPD flipping process was systematically compared to flipping undamaged nucleotides in various DNA global states and bound to photolyase enzyme. The global DNA deformation alone (without protein) significantly reduces the flipping penalty and induces a partially looped out state of the damage but not undamaged nucleotides. Bound enzyme further lowers the penalty for CPD damage flipping with a lower free energy of the flipped nucleotides in the active site compared to intra-helical state (not for undamaged DNA). Both the reduced penalty and partial looping by global DNA deformation contribute to a significantly shorter mean first passage time for CPD flipping compared to regular nucleotides which increases the repair likelihood upon short time encounter between repair enzyme and DNA.

Photo-Tethers for the (Multi-)Cyclic, Conformational Caging of Long Oligonucleotides

Intramolecular circularization of DNA oligonucleotides was accomplished by incorporation of alkyne-modified photolabile nucleosides into DNA sequences, followed by a Cu^I -catalyzed alkyne-azide cycloaddition with bis-azido linker molecules. We determined a range of ring sizes, in which the caged circular oligonucleotides exhibit superior duplex destabilizing properties. Specific binding of a full-length 90 nt C10 aptamer recognizing human Burkitt's lymphoma cells was then temporarily inhibited by locking the aptamer in a bicircularized structure. Irradiation restored the native aptamer conformation resulting in efficient cell binding and uptake. The photo-tether strategy presented here provides a robust and versatile tool for the light-activation of longer functional oligonucleotides, noteworthy without prior knowledge on the structure and the importance of specific nucleotides within a DNA aptamer.

Global deformation facilitates flipping of damaged 8-oxo-guanine and guanine in DNA

Oxidation of guanine (Gua) to form 7,8-dihydro-8-oxoguanine (8oxoG) is a frequent mutagenic DNA lesion. DNA repair glycosylases such as the bacterial MutM can effciently recognize and eliminate the 8oxoG damage by base excision. The base excision requires a 8oxoG looping out (flipping) from an intrahelical base paired to an extrahelical state where the damaged base is in the enzyme active site. It is still unclear how the damage is identified and flipped from an energetically stable stacked and paired state without any external energy source. Free energy simulations have been employed to study the flipping process for globally deformed DNA conformational states. DNA deformations were generated by systematically untwisting the DNA to mimic its conformation in repair enzyme encounter complex. The simulations indicate that global DNA untwisting deformation toward the enzyme bound form alone (without protein) significantly reduces the penalty for damage flipping to about half of the penalty observed in regular DNA. The finding offers a mechanistic explanation how binding free energy that is transformed to binding induced DNA deformation facilitates flipping and helps to rapidly detect a damaged base. It is likely of general relevance since repair enzyme binding frequently results in significant deformation of the target DNA.

Anisotropy of B-DNA Groove Bending

DNA bending is critical for DNA packaging, recognition, and repair, and occurs toward either the major or the minor groove. The anisotropy of B-DNA groove bending was quantified for eight DNA sequences by free energy simulations employing a novel reaction coordinate. The simulations show that bending toward the major groove is preferred for non-A-tracts while the A-tract has a high tendency of bending toward the minor groove. Persistence lengths were generally larger for bending toward the minor groove, which is thought to originate from differences in groove hydration. While this difference in stiffness is one of the factors determining the overall preference of bending direction, the dominant contribution is shown to be a free energy offset between major and minor groove bending. The data suggests that, for the A-tract, this offset is largely determined by inherent structural properties, while differences in groove hydration play a large role for non-A-tracts. By quantifying the energetics of DNA groove bending and rationalizing the origins of the anisotropy, the calculations provide important new insights into a key biological process.

Rigidity of melting DNA

The temperature dependence of DNA flexibility is studied in the presence of stretching and unzipping forces. Two classes of models are considered. In one case the origin of elasticity is entropic due to the polymeric correlations, and in the other the double-stranded DNA is taken to have an intrinsic rigidity for bending. In both cases single strands are completely flexible. The change in the elastic constant for the flexible case due to thermally generated bubbles is obtained exactly. For the case of intrinsic rigidity, the elastic constant is found to be proportional to the square root of the bubble number fluctuation.

Free-energy landscape and characteristic forces for the initiation of DNA unzipping

DNA unzipping, the separation of its double helix into single strands, is crucial in modulating a host of genetic processes. Although the large-scale separation of double-stranded DNA has been studied with a variety of theoretical and experimental techniques, the minute details of the very first steps of unzipping are still unclear. Here, we use atomistic molecular-dynamics simulations, coarse-grained simulations, and a statistical-mechanical model to study the initiation of DNA unzipping by an external force. Calculation of the potential of mean force profiles for the initial separation of the first few terminal basepairs in a DNA oligomer revealed that forces ranging between 130 and 230 pN are needed to disrupt the first basepair, and these values are an order of magnitude larger than those needed to disrupt basepairs in partially unzipped DNA. The force peak has an echo of ∼50 pN at the distance that unzips the second basepair. We show that the high peak needed to initiate unzipping derives from a free-energy basin that is distinct from the basins of subsequent basepairs because of entropic contributions, and we highlight the microscopic origin of the peak. To our knowledge, our results suggest a new window of exploration for single-molecule experiments.

The DNA glycosylase AlkD uses a non-base-flipping mechanism to excise bulky lesions

Threats to genomic integrity arising from DNA damage are mitigated by DNA glycosylases, which initiate the base excision repair pathway by locating and excising aberrant nucleobases. How these enzymes find small modifications within the genome is a current area of intensive research. A hallmark of these and other DNA repair enzymes is their use of base flipping to sequester modified nucleotides from the DNA helix and into an active site pocket. Consequently, base flipping is generally regarded as an essential aspect of lesion recognition and a necessary precursor to base excision. Here we present the first, to our knowledge, DNA glycosylase mechanism that does not require base flipping for either binding or catalysis. Using the DNA glycosylase AlkD from Bacillus cereus, we crystallographically monitored excision of an alkylpurine substrate as a function of time, and reconstructed the steps along the reaction coordinate through structures representing substrate, intermediate and product complexes. Instead of directly interacting with the damaged nucleobase, AlkD recognizes aberrant base pairs through interactions with the phosphoribose backbone, while the lesion remains stacked in the DNA duplex. Quantum mechanical calculations revealed that these contacts include catalytic charge-dipole and CH-π interactions that preferentially stabilize the transition state. We show in vitro and in vivo how this unique means of recognition and catalysis enables AlkD to repair large adducts formed by yatakemycin, a member of the duocarmycin family of antimicrobial natural products exploited in bacterial warfare and chemotherapeutic trials. Bulky adducts of this or any type are not excised by DNA glycosylases that use a traditional base-flipping mechanism. Hence, these findings represent a new model for DNA repair and provide insights into catalysis of base excision.

Lesion search and recognition by thymine DNA glycosylase revealed by single molecule imaging

The ability of DNA glycosylases to rapidly and efficiently detect lesions among a vast excess of nondamaged DNA bases is vitally important in base excision repair (BER). Here, we use single molecule imaging by atomic force microscopy (AFM) supported by a 2-aminopurine fluorescence base flipping assay to study damage search by human thymine DNA glycosylase (hTDG), which initiates BER of mutagenic and cytotoxic G:T and G:U mispairs in DNA. Our data reveal an equilibrium between two conformational states of hTDG–DNA complexes, assigned as search complex (SC) and interrogation complex (IC), both at target lesions and undamaged DNA sites. Notably, for both hTDG and a second glycosylase, hOGG1, which recognizes structurally different 8-oxoguanine lesions, the conformation of the DNA in the SC mirrors innate structural properties of their respective target sites. In the IC, the DNA is sharply bent, as seen in crystal structures of hTDG lesion recognition complexes, which likely supports the base flipping required for lesion identification. Our results support a potentially general concept of sculpting of glycosylases to their targets, allowing them to exploit the energetic cost of DNA bending for initial lesion sensing, coupled with continuous (extrahelical) base interrogation during lesion search by DNA glycosylases.

Molecular modeling of nucleic Acid structure: electrostatics and solvation

This unit presents an overview of computer simulation techniques as applied to nucleic acid systems, ranging from simple in vacuo molecular modeling techniques to more complete all-atom molecular dynamics treatments that include an explicit representation of the environment. The third in a series of four units, this unit focuses on critical issues in solvation and the treatment of electrostatics. UNITS 7.5 & 7.8 introduced the modeling of nucleic acid structure at the molecular level. This included a discussion of how to generate an initial model, how to evaluate the utility or reliability of a given model, and ultimately how to manipulate this model to better understand its structure, dynamics, and interactions. Subject to an appropriate representation of the energy, such as a specifically parameterized empirical force field, the techniques of minimization and Monte Carlo simulation, as well as molecular dynamics (MD) methods, were introduced as a way of sampling conformational space for a better understanding of the relevance of a given model. This discussion highlighted the major limitations with modeling in general. When sampling conformational space effectively, difficult issues are encountered, such as multiple minima or conformational sampling problems, and accurately representing the underlying energy of interaction. In order to provide a realistic model of the underlying energetics for nucleic acids in their native environments, it is crucial to include some representation of solvation (by water) and also to properly treat the electrostatic interactions. These subjects are discussed in detail in this unit.

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.

Modified Poisson–Boltzmann equations for characterizing biomolecular solvation

Methods for computing electrostatic interactions often account implicitly for the solvent, due to the much smaller number of degrees of freedom involved. In the Poisson–Boltzmann (PB) approach the electrostatic potential is obtained by solving the Poisson–Boltzmann equation (PBE), where the solvent region is modeled as a homogeneous medium with a high dielectric constant. PB however is not exempt of problems. It does not take into account for example the sizes of the ions in the atmosphere surrounding the solute, nor does it take into account the inhomogeneous dielectric response of water due to the presence of a highly charged surface. In this paper we review two major modifications of PB that circumvent these problems, namely the size-modified PB (SMPB) equation and the Dipolar Poisson–Boltzmann Langevin (DPBL) model. In SMPB, steric effects between ions are accounted for with a lattice gas model. In DPBL, the solvent region is no longer modeled as a homogeneous dielectric media but rather as an assembly of self-orienting interacting dipoles of variable density. This model results in a dielectric profile that transits smoothly from the solute to the solvent region as well as in a variable solvent density that depends on the charges of the solute. We show successful applications of the DPBL formalism to computing the solvation free energies of isolated ions in water. Further developments of more accurately modified PB models are discussed.

A simple contact mapping algorithm for identifying potential peptide mimetics in protein-protein interaction partners

A simple, static contact mapping algorithm has been developed as a first step at identifying potential peptide biomimetics from protein interaction partner structure files. This rapid and simple mapping algorithm, “OpenContact” provides screened or parsed protein interaction files based on specified criteria for interatomic separation distances and interatomic potential interactions. The algorithm, which uses all-atom Amber03 force field models, was blindly tested on several unrelated cases from the literature where potential peptide mimetics have been experimentally developed to varying degrees of success. In all cases, the screening algorithm efficiently predicted proposed or potential peptide biomimetics, or close variations thereof, and provided complete atom-atom interaction data necessary for further detailed analysis and drug development. In addition, we used the static parsing/mapping method to develop a peptide mimetic to the cancer protein target, epidermal growth factor receptor. In this case, secondary, loop structure for the peptide was indicated from the intra-protein mapping, and the peptide was subsequently synthesized and shown to exhibit successful binding to the target protein. The case studies, which all involved experimental peptide drug advancement, illustrate many of the challenges associated with the development of peptide biomimetics, in general. Proteins 2014; 82:2253–2262. © 2014 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.

Reduced atomic pair-interaction design (RAPID) model for simulations of proteins

Increasingly, theoretical studies of proteins focus on large systems. This trend demands the development of computational models that are fast, to overcome the growing complexity, and accurate, to capture the physically relevant features. To address this demand, we introduce a protein model that uses all-atom architecture to ensure the highest level of chemical detail while employing effective pair potentials to represent the effect of solvent to achieve the maximum speed. The effective potentials are derived for amino acid residues based on the condition that the solvent-free model matches the relevant pair-distribution functions observed in explicit solvent simulations. As a test, the model is applied to alanine polypeptides. For the chain with 10 amino acid residues, the model is found to reproduce properly the native state and its population. Small discrepancies are observed for other folding properties and can be attributed to the approximations inherent in the model. The transferability of the generated effective potentials is investigated in simulations of a longer peptide with 25 residues. A minimal set of potentials is identified that leads to qualitatively correct results in comparison with the explicit solvent simulations. Further tests, conducted for multiple peptide chains, show that the transferable model correctly reproduces the experimentally observed tendency of polyalanines to aggregate into β-sheets more strongly with the growing length of the peptide chain. Taken together, the reported results suggest that the proposed model could be used to succesfully simulate folding and aggregation of small peptides in atomic detail. Further tests are needed to assess the strengths and limitations of the model more thoroughly.

Binding of nucleoid-associated protein fis to DNA is regulated by DNA breathing dynamics

Physicochemical properties of DNA, such as shape, affect protein-DNA recognition. However, the properties of DNA that are most relevant for predicting the binding sites of particular transcription factors (TFs) or classes of TFs have yet to be fully understood. Here, using a model that accurately captures the melting behavior and breathing dynamics (spontaneous local openings of the double helix) of double-stranded DNA, we simulated the dynamics of known binding sites of the TF and nucleoid-associated protein Fis in Escherichia coli. Our study involves simulations of breathing dynamics, analysis of large published in vitro and genomic datasets, and targeted experimental tests of our predictions. Our simulation results and available in vitro binding data indicate a strong correlation between DNA breathing dynamics and Fis binding. Indeed, we can define an average DNA breathing profile that is characteristic of Fis binding sites. This profile is significantly enriched among the identified in vivo E. coli Fis binding sites. To test our understanding of how Fis binding is influenced by DNA breathing dynamics, we designed base-pair substitutions, mismatch, and methylation modifications of DNA regions that are known to interact (or not interact) with Fis. The goal in each case was to make the local DNA breathing dynamics either closer to or farther from the breathing profile characteristic of a strong Fis binding site. For the modified DNA segments, we found that Fis-DNA binding, as assessed by gel-shift assay, changed in accordance with our expectations. We conclude that Fis binding is associated with DNA breathing dynamics, which in turn may be regulated by various nucleotide modifications.

Cellular transcription factors (TFs) are proteins that regulate gene expression, and thereby cellular activity and fate, by binding to specific DNA segments. The physicochemical determinants of protein-DNA binding specificity are not completely understood. Here, we report that the propensity of transient opening and re-closing of the double helix, resulting from thermal fluctuations, aka “DNA breathing” or “DNA bubbles,” can be associated with binding affinity in the case of Fis, a well-studied nucleoid-associated protein in Escherichia coli. We found that a particular breathing profile is characteristic of high-affinity Fis binding sites and that DNA fragments known to bind Fis in vivo are statistically enriched for this profile. Furthermore, we used simulations of DNA breathing dynamics to guide design of gel-shift experiments aimed at testing the idea that local breathing influences Fis binding. As a result, we show that via nucleotide modifications but without modifying nucleotides that directly contact Fis, we were able to transform a low-affinity Fis binding site into a high-affinity site and vice versa. The nucleotide modifications were designed only based on DNA breathing simulations. Our study suggests that strong Fis-DNA binding depends on DNA breathing - a novel physicochemical characteristic that could be used for prediction and rational design of TF binding sites.

Effect of Bending Rigidity on the Knotting of a Polymer under Tension

A coarse-grained computational model is used to investigate how the bending rigidity of a polymer under tension affects the formation of a trefoil knot. Thermodynamic integration techniques are applied to demonstrate that the free-energy cost of forming a knot has a minimum at nonzero bending rigidity. The position of the minimum exhibits a power-law dependence on the applied tension. For knotted polymers with nonuniform bending rigidity, the knots preferentially localize in the region with a bending rigidity that minimizes the free energy.

Cyclic caged morpholinos: conformationally gated probes of embryonic gene function

Feeling a bit cagey: morpholino-based antisense reagents have been caged through oligonucleotide cyclization, enabling photocontrol of gene expression in zebrafish embryos and larvae. Using these reagents, the timing of exocrine cell fate commitment in the developing pancreas has been examined.