Mechanism of polypurine tract primer generation by HIV-1 reverse transcriptase

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.

May I Help You with Your Coat? HIV-1 Capsid Uncoating and Reverse Transcription

The human immunodeficiency virus type 1 (HIV-1) capsid is a protein core formed by multiple copies of the viral capsid (CA) protein. Inside the capsid, HIV-1 harbours all the viral components required for replication, including the genomic RNA and viral enzymes reverse transcriptase (RT) and integrase (IN). Upon infection, the RT transforms the genomic RNA into a double-stranded DNA molecule that is subsequently integrated into the host chromosome by IN. For this to happen, the viral capsid must open and release the viral DNA, in a process known as uncoating. Capsid plays a key role during the initial stages of HIV-1 replication; therefore, its stability is intimately related to infection efficiency, and untimely uncoating results in reverse transcription defects. How and where uncoating takes place and its relationship with reverse transcription is not fully understood, but the recent development of novel biochemical and cellular approaches has provided unprecedented detail on these processes. In this review, we present the latest findings on the intricate link between capsid stability, reverse transcription and uncoating, the different models proposed over the years for capsid uncoating, and the role played by other cellular factors on these processes.

Structural and computational studies of HIV-1 RNA

Viruses remain a global threat to animals, plants, and humans. The type 1 human immunodeficiency virus (HIV-1) is a member of the retrovirus family and carries an RNA genome, which is reverse transcribed into viral DNA and further integrated into the host-cell DNA for viral replication and proliferation. The RNA structures from the HIV-1 genome provide valuable insights into the mechanisms underlying the viral replication cycle. Moreover, these structures serve as models for designing novel therapeutic approaches. Here, we review structural data on RNA from the HIV-1 genome as well as computational studies based on these structural data. The review is organized according to the type of structured RNA element which contributes to different steps in the viral replication cycle. This is followed by an overview of the HIV-1 transactivation response element (TAR) RNA as a model system for understanding dynamics and interactions in the viral RNA systems. The review concludes with a description of computational studies, highlighting the impact of biomolecular simulations in elucidating the mechanistic details of various steps in the HIV-1's replication cycle.

HIV-1 3'-Polypurine Tract Mutations Confer Dolutegravir Resistance by Switching to an Integration-Independent Replication Mechanism via 1-LTR Circles

Several recent studies indicate that mutations in the human immunodeficiency virus type 1 (HIV-1) 3'polypurine tract (3'PPT) motif can reduce sensitivity to the integrase inhibitor dolutegravir (DTG). Using an in vivo systematic evolution of ligands by exponential enrichment (SELEX) approach, we discovered that multiple different mutations in this viral RNA element can confer DTG resistance, suggesting that the inactivation of this critical reverse transcription element causes resistance. An analysis of the viral DNA products formed upon infection by these 3'PPT mutants revealed that they replicate without integration into the host cell genome, concomitant with an increased production of 1-LTR circles. As the replication of these virus variants is activated by the human T-lymphotropic virus 1 (HTLV-1) Tax protein, a factor that reverses epigenetic silencing of episomal HIV DNA, these data indicate that the 3'PPT-mutated viruses escape from the integrase inhibitor DTG by switching to an integration-independent replication mechanism. IMPORTANCE The integrase inhibitor DTG is a potent inhibitor of HIV replication and is currently recommended in drug regimens for people living with HIV. Whereas HIV normally escapes from antiviral drugs by the acquisition of specific mutations in the gene that encodes the targeted enzyme, mutational inactivation of the viral 3'PPT sequence, an RNA element that has a crucial role in the viral reverse transcription process, was found to allow HIV replication in the presence of DTG in cell culture experiments. While the integration of the viral DNA into the cellular genome is considered one of the hallmarks of retroviruses, including HIV, 3'PPT inactivation caused integration-independent replication, which can explain the reduced DTG sensitivity. Whether this exotic escape route can also contribute to viral escape in HIV-infected persons remains to be determined, but our results indicate that screening for 3'PPT mutations in patients that fail on DTG therapy should be considered.

Conformational Heterogeneity of RNA Stem-Loop Hairpins Bound to FUS-RNA Recognition Motif with Disordered RGG Tail Revealed by Unbiased Molecular Dynamics Simulations

RNA-protein complexes use diverse binding strategies, ranging from structurally well-defined interfaces to completely disordered regions. Experimental characterization of flexible segments is challenging and can be aided by atomistic molecular dynamics (MD) simulations. Here, we used an extended set of microsecond-scale MD trajectories (400 μs in total) to study two FUS-RNA constructs previously characterized by nuclear magnetic resonance (NMR) spectroscopy. The FUS protein contains a well-structured RNA recognition motif domain followed by a presumably disordered RGG tail that binds RNA stem-loop hairpins. Our simulations not only provide several suggestions complementing the experiments but also reveal major methodological difficulties in studies of such complex RNA-protein interfaces. Despite efforts to stabilize the binding via system-specific force-field adjustments, we have observed progressive distortions of the RNA-protein interface inconsistent with experimental data. We propose that the dynamics is so rich that its converged description is not achievable even upon stabilizing the system. Still, after careful analysis of the trajectories, we have made several suggestions regarding the binding. We identify substates in the RNA loops, which can explain the NMR data. The RGG tail localized in the minor groove remains disordered, sampling countless transient interactions with the RNA. There are long-range couplings among the different elements contributing to the recognition, which can lead to allosteric communication throughout the system. Overall, the RNA-FUS systems form dynamical ensembles that cannot be fully represented by single static structures. Thus, albeit imperfect, MD simulations represent a viable tool to investigate dynamic RNA-protein complexes.

Structural basis of R-loop recognition by the S9.6 monoclonal antibody

R-loops are ubiquitous, dynamic nucleic-acid structures that play fundamental roles in DNA replication and repair, chromatin and transcription regulation, as well as telomere maintenance. The DNA-RNA hybrid–specific S9.6 monoclonal antibody is widely used to map R-loops. Here, we report crystal structures of a S9.6 antigen-binding fragment (Fab) free and bound to a 13-bp hybrid duplex. We demonstrate that S9.6 exhibits robust selectivity in binding hybrids over double-stranded (ds) RNA and in categorically rejecting dsDNA. S9.6 asymmetrically recognizes a compact epitope of two consecutive RNA nucleotides via their 2′-hydroxyl groups and six consecutive DNA nucleotides via their backbone phosphate and deoxyribose groups. Recognition is mediated principally by aromatic and basic residues of the S9.6 heavy chain, which closely track the curvature of the hybrid minor groove. These findings reveal the molecular basis for S9.6 recognition of R-loops, detail its binding specificity, identify a new hybrid-recognition strategy, and provide a framework for S9.6 protein engineering.

The S9.6 monoclonal antibody is widely used to map R-loops genome wide. Here, Bou-Nader et al., define the nucleic acid-binding specificity of S9.6 and report its crystal structures free and bound to a hybrid, which reveal the asymmetric recognition of the RNA and DNA strands and its A-form conformation.

DARTS: An Algorithm for Domain-Associated Retrotransposon Search in Genome Assemblies

Retrotransposons comprise a substantial fraction of eukaryotic genomes, reaching the highest proportions in plants. Therefore, identification and annotation of retrotransposons is an important task in studying the regulation and evolution of plant genomes. The majority of computational tools for mining transposable elements (TEs) are designed for subsequent genome repeat masking, often leaving aside the element lineage classification and its protein domain composition. Additionally, studies focused on the diversity and evolution of a particular group of retrotransposons often require substantial customization efforts from researchers to adapt existing software to their needs. Here, we developed a computational pipeline to mine sequences of protein-coding retrotransposons based on the sequences of their conserved protein domains—DARTS (Domain-Associated Retrotransposon Search). Using the most abundant group of TEs in plants—long terminal repeat (LTR) retrotransposons (LTR-RTs)—we show that DARTS has radically higher sensitivity for LTR-RT identification compared to the widely accepted tool LTRharvest. DARTS can be easily customized for specific user needs. As a result, DARTS returns a set of structurally annotated nucleotide and amino acid sequences which can be readily used in subsequent comparative and phylogenetic analyses. DARTS may facilitate researchers interested in the discovery and detailed analysis of the diversity and evolution of retrotransposons, LTR-RTs, and other protein-coding TEs.

Role of salt-bridging interactions in recognition of viral RNA by arginine-rich peptides

Interactions between RNA molecules and proteins are critical to many cellular processes and are implicated in various diseases. The RNA-peptide complexes are good model systems to probe the recognition mechanism of RNA by proteins. In this work, we report studies on the binding-unbinding process of a helical peptide from a viral RNA element using nonequilibrium molecular dynamics simulations. We explored the existence of various dissociation pathways with distinct free-energy profiles that reveal metastable states and distinct barriers to peptide dissociation. We also report the free-energy differences for each of the four pathways to be 96.47 ± 12.63, 96.1 ± 10.95, 91.83 ± 9.81, and 92 ± 11.32 kcal/mol. Based on the free-energy analysis, we further propose the preferred pathway and the mechanism of peptide dissociation. The preferred pathway is characterized by the formation of sequential hydrogen-bonding and salt-bridging interactions between several key arginine amino acids and the viral RNA nucleotides. Specifically, we identified one arginine amino acid (R8) of the peptide to play a significant role in the recognition mechanism of the peptide by the viral RNA molecule.

HIV-1 Reverse Transcriptase Inhibition by Major Compounds in a Kenyan Multi-Herbal Composition (CareVid™): In Vitro and In Silico Contrast

CareVid is a multi-herbal product used in southwest Kenya as an immune booster and health tonic and has been anecdotally described as improving the condition of HIV-positive patients. The product is made up of roots, barks and whole plant of 14 African medicinal plants: Acacia nilotica (L.) Willd. ex Delile (currently, Vachelia nilotica (L.) P.J.H Hurter & Mabb.), Adenia gummifera (Harv.) Harms, Anthocleista grandiflora Gilg, Asparagus africanus Lam., Bersama abyssinica Fresen., Clematis hirsuta Guill. & Perr., Croton macrostachyus Hochst. ex Delile, Clutia robusta Pax (accepted as Clutia kilimandscharica Engl.), Dovyalis abyssinica (A. Rich.) Warb, Ekebergia capensis Sparm., Periploca linearifolia Quart.-Dill. & A. Rich., Plantago palmata Hook.f., Prunus africana Hook.f. Kalkman and Rhamnus prinoides L’Her. The objective of this study was to determine the major chemical constituents of CareVid solvent extracts and screen them for in vitro and in silico activity against the HIV-1 reverse transcriptase enzyme. To achieve this, CareVid was separately extracted using CH₂Cl₂, MeOH, 80% EtOH in H₂O, cold H₂O, hot H₂O and acidified H₂O (pH 1.5–3.5). The extracts were analysed using HPLC–MS equipped with UV diode array detection. HIV-1 reverse transcriptase inhibition was performed in vitro and compared to in silico HIV-1 reverse transcriptase inhibition, with the latter carried out using MOE software, placing the docking on the hydrophobic pocket in the subdomain of p66, the NNRTI pocket. The MeOH and 80% EtOH extracts showed strong in vitro HIV-1 reverse transcriptase inhibition, with an EC₅₀ of 7 μg·mL⁻¹. The major components were identified as sucrose, citric acid, ellagic acid, catechin 3-hexoside, epicatechin 3-hexoside, procyanidin B, hesperetin O-rutinoside, pellitorine, mangiferin, isomangiferin, 4-O-coumaroulquinic acid, ellagic acid, ellagic acid O-pentoside, crotepoxide, oleuropein, magnoflorine, tremulacin and an isomer of dammarane tetrol. Ellagic acid and procyanidin B inhibited the HIV-1 reverse transcription process at 15 and 3.2 µg/mL⁻¹, respectively. Docking studies did not agree with in vitro results because the best scoring ligand was crotepoxide (ΔG = −8.55 kcal/mol), followed by magnoflorine (ΔG = −8.39 kcal/mol). This study showed that CareVid has contrasting in vitro and in silico activity against HIV-1 reverse transcriptase. However, the strongest in vitro inhibitors were ellagic acid and procyanidin B.

HIV Modifies the m6A and m5C Epitranscriptomic Landscape of the Host Cell

The study of RNA modifications, today known as epitranscriptomics, is of growing interest. The N6-methyladenosine (m6A) and 5-methylcytosine (m5C) RNA modifications are abundantly present on mRNA molecules, and impact RNA interactions with other proteins or molecules, thereby affecting cellular processes, such as RNA splicing, export, stability, and translation. Recently m6A and m5C marks were found to be present on human immunodeficiency (HIV) transcripts as well and affect viral replication. Therefore, the discovery of RNA methylation provides a new layer of regulation of HIV expression and replication, and thus offers novel array of opportunities to inhibit replication. However, no study has been performed to date to investigate the impact of HIV replication on the transcript methylation level in the infected cell. We used a productive HIV infection model, consisting of the CD4+ SupT1 T cell line infected with a VSV-G pseudotyped HIVeGFP-based vector, to explore the temporal landscape of m6A and m5C epitranscriptomic marks upon HIV infection, and to compare it to mock-treated cells. Cells were collected at 12, 24, and 36 h post-infection for mRNA extraction and FACS analysis. M6A RNA modifications were investigated by methylated RNA immunoprecipitation followed by high-throughput sequencing (MeRIP-Seq). M5C RNA modifications were investigated using a bisulfite conversion approach followed by high-throughput sequencing (BS-Seq). Our data suggest that HIV infection impacted the methylation landscape of HIV-infected cells, inducing mostly increased methylation of cellular transcripts upon infection. Indeed, differential methylation (DM) analysis identified 59 m6A hypermethylated and only 2 hypomethylated transcripts and 14 m5C hypermethylated transcripts and 7 hypomethylated ones. All data and analyses are also freely accessible on an interactive web resource (http://sib-pc17.unil.ch/HIVmain.html). Furthermore, both m6A and m5C methylations were detected on viral transcripts and viral particle RNA genomes, as previously described, but additional patterns were identified. This work used differential epitranscriptomic analysis to identify novel players involved in HIV life cycle, thereby providing innovative opportunities for HIV regulation.

Distinct Retrotransposon Evolution Profile in the Genome of Rabbit (Oryctolagus cuniculus)

Although the rabbit genome has already been annotated, it is mobilome remains largely unknown. Here, multiple pipelines were used to de novo mine and annotate the mobilome in rabbit. Four families and 19 subfamilies of LINE1s, two families and nine subfamilies of SINEs, and 12 ERV families were defined in rabbit based on sequence identity, structural organization, and phylogenetic tree. The analysis of insertion age and polymerase chain reaction suggests that a number of families are very young and may remain active, such as L1B, L1D, OcuSINEA, and OcuERV1. RepeatMasker annotation revealed a distinct transposable element landscape within the genome, with approximately two million copies of SINEs, representing the greatest proportion of the genome (19.61%), followed by LINEs (15.44%), and LTRs (4.11%), respectively, considerably different from most other mammal mobilomes except hedgehog and tree shrew, in which LINEs have the highest proportion. Furthermore, a very high rate of insertion polymorphisms (>85%) for the youngest subfamily (OcuSINEA1) was identified by polymerase chain reaction. The majority of retrotransposon insertions overlapped with protein-coding regions (>80%) and lncRNA (90%) genes. Genomic distribution bias was observed for retrotransposons, with those immediately upstream (-1 kb) and downstream (1 kb) of genes significantly depleted. Local GC content in 50-kb widows had significantly negative correlations with LINE (rs=-0.996) and LTR (rs=-0.829) insertions. The current study revealed a distinct mobilome landscape in rabbit, which will assist in the elucidation of the evolution of the genome of lagomorphs, and even other mammals.

Analysis and Molecular Determinants of HIV RNase H Cleavage Specificity at the PPT/U3 Junction

HIV reverse transcriptases (RTs) convert viral genomic RNA into double-stranded DNA. During reverse transcription, polypurine tracts (PPTs) resilient to RNase H cleavage are used as primers for plus-strand DNA synthesis. Nonnucleoside RT inhibitors (NNRTIs) can interfere with the initiation of plus-strand DNA synthesis by enhancing PPT removal, while HIV RT connection subdomain mutations N348I and N348I/T369I mitigate this effect by altering RNase H cleavage specificity. Now, we demonstrate that among approved nonnucleoside RT inhibitors (NNRTIs), nevirapine and doravirine show the largest effects. The combination N348I/T369I in HIV-1_BH10 RT has a dominant effect on the RNase H cleavage specificity at the PPT/U3 site. Biochemical studies showed that wild-type HIV-1 and HIV-2 RTs were able to process efficiently and accurately all tested HIV PPT sequences. However, the cleavage accuracy at the PPT/U3 junction shown by the HIV-2_EHO RT was further improved after substituting the sequence YQEPFKNLKT of HIV-1_BH10 RT (positions 342–351) for the equivalent residues of the HIV-2 enzyme (HQGDKILKV). Our results highlight the role of β-sheets 17 and 18 and their connecting loop (residues 342–350) in the connection subdomain of the large subunit, in determining the RNase H cleavage window of HIV RTs.

Disulfide bridge cross-linking between protein and the RNA backbone as a tool to study RNase H1

The chemical cross-linking of complexes of proteins with nucleic acids is often used in structural and mechanistic studies of these oftentimes unstable and transient complexes. To date, no method has been reported for the thiol-based conjugation of proteins with an RNA backbone, mainly because of instability of the modified ribonucleic acid that is functionalized at the phosphodiester and its rapid hydrolysis. Here, we report the site-specific synthesis of stable RNA oligonucleotides with a thiol-bearing linker that was attached to the phosphodiester backbone, where the ribonucleotide at the cross-linking site was either replaced with 2'-deoxy- or 2'-fluororibonucleotide. The utility of this approach was validated in cross-linking tests with RNase H1, a model protein for RNA/DNA binding and key effector in DNA-like antisense drug therapy. Furthermore, scale-up cross-linking and purification of the complexes confirmed that the method is useful for obtaining preparations of protein-RNA/DNA complexes with purity and stability that are suitable for further biochemical and structural studies. The present approach broadens the repertoire of disulfide-based cross-linking strategies and is a novel tool for the stabilization of protein-RNA complexes in which the interaction occurs via the RNA backbone. This methodology may be broadly applicable to studies of otherwise unstable or transient complexes of proteins with RNA and RNA/DNA.

GC-Content Dependence of Elastic and Overstretching Properties of DNA:RNA Hybrid Duplexes

DNA:RNA hybrid duplex plays important roles in various biological processes. Both its structural stability and interactions with proteins are highly sequence dependent. In this study, we utilize homebuilt optical tweezers to investigate how GC contents in the sequence influence the structural and mechanical properties of DNA:RNA hybrid by measuring its contour length, elasticities, and overstretching dynamics. Our results support that the DNA:RNA hybrid adopts a conformation between the A- and B-form helix, and the GC content does not affect its structural and elastic parameters obviously when varying from 40 to 60% before the overstretching transition of DNA:RNA hybrid occurs. In the overstretching transition, however, our study unravels significant heterogeneity and strong sequence dependence, suggesting the presence of a highly dynamic competition between the two processes, namely the S-form duplex formation (nonhysteretic) and the unpeeling (hysteretic). Analyzing the components left in DNA:RNA hybrid after the overstretching transition suggests that the RNA strand is more easily unpeeled than the DNA strand, whereas an increase in the GC content from 40 to 60% can significantly reduce unpeeling. Large hysteresis is observed between the stretching and relaxation processes, which is also quantitatively correlated with the percentage of unpeeling in the DNA:RNA duplex. Increasing in both the salt concentration and GC content can effectively reduce the hysteresis with the latter being more significant. Together, our study reveals that the mechanical properties of DNA:RNA hybrid duplexes are significantly different from double-stranded DNA and double-stranded RNA, and its overstretching behavior is highly sequence dependent. These results should be taken into account in the future studies on DNA:RNA-hybrid-related functional structures and motor proteins.

RuvC uses dynamic probing of the Holliday junction to achieve sequence specificity and efficient resolution

Holliday junctions (HJs) are four-way DNA structures that occur in DNA repair by homologous recombination. Specialized nucleases, termed resolvases, remove (i.e., resolve) HJs. The bacterial protein RuvC is a canonical resolvase that introduces two symmetric cuts into the HJ. For complete resolution of the HJ, the two cuts need to be tightly coordinated. They are also specific for cognate DNA sequences. Using a combination of structural biology, biochemistry, and a computational approach, here we show that correct positioning of the substrate for cleavage requires conformational changes within the bound DNA. These changes involve rare high-energy states with protein-assisted base flipping that are readily accessible for the cognate DNA sequence but not for non-cognate sequences. These conformational changes and the relief of protein-induced structural tension of the DNA facilitate coordination between the two cuts. The unique DNA cleavage mechanism of RuvC demonstrates the importance of high-energy conformational states in nucleic acid readouts.

RNases H: Structure and mechanism

RNases H are a family of endonucleases that hydrolyze RNA residues in various nucleic acids. These enzymes are present in all branches of life, and their counterpart domains are also found in reverse transcriptases (RTs) from retroviruses and retroelements. RNases H are divided into two main classes (RNases H1 and H2 or type 1 and type 2 enzymes) with common structural features of the catalytic domain but different range of substrates for enzymatic cleavage. Additionally, a third class is found in some Archaea and bacteria. Besides distinct cellular functions specific for each type of RNases H, this family of proteins is generally involved in the maintenance of genome stability with overlapping and cooperative role in removal of R-loops thus preventing their accumulation. Extensive biochemical and structural studies of RNases H provided not only a comprehensive and complete picture of their mechanism but also revealed key basic principles of nucleic acid recognition and processing. RNase H1 is present in prokaryotes and eukaryotes and cleaves RNA in RNA/DNA hybrids. Its main function is hybrid removal, notably in the context of R-loops. RNase H2, which is also present in all branches of life, can play a similar role but it also has a specialized function in the cleavage of single ribonucleotides embedded in the DNA. RNase H3 is present in Archaea and bacteria and is closely related to RNase H2 in sequence and structure but has RNase H1-like biochemical properties. This review summarizes the mechanisms of substrate recognition and enzymatic cleavage by different classes of RNases H with particular insights into structural features of nucleic acid binding, specificity towards RNA and/or DNA strands and catalysis.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.