PubMed Health⌕ Search

Biomedical subjects

Nilesh K Banavali

Publications and source records attributed to Nilesh K Banavali.

7 recordsLinked to original sources

Conserved patterns in backbone torsional changes allow for single base flipping from duplex DNA with minimal distortion of the double helix.

Base flipping is a structural mechanism common to many DNA processing and repair enzymes. Changes in the local backbone torsions that occur during base flipping and the effect of environment on their behavior are of particular interest in understanding different base flipping mechanisms. In the present study, structures sampled during umbrella sampling molecular dynamics (MD) simulations of base flipping in aqueous and protein-bound environments, carried out with two different MD simulation strategies, are analyzed to find the most significant phosphodiester backbone distortions in the vicinity of the flipping base. Torsional sampling on the 5' side of the flipping base during flipping through the major groove shows similarities to the torsional sampling on the 3' side during flipping through the minor groove and vice versa. In differing environments, this behavior varies only marginally. These compensating torsional changes in the DNA backbone on 5' and 3' sides of the flipping base limit overall distortion of the DNA double helix during single base flipping. Rotameric intermediate states observed during base flipping are identified and postulated to be metastable states implicated in both large-scale structural changes and functional effects of chemical modifications in DNA.

DNA↗

Free energy landscape of A-DNA to B-DNA conversion in aqueous solution.

The interconversion between the well-characterized A- and B-forms of DNA is a structural transition for which the intermediate states and the free energy difference between the two endpoints are not known precisely. In the present study, the difference between the Root Mean Square Distance (RMSD) from canonical A-form and B-form DNA is used as an order parameter to characterize this free energy difference using umbrella sampling molecular dynamics (MD) simulations with explicit solvent. The constraint imposed along this order parameter allows relatively unrestricted evolution of the intermediate structures away from both canonical A- and B-forms. The free energy difference between the A- and B-forms for the hexamer DNA sequence CTCGAG in aqueous solution is conservatively estimated to be at least 2.8 kcal/mol. A continuum of intermediate structures with no well-defined local minima links the two forms. The absence of any major barriers in the free energy surface is consistent with spontaneous conversion of the A-form DNA to B-form DNA in unconstrained simulations. The extensive sampling in the MD simulations (>0.1 mus) also allowed quantitative energetic characterization of local backbone conformational variables such as sugar pseudorotation angles and BI/BII state equilibria and their dependence on base identity. The absolute minimum in the calculated free energy profile corresponds closely to the crystal structure of the hexamer sequence, indicating that the present method has the potential to identify the most stable state for an arbitrary DNA sequence in water.

Computer Simulation↗

The N-terminal end of the catalytic domain of SRC kinase Hck is a conformational switch implicated in long-range allosteric regulation.

Signal transduction in cell growth and proliferation involves regulation of kinases through long-range allostery between remote protein regions. Molecular dynamics free energy calculations are used to clarify the coupling between the catalytic domain of Src kinase Hck and its N-terminal end connecting to the regulatory SH2 and SH3 modules. The N-terminal end is stable in the orientation required for the regulatory modules to remain properly bound only in the inactive catalytic domain. In the active catalytic domain, the N-terminal end prefers a different conformation consistent with dissociation of the regulatory modules. The free energy surface shows that the N-terminal end acts as a reversible two-state conformational switch coupling the catalytic domain to the regulatory modules. Structural analogy with insulin receptor kinase and c-Src suggests that such reversible conformational switching in a critical hinge region could be a common mechanism in long-range allosteric regulation of protein kinase activity.

Allosteric Regulation↗

Caught in the act: visualization of an intermediate in the DNA base-flipping pathway induced by HhaI methyltransferase.

Rotation of a DNA or RNA nucleotide out of the double helix and into a protein pocket ('base flipping') is a mechanistic feature common to some DNA/RNA-binding proteins. Here, we report the structure of HhaI methyltransferase in complex with DNA containing a south-constrained abasic carbocyclic sugar at the target site in the presence of the methyl donor byproduct AdoHcy. Unexpectedly, the locked south pseudosugar appears to be trapped in the middle of the flipping pathway via the DNA major groove, held in place primarily through Van der Waals contacts with a set of invariant amino acids. Molecular dynamics simulations indicate that the structural stabilization observed with the south-constrained pseudosugar will not occur with a north-constrained pseudosugar, which explains its lowered binding affinity. Moreover, comparison of structural transitions of the sugar and phosphodiester backbone observed during computational studies of base flipping in the M.HhaI-DNA-AdoHcy ternary complex indicate that the south-constrained pseudosugar induces a conformation on the phosphodiester backbone that corresponds to that of a discrete intermediate of the base-flipping pathway. As previous crystal structures of M.HhaI ternary complex with DNA displayed the flipped sugar moiety in the antipodal north conformation, we suggest that conversion of the sugar pucker from south to north beyond the middle of the pathway is an essential part of the mechanism through which flipping must proceed to reach its final destination. We also discuss the possibility of the south-constrained pseudosugar mimicking a transition state in the phosphodiester and sugar moieties that occurs during DNA base flipping in the presence of M.HhaI.

Binding Sites↗

Protein-facilitated base flipping in DNA by cytosine-5-methyltransferase.

DNA methylation, various DNA repair mechanisms, and possibly early events in the opening of DNA as required for transcription and replication are initiated by flipping of a DNA base out of the DNA double helix. The energetics and structural mechanism of base flipping in the presence of the DNA-processing enzyme, cytosine 5-methyltransferase from HhaI (M.HhaI), were obtained through molecular dynamics based upon free-energy calculations. Free-energy profiles for base flipping show that, when in the closed conformation, M.HhaI lowers the free-energy barrier to flipping by 17 kcalmol and stabilizes the fully flipped state. Flipping is shown to occur via the major groove of the DNA. Structural analysis indicates that flipping is facilitated by destabilization of the DNA double-helical structure and substitution of DNA base-pairing and base-stacking interactions with DNA-protein interactions. The fully flipped state is stabilized by DNA-protein interactions that are enhanced upon binding of coenzyme. This study represents an atomic detail description of the mechanism by which a protein facilitates specific structural distortion in DNA.

Base Composition↗

Free energy and structural pathways of base flipping in a DNA GCGC containing sequence.

Structural distortions of DNA are essential for its biological function due to the genetic information of DNA not being physically accessible in the duplex state. Base flipping is one of the simplest structural distortions of DNA and may represent an initial event in strand separation required to access the genetic code. Flipping is also utilized by DNA-modifying and repair enzymes to access specific bases. It is typically thought that base flipping (or base-pair opening) occurs via the major groove whereas minor groove flipping is only possible when mediated by DNA-binding proteins. Here, umbrella sampling with a novel center-of-mass pseudodihedral reaction coordinate was used to calculate the individual potentials of mean force (PMF) for flipping of the Watson-Crick (WC) paired C and G bases in the CCATGCGCTGAC DNA dodecamer. The novel reaction coordinate allowed explicit investigation of the complete flipping process via both the minor and major groove pathways. The minor and major groove barriers to flipping are similar for C base flipping while the major groove barrier is slightly lower for G base flipping. Minor groove flipping requires distortion of the WC partner while the flipping base pulls away from its partner during major groove flipping. The flipped states are represented by relatively flat free energy surfaces, with a small, local minimum observed for the flipped G base. Conserved patterns of phosphodiester backbone dihedral distortions during flipping indicate their essential role in the flipping process. During flipping, the target base tracks along the respective grooves, leading to hydrogen-bonding interactions with neighboring base-pairs. Such hydrogen-bonding interactions with the neighboring sequence suggest a novel mechanism of sequence dependence in DNA dynamics.

Base Pairing↗

DNA bending induced by carbocyclic sugar analogs constrained to the north conformation.

DNA bending caused by introduction of carbocyclic sugars constrained to the north conformation was studied, using explicit solvent molecular dynamic (MD) simulations. The native Drew-Dickerson (DD) dodecamer and its three modifications containing north carbocyclic sugars in the 7th (T7*), 8th (T8*) or both 7th and 8th (T7T8*) nucleotide positions were examined. Introduction of the carbocyclic sugar results in A-form conformations for the alpha, beta, chi, zeta, and sugar pucker backbone parameters in the modified nucleotides. Increased steric repulsion between the sugar and its parent base in the modified oligonucleotides impacts the roll and cup dinucleotide step parameters, increasing the bending of the oligomer axis. Increased buckling of the substituted nucleotides disrupts the usual stabilizing base stacking interactions. The level of overall bending depends on the number and position of carbocyclic sugars introduced in the DNA sequence. Single sugar substitutions are unable to induce substantial bending due to the neighboring unmodified nucleotides counterbalancing the distortion. Significant bending can, however, be induced by two consecutive north sugars (T7T8*), which is in agreement with experimental results. The modified oligomers populate a wide range of bend angles, indicating that they maintain flexibility in the bent state. The present results suggest that insertion of carbocyclic sugars into DNA or RNA duplexes can be used to engineer bending of the duplexes without impacting the electrostatic or chemical properties of the phosphodiester backbone, thereby serving as excellent tools for experimental elucidation of nucleic acid structure-function relationships.

Base Pairing↗