A theoretical study of the interaction of guanine and cytosine with specific amino acid side chains

Molecular contacts in the Cren7-DNA complex: A quantitative investigation for electrostatic interaction

The molecular association of proteins with nucleic acids leading to the formation of macromolecular complexes is a crucial step in several biological processes. Stabilization of these complexes involves electrostatic interactions between ion pairs (salt bridges) of nucleic acid phosphates and protein side chains. The crenarchaeal DNA binding protein, Cren7 plays a key role in the regulation of chromosomal structure and gene expression in eukaryotic extremophiles. However, the molecular contacts that occur at the interface of protein-DNA complexes and their contribution to the electrostatic interaction have not been fully elucidated. This work presents a quantitative description of the mechanism of the electrostatic interaction between the protein and DNA. We have identified a few residues located at the Cren7-DNA interface that could potentially be responsible for the interaction. Structural studies using circular dichroism indicate mutation of these surface residues minimally affect their structure and stability. The binding affinity of these mutants for the DNA duplexes was examined from reverse titration, biolayer interferometry, and fluorescence anisotropy measurements with calf thymus DNA, polynucleotides, and small DNA oligonucleotides. The resulting kinetic parameters highlight a difference in electrostatic interactions potentials exhibited by residues positioned at different locations of the protein-DNA interface. Computational studies attribute this difference to their surrounding atmosphere and energetic stabilization parameters. The biophysical approach described here can be extended for other proteins that play a crucial role in DNA bending and compaction, to properly evaluate the role of specific residues on the mechanisms of DNA binding.

Genome-wide interrogation of gene functions through base editor screens empowered by barcoded sgRNAs

Canonical CRISPR-knockout (KO) screens rely on Cas9-induced DNA double-strand breaks (DSBs) to generate targeted gene KOs. These methodologies may yield distorted results because DSB-associated effects are often falsely assumed to be consequences of gene perturbation itself, especially when high copy-number sites are targeted. In the present study, we report a DSB-independent, genome-wide CRISPR screening method, termed iBARed cytosine base editing-mediated gene KO (BARBEKO). This method leverages CRISPR cytosine base editors for genome-scale KO screens by perturbing gene start codons or splice sites, or by introducing premature termination codons. Furthermore, it is integrated with iBAR, a strategy we devised for improving screening quality and efficiency. By constructing such a cell library through lentiviral infection at a high multiplicity of infection (up to 10), we achieved efficient and accurate screening results with substantially reduced starting cells. More importantly, in comparison with Cas9-mediated fitness screens, BARBEKO screens are no longer affected by DNA cleavage-induced cytotoxicity in HeLa-, K562- or DSB-sensitive retinal pigmented epithelial 1 cells. We anticipate that BARBEKO offers a valuable tool to complement the current CRISPR-KO screens in various settings.

Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy

We present an overview of the SIBFA polarizable molecular mechanics procedure, which is formulated and calibrated on the basis of quantum chemistry (QC). It embodies nonclassical effects such as electrostatic penetration, exchange-polarization, and charge transfer. We address the issues of anisotropy, nonadditivity, and transferability by performing parallel QC computations on multimolecular complexes. These encompass multiply H-bonded complexes and polycoordinated complexes of divalent cations. Recent applications to the docking of inhibitors to Zn-metalloproteins are presented next, namely metallo-beta-lactamase, phosphomannoisomerase, and the nucleocapsid of the HIV-1 retrovirus. Finally, toward third-generation intermolecular potentials based on density fitting, we present the development of a novel methodology, the Gaussian electrostatic model (GEM), which relies on ab initio-derived fragment electron densities to compute the components of the total interaction energy. As GEM offers the possibility of a continuous electrostatic model going from distributed multipoles to densities, it allows an inclusion of short-range quantum effects in the molecular mechanics energies. The perspectives of an integrated SIBFA/GEM/QM procedure are discussed.

Experimental and theoretical study of energetics of complex formation between nucleic acid bases and the bases with amide group

A number of nucleic acid base pairs and complexes between the bases and the amide group of acrylamide have been studied experimentally by using mass spectrometry and theoretically by the method of atom-atom potential function calculations. It has been found from temperature dependencies of peak intensities in mass spectra of m2.2.9(3) Gua.m1Ura, m9 Ade.m1Cyt, m2.2.9(3) Gua.m1Gua.m1Cyt pairs that enthalpy values, delta H, of the complex formation are equal to 14.2 +/- 1.1, 13.5 +/- 1.3 and 16.4 +/- 1.4 kcal/M, respectively, and those of acrylamide with m1.3(2) Ura and m1Thy corresponds to 9.7 +/- 1.0 and 6.8 +/- 0.6 kcal/M. There is a good agreement of the experimental data with calculations when taking into account both the amino-oxo and the amino-hydroxy tautomeric forms of guanine. A combined use of the data allows us to determine the energy, the modes of interaction and the structure of the complexes. The results are discussed in connection with the modelling of molecular structure of biopolymers by the method of classical potential functions, protein-nucleic acids recognition and fidelity of nucleic acids biosynthesis.

Theoretical design of a bistetrapeptide derivative of mitoxantrone targeted towards the double-stranded hexanucleotide sequence d(GGCGCC)2

The hexanucleotide d(GGCGCC)2 is encountered in recurrent fashion within transcriptional activating sequences in retroviruses and protooncogenes. Our first theoretical design of novel oligopeptide derivatives of mitoxantrone, MTX (1), had enabled us to predict derivatives depsiGly-Lys(L) and depsiGly-Gly-Orn(D) to preferentially target the tetrameric core d(CCGG)2. Owing to the crucial importance of hexamer d(GGCGCC)2, we have attempted to extend the realm of our approach by now targeting this specific hexanucleotide. For that purpose, we undertook the design of further oligopeptide derivatives of MTX, in which each arm was identically amidated (rather than esterified as in (1)) by tri- or tetrapeptides of varying sequences and individual residue configurations. The binding affinities of these derivatives to the palindromic sequences d(GGCGCC)2, d(CGCGCG)2, d(GCCGGC)2 and d(CCCGGG)2, were compared by energy-minimization. We report here the results obtained with the most promising derivative, having the sequence Arg(L)-Gly-Val(L)-Glu(L), and displaying a considerable energy preference for d(GGCGCC)2 over the other candidate hexameric sites (referred to as I). In the corresponding complexes, the two arms are in two mutually antiparallel directions in the major groove, and adopt a beta-sheet like arrangement stabilized by two H-bonds involving the carbonyl and amide groups of the Gly residues. Each Arg side chain on a given arm chelates O6 and N7 atoms of G1, G2/G1', G2' with its imino and cis amino hydrogen, and is simultaneously bound through two amino hydrogens in a bidentate interaction with the Glu residue.

Associate hydrations and energies of nucleotide bases as revealed by low-temperature field ionization mass-spectrometric data

The formation of water clusters, polyhydrates of nucleotide bases and their associates during simultaneous condensation of water and base molecules in vacuo onto a surface of a needle emitter cooled to 170 K was studied by field ionization mass spectrometry. It was found that different emitter temperatures are characterized by a specific distribution of intensities of cluster currents, depending on the number of water molecules in clusters. These distributions correlate with structural peculiarities and the relative energetics of formation of water clusters, polyhydrates of nucleotide bases and their associates at low temperature. The features observed in mass spectra for clusters m9Ade (H2O)5, m1Ura (H2O)4 and m9Ade m1Ura (H2O)2 are treated as a result of formation of energetically favorable structures stabilized by H-bonded bridges of water molecules. The relative association constants and formation enthalpies of the noncomplementary pairs Ade Cyt, Gua Ura and the associates which model the aminoacid-base complexes m1Ura Gln and m1.3(2)Thy Gln were determined from the temperature dependencies of the intensities of mass spectra peaks in the range 290-320 K.

Determinants of EcoRI endonuclease sequence discrimination

The arginine at position 200 of EcoRI endonuclease is thought to make two hydrogen bonds to the guanine of the sequence GAATTC and thus be an important determinant of sequence discrimination. Arg-200 was replaced by each of the other 19 naturally occurring amino acids, and the mutant endonucleases were assessed for activities in vivo and in vitro. The mutant endonuclease with lysine at position 200 exhibits the most in vivo activity of all the position 200 mutants, although the in vitro activity is less than 1/100th of wild-type activity. Five other mutants show more drastically reduced levels of in vivo activity (Cys, Pro, Val, Ser, and Trp). The Cys, Val, and Ser mutant enzymes appear to have in vivo activity which is specific for the wild-type canonical site despite the loss of hydrogen bonding potential at position 200. The Pro and Trp mutants retain in vivo activity which is independent of the presence of the EcoRI methylase. In crude cell lysates, only the Cys mutant shows a very low level of in vitro activity. None of the mutant enzymes show a preference for alternative sites in assays in vitro. The implications of these results are discussed.

CTP can replace GTP in reactions catalyzed by eukaryotic peptide elongation factor 1

In several reactions catalyzed by highly purified peptide elongation factor 1 from rabbit reticulocytes, GTP may be fully replaced by CTP but not by ATP or UTP. This holds true for the factor-dependent binding of aminoacyl-tRNA to ribosomes, GTPase activity, GTP-dependent autophosphorylation of the factor protein and binding of cholesteryl 14-methylhexadecanoate by the factor.