Recognition of 5'-YpG-3' sequences by coupled stacking/hydrogen bonding interactions with amino acid residues

Molecular mechanism of specific DNA sequence recognition by NRF1

Nuclear respiratory factor 1 (NRF1) regulates the expression of genes that are vital for mitochondrial biogenesis, respiration, and various other cellular processes. While NRF1 has been reported to bind specifically to GC-rich promoters as a homodimer, the precise molecular mechanism governing its recognition of target gene promoters has remained elusive. To unravel the recognition mechanism, we have determined the crystal structure of the NRF1 homodimer bound to an ATGCGCATGCGCAT dsDNA. In this complex, NRF1 utilizes a flexible linker to connect its dimerization domain (DD) and DNA binding domain (DBD). This configuration allows one NRF1 monomer to adopt a U-turn conformation, facilitating the homodimer to specifically bind to the two TGCGC motifs in the GCGCATGCGC consensus sequence from opposite directions. Strikingly, while the NRF1 DBD alone could also bind to the half-site (TGCGC) DNA of the consensus sequence, the cooperativity between DD and DBD is essential for the binding of the intact GCGCATGCGC sequence and the transcriptional activity of NRF1. Taken together, our results elucidate the molecular mechanism by which NRF1 recognizes specific DNA sequences in the promoters to regulate gene expression.

Ned Seeman and the prediction of amino acid-basepair motifs mediating protein-nucleic acid recognition

Fifty years ago, the first atomic-resolution structure of a nucleic acid double helix, the mini-duplex (ApU)2, revealed details of basepair geometry, stacking, sugar conformation, and backbone torsion angles, thereby superseding earlier models based on x-ray fiber diffraction, including the original DNA double helix proposed by Watson and Crick. Just 3 years later, in 1976, Ned Seeman, John Rosenberg, and Alex Rich leapt from their structures of mini-duplexes and H-bonding motifs between bases in small-molecule structures and transfer RNA to predicting how proteins could sequence specifically recognize double helix nucleic acids. They proposed interactions between amino acid side chains and nucleobases mediated by two hydrogen bonds in the major or minor grooves. One of these, the arginine-guanine pair, emerged as the most favored amino acid-base interaction in experimental structures of protein-nucleic acid complexes determined since 1986. In this brief review we revisit the pioneering work by Seeman et al. and discuss the importance of the arginine-guanine pairing motif.

Structural basis of DNA methylation-dependent site selectivity of the Epstein-Barr virus lytic switch protein ZEBRA/Zta/BZLF1

In infected cells, Epstein-Barr virus (EBV) alternates between latency and lytic replication. The viral bZIP transcription factor ZEBRA (Zta, BZLF1) regulates this cycle by binding to two classes of ZEBRA response elements (ZREs): CpG-free motifs resembling the consensus AP-1 site recognized by cellular bZIP proteins and CpG-containing motifs that are selectively bound by ZEBRA upon cytosine methylation. We report structural and mutational analysis of ZEBRA bound to a CpG-methylated ZRE (meZRE) from a viral lytic promoter. ZEBRA recognizes the CpG methylation marks through a ZEBRA-specific serine and a methylcytosine-arginine-guanine triad resembling that found in canonical methyl-CpG binding proteins. ZEBRA preferentially binds the meZRE over the AP-1 site but mutating the ZEBRA-specific serine to alanine inverts this selectivity and abrogates viral replication. Our findings elucidate a DNA methylation-dependent switch in ZEBRA's transactivation function that enables ZEBRA to bind AP-1 sites and promote viral latency early during infection and subsequently, under appropriate conditions, to trigger EBV lytic replication by binding meZREs.

Structures of Neisseria gonorrhoeae MtrR-operator complexes reveal molecular mechanisms of DNA recognition and antibiotic resistance-conferring clinical mutations

Mutations within the mtrR gene are commonly found amongst multidrug resistant clinical isolates of Neisseria gonorrhoeae, which has been labelled a superbug by the Centers for Disease Control and Prevention. These mutations appear to contribute to antibiotic resistance by interfering with the ability of MtrR to bind to and repress expression of its target genes, which include the mtrCDE multidrug efflux transporter genes and the rpoH oxidative stress response sigma factor gene. However, the DNA-recognition mechanism of MtrR and the consensus sequence within these operators to which MtrR binds has remained unknown. In this work, we report the crystal structures of MtrR bound to the mtrCDE and rpoH operators, which reveal a conserved, but degenerate, DNA consensus binding site 5'-MCRTRCRN4YGYAYGK-3'. We complement our structural data with a comprehensive mutational analysis of key MtrR-DNA contacts to reveal their importance for MtrR-DNA binding both in vitro and in vivo. Furthermore, we model and generate common clinical mutations of MtrR to provide plausible biochemical explanations for the contribution of these mutations to multidrug resistance in N. gonorrhoeae. Collectively, our findings unveil key biological mechanisms underlying the global stress responses of N. gonorrhoeae.

A TetR-family transcription factor regulates fatty acid metabolism in the archaeal model organism Sulfolobus acidocaldarius

Fatty acid metabolism and its regulation are known to play important roles in bacteria and eukaryotes. By contrast, although certain archaea appear to metabolize fatty acids, the regulation of the underlying pathways in these organisms remains unclear. Here, we show that a TetR-family transcriptional regulator (FadR_Sa) is involved in regulation of fatty acid metabolism in the crenarchaeon Sulfolobus acidocaldarius. Functional and structural analyses show that FadR_Sa binds to DNA at semi-palindromic recognition sites in two distinct stoichiometric binding modes depending on the operator sequence. Genome-wide transcriptomic and chromatin immunoprecipitation analyses demonstrate that the protein binds to only four genomic sites, acting as a repressor of a 30-kb gene cluster comprising 23 open reading frames encoding lipases and β-oxidation enzymes. Fatty acyl-CoA molecules cause dissociation of FadR_Sa binding by inducing conformational changes in the protein. Our results indicate that, despite its similarity in overall structure to bacterial TetR-family FadR regulators, FadR_Sa displays a different acyl-CoA binding mode and a distinct regulatory mechanism.

Certain archaea appear to metabolize fatty acids, but the regulation of these pathways is unclear. Here, Wang et al. provide genetic, functional and structural evidence supporting that a TetR-family transcriptional regulator is involved in regulation of fatty acid metabolism in Sulfolobus acidocaldarius.

Structural basis of transcriptional regulation by the HigA antitoxin

Bacterial toxin-antitoxin systems are important factors implicated in growth inhibition and plasmid maintenance. Type II toxin-antitoxin pairs are regulated at the transcriptional level by the antitoxin itself. Here, we examined how the HigA antitoxin regulates the expression of the Proteus vulgaris higBA toxin-antitoxin operon from the Rts1 plasmid. The HigBA complex adopts a unique architecture suggesting differences in its regulation as compared to classical type II toxin-antitoxin systems. We find that the C-terminus of the HigA antitoxin is required for dimerization and transcriptional repression. Further, the HigA structure reveals that the C terminus is ordered and does not transition between disorder-to-order states upon toxin binding. HigA residue Arg40 recognizes a TpG dinucleotide in higO2, an evolutionary conserved mode of recognition among prokaryotic and eukaryotic transcription factors. Comparison of the HigBA and HigA-higO2 structures reveals the distance between helix-turn-helix motifs of each HigA monomer increases by ~4 Å in order to bind to higO2. Consistent with these data, HigBA binding to each operator is twofold less tight than HigA alone. Together, these data show the HigB toxin does not act as a co-repressor suggesting potential novel regulation in this toxin-antitoxin system.

The MerR-like protein BldC binds DNA direct repeats as cooperative multimers to regulate Streptomyces development

Streptomycetes are notable for their complex life cycle and production of most clinically important antibiotics. A key factor that controls entry into development and the onset of antibiotic production is the 68-residue protein, BldC. BldC is a putative DNA-binding protein related to MerR regulators, but lacks coiled-coil dimerization and effector-binding domains characteristic of classical MerR proteins. Hence, the molecular function of the protein has been unclear. Here we show that BldC is indeed a DNA-binding protein and controls a regulon that includes other key developmental regulators. Intriguingly, BldC DNA-binding sites vary significantly in length. Our BldC-DNA structures explain this DNA-binding capability by revealing that BldC utilizes a DNA-binding mode distinct from MerR and other known regulators, involving asymmetric head-to-tail oligomerization on DNA direct repeats that results in dramatic DNA distortion. Notably, BldC-like proteins radiate throughout eubacteria, establishing BldC as the founding member of a new structural family of regulators.

BldC regulates the onset of differentiation in Streptomycetes by a yet unknown molecular mechanism. Using a combination of structural, biochemical and in vivo approaches, the authors show that BldC controls the transcription of several developmental regulators and unravel its DNA binding mode.

High-affinity and selective detection of pyrophosphate in water by a resorcinarene salt receptor

Pyrophosphate (PPi) is a byproduct of DNA and RNA synthesis, and abnormal levels are indicative of disease. We report the high-affinity binding of PPi in water by N-alkyl ammonium resorcinarene chloride receptors. Experimental analysis using ¹H and ³¹P NMR, isothermal titration calorimetry, mass spectrometry, and UV-vis spectroscopy all support exceptional selectivity of these systems for PPi in water. The measured affinity of K₁ = 1.60 × 10⁷ M^-1 for PPi is three orders of magnitude larger than that observed for binding to another phosphate, ATP. This exceptional anion-binding affinity in water is explored through a detailed density functional theory computational study. These systems provide a promising avenue for the development of future innovative medical diagnostic tools.

Direct Determination of Site-Specific Noncovalent Interaction Strengths of Proteins from NMR-Derived Fast Side Chain Motional Parameters

A novel approach to accurately determine residue-specific noncovalent interaction strengths (ξ) of proteins from NMR-measured fast side chain motional parameters (O_axis²) is presented. By probing the environmental sensitivity of side chain conformational energy surfaces of individual residues of a diverse set of proteins, the microscopic connections between ξ, O_axis², conformational entropy (S_conf), conformational barriers, and rotamer stabilities established here are found to be universal among proteins. The results reveal that side chain flexibility and conformational entropy of each residue decrease with increasing ξ and that for each residue type there exists a critical range of ξ, determined primarily by the mean side chain conformational barriers, within which flexibility of any residue can be reversibly tuned from highly flexible (with O_axis² ∼ 0) to highly restricted (with O_axis² ∼ 1) by increasing ξ by ∼3 kcal/mol. Beyond this critical range of ξ, both side chain flexibility and conformational entropy are insensitive to ξ. The interrelationships between conformational dynamics, conformational entropy, and noncovalent interactions of protein side chains established here open up new avenues to probe perturbation-induced (for example, ligand-binding, temperature, pressure) changes in fast side chain dynamics and thermodynamics of proteins by comparing their conformational energy surfaces in the native and perturbed states.

Molecular insights into DNA binding and anchoring by the Bacillus subtilis sporulation kinetochore-like RacA protein

During Bacillus subtilis sporulation, segregating sister chromosomes are anchored to cell poles and the chromosome is remodeled into an elongated structure called the axial filament. Data indicate that a developmentally regulated protein called RacA is involved in these functions. To gain insight into how RacA performs these diverse processes we performed a battery of structural and biochemical analyses. These studies show that RacA contains an N-terminal winged-helix-turn-helix module connected by a disordered region to a predicted coiled-coil domain. Structures capture RacA binding the DNA using distinct protein-protein interfaces and employing adjustable DNA docking modes. This unique DNA binding mechanism indicates how RacA can both specifically recognize its GC-rich centromere and also non-specifically bind the DNA. Adjacent RacA molecules within the protein-DNA structure interact leading to DNA compaction, suggesting a mechanism for axial filament formation. We also show that the RacA C-domain coiled coil directly contacts the coiled coil region of the polar protein DivIVA, which anchors RacA and hence the chromosome to the pole. Thus, our combined data reveal unique DNA binding properties by RacA and provide insight into the DNA remodeling and polar anchorage functions of the protein.

An Evolving Insight into Chiral H-Bond Catalyzed Aza-Henry Reactions: A Cooperative Role for Noncovalent Attractive Interactions Unveiled by Density Functional Theory

The role of cooperative effects arising from noncovalent attractive interactions as a vital factor governing stereoinduction in chiral H-bond catalyzed aza-Henry reactions is reported. Supporting this finding were density functional theory (DFT) calculations which revealed a shape and size dependency existed between the catalyst and substrates that when matched lead to high enantioselectivity, as reflected by favorable activation parameters. Associated with optimal catalyst and substrate pairing were a closed catalytic binding pocket and a synclinal orientation of the substrates that reinforced favorable stereoelectronic effects and dispersive type forces. Meanwhile, unfavorable steric interactions were found to be a dominant effect controlling diastereoselection.

Density Functionals for Noncovalent Interaction Energies of Biological Importance

Forty density functionals and one wavefunction method are assessed against a recently published database of accurate noncovalent interaction energies of biological importance. The comparison shows that two newly developed density functional theory (DFT) methods, PWB6K and M05-2X, give the best performance for this benchmark database of 22 noncovalent complexes, including both hydrogen-bonding and dispersion-dominated complexes. In contrast, the more popular B3LYP and PBEh functionals fail to describe the interactions in the dispersion-dominated complexes. The local spin density approximation and BHandH functionals give good performance for dispersion-dominated interactions at the expense of a large error for hydrogen bonding. PWB6K and M05-2X constitute a new generation of DFT methods based on simultaneously optimized exchange and correlation functionals that include kinetic energy density in both the exchange and correlation functional, and the present study confirms that they have greatly improved performance for noncovalent interactions as compared to previous DFT methods. We interpret this as being due to an improved treatment of medium-range correlation effects by the exchange-correlation functional. We recommend the PWB6K and M05-2X methods for investigating large biological systems and soft materials.

Structures of regulatory machinery reveal novel molecular mechanisms controlling B. subtilis nitrogen homeostasis

In Bacillus subtilis, nitrogen homeostasis is controlled by a unique circuitry composed of the regulator TnrA and the repressor GlnR. Here, Schumacher et al. describe a comprehensive molecular dissection of this network that reveals novel mechanisms, including oligomeric transformations, by which their inducible signal transduction domains are employed to provide a readout of nitrogen levels.

All cells must sense and adapt to changing nutrient availability. However, detailed molecular mechanisms coordinating such regulatory pathways remain poorly understood. In Bacillus subtilis, nitrogen homeostasis is controlled by a unique circuitry composed of the regulator TnrA, which is deactivated by feedback-inhibited glutamine synthetase (GS) during nitrogen excess and stabilized by GlnK upon nitrogen depletion, and the repressor GlnR. Here we describe a complete molecular dissection of this network. TnrA and GlnR, the global nitrogen homeostatic transcription regulators, are revealed as founders of a new structural family of dimeric DNA-binding proteins with C-terminal, flexible, effector-binding sensors that modulate their dimerization. Remarkably, the TnrA sensor domains insert into GS intersubunit catalytic pores, destabilizing the TnrA dimer and causing an unprecedented GS dodecamer-to-tetradecamer conversion, which concomitantly deactivates GS. In contrast, each subunit of the GlnK trimer “templates” active TnrA dimers. Unlike TnrA, GlnR sensors mediate an autoinhibitory dimer-destabilizing interaction alleviated by GS, which acts as a GlnR chaperone. Thus, these studies unveil heretofore unseen mechanisms by which inducible sensor domains drive metabolic reprograming in the model Gram-positive bacterium B. subtilis.

Molecular dynamics study of the recognition of dimethylated CpG sites by MBD1 protein

The cell is able to regulate which genes to express via chemical marks on the DNA and on the histone proteins. In all vertebrates, the modification on the DNA is methylation at position 5 of the two cytosines present in the dinucleotide sequence CpG. The information encoded by these chemical marks on the DNA is processed by a family of protein factors containing a conserved methyl-CpG binding domain (MBD). Essential to their function, the MBD proteins are able to bind DNA containing dimethylated CpG sites, whereas binding to unmethylated sites is not observed. In this paper, we perform molecular dynamics simulations to investigate the mechanism by which the mCpG binding domain of MBD1 is able to bind specifically dimethylated CpG sites. We find that the binding affinity of MBD1 to a DNA containing dimethylated CpG site is stronger by 26.4 kJ/mol relative to binding the same DNA but with an unmethylated CpG site. The contribution of each of the methyl groups to the change in free energy is very similar and additive. Therefore, this binding affinity (to a dimethylated DNA) is halved when considered relative to binding a hemimethylated DNA, a result that is also supported by experimental observations. Despite their equal contributions, the two methyl groups are recognized differently by MBD1. In one case, demethylation induces conformational changes in which the hydrophobic patch formed by the conserved residues Val20, Arg22, and Tyr34 moves away from the (methyl)cytosine, weakening the DNA-protein interactions. This is accompanied by an intrusion of a bulk water into the binding site at the protein-DNA interface. As a consequence, there is a reduction and rearrangements of the protein-DNA hydrogen bonds including a loss of a crucial hydrogen bond between Tyr34 and the (methyl)cytosine. The methylcytosine on the opposite strand is recognized by conformational changes of the surrounding conserved hydrophobic residues, Arg44 and Ser45, in which Arg44 participate in the 5mC-Arg-G triad. More specifically, the hydrogens of the methyl group form weak hydrogen bonds with the guanidino group and backbone carbonyl of the conserved Arg44, interactions that are absent when the cytosine is unmethylated. The results presented in this paper contribute to our knowledge of the different ways the chemical mark on the DNA is recognized by the epigenetic machinery.

Mechanistic basis of plasmid-specific DNA binding of the F plasmid regulatory protein, TraM

The conjugative transfer of bacterial F plasmids relies on TraM, a plasmid-encoded protein that recognizes multiple DNA sites to recruit the plasmid to the conjugative pore. In spite of the high degree of amino acid sequence conservation between TraM proteins, many of these proteins have markedly different DNA binding specificities that ensure the selective recruitment of a plasmid to its cognate pore. Here we present the structure of F TraM RHH (ribbon-helix-helix) domain bound to its sbmA site. The structure indicates that a pair of TraM tetramers cooperatively binds an underwound sbmA site containing 12 base pairs per turn. The sbmA is composed of 4 copies of a 5-base-pair motif, each of which is recognized by an RHH domain. The structure reveals that a single conservative amino acid difference in the RHH β-ribbon between F and pED208 TraM changes its specificity for its cognate 5-base-pair sequence motif. Specificity is also dictated by the positioning of 2-base-pair spacer elements within sbmA; in F sbmA, the spacers are positioned between motifs 1 and 2 and between motifs 3 and 4, whereas in pED208 sbmA, there is a single spacer between motifs 2 and 3. We also demonstrate that a pair of F TraM tetramers can cooperatively bind its sbmC site with an affinity similar to that of sbmA in spite of a lack of sequence similarity between these DNA elements. These results provide a basis for the prediction of the DNA binding properties of the family of TraM proteins.

Kinetics, structure, and mechanism of 8-Oxo-7,8-dihydro-2'-deoxyguanosine bypass by human DNA polymerase η

DNA damage incurred by a multitude of endogenous and exogenous factors constitutes an inevitable challenge for the replication machinery. Cells rely on various mechanisms to either remove lesions or bypass them in a more or less error-prone fashion. The latter pathway involves the Y-family polymerases that catalyze trans-lesion synthesis across sites of damaged DNA. 7,8-Dihydro-8-oxo-2'-deoxyguanosine (8-oxoG) is a major lesion that is a consequence of oxidative stress and is associated with cancer, aging, hepatitis, and infertility. We have used steady-state and transient-state kinetics in conjunction with mass spectrometry to analyze in vitro bypass of 8-oxoG by human DNA polymerase η (hpol η). Unlike the high fidelity polymerases that show preferential insertion of A opposite 8-oxoG, hpol η is capable of bypassing 8-oxoG in a mostly error-free fashion, thus preventing GC→AT transversion mutations. Crystal structures of ternary hpol η-DNA complexes and incoming dCTP, dATP, or dGTP opposite 8-oxoG reveal that an arginine from the finger domain assumes a key role in avoiding formation of the nascent 8-oxoG:A pair. That hpol η discriminates against dATP exclusively at the insertion stage is confirmed by structures of ternary complexes that allow visualization of the extension step. These structures with G:dCTP following either 8-oxoG:C or 8-oxoG:A pairs exhibit virtually identical active site conformations. Our combined data provide a detailed understanding of hpol η bypass of the most common oxidative DNA lesion.

Cations form sequence selective motifs within DNA grooves via a combination of cation-pi and ion-dipole/hydrogen bond interactions

The fine conformational subtleties of DNA structure modulate many fundamental cellular processes including gene activation/repression, cellular division, and DNA repair. Most of these cellular processes rely on the conformational heterogeneity of specific DNA sequences. Factors including those structural characteristics inherent in the particular base sequence as well as those induced through interaction with solvent components combine to produce fine DNA structural variation including helical flexibility and conformation. Cation-pi interactions between solvent cations or their first hydration shell waters and the faces of DNA bases form sequence selectively and contribute to DNA structural heterogeneity. In this paper, we detect and characterize the binding patterns found in cation-pi interactions between solvent cations and DNA bases in a set of high resolution x-ray crystal structures. Specifically, we found that monovalent cations (Tl⁺) and the polarized first hydration shell waters of divalent cations (Mg²⁺, Ca²⁺) form cation-pi interactions with DNA bases stabilizing unstacked conformations. When these cation-pi interactions are combined with electrostatic interactions a pattern of specific binding motifs is formed within the grooves.

A common mode of recognition for methylated CpG

Much is known about vertebrate DNA methylation, however it is not known how methylated CpG within particular sequences is recognized. Two recent structures of C2H2 zinc finger (ZnF) proteins in complex with methylated DNA reveal a common recognition mode for 5-methylcytosine (5mC) that involves a 5mC-Arg-G triad. In the two ZnF proteins, an arginine that precedes the first Zn-binding histidine (RH motif) can interact with a 5mCpG or TpG dinucleotide. Among a family of >300 human Krüppel-associated box (KRAB) domain containing ZnF proteins examined, two-thirds contained at least one ZnF that included an RH motif. We propose that the RH-ZnF motifs provide specificity for 5mCpG, whereas the neighboring Zn fingers recognize the surrounding DNA sequence context.

The Sum1/Ndt80 transcriptional switch and commitment to meiosis in Saccharomyces cerevisiae

Cells encounter numerous signals during the development of an organism that induce division, differentiation, and apoptosis. These signals need to be present for defined intervals in order to induce stable changes in the cellular phenotype. The point after which an inducing signal is no longer needed for completion of a differentiation program can be termed the "commitment point." Meiotic development in the yeast Saccharomyces cerevisiae (sporulation) provides a model system to study commitment. Similar to differentiation programs in multicellular organisms, the sporulation program in yeast is regulated by a transcriptional cascade that produces early, middle, and late sets of sporulation-specific transcripts. Although critical meiosis-specific events occur as early genes are expressed, commitment does not take place until middle genes are induced. Middle promoters are activated by the Ndt80 transcription factor, which is produced and activated shortly before most middle genes are expressed. In this article, I discuss the connection between Ndt80 and meiotic commitment. A transcriptional regulatory pathway makes NDT80 transcription contingent on the prior expression of early genes. Once Ndt80 is produced, the recombination (pachytene) checkpoint prevents activation of the Ndt80 protein. Upon activation, Ndt80 triggers a positive autoregulatory loop that leads to the induction of genes that promote exit from prophase, the meiotic divisions, and spore formation. The pathway is controlled by multiple feed-forward loops that give switch-like properties to the commitment transition. The conservation of regulatory components of the meiotic commitment pathway and the recently reported ability of Ndt80 to increase replicative life span are discussed.

Optimization of the CHARMM additive force field for DNA: Improved treatment of the BI/BII conformational equilibrium

The B-form of DNA can populate two different backbone conformations: BI and BII, defined by the difference between the torsion angles ε and ζ (BI = ε-ζ < 0 and BII = ε-ζ > 0). BI is the most populated state, but the population of the BII state, which is sequence dependent, is significant and accumulating evidence shows that BII affects the overall structure of DNA, and thus influences protein-DNA recognition. This work presents a reparametrization of the CHARMM27 additive nucleic acid force field to increase the sampling of the BII form in MD simulations of DNA. In addition, minor modifications of sugar puckering were introduced to facilitate sampling of the A form of DNA under the appropriate environmental conditions. Parameter optimization was guided by quantum mechanical data on model compounds, followed by calculations on several DNA duplexes in the condensed phase. The selected optimized parameters were then validated against a number of DNA duplexes, with the most extensive tests performed on the EcoRI dodecamer, including comparative calculations using the Amber Parm99bsc0 force field. The new CHARMM model better reproduces experimentally observed sampling of the BII conformation, including sampling as a function of sequence. In addition, the model reproduces the A form of the 1ZF1 duplex in 75 % ethanol, and yields a stable Z-DNA conformation of duplex (GTACGTAC) in its crystal environment. The resulting model, in combination with a recent reoptimization of the CHARMM27 force field for RNA, will be referred to as CHARMM36.

Energies and physicochemical properties of cation-π interactions in biological structures

The cation-π interactions occur frequently within or between proteins due to six (Phe, Tyr, Trp, Arg, Lys, and His) of the twenty natural amino acids potentially interacting with metallic cations via these interactions. In this study, quantum chemical calculations and molecular orbital (MO) theory are used to study the energies and properties of cation-π interactions in biological structures. The cation-π interactions of H⁺ and Li⁺ are similar to hydrogen bonds and lithium bonds, respectively, in which the small, naked cations H⁺ and Li⁺ are buried deep within the π-electron density of aromatic molecules, forming stable cation-π bonds that are much stronger than the cation-π interactions of other alkali metal cations. The cation-π interactions of metallic cations with atomic masses greater than that of Li⁺ arise mainly from the coordinate bond comprising empty valence atomic orbitals (AOs) of metallic cations and π-MOs of aromatic molecules, though electrostatic interactions may also contribute to the cation-π interaction. The binding strength of cation-π interactions is determined by the charge and types of AOs in the metallic cations. Cation-π interaction energies are distance- and orientation-dependent; energies decrease with the distance (r) and the orientation angle (θ). In solution, the cation-π energies decrease with the increase of the dielectric constant (ɛ) of the solvent; however, solvation has less influence on the H⁺-π and H₃O⁺-π interactions than on interactions with other cations. The conclusions from this study provide useful theoretical insights into the nature of cation-π interactions and may contribute to the development of better force field parameters for describing the molecular dynamics of cation-π interactions within and between proteins.

Recognition of methylated DNA through methyl-CpG binding domain proteins

DNA methylation is a key regulatory control route in epigenetics, involving gene silencing and chromosome inactivation. It has been recognized that methyl-CpG binding domain (MBD) proteins play an important role in interpreting the genetic information encoded by methylated DNA (mDNA). Although the function of MBD proteins has attracted considerable attention and is well characterized, the mechanism underlying mDNA recognition by MBD proteins is still poorly understood. In this article, we demonstrate that the methyl-CpG dinucleotides are recognized at the MBD–mDNA interface by two MBD arginines through an interplay of hydrogen bonding and cation-π interaction. Through molecular dynamics and quantum-chemistry calculations we investigate the methyl-cytosine recognition process and demonstrate that methylation enhances MBD–mDNA binding by increasing the hydrophobic interfacial area and by strengthening the interaction between mDNA and MBD proteins. Free-energy perturbation calculations also show that methylation yields favorable contribution to the binding free energy for MBD–mDNA complex.

Empirical formulation and parameterization of cation-π interactions for protein modeling

Cation-π interaction is comparable and as important as other main molecular interaction types, such as hydrogen bond, electrostatic interaction, van der Waals interaction, and hydrophobic interaction. Cation-π interactions frequently occur in protein structures, because six (Phe, Tyr, Trp, Arg, Lys, and His) of 20 natural amino acids and all metallic cations could be involved in cation-π interaction. Cation-π interactions arise from complex physicochemical nature and possess unique interaction behaviors, which cannot be modeled and evaluated by existing empirical equations and force field parameters that are widely used in the molecular dynamics. In this study, the authors present an empirical approach for cation-π interaction energy calculations in protein interactions. The accurate cation-π interaction energies of aromatic amino acids (Phe, Tyr, and Try) with protonated amino acids (Arg and Lys) and metallic cations (Li(+), Na(+), K(+), and Ca(2+)) are calculated using B3LYP/6-311+G(d,p) method as the benchmark for the empirical formulization and parameterization. Then, the empirical equations are built and the parameters are optimized based on the benchmark calculations. The cation-π interactions are distance and orientation dependent. Correspondingly, the empirical equations of cation-π interactions are functions of two variables, the distance r and the orientation angle θ. Two types of empirical equations of cation-π interactions are proposed. One is a modified distance and orientation dependent Lennard-Jones equation. The second is a polynomial function of two variables r and θ. The amino acid-based empirical equations and parameters provide simple and useful tools for evaluations of cation-π interaction energies in protein interactions.

Putting anion-π interactions into perspective

Supramolecular chemistry is a field of scientific exploration that probes the relationship between molecular structure and function. It is the chemistry of the noncovalent bond, which forms the basis of highly specific recognition, transport, and regulation events that actuate biological processes. The classic design principles of supramolecular chemistry include strong, directional interactions like hydrogen bonding, halogen bonding, and cation-π complexation, as well as less directional forces like ion pairing, π-π, solvophobic, and van der Waals potentials. In recent years, the anion-π interaction (an attractive force between an electron-deficient aromatic π system and an anion) has been recognized as a hitherto unexplored noncovalent bond, the nature of which has been interpreted through both experimental and theoretical investigations. The design of selective anion receptors and channels based on this interaction represent important advances in the field of supramolecular chemistry. The objectives of this Review are 1) to discuss current thinking on the nature of this interaction, 2) to survey key experimental work in which anion-π bonding is demonstrated, and 3) to provide insights into the directional nature of anion-π contact in X-ray crystal structures.

Structures of carbon catabolite protein A-(HPr-Ser46-P) bound to diverse catabolite response element sites reveal the basis for high-affinity binding to degenerate DNA operators

In Gram-positive bacteria, carbon catabolite protein A (CcpA) is the master regulator of carbon catabolite control, which ensures optimal energy usage under diverse conditions. Unlike other LacI-GalR proteins, CcpA is activated for DNA binding by first forming a complex with the phosphoprotein HPr-Ser46-P. Bacillus subtilis CcpA functions as both a transcription repressor and activator and binds to more than 50 operators called catabolite response elements (cres). These sites are highly degenerate with the consensus, WTGNNARCGNWWWCAW. How CcpA–(HPr-Ser46-P) binds such diverse sequences is unclear. To gain insight into this question, we solved the structures of the CcpA–(HPr-Ser46-P) complex bound to three different operators, the synthetic (syn) cre, ackA2 cre and gntR-down cre. Strikingly, the structures show that the CcpA-bound operators display different bend angles, ranging from 31° to 56°. These differences are accommodated by a flexible linkage between the CcpA helix-turn-helix-loop-helix motif and hinge helices, which allows independent docking of these DNA-binding modules. This flexibility coupled with an abundance of non-polar residues capable of non-specific nucleobase interactions permits CcpA–(HPr-Ser46-P) to bind diverse operators. Indeed, biochemical data show that CcpA–(HPr-Ser46-P) binds the three cre sites with similar affinities. Thus, the data reveal properties that license this protein to function as a global transcription regulator.

A generalized conformational energy function of DNA derived from molecular dynamics simulations

Proteins recognize DNA sequences by two different mechanisms. The first is direct readout, in which recognition is mediated by direct interactions between the protein and the DNA bases. The second is indirect readout, which is caused by the dependence of conformation and the deformability of the DNA structure on the sequence. Various energy functions have been proposed to evaluate the contribution of indirect readout to the free-energy changes in complex formations. We developed a new generalized energy function to estimate the dependence of the deformability of DNA on the sequence. This function was derived from molecular dynamics simulations previously conducted on B-DNA dodecamers, each of which had one possible tetramer sequence embedded at its center. By taking the logarithm of the probability distribution function (PDF) for the base-step parameters of the central base-pair step of the tetramer, its ability to distinguish the native sequence from random ones was superior to that with the previous method that approximated the energy function in harmonic form. From a comparison of the energy profiles calculated with these two methods, we found that the harmonic approximation caused significant errors in the conformational energies of the tetramers that adopted multiple stable conformations.

Investigation of transcription factor Ndt80 affinity differences for wild type and mutant DNA: a molecular dynamics study

Molecular dynamics simulations and free energy calculations have been performed on the transcription factor Ndt80 either in complex with the native DNA sequence or with a mutant DNA with a switched central base pair, C5-G5' to G5-C5'. This mutant has been shown to have a 100-fold decrease in binding affinity of Ndt80, and in this study we explain this both structurally and energetically. The major interactions between the protein and the DNA were maintained in the simulations, apart from around the mutation site. The crystal structure of the Ndt80-DNA complex revealed that R177 makes a base specific bidentate interaction with G5' which is part of a conserved 5'-YpG-3' step. In the simulation with the mutant DNA, the side chain of R177 changes conformation and makes three new stable hydrogen bonds to the DNA backbone. This in turn induces a conformational change in the DNA backbone of the T6'-G5' step from the unusual BII state to the canonical BI state. The affinity difference for the protein-DNA complex with the native DNA compared with the mutant DNA is only about 3 kcal/mol. The free energy calculations of the base pair switch indicated a larger difference than what was found experimentally, about 7.7 kcal/mol, but this is explained in structural terms using the 10 ns simulations of the solvated complexes and the rearrangement of the R177 side chain.

Sequence-dependent DNA deformability studied using molecular dynamics simulations

Proteins recognize specific DNA sequences not only through direct contact between amino acids and bases, but also indirectly based on the sequence-dependent conformation and deformability of the DNA (indirect readout). We used molecular dynamics simulations to analyze the sequence-dependent DNA conformations of all 136 possible tetrameric sequences sandwiched between CGCG sequences. The deformability of dimeric steps obtained by the simulations is consistent with that by the crystal structures. The simulation results further showed that the conformation and deformability of the tetramers can highly depend on the flanking base pairs. The conformations of xATx tetramers show the most rigidity and are not affected by the flanking base pairs and the xYRx show by contrast the greatest flexibility and change their conformations depending on the base pairs at both ends, suggesting tetramers with the same central dimer can show different deformabilities. These results suggest that analysis of dimeric steps alone may overlook some conformational features of DNA and provide insight into the mechanism of indirect readout during protein-DNA recognition. Moreover, the sequence dependence of DNA conformation and deformability may be used to estimate the contribution of indirect readout to the specificity of protein-DNA recognition as well as nucleosome positioning and large-scale behavior of nucleic acids.

Structure of Aart, a designed six-finger zinc finger peptide, bound to DNA

Cys2-His2 zinc fingers are one of the most common types of DNA-binding domains. Modifications to zinc-finger binding specificity have recently enabled custom DNA-binding proteins to be designed to a wide array of target sequences. We present here a 1.96 A structure of Aart, a designed six-zinc finger protein, bound to a consensus DNA target site. This is the first structure of a designed protein with six fingers, and was intended to provide insights into the unusual affinity and specificity characteristics of this protein. Most protein-DNA contacts were found to be consistent with expectations, while others were unanticipated or insufficient to explain specificity. Several were unexpectedly mediated by glycerol, water molecules or amino acid-base stacking interactions. These results challenge some conventional concepts of recognition, particularly the finding that triplets containing 5'A, C, or T are typically not specified by direct interaction with the amino acid in position 6 of the recognition helix.

Principles of protein-DNA recognition revealed in the structural analysis of Ndt80-MSE DNA complexes

The Saccharomyces cerevisiae transcription factor Ndt80 selectively binds a DNA consensus sequence (the middle sporulation element [MSE]) to activate gene expression after the successful completion of meiotic recombination. Here we report the X-ray crystal structures of Ndt80 bound to ten distinct MSE variants. Comparison of these structures with the structure of Ndt80 bound to a consensus MSE reveals structural principles that determine the DNA binding specificity of this transcription factor. The 5' GC-rich end of the MSE contains distinct 5'-YpG-3' steps that are recognized by arginine side chains through a combination of hydrogen bonding and cation-pi interactions. The 3' AT-rich region is recognized via minor groove contacts that sterically exclude the N2 atom of GC base pairs. The conformation of the AT-rich region is fixed by interactions with the protein that favor recognition of poly(A)-poly(T) versus mixed AT sequences through an avoidance of major groove steric clashes at 5'-ApT-3' steps.

Identification and characterization of the methyl arginines in the fragile X mental retardation protein Fmrp

Fragile X syndrome is the most common form of inherited mental retardation and is caused by the absence of expression of the FMR1 gene. The protein encoded by this gene, Fmrp, is an RNA-binding protein that binds a subset of mRNAs and regulates their translation, leading to normal cognitive function. Although the association with RNAs is well established, it is still unknown how Fmrp finds and assembles with its RNA cargoes and how these activities are regulated. We show here that Fmrp is post-translationally methylated, primarily on its arginine-glycine-glycine box. We identify the four arginines that are methylated and show that cellular Fmrp is monomethylated and asymmetrically dimethylated. We also show that the autosomal paralog Fxr1 and the Drosophila ortholog dFmr1 are methylated post-translationally. Recombinant protein arginine methyl transferase 1 (PRMT1) methylates Fmrp on the same arginines in vitro as in cells. In vitro methylation of Fmrp results in reduced binding to the minimal RNA sequence sc1, which encodes a stem loop G-quartet structure. Our data identify an additional mechanism, arginine methylation, for modifying Fmrp function and suggest that methylation occurs to limit or modulate RNA binding by Fmrp.

Protein-DNA Recognition Patterns and Predictions

Structural data on protein-DNA complexes provide clues for understanding the mechanism of protein-DNA recognition. Although the structures of a large number of protein-DNA complexes are known, the mechanisms underlying their specific binding are still only poorly understood. Analysis of these structures has shown that there is no simple one-to-one correspondence between bases and amino acids within protein-DNA complexes; nevertheless, the observed patterns of interaction carry important information on the mechanisms of protein-DNA recognition. In this review, we show how the patterns of interaction, either observed in known structures or derived from computer simulations, confer recognition specificity, and how they can be used to examine the relationship between structure and specificity and to predict target DNA sequences used by regulatory proteins.

Multicoefficient Extrapolated Density Functional Theory Studies of π···π Interactions: The Benzene Dimer

We report tests of new (2005) and established (1999-2003) multilevel methods against essentially converged benchmark results for nonbonded interactions in benzene dimers. We found that the newly developed multicoefficient extrapolated density functional theory (DFT) methods (which combine DFT with correlated wave function methods) give better performance than multilevel methods such as G3SX, G3SX(MP3), and CBS-QB3 that are based purely on wave function theory (WFT); furthermore, they have a lower computational cost. We conclude that our empirical approach for combining WFT methods with DFT methods is a very efficient and effective way for describing not only covalent interactions (as shown previously) but also nonbonded interactions.

Structural analysis of cation-pi interactions in DNA binding proteins

Cation-pi interactions play an important role in the stability of protein structures. In this work, we have analyzed the influence of cation-pi interactions in DNA binding proteins. We observed cation-pi interactions in 45 out of 62 DNA binding proteins and there is no significant correlation between the number of amino acid residues and number of cation-pi interactions. These interactions are mainly formed by long-range contacts, and the role of short and medium-range contacts is minimal. The preference of Arg is higher than Lys to form cation-pi interactions. The pair-wise cation-pi interaction energy between aromatic and positively charged residues shows that Arg-Tyr energy is the strongest among the possible six pairs. The structural analysis of cation-pi interaction forming residues shows that Lys, Trp, and Tyr prefer to be in the binding site of protein-DNA complexes. Further, the accessible surface areas of cation-pi interaction forming cationic residues are significantly less than that of other residues. The preference of cation-pi interaction forming residues in different secondary structures shows that Lys prefers to be in strand and Phe prefers to be in turn regions. The results obtained in the present study will be useful in understanding the contribution of cation-pi interactions to the stability and specificity of protein-DNA complexes.

Cation-π interactions in protein-protein interfaces

Arginine is an abundant residue in protein-protein interfaces. The importance of this residue relates to the versatility of its side chain in intermolecular interactions. Different classes of protein-protein interfaces were surveyed for cation-pi interactions. Approximately half of the protein complexes and one-third of the homodimers analyzed were found to contain at least one intermolecular cation-pi pair. Interactions between arginine and tyrosine were found to be the most abundant. The electrostatic interaction energy was calculated to be approximately 3 kcal/mol, on average. A distance-based search of guanidinium:aromatic interactions was also performed using the Macromolecular Structure Database (MSD). This search revealed that half of the guanidinium:aromatic pairs pack in a coplanar manner. Furthermore, it was found that the cationic group of the cation-pi pair is frequently involved in intermolecular hydrogen bonds. In this manner the arginine side chain can participate in multiple interactions, providing a mechanism for inter-protein specificity. Thus, the cation-pi interaction is established as an important contributor to protein-protein interfaces.

How well can new-generation density functional methods describe stacking interactions in biological systems?

We compare the performance of four recently developed DFT methods (MPW1B95, MPWB1K, PW6B95, and PWB6K) and two previous, generally successful DFT methods (B3LYP and B97-1) for the calculation of stacking interactions in six nucleic acid bases complexes and five amino acid pairs and for the calculation of hydrogen bonding interactions in two Watson-Crick type base pairs. We found that the four newly developed DFT methods give reasonable results for the stacking interactions in the DNA base pairs and amino acid pairs, whereas the previous DFT methods fail to describe interactions in these stacked complexes. We conclude that the new generation of DFT methods have greatly improved performance for stacking interaction as compared to previously available methods. We recommend the PWB6K method for investigating large DNA or protein systems where stacking plays an important role.

The N6-Methyl Group of Adenine Further Increases the BI Stability of DNA Compared to C5-Methyl Groups

Methylated DNA bases are natural modifications which play an important role in protein-DNA interactions. Recent experimental and theoretical results have shown an influence of the base modification on the conformational behavior of the DNA backbone. MD simulations of four different B-DNA dodecamers (d(GC)(6), d(AT)(6), d(G(5mCG)(5)C), and d(A(T6mA)(5)T)) have been performed with the aim to examine the influence of methyl groups on the B-DNA backbone behavior. An additional control simulation of d(AU)(6) has also been performed to examine the further influence of the C5-methyl group in thymine. Methyl groups in the major groove (as in C5-methylcytosine, thymine, or N6-methyladenine) decrease the BII substate population of RpY steps. Due to methylation a clearer distinction of the BI substate stability between YpR and RpY (CpG/GpC or TpA/ApT) steps arises. A positive correlation between the BII substate population and base stacking distances is seen only for poly(GC). A methyl group added into the major groove increases mean water residence times around the purine N7 atom, which may stabilize the BI substate by improving the hydration network between the DNA backbone and the major groove. The N6-methyl group also forms a water molecule bridge between the N6 and O4 atoms, and thus further stabilizes the BI substate.

Stair motifs at protein-DNA interfaces: nonadditivity of H-bond, stacking, and cation-pi interactions

At the interface between protein and double-stranded DNA, stair motifs simultaneously involve three different types of pairwise interactions: aromatic base stacking, hydrogen bonding, and cation-pi. The relative importance of these interactions is studied in the stair motif occurring in the 1TC3 crystal structure, which involves an arginine and two stacked guanines, by means of Hartree-Fock (HF) and Møller-Plesset energy and free energy calculations, including vibrational, rotational, translational contributions, both in a vacuum and various solvents. The results obtained show an anti-cooperative tendency of the HF energy and vibrational free energy terms, and the cooperativity of the rotational, translational, and solvation free energies. Hence, the cooperativity of the stair motif interactions, in the context of protein-DNA recognition, can be viewed as arising from the environment.

Characterization of critical interactions between Ndt80 and MSE DNA defining a novel family of Ig-fold transcription factors

The Ndt80 protein of the yeast Saccharomyces cerevisiae is the founding member of a new sub-family of proteins in the Ig-fold superfamily of transcription factors. The crystal structure of Ndt80 bound to DNA shows that it makes contacts through several loops on one side of the protein that connect beta-strands which form the beta-sandwich fold common to proteins in this superfamily. However, the DNA-binding domain of Ndt80 is considerably larger than many other members of the Ig-fold superfamily and it appears to make a larger number of contacts with the DNA than these proteins. To determine the contribution of each of these contacts and to examine if the mechanism of Ndt80 DNA binding was similar to other members of the Ig-fold superfamily, amino acid substitutions were introduced at each residue that contacts the DNA and assayed for their effect on Ndt80 activity. Many of the mutations caused significant decreases in DNA-binding affinity and transcriptional activation. Several of these are in residues that are not found in other sub-families of Ig-fold proteins. These additional contacts are likely responsible for Ndt80's ability to bind DNA as a monomer while most other members require additional domains or cofactors to recognize their sites.