INVESTIGATING THE EFFECTS OF IONIC LIQUIDS ON DNA GQUADRUPLEX AND PROTEIN STRUCTURE USING MOLECULAR DYNAMICS SIMULATIONS

Nucleic acids and proteins have huge implications in biomedicine and bioengineering, however their storage instability limits their applicability and current storage protocols are expensive and globally-inaccessible. Finding an alternative biocompatible media to store nucleic acids and proteins would reduce costs and increase their applicability. Ionic liquids (ILs) are molten salt compounds that have been shown to modulate the stability and activity of nucleic acids and proteins. In this thesis, molecular modeling studies of DNA/RNA and protein structure in ILs will be discussed (Chapter 1) and this method will be used to study the IL effects on the structure on the Pu22 c-MYC DNA G-quadruplex (Chapter 2) and the azurin protein (Chapter 3). ILs have been observed to stabilize/destabilize DNA G-quadruplexes linked to cancer oncogene expression, however the structural effects of imidazolium-based ILs on G-quadruplexes remain unknown. Bioengineering of azurin is attractive for soil bioremediation, thus understanding the structural changes induced by TMG amino acid-based ILs will mediate future IL design for enhancing azurin\u27s activity. In Chapter 2, molecular dynamics (MD) simulations will elucidate the stabilizing mechanism of four imidazolium-based ILs of increasing hydrophobicity to Pu22, using the G-quadruplex stabilizer TMPyP4 as a molecular probe. In Chapter 3, conventional and replica-exchange MD simulations will provide insight into the enthalpic and entropic change induced by two TMG-AA based ILs on the folded and unfolded azurin conformations

Shigella sonnei genome sequencing and phylogenetic analysis indicate recent global dissemination from Europe

Shigella are human-adapted Escherichia coli that have gained the ability to invade the human gut mucosa and cause dysentery1,2, spreading efficiently via low-dose fecal-oral transmission3,4. Historically, S. sonnei has been predominantly responsible for dysentery in developed countries, but is now emerging as a problem in the developing world, apparently replacing the more diverse S. flexneri in areas undergoing economic development and improvements in water quality4-6. Classical approaches have shown S. sonnei is genetically conserved and clonal7. We report here whole-genome sequencing of 132 globally-distributed isolates. Our phylogenetic analysis shows that the current S. sonnei population descends from a common ancestor that existed less than 500 years ago and has diversified into several distinct lineages with unique characteristics. Our analysis suggests the majority of this diversification occurred in Europe, followed by more recent establishment of local pathogen populations in other continents predominantly due to the pandemic spread of a single, rapidly-evolving, multidrug resistant lineage

Evolution of the insecticide target Rdl in African Anopheles is driven by interspecific and interkaryotypic introgression.

The evolution of insecticide resistance mechanisms in natural populations of Anopheles malaria vectors is a major public health concern across Africa. Using genome sequence data, we study the evolution of resistance mutations in the resistance to dieldrin locus (Rdl), a GABA receptor targeted by several insecticides, but most notably by the long-discontinued cyclodiene, dieldrin. The two Rdl resistance mutations (296G and 296S) spread across West and Central African Anopheles via two independent hard selective sweeps that included likely compensatory nearby mutations, and were followed by a rare combination of introgression across species (from A. gambiae and A. arabiensis to A. coluzzii) and across non-concordant karyotypes of the 2La chromosomal inversion. Rdl resistance evolved in the 1950s as the first known adaptation to a large-scale insecticide-based intervention, but the evolutionary lessons from this system highlight contemporary and future dangers for management strategies designed to combat development of resistance in malaria vectors

Structure and dynamics of the pan-genome of Streptococcus pneumoniae and closely related species

Background Streptococcus pneumoniae is one of the most important causes of microbial diseases in humans. The genomes of 44 diverse strains of S. pneumoniae were analyzed and compared with strains of non-pathogenic streptococci of the Mitis group. Results Despite evidence of extensive recombination, the S. pneumoniae phylogenetic tree revealed six major lineages. With the exception of serotype 1, the tree correlated poorly with capsular serotype, geographical site of isolation and disease outcome. The distribution of dispensable genes, genes present in not all, but more than one strain, was consistent with phylogeny, although horizontal gene transfer events attenuated this correlation in the case of ancient lineages. Homologous recombination, involving short stretches of DNA, was the dominant 13 evolutionary process of the core genome of S. pneumoniae. Genetic exchange occurred both within and across the borders of the species, and S. mitis was the main reservoir of genetic diversity of S. pneumoniae. The pan-genome size of S. pneumoniae increased logarithmically with the number of strains and linearly with the number of polymorphic sites of the sampled genomes, suggesting that acquired genes accumulate proportionately to the age of clones. Most genes associated with pathogenicity were shared by all S. pneumoniae strains, but were also present in S. mitis, S. oralis and S. infantis, indicating that these genes are not sufficient to determine virulence. Conclusion Genetic exchange with related species sharing the same ecological niche is the main mechanism of evolution of S. pneumoniae. The open pan genome guarantees the species a quick and economical response to diverse environments

Outlier SNPs detect weak regional structure against a background of genetic homogeneity in the Eastern Rock Lobster, Sagmariasus verreauxi

Genetic differentiation is characteristically weak in marine species making assessments of population connectivity and structure difficult. However, the advent of genomic methods has increased genetic resolution, enabling studies to detect weak, but significant population differentiation within marine species. With an increasing number of studies employing high resolution genome-wide techniques, we are realising that the connectivity of marine populations is often complex and quantifying this complexity can provide an understanding of the processes shaping marine species genetic structure and to inform long-term, sustainable management strategies. This study aims to assess the genetic structure, connectivity, and local adaptation of the Eastern Rock Lobster (Sagmariasus verreauxi), which has a maximum pelagic larval duration of 12 months and inhabits both subtropical and temperate environments. We used 645 neutral and 15 outlier SNPs to genotype lobsters collected from the only two known breeding populations and a third episodic population—encompassing S. verreauxi's known range. Through examination of the neutral SNP panel, we detected genetic homogeneity across the three regions, which extended across the Tasman Sea encompassing both Australian and New Zealand populations. We discuss differences in neutral genetic signature of S. verreauxi and a closely related, co-distributed rock lobster, Jasus edwardsii, determining a regional pattern of genetic disparity between the species, which have largely similar life histories. Examination of the outlier SNP panel detected weak genetic differentiation between the three regions. Outlier SNPs showed promise in assigning individuals to their sampling origin and may prove useful as a management tool for species exhibiting genetic homogeneity

Differing patterns of selection and geospatial genetic diversity within two leading Plasmodium vivax candidate vaccine antigens

Although Plasmodium vivax is a leading cause of malaria around the world, only a handful of vivax antigens are being studied for vaccine development. Here, we investigated genetic signatures of selection and geospatial genetic diversity of two leading vivax vaccine antigens--Plasmodium vivax merozoite surface protein 1 (pvmsp-1) and Plasmodium vivax circumsporozoite protein (pvcsp). Using scalable next-generation sequencing, we deep-sequenced amplicons of the 42 kDa region of pvmsp-1 (n = 44) and the complete gene of pvcsp (n = 47) from Cambodian isolates. These sequences were then compared with global parasite populations obtained from GenBank. Using a combination of statistical and phylogenetic methods to assess for selection and population structure, we found strong evidence of balancing selection in the 42 kDa region of pvmsp-1, which varied significantly over the length of the gene, consistent with immune-mediated selection. In pvcsp, the highly variable central repeat region also showed patterns consistent with immune selection, which were lacking outside the repeat. The patterns of selection seen in both genes differed from their P. falciparum orthologs. In addition, we found that, similar to merozoite antigens from P. falciparum malaria, genetic diversity of pvmsp-1 sequences showed no geographic clustering, while the non-merozoite antigen, pvcsp, showed strong geographic clustering. These findings suggest that while immune selection may act on both vivax vaccine candidate antigens, the geographic distribution of genetic variability differs greatly between these two genes. The selective forces driving this diversification could lead to antigen escape and vaccine failure. Better understanding the geographic distribution of genetic variability in vaccine candidate antigens will be key to designing and implementing efficacious vaccines

Effects of Ionic Liquids on Laccase from <i>Trametes versicolor</i>

Interactions between ionic liquids and biomolecules are of great interest due to the intrinsic properties of ionic liquids and the flexibility allowed by mixing and matching cations and anions to create unique ionic liquids. A number of ionic liquid–biomolecule studies have focused on interactions with proteins, including industrially relevant enzymes. One of these, laccase from Trametes versicolor, is a naturally derived enzyme used in the breakdown of phenolic compounds in a wide variety of industries, especially useful in breakdown of lignocellulosic materials. Here, a combination of experiments and molecular dynamics (MD) simulations was used to investigate the interactions of ionic liquids with laccase. Enzyme kinetics assays indicated that ionic liquids composed of tetramethylguanidine (TMG) and either serine or threonine caused significant reduction in enzymatic activity, while kinetics was not impacted by TMG-Asp or TMG-Glu ionic liquids. Similarly, intrinsic fluorescence of laccase in the presence of TMG-Ser and TMG-Thr exhibited a shift in spectral properties consistent with structural destabilization, but again TMG-Asp and TMG-Glu had no impact. MD simulations of laccase and ABTS with/without TMG-Ser ionic liquid provided insight into the deactivation mechanism of laccase. The simulations indicated that TMG-Ser disrupts laccase’s electron transfer mechanism

Inhibition Mechanism of Anti-TB Drug SQ109: Allosteric Inhibition of TMM Translocation of Mycobacterium Tuberculosis MmpL3 Transporter

The mycolic acid transporter MmpL3 is driven by proton motive forces (PMF) and functions via an antiport mechanism. Although the crystal structures of the Mycobacterium smegmatis MmpL3 transporter alone and in complex with a trehalose monomycolate (TMM) substrate and an antituberculosis drug candidate SQ109 under Phase 2b-3 Clinical Trials are available, no water and no conformational change in MmpL3 were observed in these structures to explain SQ109’s inhibition mechanism of proton and TMM transportation. In this study, molecular dynamics simulations of both apo form and inhibitor-bound MmpL3 in an explicit membrane were used to decipher the inhibition mechanism of SQ109. In the apo system, the close-open motion of the two TM domains, likely driven by the proton translocation, drives the close-open motion of the two PD domains, presumably allowing for TMM translocation. In contrast, in the holo system, the two PD domains are locked in a closed state, and the two TM domains are locked in an off pathway wider open state due to the binding of the inhibitor. Consistent with the close-open motion of the two PD domains, TMM entry size changes in the apo system, likely loading and moving the TMM, but does not vary much in the holo system and probably impair the movement of the TMM. Furthermore, we observed that water molecules passed through the central channel of the MmpL3 transporter to the cytoplasmic side in the apo system but not in the holo system, with a mean passing time of ∼135 ns. Because water wires play an essential role in transporting protons, our findings shed light on the importance of PMF in driving the close-open motion of the two TM domains. Interestingly, the key channel residues involved in water passage display considerable overlap with conserved residues within the MmpL protein family, supporting their critical function role

Inhibition Mechanism of Anti-TB Drug SQ109: Allosteric Inhibition of TMM Translocation of Mycobacterium Tuberculosis MmpL3 Transporter

The mycolic acid transporter MmpL3 is driven by proton motive forces (PMF) and functions via an antiport mechanism. Although the crystal structures of the Mycobacterium smegmatis MmpL3 transporter alone and in complex with a trehalose monomycolate (TMM) substrate and an antituberculosis drug candidate SQ109 under Phase 2b-3 Clinical Trials are available, no water and no conformational change in MmpL3 were observed in these structures to explain SQ109’s inhibition mechanism of proton and TMM transportation. In this study, molecular dynamics simulations of both apo form and inhibitor-bound MmpL3 in an explicit membrane were used to decipher the inhibition mechanism of SQ109. In the apo system, the close-open motion of the two TM domains, likely driven by the proton translocation, drives the close-open motion of the two PD domains, presumably allowing for TMM translocation. In contrast, in the holo system, the two PD domains are locked in a closed state, and the two TM domains are locked in an off pathway wider open state due to the binding of the inhibitor. Consistent with the close-open motion of the two PD domains, TMM entry size changes in the apo system, likely loading and moving the TMM, but does not vary much in the holo system and probably impair the movement of the TMM. Furthermore, we observed that water molecules passed through the central channel of the MmpL3 transporter to the cytoplasmic side in the apo system but not in the holo system, with a mean passing time of ∼135 ns. Because water wires play an essential role in transporting protons, our findings shed light on the importance of PMF in driving the close-open motion of the two TM domains. Interestingly, the key channel residues involved in water passage display considerable overlap with conserved residues within the MmpL protein family, supporting their critical function role