A collaborative approach to improving representation in viral genomic surveillance.

The lack of routine viral genomic surveillance delayed the initial detection of SARS-CoV-2, allowing the virus to spread unfettered at the outset of the U.S. epidemic. Over subsequent months, poor surveillance enabled variants to emerge unnoticed. Against this backdrop, long-standing social and racial inequities have contributed to a greater burden of cases and deaths among minority groups. To begin to address these problems, we developed a new variant surveillance model geared toward building 'next generation' genome sequencing capacity at universities in or near rural areas and engaging the participation of their local communities. The resulting genomic surveillance network has generated more than 1,000 SARS-CoV-2 genomes to date, including the first confirmed case in northeast Louisiana of Omicron, and the first and sixth confirmed cases in Georgia of the emergent BA.2.75 and BQ.1.1 variants, respectively. In agreement with other studies, significantly higher viral gene copy numbers were observed in Delta variant samples compared to those from Omicron BA.1 variant infections, and lower copy numbers were seen in asymptomatic infections relative to symptomatic ones. Collectively, the results and outcomes from our collaborative work demonstrate that establishing genomic surveillance capacity at smaller academic institutions in rural areas and fostering relationships between academic teams and local health clinics represent a robust pathway to improve pandemic readiness.

A collaborative approach to improve representation in viral genomic surveillance.

The lack of routine viral genomic surveillance delayed the initial detection of SARS-CoV-2, allowing the virus to spread unfettered at the outset of the U.S. epidemic. Over subsequent months, poor surveillance enabled variants to emerge unnoticed. Against this backdrop, long-standing social and racial inequities have contributed to a greater burden of cases and deaths among minority groups. To begin to address these problems, we developed a new variant surveillance model geared toward building microbial genome sequencing capacity at universities in or near rural areas and engaging the participation of their local communities. The resulting genomic surveillance network has generated more than 1,000 SARS-CoV-2 genomes to date, including the first confirmed case in northeast Louisiana of Omicron, and the first and sixth confirmed cases in Georgia of the emergent BA.2.75 and BQ.1.1 variants, respectively. In agreement with other studies, significantly higher viral gene copy numbers were observed in Delta variant samples compared to those from Omicron BA.1 variant infections, and lower copy numbers were seen in asymptomatic infections relative to symptomatic ones. Collectively, the results and outcomes from our collaborative work demonstrate that establishing genomic surveillance capacity at smaller academic institutions in rural areas and fostering relationships between academic teams and local health clinics represent a robust pathway to improve pandemic readiness. Author summary: Genomic surveillance involves decoding a pathogen’s genetic code to track its spread and evolution. During the pandemic, genomic surveillance programs around the world provided valuable data to scientists, doctors, and public health officials. Knowing the complete SARS-CoV-2 genome has helped detect the emergence of new variants, including ones that are more transmissible or cause more severe disease, and has supported the development of diagnostics, vaccines, and therapeutics. The impact of genomic surveillance on public health depends on representative sampling that accurately reflects the diversity and distribution of populations, as well as rapid turnaround time from sampling to data sharing. After a slow start, SARS-CoV-2 genomic surveillance in the United States grew exponentially. Despite this, many rural regions and ethnic minorities remain poorly represented, leaving significant gaps in the data that informs public health responses. To address this problem, we formed a network of universities and clinics in Louisiana, Georgia, and Mississippi with the goal of increasing SARS-CoV-2 sequencing volume, representation, and equity. Our results demonstrate the advantages of rapidly sequencing pathogens in the same communities where the cases occur and present a model that leverages existing academic and clinical infrastructure for a powerful decentralized genomic surveillance system.

Genome Dashboards: Framework and Examples.

Genomics is a sequence-based informatics science and a three-dimensional-structure-based material science. However, in practice, most genomics researchers utilize sequence-based informatics approaches or three-dimensional-structure-based material science techniques, not both. This division is, at least in part, the result of historical developments rather than a fundamental necessity. The underlying computational tools, experimental techniques, and theoretical models were developed independently. The primary result presented here is a framework for the unification of informatics- and physics-based data associated with DNA, nucleosomes, and chromatin. The framework is based on the mathematical representation of geometrically exact rods and the generalization of DNA basepair step parameters. Data unification enables researchers to integrate computational, experimental, and theoretical approaches for the study of chromatin biology. The framework can be implemented using model-view-controller design principles, existing genome browsers, and existing molecular visualization tools. We developed a minimal, web-based genome dashboard, G-Dash-min, and applied it to two simple examples to demonstrate the usefulness of data unification and proof of concept. Genome dashboards developed using the framework and design principles presented here are extensible and customizable and are therefore more broadly applicable than the examples presented. We expect a number of purpose-specific genome dashboards to emerge as a novel means of investigating structure-function relationships for genomes that range from basepairs to entire chromosomes and for generating, validating, and testing mechanistic hypotheses.

Computational Design and Analysis of a Magic Snake.

The Magic Snake (Rubik's Snake) is a toy that was invented decades ago. It draws much less attention than Rubik's Cube, which was invented by the same professor, Erno Rubik. The number of configurations of a Magic Snake, determined by the number of discrete rotations about the elementary wedges in a typical snake, is far less than the possible configurations of a typical cube. However, a cube has only a single three-dimensional (3D) structure while the number of sterically allowed 3D conformations of the snake is unknown. Here, we demonstrate how to represent a Magic Snake as a one-dimensional (1D) sequence that can be converted into a 3D structure. We then provide two strategies for designing Magic Snakes to have specified 3D structures. The first enables the folding of a Magic Snake onto any 3D space curve. The second introduces the idea of "embedding" to expand an existing Magic Snake into a longer, more complex, self-similar Magic Snake. Collectively, these ideas allow us to rapidly list and then compute all possible 3D conformations of a Magic Snake. They also form the basis for multidimensional, multi-scale representations of chain-like structures and other slender bodies including certain types of robots, polymers, proteins, and DNA.

TMB Library of Nucleosome Simulations.

Nucleosomes are the fundamental building blocks of chromatin, the biomaterial that houses the genome in all higher organisms. A nucleosome consists of 145-147 base pairs of DNA wrapped 1.7 times around eight histones. Given a four-letter code (A, C, G, T), there are approximately 4147 or 1088 oligonucleotides that can form a nucleosome. Comparative, rather than comprehensive, studies are required. Here we introduce the TMB Library of nucleosome simulations and present a meta-analysis of over 20 µs of all atom molecular dynamics simulations representing 518 different realizations of the nucleosome. The TMB Library serves as a reference for future comparative, on-demand simulations of nucleosomes and a demonstration of iBIOMES Lite as a tool for managing a laboratory's simulation library. For every simulation, dewatered trajectories, RMSD, and DNA helical parameter data are provided through iBIOMES Lite in a Web browser and a file browser format. A novel view of nucleosomal DNA emerges from our meta-analysis of the TMB Library. DNA conformation is restricted to a specific left-handed superhelix, but the range of conformations observed for individual bases and base pairs is not more restricted nor more highly deformed than DNA free in solution. With the exception of Roll, mean DNA helical parameter values obtained from simulations of nucleosomes are largely within the range of thermal motion of DNA free in solution. The library provides evidence of DNA kinking in the nucleosome and clearly demonstrates the effects of DNA sequence on the gross structure and dynamics of nucleosomes. These effects and mispositioning of the 601 super strong nucleosome positioning sequence can be detected in short simulations (10 ns). Collectively, the results provide a basis for comparative simulation studies of nucleosomes and extend our understanding of the binding of proteins and drugs to nucleosomal DNA. The TMB Library can be found at http://dna.engr.latech.edu/~tmbshare/ .

µABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA.

We present the results of microsecond molecular dynamics simulations carried out by the ABC group of laboratories on a set of B-DNA oligomers containing the 136 distinct tetranucleotide base sequences. We demonstrate that the resulting trajectories have extensively sampled the conformational space accessible to B-DNA at room temperature. We confirm that base sequence effects depend strongly not only on the specific base pair step, but also on the specific base pairs that flank each step. Beyond sequence effects on average helical parameters and conformational fluctuations, we also identify tetranucleotide sequences that oscillate between several distinct conformational substates. By analyzing the conformation of the phosphodiester backbones, it is possible to understand for which sequences these substates will arise, and what impact they will have on specific helical parameters.

ICM Web: the interactive chromatin modeling web server.

The interactive chromatin modeling web server (ICM Web) is an interactive tool that allows users to rapidly assess nucleosome stability and fold sequences of DNA into putative chromatin templates. ICM Web takes a sequence composed of As, Cs, Gs, and Ts as input and generates (i) a nucleosome energy level diagram, (ii) coarse-grained representations of free DNA and chromatin and (iii) plots of the helical parameters (Tilt, Roll, Twist, Shift, Slide and Rise) as a function of position. The user can select from several different energy models, nucleosome structures and methods for placing nucleosomes in the energy landscape. Alternatively, if nucleosome footprints are known from experiment, ICM Web can use these positions to create a nucleosome array. The default energy model achieves a correlation coefficient of 0.7 with 100 experimentally determined values of stability and properly predicts the location of six positioned nucleosomes in the mouse mammary tumor virus (MMTV) promoter. ICM Web is suitable for interactively investigating nucleosome stability and chromatin folding for sequences up to tens of kilobases in length. No login is required to use ICM Web.

Evaluation of elastic rod models with long range interactions for predicting nucleosome stability.

The ability of a dinucleotide-step based elastic-rod model of DNA to predict nucleosome binding free energies is investigated using four available sets of elastic parameters. We compare the predicted free energies to experimental values derived from nucleosome reconstitution experiments for 84 DNA sequences. Elastic parameters (conformation and stiffnessess) obtained from MD simulations are shown to be the most reliable predictors, as compared to those obtained from analysis of base-pair step melting temperatures, or from analysis of x-ray structures. We have also studied the effect of varying the folded conformation of nucleosomal DNA by means of our Fourier - filtering knock-out and knock-in procedure. This study confirmed the above ranking of elastic parameters, and helped to reveal problems inherent in models using only a local elastic energy function. Long-range interactions were added to the elastic-rod model in an effort to improve its predictive ability. For this purpose a Debye-Huckel energy term with a single, homogenous point charge per base-pair was introduced. This term contains only three parameters, - its weight relative to the elastic energy, the Debye screening length, and a minimum sequence distance for including pairwise interactions between charges. After optimization of these parameters, our Debye-Huckel term is attractive, and yields the same level of correlation with experiment (R=0.75) as was achieved merely by varying the nucleosomal shape in the elastic-rod model. We suggest this result indicates a linker DNA - histone attraction or, possibly, entropic effects, that lead to a stabilization of a nucleosome away from the ends of DNA segments longer than 147 bp. Such effects are not accounted for by a localized elastic energy model.

A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA.

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein-DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50-100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA.

VDNA: the virtual DNA plug-in for VMD.

UNLABELLED: The DNA inter base pair step parameters (Tilt, Roll, Twist, Shift, Slide, Rise) are a standard internal coordinate representation of DNA. In the absence of bend and shear, it is relatively easy to mentally visualize how Twist and Rise generate the familiar double helix. More complex structures do not readily yield to such intuition. For this reason, we developed a plug-in for VMD that accepts a set of mathematical expressions as input and generates a coarse-grained model of DNA as output. This feature of VDNA appears to provide a unique approach to DNA modeling. Predefined expressions include: linear, sheared, bent and circular DNA, and models of the nucleosome superhelix, chromatin, thermal motion and nucleosome unwrapping. AVAILABILITY: VDNA is pre-installed in VMD, http://www.ks.uiuc.edu/Research/vmd. Updates are at http://dna.ccs.tulane.edu.

The nucleosome family: dynamic and growing.

Ever since the discovery of the nucleosome in 1974, scientists have stumbled upon discrete particles in which DNA is wrapped around histone complexes of different stoichiometries: octasomes, hexasomes, tetrasomes, "split" half-nucleosomes, and, recently, bona fide hemisomes. Do all these particles exist in vivo? Under what conditions? What is their physiological significance in the complex DNA transactions in the eukaryotic nucleus? What are their dynamics? This review summarizes research spanning more than three decades and provides a new meaning to the term "nucleosome." The nucleosome can no longer be viewed as a single static entity: rather, it is a family of particles differing in their structural and dynamic properties, leading to different functionalities.

Relaxation dynamics of nucleosomal DNA.

Recent experimental and theoretical evidence demonstrates that proteins and water in the hydration layer can follow complex stretched exponential or power law relaxation dynamics. Here, we report on a 50 ns all atom molecular dynamics (MD) simulation of the yeast nucleosome, where the interactions between DNA, histones, surrounding water and ions are explicitly included. DNA interacts with the histone core in 14 locations, approximately every 10.4 base pairs. We demonstrate that all sites of interaction exhibit anomalously slow power law relaxation, extending up to 10 ns, while fast exponential relaxation dynamics of hundreds of picoseconds applies to DNA regions outside these locations. The appearance of 1/f(alpha) noise or pink noise in DNA dynamics is ubiquitous. For histone-bound nucleotide dynamics alpha --> 1 and is a signature of complexity of the protein-DNA interactions. For control purposes two additional DNA simulations free of protein are conducted. Both utilize the same sequence of DNA, as found the in the nucleosome. In one simulation the initial conformation of the double helix is a straight B-form. In the other, the initial conformation is super helical. Neither of these simulations exhibits the variation of alpha as a function of position, the measure of power law for dynamical behavior, which we observe in the nucleosome simulation. The unique correspondence (high alpha to DNA-histone interaction sites, low alpha to free DNA sites), suggests that alpha may be an important and new quantification of protein-DNA interactions for future experiments.

Geometry of the nucleosomal DNA superhelix.

Nucleosome stability is largely an indirect measure of DNA sequence based on the material properties of DNA and the ability of a sequence to assume the required left-handed superhelical conformation. Here we focus attention only on the geometry of the superhelix and present two distinct mathematical expressions that rely on the DNA helical parameters (Shift, Slide, Rise, Tilt, Roll, Twist). One representation requires torsion for superhelix formation; the other requires shear. To compare these mathematical expressions to experimental data we develop a strategy for Fourier-filtering the helical parameters that identifies necessary and sufficient conditions to achieve a high-resolution model of the nucleosome superhelix. We apply this filtering strategy to 24 high-resolution structures of the nucleosome and demonstrate that all structures have a highly conserved distribution of Roll, Slide and Twist that involves two length scales. One length scale spans the entire length of nucleosomal DNA. The other is associated with the helix repeat. Our strategy also enables us to identify ground state or simple nucleosomes and altered nucleosome structures. These results form a basis for characterizing structural variations in the emerging family of nucleosome structures and a method for further developing structure-based models of nucleosome stability.

Molecular dynamics simulations of a nucleosome and free DNA.

All atom molecular dynamics simulations (10ns) of a nucleosome and of its 146 basepairs of DNA free in solution have been conducted. DNA helical parameters (Roll, Tilt, Twist, Shift, Slide, Rise) were extracted from each trajectory to compare the conformation, effective force constants, persistence length measures, and fluctuations of nucleosomal DNA to free DNA. The conformation of DNA in the nucleosome, as determined by helical parameters, is found to be largely within the range of thermally accessible values obtained for free DNA. DNA is found to be less flexible on the nucleosome than when free in solution, however such measures are length scale dependent. A method for disassembling and reconstructing the conformation and dynamics of the nucleosome using Fourier analysis is presented. Long length variations in the conformation of nucleosomal DNA are identified other than those associated with helix repeat. These variations are required to create a proposed tetrasome conformation or to qualitatively reconstruct the 1.75 turns of the nucleosome's superhelix. Reconstruction of free DNA using selected long wavelength variations in conformation can produce either a left-handed or a right-handed superhelix. The long wavelength variations suggest 146 basepairs is a natural length of DNA to wrap around the histone core.

Identification of important residues in metal-chelate recognition by monoclonal antibodies.

The molecular characterization of antibodies directed against metal-chelate complexes will provide important insights into the design and development of radiotherapeutic and radioimaging reagents. In this study, two monoclonal antibodies directed against different metal-chelate complexes were expressed as recombinant Fab fragments. Covalent modification and site-directed mutagenesis were employed to ascertain those residues important in antigen recognition. Antibody 5B2 was raised to a Pb(II)-loaded isothiocyanatobenzyl-diethylenetriamine pentaacetic acid (DTPA)-protein conjugate. The native antibody bound to complexes of Pb(II)-p-aminobenzyl-DTPA with an affinity of 4.6 x 10(-9) M. A monovalent Fab fragment prepared from the native protein and a bivalent recombinant fragment exhibited comparable affinities for the same Pb(II)-chelate complex, approximately 6-fold lower than that of the intact antibody. Covalent modification and molecular modeling predicted that Lys(58) in the heavy chain contacted the Pb(II)-chelate ligand. Mutational analysis supported a role for Lys(58) in ion pair or hydrogen bond formation with the carboxylate groups on the chelate. Antibody E5 was directed toward an isothiocyanatobenzyl-ethylenediamine tetraacetic acid (EDTA)-protein conjugate loaded with ionic Cd(II). The native immunoglobulin recognized Cd(II)-p-aminobenzyl-EDTA with an affinity of 8.2 x 10(-12) M. A proteolytically derived fragment and a bivalent recombinant fragment bound to the same Cd(II)-chelate complex with affinities that were comparable to that of the native antibody. Homology modeling and mutagenesis identified three residues (Trp(52) and His(96) in the heavy chain and Arg(96) in the light chain) that were important for Cd(II)-chelate recognition. His(96) likely mediates a direct ligation to the Cd(II) ion and Trp(52) appears to be involved in hydrophobic stacking with the benzyl moiety of the chelator. Arg(96) appeared to mediate an electrostatic or hydrogen bond to the chelate portion of the complex. These studies demonstrate that antibody recognition of metal-chelate haptens occurs through a limited number of molecular contacts and that these molecular interactions involve both direct ligation between the antibody and the metal ion and interactions between the antibody and the chelator.

Enzymatic E-colicins bind to their target receptor BtuB by presentation of a small binding epitope on a coiled-coil scaffold.

Toxins and viruses often initiate their attacks by binding to specific proteins on the surfaces of target cells. Bacterial toxins (e.g. bacteriocins) and viruses (bacteriophages) targeting Gram-negative bacteria typically bind to outer membrane proteins. Bacterial E-colicins target Escherichia coli by binding to the outer membrane cobalamin transporter BtuB. Colicins are tripartite molecules possessing receptor-binding, translocation, and toxin domains connected by long coiled-coil alpha-helices. Surprisingly, the crystal structure of colicin E3 does not possess a recognizable globular fold in its receptor-binding domain. We hypothesized that the binding epitope of enzymatic E-colicins is a short loop connecting the two alpha-helices that comprise the coiled-coil region and that this flanking coiled-coil region serves to present the loop in a binding-capable conformation. To test this hypothesis, we designed and synthesized a 34-residue peptide (E-peptide-1) corresponding to residues Ala366-Arg399 of the helix-loop-helix region of colicin E3. Cysteines placed near the ends of the peptide (I372C and A393C) enabled crosslinking for reduction of conformational entropy and formation of a peptide structure that would present the loop epitope. A fluorescent analog was also made for characterization of binding by measurement of fluorescence polarization. Our analysis shows the following. (i). E-peptide-1 is predominantly random coil in aqueous solution, but disulfide bond formation increases its alpha-helical content in both aqueous buffer and solvents that promote helix formation. (ii). Fluorescein-labeled E-peptide-1 binds to purified BtuB in a calcium-dependent manner with a Kd of 43.6 +/- 4.9 nm or 2370 +/- 670 nm in the presence or absence of calcium, respectively. (iii). In the presence of calcium, cyanocobalamin (CN-Cbl) displaces E-peptide-1 with a nanomolar inhibition constant (Ki = 78.9 +/- 5.6 nm). We conclude that the BtuB binding sites for cobalamins and enzymatic E-colicins are overlapping but inequivalent and that the distal loop and (possibly) the short alpha-helical flanking regions are sufficient for high affinity binding.

Mechanical model of the nucleosome and chromatin.

A theoretical framework for evaluating the approximate energy and dynamic properties associated with the folding of DNA into nucleosomes and chromatin is presented. Experimentally determined elastic constants of linear DNA and a simple fold geometry are assumed in order to derive elastic constants for extended and condensed chromatin. The model predicts the Young s modulus of extended and condensed chromatin to within an order of magnitude of experimentally determined values. Thus we demonstrate that the elastic properties of DNA are a primary determinant of the elastic properties of the higher order folded states. The derived elastic constants are used to predict the speed of propagation of small amplitude waves that excite an extension(sound), twist, bend or shear motion in each folded state. Taken together the results demonstrate that folding creates a hierarchy of time, length and energy scales.