Browsing by Author "Lemkul, Justin A."

Now showing 1 - 20 of 24
  • Thumbnail Image
    A Bioinformatics Approach to Identifying Radical SAM (S-Adenosyl-L-Methionine) Enzymes
    Gagliano, Elisa (Virginia Tech, 2020-06-03)
    Radical SAM enzymes are ancient, essential enzymes. They perform radical chemical reactions in virtually all living organisms and are involved in producing antibiotics, generating greenhouse gases, human health, and likely many other essential roles that have yet to be established. A wide variety of reactions have been characterized from this group of enzymes, including hydrogen abstractions, the transferring of methylthio groups, complex cyclization and rearrangement reactions, and others. However, many radical SAM enzymes have yet to be identified or characterized. There have been great leaps forward in the amount of enzyme sequences that are available in public databases, but experiments to investigate what chemical reactions the enzymes perform take a great deal of time. In our work, we utilize Hidden Markov Models to identify possible radical SAM enzymes and predict their possible functions through BLAST alignments and homology modelling. We also explore their distribution across the tree of life and determine how it is correlated with organism oxygen tolerances, because the core iron-sulfur cluster is oxygen sensitive. Trends in the abundances of radical SAM enzymes depending on oxygen tolerances were more apparent in prokaryotes than in eukaryotes. Although eukaryotes tend to have fewer radical SAM enzymes than prokaryotes, we were able to analyze uncharacterized radical SAM enzymes from both an aerobic eukaryote (Entamoeba histolytica) and a eukaryote capable of oxygenic photosynthesis (Gossypium barbadense), and predict the reactions they catalyze. This work sets the stage for the functional characterization of these essential yet elusive enzymes in future laboratory experiments.
  • Thumbnail Image
    Cation competition and recruitment around the c-kit1 G-quadruplex using polarizable simulations
    Salsbury, Alexa M.; Lemkul, Justin A. (2021-06-01)
    Nucleic acid-ion interactions are fundamentally important to the physical, energetic, and conformational properties of DNA and RNA. These interactions help fold and stabilize highly ordered secondary and tertiary structures, such as G-quadruplexes (GQs), which are functionally relevant in telomeres, replication initiation sites, and promoter sequences. The c-kit protooncogene encodes for a receptor tyrosine kinase and is linked to gastrointestinal stromal tumors, mast cell disease, and leukemia. This gene contains three unique GQ-forming sequences that have proposed antagonistic effects on gene expression. The dominant GQ, denoted c-kit1, has been shown to decrease expression of c-kit transcripts, making the c-kit1GQa promising drug target. Toward disease intervention, more information is needed regarding its conformational dynamics and ion binding properties. Therefore, we performed molecular dynamics simulations of the c-kit1 GQ with K+, Na+, Li+, and mixed salt solutions using the Drude-2017 polarizable force field. We evaluated GQ structure, ion sampling, core energetics, ion dehydration and binding, and ion competition and found that each analysis supported the known GQ-ion specificity trend (K+ > Na+ > Li+). We also found that K+ ions coordinate in the tetrad core antiprismatically, whereas Na+ and Li+ align coplanar to guanine tetrads, partially because of their attraction to surrounding water. Further, we showed that K+ occupancy is higher around the c-kit1 GQ and its nucleobases than Na+ and Li+, which tend to interact with backbone and sugar moieties. Finally, we showed that K+ binding to the c-kit1GQ is faster and more frequent than Na+ and Li+. Such descriptions of GQ-ion dynamics suggest the rate of dehydration as the dominant factor for preference of K+ by DNA GQs and provide insight into noncanonical nucleic acids for which little experimental data exist.
  • Thumbnail Image
    CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field
    Kognole, Abhishek A.; Lee, Jumin; Park, Sang-Jun; Jo, Sunhwan; Chatterjee, Payal; Lemkul, Justin A.; Huang, Jing; MacKerell, Alexander D.; Im, Wonpil (Wiley, 2021-12-07)
    Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. In addition, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. The capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.
  • Thumbnail Image
    charmm2gmx: An Automated Method to Port the CHARMM Additive Force Field to GROMACS
    Wacha, Andras F.; Lemkul, Justin A. (American Chemical Society, 2023-07-03)
    CHARMM is one of the most widely used biomolecular force fields. Although developed in close connection with a dedicated molecular simulation engine of the same name, it is also usable with other codes. GROMACS is a well-established, highly optimized, and multipurpose software for molecular dynamics, versatile enough to accommodate many different force field potential functions and the associated algorithms. Due to conceptional differences related to software design and the large amount of numeric data inherent to residue topologies and parameter sets, conversion from one software format to another is not straightforward. Here, we present an automated and validated means to port the CHARMM force field to a format read by the GROMACS engine, harmonizing the different capabilities of the two codes in a self-documenting and reproducible way with a bare minimum of user interaction required. Being based entirely on the upstream data files, the presented approach does not involve any hard-coded data, in contrast with previous attempts to solve the same problem. The heuristic approach used for perceiving the local internal geometry is directly applicable for analogous transformations of other force fields.
  • Thumbnail Image
    Complete Genome Sequence of Fusobacterium necrophorum subsp. necrophorum ATCC 25286
    Umaña, Ariana; Lemkul, Justin A.; Slade, Daniel J. (American Society for Microbiology, 2019-02-21)
    Fusobacterium necrophorum is a pathogenic Gram-negative, anaerobic bacterium. In this study, we present the first complete genome sequence of Fusobacterium necrophorum subsp. necrophorum ATCC 25286. These data provide a critical advancement in our understanding of virulence factors that could contribute to F. necrophorum pathogenesis in both human and livestock infections.
  • Thumbnail Image
    CT610: A Mn-Dependent Self-Sacrificing Oxygenase in p-Aminobenzoate Biosynthesis in Chlamydia trachomatis
    Wooldridge, Rowan Scott (Virginia Tech, 2022-06-09)
    Folate is an essential cofactor required for several processes including DNA and amino acid biosynthesis. Folate molecules are made up of three parts: a pteridine ring, p-aminobenzoate (pABA), and a variable number of glutamate residues. Chlamydia trachomatis synthesizes folate de novo; however, several genes encoding enzymes required for the canonical folate biosynthesis pathway are missing, including pabA/B and pabC, which are normally required for pABA biosynthesis from chorismate. Previous studies have found that a single gene in C. trachomatis, CT610, functionally replaces the canonical pABA biosynthesis genes. Interestingly, CT610 does not use chorismate as a substrate. Instead, the CT610-route for pABA biosynthesis incorporates isotopically labeled tyrosine into the synthesized pABA molecule. However, in vitro experiments revealed that CT610 produces pABA without any added substrates (including tyrosine) in the presence of a reducing agent and molecular oxygen. CT610 shares low sequence similarity to non-heme diiron oxygenases and the previously solved crystal structure revealed a diiron active site. Taken together, CT610 is proposed to be a novel self-sacrificing enzyme that uses one of its active site tyrosine residues as a precursor to pABA in a reaction that requires O2 and a reduced metallocofactor. Here, we discuss our progress towards understanding CT610-catalyzed pABA synthesis. Upon investigation of the pABA production and oxygenase activities of several active site tyrosine to phenylalanine variants, we found that Y27 and/or Y43 are the most likely precursors to the resulting pABA molecule. Further, activity was nearly completely abolished with a K152R variant, suggesting that this conserved lysine may be the required amino group donor. We also developed an in vitro Fe(II) reconstitution procedure, where the reconstituted enzyme exhibited a drastic increase in oxygenase activity but, surprisingly, a significant decrease in pABA synthase activity. Interestingly, a significant increase in pABA synthase activity was observed when the enzyme was reconstituted with manganese as opposed to iron, suggesting that the diiron active site of this enzyme might not be directly involved in CT610-dependent production of pABA and instead Mn may be the actual cofactor. Finally, we show that two 18O atoms from molecular oxygen are incorporated into the pABA molecule when synthesized by Mn-reconstituted CT610, providing further evidence for the oxygenase activity of CT610 and supporting our proposed mechanism that involves two monooxygenase reactions.
  • Thumbnail Image
    Differences in Conformational Sampling and Intrinsic Electric Fields Drive Ion Binding in Telomeric and TERRA G-Quadruplexes
    Poleto, Marcelo D.; Lemkul, Justin A. (American Chemical Society, 2023-10-17)
    The formation of G-quadruplexes (GQs) occurs in guanine-rich sequences of DNA and RNA, producing highly stable and structurally diverse noncanonical nucleic acid structures. GQs play crucial roles in regulating transcription, translation, and replication and maintaining the genome, among others; thus, changes to their structures can lead to diseases such as cancer. Previous studies using polarizable molecular dynamics simulations have shown differences in ion binding properties between telomeric and telomeric repeat-containing RNA GQs despite architectural similarities. Here, we used volume-based metadynamics and repulsive potential simulations in conjunction with polarizable force fields to quantify the impact of ion binding on the GQ dynamics and ion binding free energies. Furthermore, we describe how GQs exert electric fields on their surroundings to link dynamics with variations in the electronic structure. Our findings provide new insights into the energetic, physical, and conformational properties of GQs and expose subtle but important differences between DNA and RNA GQs with the same fold.
  • Thumbnail Image
    Effects of Familial Alzheimer's Disease Mutations on the Folding Free Energy and Dipole-Dipole Interactions of the Amyloid β-Peptide
    Davidson, Darcy S.; Kraus, Joshua A.; Montgomery, Julia M.; Lemkul, Justin A. (American Chemical Society, 2022-10-06)
    Familial Alzheimer's disease (FAD) mutations of the amyloid β-peptide (Aβ) are known to lead to early onset and more aggressive Alzheimer's disease. FAD mutations such as "Iowa" (D23N), "Arctic" (E22G), "Italian" (E22K), and "Dutch" (E22Q) have been shown to accelerate Aβ aggregation relative to the wild-type (WT). The mechanism by which these mutations facilitate increased aggregation is unknown, but each mutation results in a change in the net charge of the peptide. Previous studies have used nonpolarizable force fields to study Aβ, providing some insight into how this protein unfolds. However, nonpolarizable force fields have fixed charges that lack the ability to redistribute in response to changes in local electric fields. Here, we performed polarizable molecular dynamics simulations on the full-length Aβ42of WT and FAD mutations and calculated folding free energies of the Aβ15-27fragment via umbrella sampling. By studying both the full-length Aβ42and a fragment containing mutations and the central hydrophobic cluster (residues 17-21), we were able to systematically study how these FAD mutations impact secondary and tertiary structure and the thermodynamics of folding. Electrostatic interactions, including those between permanent and induced dipoles, affected side-chain properties, salt bridges, and solvent interactions. The FAD mutations resulted in shifts in the electronic structure and solvent accessibility at the central hydrophobic cluster and the hydrophobic C-terminal region. Using umbrella sampling, we found that the folding of the WT and E22 mutants is enthalpically driven, whereas the D23N mutant is entropically driven, arising from a different unfolding pathway and peptide-bond dipole response. Together, the unbiased, full-length, and umbrella sampling simulations of fragments reveal that the FAD mutations perturb nearby residues and others in hydrophobic regions to potentially alter solubility. These results highlight the role electronic polarizability plays in amyloid misfolding and the role of heterogeneous microenvironments that arise as conformational change takes place.
  • Electronic Polarization at the Interface between the p53 Transactivation Domain and Two Binding Partners
    Corrigan, Alexsandra N.; Lemkul, Justin A. (American Chemical Society, 2022-07-07)
    Intrinsically disordered proteins (IDPs) are an abundant class of highly charged proteins that participate in numerous crucial biological processes, often in regulatory roles. IDPs do not have one major free energy minimum with a dominant structure, instead existing as conformational ensembles of multiple semistable conformations. p53 is a prototypical protein with disordered regions and binds to many structurally diverse partners, making it a useful model for exploring the role of electrostatic interactions at IDP binding interfaces. In this study, we used the Drude-2019 force field to simulate the p53 transactivation domain with two protein partners to probe the role of electrostatic interactions in IDP protein-protein interactions. We found that the Drude-2019 polarizable force field reasonably reproduced experimental chemical shifts of the p53 transactivation domain (TAD) in one complex for which these data are available. We also found that the proteins in these complexes displayed dipole response at specific residues of each protein and that residues primarily involved in binding showed a large percent change in dipole moment between the unbound and complexed states. Probing the role of electrostatic interactions in IDP binding can allow us greater fundamental understanding of these interactions and may help with targeting p53 or its partners for drug design.
  • Thumbnail Image
    Electrostatic properties at the interface of p53 Transactivation domain binding
    Corrigan, Alexsandra Nikol (Virginia Tech, 2021-05-25)
    Intrinsically disordered proteins (IDPs) are an abundant class of proteins and protein regions which rapidly change between multiple structures without an equilibrium position. IDPs exist as a series of conformational ensembles of semi-stable conformations that can be adopted based on a hilly landscape of shallow free energy minima. Disordered sequences share characteristic features differentiating them from globular proteins, including low sequence complexity, low occurrence of hydrophobic residues, high polar and charged residue content, and high flexibility. IDPs are commonly involved in regulation in the cell, and frequently function as, or interact with, hub proteins in protein-protein interaction networks, making them an important class of macromolecules for understanding regulatory and other processes. Given their functional importance, these proteins are widely studied. Many analytical techniques are used, though rapid conformational sampling by IDPs makes it difficult to detect many states with NMR or other techniques. Computational approaches such as molecular dynamics are increasingly used to probe the binding and conformational sampling of these proteins, allowing for observation of factors that cannot be observed with traditional analytical methods such as NMR, such as differing conformational ensembles and the dipoles of individual residues. Here, we studied the role of electrostatic interactions in IDP protein-protein interaction using molecular dynamics simulations performed with the Drude-2019 force field (FF), a polarizable model that allows for more accurate representation of electrostatics, an important factor for highly charged systems like IDPs. For this project, a prototypical protein with disordered regions, p53, was simulated with two protein partners, the nuclear coactivator domain of the CREB binding protein (CBP), and the E3 ubiquitin-protein ligase mouse double minute 2 (MDM2). p53 is widely studied, and the p53 transactivation domain (TAD) is disordered and binds to many structurally diverse partners, making this protein domain a useful model for probing the role of electrostatic interactions formed by IDPs at protein-protein binding interfaces. We found that the Drude-2019 FF allows for simulation of the p53 TAD with Cα chemical shifts comparable to those observed with NMR, supporting that the Drude-2019 FF performs well in simulating IDPs. We observed large relative change in sidechain dipole moments when comparing the p53 TAD alone and when bound to either CBP or MDM2. We observed that aliphatic and aromatic amino acids experienced the largest relative change in sidechain dipole moments, and that there is sensitivity to binding shown in this dipole response. The largest percent changes in sidechain dipole moment were found to localize at and around the binding interface. Understanding the binding interactions of IDPs at a fundamental level, including the role of electrostatic interactions, may help with targeting IDPs or their partners for drug design.
  • Thumbnail Image
    Exploring the Forces Underlying the Dynamics and Energetics of G-quadruplexes with Polarizable Molecular Dynamics Simulations
    Salsbury, Alexa Marie (Virginia Tech, 2021-05-24)
    G-quadruplexes (GQs) are highly stable noncanonical nucleic acid structures that form in the DNA of human cells and play fundamental roles in maintaining genomic stability and regulating gene expression. These unique structures exert broad influence over biologically important processes and can modulate cell survival and human health. In fact, mutations, hyper-stability, and dissociation of GQs are implicated in neurodegenerative disease, mental retardation, premature-aging conditions, and various cancers. As such, GQs are novel drug targets. GQ-targeting therapeutics are developed to influence the folding and genetic interactions of GQs that are implicated in diseased states. To do so requires a greater understanding of GQ structure and dynamics and molecular dynamics (MD) simulations are well suited to provide these fundamental insights. Previous MD simulations of GQs have provided limited information due to inaccuracies in their models, namely the nonpolarizable nature of their force fields (FFs). The cutting-edge Drude polarizable FF models electronic degrees of freedom, allowing charge distribution to change in response to its environment. This is an important component for modeling ion-ion and ion-DNA interactions and can influence the overall stability of GQ structures. The work herein employs the Drude polarizable FF to 1) describe the role of electronic structure on the dynamics and folded stability of GQs, 2) determine the impact of ion interaction on GQ stability, and 3) characterize the role of G-hairpin motifs in GQ intermediates. Such fundamental investigations will help clarify GQs role in healthy and diseased states and transform our understanding of noncanonical DNA, improving human health, therapeutic design, and fundamental science.
  • Thumbnail Image
    FusoPortal: an Interactive Repository of Hybrid MinIONSequenced Fusobacterium Genomes Improves Gene Identification and Characterization
    Sanders, Blake E.; Umaña, Ariana; Lemkul, Justin A.; Slade, Daniel J. (American Society for Microbiology, 2018-07-05)
    Here we present FusoPortal, an interactive repository of Fusobacterium genomes that were sequenced using a hybrid MinION long-read sequencing pipeline, followed by assembly and annotation using a diverse portfolio of predominantly open-source software. Significant efforts were made to provide genomic and bioinformatic data as downloadable files, including raw sequencing reads, genome maps, gene annotations, protein functional analysis and classifications, and a custom BLAST server for FusoPortal genomes. FusoPortal has been initiated with eight complete genomes, of which seven were previously only drafts that ranged from 24 to 67 contigs. We have showcased that the genomes in FusoPortal provide accurate open reading frame annotations and have corrected a number of large (3-kb) genes that were previously misannotated due to contig boundaries. In summary, FusoPortal (http://fusoportal.org) is the first database of MinION-sequenced and completely assembled Fusobacterium genomes, and this central Fusobacterium genomic and bioinformatic resource will aid the scientific community in developing a deeper understanding of how this human pathogen contributes to an array of diseases, including periodontitis and colorectal cancer.
  • Thumbnail Image
    Impact of Electronic Polarization on Preformed, β-Strand Rich Homogenous and Heterogeneous Amyloid Oligomers
    King, Kelsie M.; Sharp, Amanda K.; Davidson, Darcy S.; Brown, Anne M.; Lemkul, Justin A. (World Scientific, 2021-12-29)
    Amyloids are a subset of intrinsically disordered proteins (IDPs) that self-assemble into cross-𝛽 oligomers and fibrils. The structural plasticity of amyloids leads to sampling of metastable, low-molecular-weight oligomers that contribute to cytotoxicity. Of interest are amyloid-𝛽 (A𝛽 and islet amyloid polypeptide (IAPP), which are involved in the pathology of Alzheimer’s disease and Type 2 diabetes mellitus, respectively. In addition to forming homogenous oligomers and fibrils, these species have been found to cross-aggregate in heterogeneous structures. Biophysical properties, including electronic effects, that are unique or conserved between homogenous and heterogeneous amyloids oligomers are thus far unexplored. Here, we simulated homogenous and heterogeneous amyloid oligomers of A𝛽16−22 and IAPP20−29 fragments using the Drude oscillator model to investigate the impact of electronic polarization on the structural morphology and stability of preformed hexamers. Upon simulation of preformed, 𝛽-strand rich oligomers with Drude, structural rearrangement occurred causing some loss of 𝛽-strand structure in favor of random coil content for all oligomers. Homogenous A𝛽16−22 was the most stable system, deriving stability from low polarization in hydrophobic residues and through salt bridge formation. Changes in polarization were observed primarily for A𝛽16−22 residues in heterogeneous cross-amyloid systems, displaying a decrease in charged residue dipole moments and an increase in hydrophobic sidechain dipole moments. This work is the first study utilizing the Drude-2019 force field with amyloid oligomers, providing insight into the impact of electronic effects on oligomer structure and highlighting the importance of different microenvironments on amyloid oligomer stability.
  • An in silico approach to determine inter-subunit affinities in human septin complexes
    Grupp, Benjamin; Lemkul, Justin A.; Gronemeyer, Thomas (Wiley, 2023-03-09)
    The septins are a conserved family of filament-forming guanine nucleotide binding proteins, often named the fourth component of the cytoskeleton. Correctly assembled septin structures are required for essential intracellular processes such as cytokinesis, vesicular transport, polarity establishment, and cellular adhesion. Structurally, septins belong to the P-Loop NTPases but they do not mediate signals to effectors through GTP binding and hydrolysis. GTP binding and hydrolysis are believed to contribute to septin complex integrity, but biochemical approaches addressing this topic are hampered by the stability of septin complexes after recombinant expression and the lack of nucleotide-depleted complexes. To overcome this limitation, we used a molecular dynamics-based approach to determine inter-subunit binding free energies in available human septin dimer structures and in their apo forms, which we generated in silico. The nucleotide in the GTPase active subunits SEPT2 and SEPT7, but not in SEPT6, was identified as a stabilizing element in the G interface. Removal of GDP from SEPT2 and SEPT7 results in flipping of a conserved Arg residue and disruption of an extensive hydrogen bond network in the septin unique element, concomitant with a decreased inter-subunit affinity. Based on these findings we propose a singular “lock-hydrolysis” mechanism stabilizing human septin filaments.
  • Thumbnail Image
    Integration of experimental data and use of automated fitting methods in developing protein force fields
    Polêto, Marcelo D.; Lemkul, Justin A. (Springer Nature, 2022-03-18)
    The development of accurate protein force fields has been the cornerstone of molecular simulations for the past 50 years. During this period, many lessons have been learned regarding the use of experimental target data and parameter fitting procedures. Here, we review recent advances in protein force field development. We discuss the recent emergence of polarizable force fields and the role of electronic polarization and areas in which additive force fields fall short. The use of automated fitting methods and the inclusion of additional experimental solution data during parametrization is discussed as a means to highlight possible routes to improve the accuracy of force fields even further.
  • Thumbnail Image
    Investigating the Electrostatic Properties and Dynamics of Amyloidogenic Proteins with Polarizable Molecular Dynamics Simulations
    Davidson, Darcy Shanley (Virginia Tech, 2022-04-14)
    Amyloidogenic diseases, such as Alzheimer's disease (AD) and Type II Diabetes (T2D), are characterized by the accumulation of amyloid aggregates. Despite having very different amino-acid sequences, the underlying amyloidogenic proteins form similar supramolecular fibril structures that are highly stable and resistant to physical and chemical denaturation. AD is characterized by two toxic lesions: extracellular amyloid β-peptide (Aβ) plaques and intracellular neurofibrillary tangles composed of microtubule-associated protein tau. Similarly, a feature of T2D is the deposition of islet amyloid polypeptide (IAPP) aggregates in and around the pancreas. The mechanisms by which Aβ, tau, and IAPP aggregate, and cause cell death is unknown; thus, gaining greater insight into the stabilizing forces and initial unfolding events is crucial to our understanding of these amyloidogenic diseases. This work uses molecular dynamics (MD) simulations to study the secondary, tertiary, and quaternary structure of Aβ, tau, and IAPP. Specifically, this work used the Drude polarizable force field (FF), which explicitly represents electronic polarization allowing charge distributions to change in response to perturbations in local electric fields. This model allows us to describe the role charge plays on protein folding and stability and how perturbations to the charge state drive pathology. Studies were conducted to address the following questions: 1) What are the stabilizing forces of fibril and oligomeric structures? 2) How do charge-altering mutations modulate the conformational ensemble and thermodynamic properties of Aβ? 3) How do charge-altering post-translational modifications of Aβ and tau modulate changes in the conformational ensembles? These studies establish that shifts in local microenvironments play a role in fibril and oligomer stability. Furthermore, these studies found that changes in protein sequence and charge are sufficient to disrupt and change the secondary and tertiary structure of these amyloidogenic proteins. Overall, this dissertation describes how charge modulates protein unfolding and characterizes the mechanism of those changes. In the long term, this work will help in the development of therapeutics that can target these changes to prevent protein aggregation that leads to cell death.
  • Thumbnail Image
    Ion-Dependent Conformational Plasticity of Telomeric G-Hairpins and G-Quadruplexes
    Salsbury, Alexa M.; Michel, Haley M.; Lemkul, Justin A. (American Chemical Society, 2022-07-12)
    Telomeric DNA is guanine-rich and can adopt structures such as G-quadruplexes (GQs) and G-hairpins. Telomeric GQs influence genome stability and telomerase activity, making understanding of enzyme-GQ interactions and dynamics important for potential drug design. GQs have a characteristic tetrad core, which is connected by loop regions. Within this architecture are G-hairpins, fold back motifs that are thought to represent the first intermediate in GQ folding. To better understand the relationship between G-hairpin motifs and GQs, we performed polarizable simulations of a two-tetrad telomeric GQ and an isolated SC11 telomeric G-hairpin. The telomeric GQ contains a G-triad, which functions as part of the tetrad core or linker regions, depending on local conformational change. This triad and another motif below the tetrad core frequently bound ions and may represent druggable sites. Further, we observed the unbiased formation of a G-triad and a G-tetrad in simulations of the SC11 G-hairpin and found that cations can be partially hydrated while facilitating the formation of these motifs. Finally, we demonstrated that K+ ions form specific interactions with guanine bases, while Na+ ions interact nonspecifically with bases in the structure. Together, these simulations provide new insights into the influence of ions on GQs, G-hairpins, and G-triad motifs.
  • Thumbnail Image
    Molecular Modeling of the Amyloid β-Peptide: Understanding the Mechanism of Alzheimer's Disease and the Potential for Therapeutic Intervention
    Lemkul, Justin A. (Virginia Tech, 2012-03-07)
    Alzheimer's disease is the leading cause of senile dementia in the elderly, and as life expectancy increases across the globe, incidence of the disease is continually increasing. Current estimates place the number of cases at 25-30 million worldwide, with more than 5.4 million of these occurring in the United States. While the exact cause of the disease remains a mystery, it has become clear that the amyloid β-peptide (Aβ) is central to disease pathogenesis. The aggregation and deposition of this peptide in the brain is known to give rise to the hallmark lesions associated with Alzheimer's disease, but its exact mechanism of toxicity remains largely uncharacterized. Molecular dynamics (MD) simulations have achieved great success in exploring molecular events with atomic resolution, predicting and explaining phenomena that are otherwise obscured from even the most sensitive experimental techniques. Due to the difficulty of obtaining high-quality structural data of Aβ and its toxic assemblies, MD simulations can be an especially useful tool in understanding the progression of Alzheimer's disease on a molecular level. The work contained herein describes the interactions of Aβ monomers and oligomers with lipid bilayers to understand the mechanism by which Aβ exerts its toxicity. Also explored is the mechanism by which flavonoid antioxidants may prevent Aβ self-association and destabilize toxic aggregates, providing insight into the chemical features that give rise to this therapeutic effect.
  • Thumbnail Image
    Phosphorylation of PPAR gamma Affects the Collective Motions of the PPAR gamma-RXR alpha-DNA Complex
    Lemkul, Justin A.; Lewis, Stephanie N.; Bassaganya-Riera, Josep; Bevan, David R. (PLOS, 2015-05-08)
    Peroxisome-proliferator activated receptor-γ (PPARγ) is a nuclear hormone receptor that forms a heterodimeric complex with retinoid X receptor-α (RXRα) to regulate transcription of genes involved in fatty acid storage and glucose metabolism. PPARγ is a target for pharmaceutical intervention in type 2 diabetes, and insight into interactions between PPARγ, RXRα, and DNA is of interest in understanding the function and regulation of this complex. Phosphorylation of PPARγ by cyclin-dependent kinase 5 (Cdk5) has been shown to dysregulate the expression of metabolic regulation genes, an effect that is counteracted by PPARγ ligands. We applied molecular dynamics (MD) simulations to study the relationship between the ligand-binding domains of PPARγ and RXRα with their respective DNA-binding domains. Our results reveal that phosphorylation alters collective motions within the PPARγ-RXRα complex that affect the LBD-LBD dimerization interface and the AF-2 coactivator binding region of PPARγ.
  • Thumbnail Image
    Polarizable Simulations of the bcl-2 DNA G-Quadruplex and FMRP RNA G-Quadruplex:Duplex Junction Binding Protein
    Ratnasinghe, Brian Damith (Virginia Tech, 2021-06-03)
    A G-quadruplex (GQ) is a type of noncanonical nucleic acid structure that can form in regions of nucleic acids rich in guanine nucleotides. The guanine bases form a square planar conformation via Hoogsteen hydrogen bonding. These stacked tetrads have inward-facing carbonyl oxygens, facilitating the coordination of ions. Improper GQ conformations can lead to improper regulation of gene expression, potentially resulting in genetic diseases or cancer. Here, we performed molecular dynamics simulations using the Drude polarizable force field (FF) to gain insight into factors contributing to the stability of two GQs. One is the bcl-2 promoter region GQ, which is implicated in several types of cancer including B-cell lymphoma, and the second is the sc1 RNA GQ, which binds to the Fragile-X Mental Retardation Protein (FMRP) and is implicated in the development of Fragile X Syndrome (FXS). Aberrant bcl-2 GQ conformations result in increased production of the BCL2 protein, which is an apoptosis inhibitor. As such, we aim to characterize the factors stabilizing the GQ for future small-molecule development to prevent apoptosis inhibition and therefore cancer. The FMRP protein functions as a regulator of sc1 conformation to control the translation of proteins required for frontal lobe development. FXS arises from a nonsense mutation that causes the deletion of the C-terminal region of FMRP, rendering it non-function. Therefore, we aim to simulate sc1 when FMRP is bound as well as unbound to provide insight into the types of interactions that must be maintained and therefore mimicked by a small molecule drug.