The mechanism of covalent inhibition of LAR phosphatase by illudalic acid.

Leukocyte antigen-related (LAR) phosphatase is a receptor-type protein tyrosine phosphatase involved in cellular signaling and associated with human disease including cancer and metabolic disorders. Selective inhibition of LAR phosphatase activity by well characterized and well validated small molecules would provide key insights into the roles of LAR phosphatase in health and disease, but identifying selective inhibitors of LAR phosphatase activity has been challenging. Recently, we described potent and selective inhibition of LAR phosphatase activity by the fungal natural product illudalic acid. Here we provide a detailed biochemical characterization of the adduct formed between LAR phosphatase and illudalic acid. A mass spectrometric analysis indicates that two cysteine residues are covalently labeled by illudalic acid and a related analog. Mutational analysis supports the hypothesis that inhibition of LAR phosphatase activity is due primarily to the adduct with the catalytic cysteine residue. A computational study suggests potential interactions between the illudalic acid moiety and the enzyme active site. Taken together, these data offer novel insights into the mechanism of inhibition of LAR phosphatase activity by illudalic acid.

Systematic Investigation of Tether Length and Phosphorus Configuration in Backbone Constrained Macrocyclic Nucleic Acids to Modulate Binding Kinetics for RNA.

We recently described a chemical strategy to pre-organize a trinucleotide subunit in a conformation suitable for Watson-Crick base pairing for modulating the binding kinetics of single-stranded oligonucleotides (ONs) using bis-phosphonate esters bridging hydrocarbon tethers to provide 11- and 15-membered macrocyclic analogues. In this manuscript, we describe the synthesis of all eight P-stereoisomers of macrocyclic 12-, 13-, 14-, and 16-membered hydrocarbon-bridged nucleotide trimers, their incorporation into ONs, and biophysical characterization of the modified ONs. The size of the macrocyclic tether and configuration at phosphorus had profound effects on hybridization kinetics. ONs containing 12- and 13-membered rings exhibited faster on-rates (up to 5-fold) and off-rates (up to 161-fold). In contrast, ONs using the larger ring size macrocycles generally exhibited smaller changes in binding kinetics relative to unmodified DNA. Interestingly, several of the analogues retained significant binding affinity for RNA based on their dissociation constants, despite being modestly destabilizing in the thermal denaturation experiments, highlighting the potential utility of measuring dissociation constants versus duplex thermal stability when evaluating novel nucleic acid analogues. Overall, our results provide additional insights into the ability of backbone-constrained macrocyclic nucleic acid analogues to modulate hybridization kinetics of modified ONs with RNA.

Benchmarking the Drude Polarizable Force Field Using the r(GACC) Tetranucleotide.

Polarizable force fields, in particular, the Drude polarizable force field (Drude FF), may hold the key to more accurately modeling biomolecules with molecular dynamics simulations by explicitly accounting for atomic polarizability. Previous work has shown promising results in simulating duplex nucleic acids and protein structures with excellent agreement with experimental values. However, benchmarking the Drude polarizable force field with highly flexible, single-stranded structures has yet to be achieved. In this work, the r(GACC) tetranucleotide is simulated over a multimicrosecond time scale, starting with various different initial conformations. Despite the starting conformation, including starting from the expected dominant A-form major conformation, the experimental structural distribution is not matched. In fact, the major NMR conformation is never resampled. Instead, the r(GACC) tetranucleotide becomes stabilized in anomalous structures that are inconsistent with the NMR data and that favor base-pairing and electrostatic interactions over base stacking. These structures are maintained for lengthy time scales (>1 µs) themselves, suggesting a misbalance of forces in the Drude polarizable force field itself. This model system is suggestive of the fact that currently the Drude polarizable force field does not appear to produce the sensitive balance of forces required to accurately model other single-stranded or noncanonical RNA structures.

AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

Ethidium bromide interactions with DNA: an exploration of a classic DNA-ligand complex with unbiased molecular dynamics simulations.

Visualization of double stranded DNA in gels with the binding of the fluorescent dye ethidium bromide has been a basic experimental technique in any molecular biology laboratory for >40 years. The interaction between ethidium and double stranded DNA has been observed to be an intercalation between base pairs with strong experimental evidence. This presents a unique opportunity for computational chemistry and biomolecular simulation techniques to benchmark and assess their models in order to see if the theory can reproduce experiments and ultimately provide new insights. We present molecular dynamics simulations of the interaction of ethidium with two different double stranded DNA models. The first model system is the classic sequence d(CGCGAATTCGCG)2 also known as the Drew-Dickerson dodecamer. We found that the ethidium ligand binds mainly stacked on, or intercalated between, the terminal base pairs of the DNA with little to no interaction with the inner base pairs. As the intercalation at the terminal CpG steps is relatively rapid, the resultant DNA unwinding, rigidification, and increased stability of the internal base pair steps inhibits further intercalation. In order to reduce these interactions and to provide a larger groove space, a second 18-mer DNA duplex system with the sequence d(GCATGAACGAACGAACGC) was tested. We computed molecular dynamics simulations for 20 independent replicas with this sequence, each with ∼27 µs of sampling time. Results show several spontaneous intercalation and base-pair eversion events that are consistent with experimental observations. The present work suggests that extended MD simulations with modern DNA force fields and optimized simulation codes are allowing the ability to reproduce unbiased intercalation events that we were not able to previously reach due to limits in computing power and the lack of extensively tested force fields and analysis tools.

Riboflavin Stabilizes Abasic, Oxidized G-Quadruplex Structures.

The G-quadruplex is a noncanonical fold of DNA commonly found at telomeres and within gene promoter regions of the genome. These guanine-rich sequences are highly susceptible to damages such as base oxidation and depurination, leading to abasic sites. In the present work, we address whether a vacancy, such as an abasic site, in a G-quadruplex serves as a specific ligand recognition site. When the G-tetrad is all guanines, the vacant (abasic) site is recognized and bound by free guanine nucleobase. However, we aim to understand whether the preference for a specific ligand recognition changes with the presence of a guanine oxidation product 8-oxo-7,8-dihydroguanine (OG) adjacent to the vacancy in the tetrad. Using molecular dynamics simulation, circular dichroism, and nuclear magnetic resonance, we examined the ability for riboflavin to stabilize abasic site-containing G-quadruplex structures. Through structural and free energy binding analysis, we observe riboflavin's ability to stabilize an abasic site-containing G-quadruplex only in the presence of an adjacent OG-modified base. Further, when compared to simulation with the vacancy filled by free guanine, we observe that the free guanine nucleobase is pushed outside of the tetrad by OG to interact with other parts of the structure, including loop residues. These results support the preference of riboflavin over free guanine to fill an OG-adjacent G-quadruplex abasic vacancy.

Backbone Hydrocarbon-Constrained Nucleic Acids Modulate Hybridization Kinetics for RNA.

The binding affinity of therapeutic oligonucleotides (ONs) for their cognate RNA is determined by the rates of association (ka) and dissociation (kd). Single-stranded ONs are highly flexible and can adopt multiple conformations in solution, some of which may not be conducive for hybridization. We investigated if restricting rotation around the sugar-phosphate backbone, by tethering two adjacent backbone phosphonate esters using hydrocarbon bridges, can modulate hybridization kinetics of the modified ONs for complementary RNA. Given the large number of possible analogues with different tether lengths and configurations at the phosphorus atoms, we employed molecular dynamic simulations to optimize the size of the hydrocarbon bridge to guide the synthetic efforts. The backbone-constrained nucleotide trimers with stereodefined configurations at the contiguous backbone phosphorus atoms were assembled using a ring-closing metathesis reaction, then incorporated into oligonucleotides by an in situ synthesis of the phosphoramidites followed by coupling to solid supports. Evaluation of the modified oligonucleotides revealed that 15-membered macrocyclic-constrained analogues displayed similar or slightly improved on-rates but significantly increased off-rates compared to unmodified DNA ONs, resulting in reduced duplex stability. In contrast, LNA ONs with conformationally preorganized furanose rings showed similar on-rates to DNA ONs but very slow off-rates, resulting in net improvement in duplex stability. Furthermore, the experimental data generally supported the molecular dynamics simulation results, suggesting that this strategy can be used as a predictive tool for designing the next generation of constrained backbone ON analogues with improved hybridization properties.

Peptoid Residues Make Diverse, Hyperstable Collagen Triple-Helices.

As the only ribosomally encoded N-substituted amino acid, proline promotes distinct secondary protein structures. The high proline content in collagen, the most abundant protein in the human body, is crucial to forming its hallmark structure: the triple-helix. For over five decades, proline has been considered compulsory for synthetic designs aimed at recapitulating collagen's structure and properties. Here we describe that N-substituted glycines (N-glys), also known as peptoid residues, exhibit a general triple-helical propensity similar to or greater than proline, enabling synthesis of stable triple-helical collagen mimetic peptides (CMPs) with unprecedented side chain diversity. Supported by atomic-resolution crystal structures as well as circular dichroism and computational characterizations spanning over 30 N-gly-containing CMPs, we discovered that N-glys stabilize the triple-helix primarily by sterically preorganizing individual chains into the polyproline-II helix. We demonstrated that N-glys with exotic side chains including a "click"-able alkyne and a photosensitive side chain enable CMPs for functional applications including the spatiotemporal control of cell adhesion and migration. The structural principles uncovered in this study open up opportunities for a new generation of collagen-mimetic therapeutics and materials.

Molecular basis for the broad substrate selectivity of a peptide prenyltransferase.

The cyanobactin prenyltransferases catalyze a series of known or unprecedented reactions on millions of different substrates, with no easily observable recognition motif and exquisite regioselectivity. Here we define the basis of broad substrate tolerance for the otherwise uncharacterized TruF family. We determined the structures of the Tyr-prenylating enzyme PagF, in complex with an isoprenoid donor analog and a panel of linear and macrocyclic peptide substrates. Unexpectedly, the structures reveal a truncated barrel fold, wherein binding of large peptide substrates is necessary to complete a solvent-exposed hydrophobic pocket to form the catalytically competent active site. Kinetic, mutational, chemical, and computational analyses revealed the structural basis of selectivity, showing a small motif within peptide substrates that is sufficient for recognition by the enzyme. Attaching this 2-residue motif to two random peptides results in their isoprenylation by PagF, demonstrating utility as a general biocatalytic platform for modifications on any peptide substrate.

Parallelization of CPPTRAJ enables large scale analysis of molecular dynamics trajectory data.

Advances in biomolecular simulation methods and access to large scale computer resources have led to a massive increase in the amount of data generated. The key enablers have been optimization and parallelization of the simulation codes. However, much of the software used to analyze trajectory data from these simulations is still run in serial, or in some cases many threads via shared memory. Here, we describe the addition of multiple levels of parallel trajectory processing to the molecular dynamics simulation analysis software CPPTRAJ. In addition to the existing OpenMP shared-memory parallelism, CPPTRAJ now has two additional levels of message passing (MPI) parallelism involving both across-trajectory processing and across-ensemble processing. All three levels of parallelism can be simultaneously active, leading to significant speed ups in data analysis of large datasets on the NCSA Blue Waters supercomputer by better leveraging the many available nodes and its parallel file system. © 2018 Wiley Periodicals, Inc.

Mg2+ Binding Promotes SLV as a Scaffold in Varkud Satellite Ribozyme SLI-SLV Kissing Loop Junction.

Though the structure of the substrate stem loop I (SLI)-stem loop V (SLV) kissing loop junction of the Varkud Satellite ribozyme has been experimentally characterized, the dynamics of this Mg2+-dependent loop-loop interaction have been elusive. Specifically, each hairpin loop contains a U-turn motif, but only SLV shows a conformational shift triggered by Mg2+ ion association. Here, we use molecular dynamics simulations to analyze the binding and dynamics of this kissing loop junction. We show that SLV acts as a scaffold, providing stability to the junction. Mg2+ ions associate with SLV when it is part of the junction in a manner similar to when it is unbound, but there is no specificity in Mg2+ binding for the SLI loop. This suggests that the entropic penalty of ordering the larger SLI is too high, allowing SLV to act as a scaffold for multiple substrate loop sequences.

Probing the influence of hypermodified residues within the tRNA3(Lys) anticodon stem loop interacting with the A-loop primer sequence from HIV-1.

Replication of the HIV-1 virus requires reverse transcription of the viral RNA genome, a process that is specifically initiated by human tRNA3(Lys) packaged within the infectious virion. The primary binding site for the tRNA involves the 3' 18 nucleotides with an additional interaction between an adenine rich loop (A-loop) in the template and the anticodon stem-loop region of the tRNA3(Lys). The loop of the tRNA primer contains two hypermodified base residues and a pseudouridine that are required for a proper binding and activity. Here, we investigate the influence on the structure, dynamics and binding stability of the three modified residues (mnm(5)s(2)U34, t(6)A37 and Ψ39) using extensive molecular dynamics and Quantum Theory of Atoms in Molecules (QTAIM) analysis. Consistent with experiment, the results suggest that the three modified residues are required for faithful binding. Residues mnm(5)s(2)U34 and Ψ39 have a major influence in stabilizing the anticodon loop whereas mnm(5)s(2)U34 and t(6)A37 appear to stabilize the formation of the complex of tRNA3(Lys) with the HIV-1 A-loop.

Highly sampled tetranucleotide and tetraloop motifs enable evaluation of common RNA force fields.

Recent modifications and improvements to standard nucleic acid force fields have attempted to fix problems and issues that have been observed as longer timescale simulations have become routine. Although previous work has shown the ability to fold the UUCG stem-loop structure, until now no group has attempted to quantify the performance of current force fields using highly converged structural populations of the tetraloop conformational ensemble. In this study, we report the use of multiple independent sets of multidimensional replica exchange molecular dynamics (M-REMD) simulations with different initial conditions to generate well-converged conformational ensembles for the tetranucleotides r(GACC) and r(CCCC), as well as the larger UUCG tetraloop motif. By generating what is to our knowledge the most complete RNA structure ensembles reported to date for these systems, we remove the coupling between force field errors and errors due to incomplete sampling, providing a comprehensive comparison between current top-performing MD force fields for RNA. Of the RNA force fields tested in this study, none demonstrate the ability to correctly identify the most thermodynamically stable structure for all three systems. We discuss the deficiencies present in each potential function and suggest areas where improvements can be made. The results imply that although "short" (nsec-µsec timescale) simulations may stay close to their respective experimental structures and may well reproduce experimental observables, inevitably the current force fields will populate alternative incorrect structures that are more stable than those observed via experiment.

Application of Thiol-yne/Thiol-ene Reactions for Peptide and Protein Macrocyclizations.

The application of thiol-yne/thiol-ene reactions to synthesize mono- and bicyclic-stapled peptides and proteins is reported. First, a thiol-ene-based peptide-stapling method in aqueous conditions was developed. This method enabled the efficient stapling of recombinantly expressed coil-coiled proteins. The resulting stapled protein demonstrated higher stability in its secondary structure than the unstapled version. Furthermore, a thiol-yne coupling was performed by using an α,ω-diyne to react with two cysteine residues to synthesize a stapled peptide with two vinyl sulfide groups. The stapled peptide could further react with another biscysteine peptide to yield a bicyclic stapled peptide with enhanced properties. For example, the cell permeability of a stapled peptide was further increased by appending an oligoarginine cell-penetrating peptide. The robustness and versatility of thiol-yne/thiol-ene reactions that can be applied to both synthetic and expressed peptides and proteins were demonstrated.

Influence of BII Backbone Substates on DNA Twist: A Unified View and Comparison of Simulation and Experiment for All 136 Distinct Tetranucleotide Sequences.

Reliable representation of the B-DNA base-pair step twist is one of the crucial requirements for theoretical modeling of DNA supercoiling and other biologically relevant phenomena in B-DNA. It has long been suspected that the twist is inaccurately described by current empirical force fields. Unfortunately, comparison of simulation results with experiments is not straightforward because of the presence of BII backbone substates, whose populations may differ in experimental and simulation ensembles. In this work, we provide a comprehensive view of the effect of BII substates on the overall B-DNA helix twist and show how to reliably compare twist values from experiment and simulation in two scenarios. First, for longer DNA segments freely moving in solution, we show that sequence-averaged twists of different BI/BII ensembles can be compared directly because of approximate cancellation of the opposing BII effects. Second, for sequence-specific data, such as a particular base-pair step or tetranucleotide twist, can be compared only for a clearly defined BI/BII backbone conformation. For the purpose of force field testing, we designed a compact set of fourteen 22-base-pair B-DNA duplexes (Set 14) containing all 136 distinct tetranucleotide sequences and carried out a total of 84 µs of molecular dynamics simulations, primarily with the OL15 force field. Our results show that the ff99bsc0εζOL1χOL4, parmbsc1, and OL15 force fields model the B-DNA helical twist in good agreement with X-ray and minicircle ligation experiments. The comprehensive understanding obtained regarding the effect of BII substates on the base-pair step geometry should aid meaningful comparisons of various conformational ensembles in future research.

Intercalation processes of copper complexes in DNA.

The family of anticancer complexes that include the transition metal copper known as Casiopeínas® shows promising results. Two of these complexes are currently in clinical trials. The interaction of these compounds with DNA has been observed experimentally and several hypotheses regarding the mechanism of action have been developed, and these include the generation of reactive oxygen species, phosphate hydrolysis and/or base-pair intercalation. To advance in the understanding on how these ligands interact with DNA, we present a molecular dynamics study of 21 Casiopeínas with a DNA dodecamer using 10 µs of simulation time for each compound. All the complexes were manually inserted into the minor groove as the starting point of the simulations. The binding energy of each complex and the observed representative type of interaction between the ligand and the DNA is reported. With this extended sampling time, we found that four of the compounds spontaneously flipped open a base pair and moved inside the resulting cavity and four compounds formed stacking interactions with the terminal base pairs. The complexes that formed the intercalation pocket led to more stable interactions.

Convergence and reproducibility in molecular dynamics simulations of the DNA duplex d(GCACGAACGAACGAACGC).

BACKGROUND: The structure and dynamics of DNA are critically related to its function. Molecular dynamics simulations augment experiment by providing detailed information about the atomic motions. However, to date the simulations have not been long enough for convergence of the dynamics and structural properties of DNA. METHODS: Molecular dynamics simulations performed with AMBER using the ff99SB force field with the parmbsc0 modifications, including ensembles of independent simulations, were compared to long timescale molecular dynamics performed with the specialized Anton MD engine on the B-DNA structure d(GCACGAACGAACGAACGC). To assess convergence, the decay of the average RMSD values over longer and longer time intervals was evaluated in addition to assessing convergence of the dynamics via the Kullback-Leibler divergence of principal component projection histograms. RESULTS: These molecular dynamics simulations-including one of the longest simulations of DNA published to date at ~44µs-surprisingly suggest that the structure and dynamics of the DNA helix, neglecting the terminal base pairs, are essentially fully converged on the ~1-5µs timescale. CONCLUSIONS: We can now reproducibly converge the structure and dynamics of B-DNA helices, omitting the terminal base pairs, on the µs time scale with both the AMBER and CHARMM C36 nucleic acid force fields. Results from independent ensembles of simulations starting from different initial conditions, when aggregated, match the results from long timescale simulations on the specialized Anton MD engine. GENERAL SIGNIFICANCE: With access to large-scale GPU resources or the specialized MD engine "Anton" it is possible for a variety of molecular systems to reproducibly and reliably converge the conformational ensemble of sampled structures. This article is part of a Special Issue entitled: Recent developments of molecular dynamics.

Using Wavelet Analysis To Assist in Identification of Significant Events in Molecular Dynamics Simulations.

Long time scale molecular dynamics (MD) simulations of biological systems are becoming increasingly commonplace due to the availability of both large-scale computational resources and significant advances in the underlying simulation methodologies. Therefore, it is useful to investigate and develop data mining and analysis techniques to quickly and efficiently extract the biologically relevant information from the incredible amount of generated data. Wavelet analysis (WA) is a technique that can quickly reveal significant motions during an MD simulation. Here, the application of WA on well-converged long time scale (tens of µs) simulations of a DNA helix is described. We show how WA combined with a simple clustering method can be used to identify both the physical and temporal locations of events with significant motion in MD trajectories. We also show that WA can not only distinguish and quantify the locations and time scales of significant motions, but by changing the maximum time scale of WA a more complete characterization of these motions can be obtained. This allows motions of different time scales to be identified or ignored as desired.

Efficient treatment of induced dipoles.

Most existing treatments of induced dipoles in polarizable molecular mechanics force field calculations use either the self-consistent variational method, which is solved iteratively, or the "direct" approximation that is non-iterative as a result of neglecting coupling between induced dipoles. The variational method is usually implemented using assumptions that are only strictly valid under tight convergence of the induced dipoles, which can be computationally demanding to enforce. In this work, we discuss the nature of the errors that result from insufficient convergence and suggest a strategy that avoids such problems. Using perturbation theory to reintroduce the mutual coupling into the direct algorithm, we present a computationally efficient method that combines the precision of the direct approach with the accuracy of the variational approach. By analyzing the convergence of this perturbation series, we derive a simple extrapolation formula that delivers a very accurate approximation to the infinite order solution at the cost of only a few iterations. We refer to the new method as extrapolated perturbation theory. Finally, we draw connections to our previously published permanent multipole algorithm to develop an efficient implementation of the electric field and Thole terms and also derive some necessary, but not sufficient, criteria that force field parameters must obey.

Structural dynamics of possible late-stage intermediates in folding of quadruplex DNA studied by molecular simulations.

Explicit solvent molecular dynamics simulations have been used to complement preceding experimental and computational studies of folding of guanine quadruplexes (G-DNA). We initiate early stages of unfolding of several G-DNAs by simulating them under no-salt conditions and then try to fold them back using standard excess salt simulations. There is a significant difference between G-DNAs with all-anti parallel stranded stems and those with stems containing mixtures of syn and anti guanosines. The most natural rearrangement for all-anti stems is a vertical mutual slippage of the strands. This leads to stems with reduced numbers of tetrads during unfolding and a reduction of strand slippage during refolding. The presence of syn nucleotides prevents mutual strand slippage; therefore, the antiparallel and hybrid quadruplexes initiate unfolding via separation of the individual strands. The simulations confirm the capability of G-DNA molecules to adopt numerous stable locally and globally misfolded structures. The key point for a proper individual folding attempt appears to be correct prior distribution of syn and anti nucleotides in all four G-strands. The results suggest that at the level of individual molecules, G-DNA folding is an extremely multi-pathway process that is slowed by numerous misfolding arrangements stabilized on highly variable timescales.