1
|
Stacking geometry between two sheared Watson-Crick basepairs: Computational chemistry and bioinformatics based prediction. Biochim Biophys Acta Gen Subj 2020; 1864:129600. [PMID: 32179130 DOI: 10.1016/j.bbagen.2020.129600] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2019] [Revised: 03/05/2020] [Accepted: 03/11/2020] [Indexed: 11/21/2022]
Abstract
BACKGROUND Molecular modeling of RNA double helices is possible using most probable values of basepair parameters obtained from crystal structure database. The A:A w:wC non-canonical basepair, involving Watson-Crick edges of two Adenines in cis orientation, appears quite frequently in database. Bimodal distribution of its Shear, due to two different H-bonding schemes, introduces the confusion in assigning most the probable value. Its effect is pronounced when the A:A w:wC basepair stacks on Sheared wobble G:U W:WC basepairs. METHODS We employed molecular dynamics simulations of three possible double helices with GAG, UAG and GAU sequence motifs at their centers and quantum chemical calculation for non-canonical A:A w:wC basepair stacked on G:U W:WC basepair. RESULTS We noticed stable structures of GAG motif with specifically negative Shear of the A:A basepair but stabilities of the other motifs were not found with A:A w:wC basepairing. Hybrid DFT-D and MP2 stacking energy analyses on dinucleotide step sequences, A:A w:wC::G:U W:WC and A:A w:wC::U:G W:WC reveal that viable orientation of A:A::G:U prefers one of the H-bonding modes with negative Shear, supported by crystal structure database. The A:A::U:G dinucleotide, however, prefers structure with only positive Shear. CONCLUSIONS The quantum chemical calculations explain why MD simulations of GAG sequence motif only appear stable. In the cases of the GAU and UAG motifs "tug of war" situation between positive and negative Shears of A:A w:wC basepair induces conformational plasticity. GENERAL SIGNIFICANCE We have projected comprehensive reason behind the promiscuous nature of A:A w:wC basepair which brings occasional structural plasticity.
Collapse
|
2
|
Mondal M, Bhattacharyya D, Gao YQ. Structural properties and influence of solvent on the stability of telomeric four-stranded i-motif DNA. Phys Chem Chem Phys 2019; 21:21549-21560. [PMID: 31536074 DOI: 10.1039/c9cp03253c] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Repetitive cytosine rich i-motif forming sequences are abundant in the telomere, centromere and promoters of several oncogenes and in some instances are known to regulate transcription and gene expression. The in vivo existence of i-motif structures demands further insight into the factors affecting their formation and stability and development of better understanding of their gene regulatory functions. Most prior studies characterizing the conformational dynamics of i-motifs are based on i-motif forming synthetic constructs. Here, we present a systematic study on the stability and structural properties of biologically relevant i-motifs of telomeric and centromeric repeat fragments. Our results based on molecular dynamics simulations and quantum chemical calculations indicate that along with base pairing interactions within the i-motif core the overall folded conformation is associated with the stable C-HO sugar "zippers" in the narrow grooves and structured water molecules along the wide grooves. The stacked geometry of the hemi-protonated cytosine pairs within the i-motif core is mainly governed by the repulsive base stacking interaction. The loop sequence can affect the structural dynamics of the i-motif by altering the loop motion and backbone conformation. Overall this study provides microscopic insight into the i-motif structure that will be helpful to understand the structural aspect of mechanisms of gene regulation by i-motif DNA.
Collapse
Affiliation(s)
- Manas Mondal
- Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Beijing National Laboratory for Molecular Sciences, Peking University, Beijing 100871, China.
| | | | | |
Collapse
|
3
|
Seelam PP, Mitra A, Sharma P. Pairing interactions between nucleobases and ligands in aptamer:ligand complexes of riboswitches: crystal structure analysis, classification, optimal structures, and accurate interaction energies. RNA (NEW YORK, N.Y.) 2019; 25:1274-1290. [PMID: 31315914 PMCID: PMC6800475 DOI: 10.1261/rna.071530.119] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/12/2019] [Accepted: 07/12/2019] [Indexed: 06/10/2023]
Abstract
In the present work, 67 crystal structures of the aptamer domains of RNA riboswitches are chosen for analysis of the structure and strength of hydrogen bonding (pairing) interactions between nucleobases constituting the aptamer binding pockets and the bound ligands. A total of 80 unique base:ligand hydrogen-bonded pairs containing at least two hydrogen bonds were identified through visual inspection. Classification of these contacts in terms of the interacting edge of the aptamer nucleobase revealed that interactions involving the Watson-Crick edge are the most common, followed by the sugar edge of purines and the Hoogsteen edge of uracil. Alternatively, classification in terms of the chemical constitution of the ligand yields five unique classes of base:ligand pairs: base:base, base:amino acid, base:sugar, base:phosphate, and base:other. Further, quantum mechanical (QM) geometry optimizations revealed that 67 out of 80 pairs exhibit stable geometries and optimal deviations from their macromolecular crystal occurrences. This indicates that these contacts are well-defined RNA aptamer:ligand interaction motifs. QM calculated interaction energies of base:ligand pairs reveal a rich hydrogen bonding landscape, ranging from weak interactions (base:other, -3 kcal/mol) to strong (base:phosphate, -48 kcal/mol) contacts. The analysis was further extended to study the biological importance of base:ligand interactions in the binding pocket of the tetrahydrofolate riboswitch and thiamine pyrophosphate riboswitch. Overall, our study helps in understanding the structural and energetic features of base:ligand pairs in riboswitches, which could aid in developing meaningful hypotheses in the context of RNA:ligand recognition. This can, in turn, contribute toward current efforts to develop antimicrobials that target RNAs.
Collapse
Affiliation(s)
- Preethi P Seelam
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| |
Collapse
|
4
|
Halder A, Data D, Seelam PP, Bhattacharyya D, Mitra A. Estimating Strengths of Individual Hydrogen Bonds in RNA Base Pairs: Toward a Consensus between Different Computational Approaches. ACS OMEGA 2019; 4:7354-7368. [PMID: 31459834 PMCID: PMC6648064 DOI: 10.1021/acsomega.8b03689] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/31/2018] [Accepted: 04/12/2019] [Indexed: 06/10/2023]
Abstract
Noncoding RNA molecules are composed of a large variety of noncanonical base pairs that shape up their functionally competent folded structures. Each base pair is composed of at least two interbase hydrogen bonds (H-bonds). It is expected that the characteristic geometry and stability of different noncanonical base pairs are determined collectively by the properties of these interbase H-bonds. We have studied the ground-state electronic properties [using density functional theory (DFT) and DFT-D3-based methods] of all the 118 normal base pairs and 36 modified base pairs, belonging to 12 different geometric families (cis and trans of WW, WH, HH, WS, HS, and SS) that occur in a nonredundant set of high-resolution RNA crystal structures. Having addressed some of the limitations of the earlier approaches, we provide here a comprehensive compilation of the average energies of different types of interbase H-bonds (E HB). We have also characterized each interbase H-bond using 13 different parameters that describe its geometry, charge distribution at its bond critical point (BCP), and n → σ*-type charge transfer from filled π orbitals of the H-bond acceptor to the empty antibonding orbital of the H-bond donor. On the basis of the extent of their linear correlation with the H-bonding energy, we have shortlisted five parameters to model linear equations for predicting E HB values. They are (i) electron density at the BCP: ρ, (ii) its Laplacian: ∇2ρ, (iii) stabilization energy due to n → σ*-type charge transfer: E(2), (iv) donor-hydrogen distance, and (v) hydrogen-acceptor distance. We have performed single variable and multivariable linear regression analysis over the normal base pairs and have modeled sets of linear relationships between these five parameters and E HB. Performance testing of our model over the set of modified base pairs shows promising results, at least for the moderately strong H-bonds.
Collapse
Affiliation(s)
- Antarip Halder
- Center
for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology
(IIIT-H), Gachibowli, Hyderabad 500032, India
| | - Dhruv Data
- Center
for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology
(IIIT-H), Gachibowli, Hyderabad 500032, India
| | - Preethi P. Seelam
- Center
for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology
(IIIT-H), Gachibowli, Hyderabad 500032, India
| | - Dhananjay Bhattacharyya
- Computational
Science Division, Saha Institute of Nuclear
Physics(SINP), 1/AF,
Bidhannagar, Kolkata 700064, India
| | - Abhijit Mitra
- Center
for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology
(IIIT-H), Gachibowli, Hyderabad 500032, India
| |
Collapse
|
5
|
Halder A, Vemuri S, Roy R, Katuri J, Bhattacharyya D, Mitra A. Evidence for Hidden Involvement of N3-Protonated Guanine in RNA Structure and Function. ACS OMEGA 2019; 4:699-709. [PMID: 30775644 PMCID: PMC6372247 DOI: 10.1021/acsomega.8b02908] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2018] [Accepted: 12/25/2018] [Indexed: 05/05/2023]
Abstract
Charged nucleobases have been found to occur in several known RNA molecules and are considered essential for their structure and function. The mechanism of their involvement is however not yet fully understood. Revelation of the role of N7-protonated guanine, in modulating the geometry and stability of noncanonical base pairs formed through its unprotonated edges [Watson-Crick (WC) and sugar], has triggered the need to evaluate the feasibility of similar roles of other protonated nucleobases [Halder et al., Phys Chem Chem Phys, 2015, 17, 26249]. In this context, N3 protonation of guanine makes an interesting case as its influence on the charge distribution of the WC edge is similar to that of N7 protonation, though its thermodynamic cost of protonation is significantly higher. In this work, we have carried out structural bioinformatics analyses and quantum mechanics-based calculations to show that N3 protonation of guanine may take place in a cellular environment, at least in the G:C W:W Trans and G:G W:H Cis base pairs. Our results provide a reasonable starting point for future investigations in order to address the larger mechanistic question.
Collapse
Affiliation(s)
- Antarip Halder
- Center
for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology,
Hyderabad (IIIT-H), Gachibowli, Hyderabad 500032, Telangana, India
| | - Saurabh Vemuri
- Center
for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology,
Hyderabad (IIIT-H), Gachibowli, Hyderabad 500032, Telangana, India
| | - Rohit Roy
- Center
for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology,
Hyderabad (IIIT-H), Gachibowli, Hyderabad 500032, Telangana, India
| | - Jayanth Katuri
- Center
for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology,
Hyderabad (IIIT-H), Gachibowli, Hyderabad 500032, Telangana, India
| | - Dhananjay Bhattacharyya
- Computational
Science Division, Saha Institute of Nuclear
Physics (SINP), 1/AF,
Bidhannagar, Kolkata 700064, India
| | - Abhijit Mitra
- Center
for Computational Natural Sciences and Bioinformatics (CCNSB), International Institute of Information Technology,
Hyderabad (IIIT-H), Gachibowli, Hyderabad 500032, Telangana, India
| |
Collapse
|
6
|
Halder A, Roy R, Bhattacharyya D, Mitra A. Consequences of Mg2+ binding on the geometry and stability of RNA base pairs. Phys Chem Chem Phys 2018; 20:21934-21948. [DOI: 10.1039/c8cp03602k] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Quantum chemical calculations reveal the role of magnesium in stabilizing the geometries of intrinsically unstable RNA base pairs.
Collapse
Affiliation(s)
- Antarip Halder
- Center for Computational Natural Sciences and Bioinformatics (CCNSB)
- International Institute of Information Technology (IIIT-H)
- Hyderabad 500032
- India
| | - Rohit Roy
- Center for Computational Natural Sciences and Bioinformatics (CCNSB)
- International Institute of Information Technology (IIIT-H)
- Hyderabad 500032
- India
| | | | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics (CCNSB)
- International Institute of Information Technology (IIIT-H)
- Hyderabad 500032
- India
| |
Collapse
|
7
|
Effect of single-residue bulges on RNA double-helical structures: crystallographic database analysis and molecular dynamics simulation studies. J Mol Model 2017; 23:311. [DOI: 10.1007/s00894-017-3480-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2017] [Accepted: 09/19/2017] [Indexed: 11/26/2022]
|
8
|
Seelam PP, Sharma P, Mitra A. Structural landscape of base pairs containing post-transcriptional modifications in RNA. RNA (NEW YORK, N.Y.) 2017; 23:847-859. [PMID: 28341704 PMCID: PMC5435857 DOI: 10.1261/rna.060749.117] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/09/2017] [Accepted: 03/23/2017] [Indexed: 05/20/2023]
Abstract
Base pairs involving post-transcriptionally modified nucleobases are believed to play important roles in a wide variety of functional RNAs. Here we present our attempts toward understanding the structural and functional role of naturally occurring modified base pairs using a combination of X-ray crystal structure database analysis, sequence analysis, and advanced quantum chemical methods. Our bioinformatics analysis reveals that despite their presence in all major secondary structural elements, modified base pairs are most prevalent in tRNA crystal structures and most commonly involve guanine or uridine modifications. Further, analysis of tRNA sequences reveals additional examples of modified base pairs at structurally conserved tRNA regions and highlights the conservation patterns of these base pairs in three domains of life. Comparison of structures and binding energies of modified base pairs with their unmodified counterparts, using quantum chemical methods, allowed us to classify the base modifications in terms of the nature of their electronic structure effects on base-pairing. Analysis of specific structural contexts of modified base pairs in RNA crystal structures revealed several interesting scenarios, including those at the tRNA:rRNA interface, antibiotic-binding sites on the ribosome, and the three-way junctions within tRNA. These scenarios, when analyzed in the context of available experimental data, allowed us to correlate the occurrence and strength of modified base pairs with their specific functional roles. Overall, our study highlights the structural importance of modified base pairs in RNA and points toward the need for greater appreciation of the role of modified bases and their interactions, in the context of many biological processes involving RNA.
Collapse
Affiliation(s)
- Preethi P Seelam
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology Hyderabad (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology Hyderabad (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| |
Collapse
|
9
|
RNAHelix: computational modeling of nucleic acid structures with Watson–Crick and non-canonical base pairs. J Comput Aided Mol Des 2017; 31:219-235. [DOI: 10.1007/s10822-016-0007-0] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2016] [Accepted: 12/29/2016] [Indexed: 12/11/2022]
|
10
|
Maiti S, Bhattacharyya D. Stacking interactions involving non-Watson–Crick basepairs: dispersion corrected density functional theory studies. Phys Chem Chem Phys 2017; 19:28718-28730. [DOI: 10.1039/c7cp04904h] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Stacking interactions between a non Watson–Crick G:A S:HT basepair and C:G basepair is predicted in terms of roll, twist and slide basepair step parameters using DFT-D augmented with coarse-grain energy penalty for sugar–phosphate backbone.
Collapse
Affiliation(s)
- Satyabrata Maiti
- Computational Science Division
- Saha Institute of Nuclear Physics
- Kolkata 700064
- India
- Homi Bhaba National Institute
| | - Dhananjay Bhattacharyya
- Computational Science Division
- Saha Institute of Nuclear Physics
- Kolkata 700064
- India
- Homi Bhaba National Institute
| |
Collapse
|
11
|
S. P. P, Sharma P, Mitra A. Higher order structures involving post transcriptionally modified nucleobases in RNA. RSC Adv 2017. [DOI: 10.1039/c7ra05284g] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Quantum chemical studies are carried out to understand the structures and stabilities of higher order structures involving post-transcriptionally modified nucleobases in RNA.
Collapse
Affiliation(s)
- Preethi S. P.
- Center for Computational Natural Sciences and Bioinformatics
- International Institute of Information Technology Hyderabad (IIIT-H)
- Hyderabad
- India
| | - Purshotam Sharma
- Computational Biochemistry Laboratory
- Department of Chemistry and Centre for Advanced Studies in Chemistry
- Panjab University
- Chandigarh
- India
| | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics
- International Institute of Information Technology Hyderabad (IIIT-H)
- Hyderabad
- India
| |
Collapse
|
12
|
Halder A, Bhattacharya S, Datta A, Bhattacharyya D, Mitra A. The role of N7 protonation of guanine in determining the structure, stability and function of RNA base pairs. Phys Chem Chem Phys 2015; 17:26249-63. [DOI: 10.1039/c5cp04894j] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Ab initio computations and bioinformatics studies reveal that stabilization of some important RNA structural motifs might involve N7 protonation of guanine.
Collapse
Affiliation(s)
- Antarip Halder
- Center for Computational Natural Sciences and Bioinformatics (CCNSB)
- International Institute of Information Technology (IIIT-H)
- Hyderabad 500032
- India
| | - Sohini Bhattacharya
- Center for Computational Natural Sciences and Bioinformatics (CCNSB)
- International Institute of Information Technology (IIIT-H)
- Hyderabad 500032
- India
| | - Ayan Datta
- Department of Spectroscopy
- Indian Association for the Cultivation of Science
- Kolkata 700032
- India
| | | | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics (CCNSB)
- International Institute of Information Technology (IIIT-H)
- Hyderabad 500032
- India
| |
Collapse
|