Computational investigations of polymerase enzymes: Structure, function, inhibition, and biotechnology

AbstractDNA and RNA polymerases (Pols) are central to life, health, and biotechnology because they allow the flow of genetic information in biological systems. Importantly, Pol function and (de)regulation are linked to human diseases, notably cancer (DNA Pols) and viral infections (RNA Pols) such as COVID‐19. In addition, Pols are used in various applications such as synthesis of artificial genetic polymers and DNA amplification in molecular biology, medicine, and forensic analysis. Because of all of this, the field of Pols is an intense research area, in which computational studies contribute to elucidating experimentally inaccessible atomistic details of Pol function. In detail, Pols catalyze the replication, transcription, and repair of nucleic acids through the addition, via a nucleotidyl transfer reaction, of a nucleotide to the 3′‐end of the growing nucleic acid strand. Here, we analyze how computational methods, including force‐field‐based molecular dynamics, quantum mechanics/molecular mechanics, and free energy simulations, have advanced our understanding of Pols. We examine the complex interaction of chemical and physical events during Pol catalysis, like metal‐aided enzymatic reactions for nucleotide addition and large conformational rearrangements for substrate selection and binding. We also discuss the role of computational approaches in understanding the origin of Pol fidelity—the ability of Pols to incorporate the correct nucleotide that forms a Watson–Crick base pair with the base of the template nucleic acid strand. Finally, we explore how computations can accelerate the discovery of Pol‐targeting drugs and engineering of artificial Pols for synthetic and biotechnological applications.This article is categorized under: Structure and Mechanism > Reaction Mechanisms and Catalysis Structure and Mechanism > Computational Biochemistry and Biophysics Software > Molecular Modelin

Molecular Dynamics Study of the Opening Mechanism for DNA Polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although Xray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme:DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics

Structural Mechanism of Substrate Specificity In Human Cytidine Deaminase Family APOBEC3s

APOBEC3s (A3s) are a family of human cytidine deaminases that play important roles in both innate immunity and cancer. A3s protect host cells against retroviruses and retrotransposons by deaminating cytosine to uracil on foreign pathogenic genomes. However, when mis-regulated, A3s can cause heterogeneities in host genome and thus promote cancer and the development of therapeutic resistance. The family consists of seven members with either one (A3A, A3C and A3H) or two zinc-binding domains (A3B, A3D, A3D and A3G). Despite overall similarity, A3 proteins have distinct deamination activity and substrate specificity. Over the past years, several crystal and NMR structures of apo A3s and DNA/RNA-bound A3s have been determined. These structures have suggested the importance of the loops around the active site for nucleotide specificity and binding. However, the structural mechanism underlying A3 activity and substrate specificity requires further examination. Using a combination of computational molecular modeling and parallel molecular dynamics (pMD) simulations followed by experimental verifications, I investigated the roles of active site residues and surrounding loops in determining the substrate specificity and RNA versus DNA binding among A3s. Starting with A3B, I revealed the structural basis and gatekeeper residue for DNA binding. I also identified a unique auto-inhibited conformation in A3B that restricts access to the active site and may underlie lower catalytic activity compared to the highly similar A3A. Besides, I investigated the structural mechanism of substrate specificity and ssDNA binding conformation in A3s. I found an interdependence between substrate conformation and specificity. Specifically, the linear DNA conformation helps accommodate CC dinucleotide motif while the U-shaped conformation prefers TC. I also identified the molecular mechanisms of substrate sequence specificity at -1’ and -2’ positions. Characterization of substrate binding to A3A revealed that intra-DNA interactions may be responsible for the specificity in A3A. Finally, I investigated the structural mechanism for exclusion of RNA from A3G catalytic activity using similar methods. Overall, the comprehensive analysis of A3s in this thesis shed light into the structural mechanism of substrate specificity and broaden the understanding of molecular interactions underlying the biological function of these enzymes. These results have implications for designing specific A3 inhibitors as well as base editing systems for gene therapy

Modelling polyketide synthases and related macromolecular complexes

Polyketide synthases (PKS) are enzyme complexes that synthesise many natural products of medicinal interest, notably a large number of antibiotics. The present work investigated the mupirocin biosynthesis system, comparing it with similar pathways such as thiomarinol and kalimantacin. The focus was on the structural modelling of the protein complexes involved in antibiotic synthesis, via molecular simulation and the analysis of structural and sequence data. Structural docking of acyl carrier proteins (ACP) cognate for an HMG-CoA synthase orthologue responsible for β-methylation (MupH) identified key residues involved in the recognitions specificity of the interacting partners, further supported by mutagenesis experiments, which thus allows prediction of β-methylation sites in PKS. Moreover, complementation and mutagenesis experiments performed on MupH homologs from kalimantacin and thiomarinol systems suggests specificity between the ACP:HCS proteins in the β-branching suggesting the possibility of engineering multiple specific β-branching modifications into the same pathway. Molecular dynamics simulations of ACPs from the mupirocin cluster revealed that the PKS ACPs form a cavity upon the attachment of the phosphopantetheine and acyl chains similar to what is seen in the fatty acid synthase ACPs and provide a better understanding of the structure function relationship in these small proteins. Molecular docking of the putative cognate substrate with the ketosynthase (KS) homo dimer of module 5 of the MmpA in the mupirocin pathway revealed a loop that may control specificity for the α-hydroxylated substrate and mutagenesis experiments support this proposition

Computational Modeling of Protein Kinases: Molecular Basis for Inhibition and Catalysis

Protein kinases catalyze protein phosphorylation reactions, i.e. the transfer of the γ-phosphoryl group of ATP to tyrosine, serine and threonine residues of protein substrates. This phosphorylation plays an important role in regulating various cellular processes. Deregulation of many kinases is directly linked to cancer development and the protein kinase family is one of the most important targets in current cancer therapy regimens. This relevance to disease has stimulated intensive efforts in the biomedical research community to understand their catalytic mechanisms, discern their cellular functions, and discover inhibitors. With the advantage of being able to simultaneously define structural as well as dynamic properties for complex systems, computational studies at the atomic level has been recognized as a powerful complement to experimental studies. In this work, we employed a suite of computational and molecular simulation methods to (1) explore the catalytic mechanism of a particular protein kinase, namely, epidermal growth factor receptor (EGFR); (2) study the interaction between EGFR and one of its inhibitors, namely erlotinib (Tarceva); (3) discern the effects of molecular alterations (somatic mutations) of EGFR to differential downstream signaling response; and (4) model the interactions of a novel class of kinase inhibitors with a common ruthenium based organometallic scaffold with different protein kinases. Our simulations established some important molecular rules in operation in the contexts of inhibitor-binding, substrate-recognition, catalytic landscapes, and signaling in the EGFR tyrosine kinase. Our results also shed insights on the mechanisms of inhibition and phosphorylation commonly employed by many kinases

“Gate-keeper” Residues and Active-Site Rearrangements in DNA Polymerase μ Help Discriminate Non-cognate Nucleotides

<div><p>Incorporating the cognate instead of non-cognate substrates is crucial for DNA polymerase function. Here we analyze molecular dynamics simulations of DNA polymerase μ (pol μ) bound to different non-cognate incoming nucleotides including A:dCTP, A:dGTP, A(syn):dGTP, A:dATP, A(syn):dATP, T:dCTP, and T:dGTP to study the structure-function relationships involved with aberrant base pairs in the conformational pathway; while a pol μ complex with the A:dTTP base pair is available, no solved non-cognate structures are available. We observe distinct differences of the non-cognate systems compared to the cognate system. Specifically, the motions of active-site residue His329 and Asp330 distort the active site, and Trp436, Gln440, Glu443 and Arg444 tend to tighten the nucleotide-binding pocket when non-cognate nucleotides are bound; the latter effect may further lead to an altered electrostatic potential within the active site. That most of these “gate-keeper” residues are located farther apart from the upstream primer in pol μ, compared to other X family members, also suggests an interesting relation to pol μ's ability to incorporate nucleotides when the upstream primer is not paired. By examining the correlated motions within pol μ complexes, we also observe different patterns of correlations between non-cognate systems and the cognate system, especially decreased interactions between the incoming nucleotides and the nucleotide-binding pocket. Altered correlated motions in non-cognate systems agree with our recently proposed hybrid conformational selection/induced-fit models. Taken together, our studies propose the following order for difficulty of non-cognate system insertions by pol μ: T:dGTP</p></div

Molecular Targets of CNS Tumors

Molecular Targets of CNS Tumors is a selected review of Central Nervous System (CNS) tumors with particular emphasis on signaling pathway of the most common CNS tumor types. To develop drugs which specifically attack the cancer cells requires an understanding of the distinct characteristics of those cells. Additional detailed information is provided on selected signal pathways in CNS tumors

Psr1p interacts with SUN/sad1p and EB1/mal3p to establish the bipolar spindle

Regular Abstracts - Sunday Poster Presentations: no. 382During mitosis, interpolar microtubules from two spindle pole bodies (SPBs) interdigitate to create an antiparallel microtubule array for accommodating numerous regulatory proteins. Among these proteins, the kinesin-5 cut7p/Eg5 is the key player responsible for sliding apart antiparallel microtubules and thus helps in establishing the bipolar spindle. At the onset of mitosis, two SPBs are adjacent to one another with most microtubules running nearly parallel toward the nuclear envelope, creating an unfavorable microtubule configuration for the kinesin-5 kinesins. Therefore, how the cell organizes the antiparallel microtubule array in the first place at mitotic onset remains enigmatic. Here, we show that a novel protein psrp1p localizes to the SPB and plays a key role in organizing the antiparallel microtubule array. The absence of psr1+ leads to a transient monopolar spindle and massive chromosome loss. Further functional characterization demonstrates that psr1p is recruited to the SPB through interaction with the conserved SUN protein sad1p and that psr1p physically interacts with the conserved microtubule plus tip protein mal3p/EB1. These results suggest a model that psr1p serves as a linking protein between sad1p/SUN and mal3p/EB1 to allow microtubule plus ends to be coupled to the SPBs for organization of an antiparallel microtubule array. Thus, we conclude that psr1p is involved in organizing the antiparallel microtubule array in the first place at mitosis onset by interaction with SUN/sad1p and EB1/mal3p, thereby establishing the bipolar spindle.postprin