Sitansh Sharma

Binding of Gold Nanoclusters with Size-Expanded DNA Bases: A Computational Study of Structural and Electronic Properties

Binding of gold nanoclusters with size-expanded DNA bases, xA, xC, xG, and xT, is studied using quantum chemical methods. Geometries of the neutral xA-Au6, xC-Au6, xG-Au6, and xT-Au6 complexes were fully optimized using the B3LYP density functional method (DFT). The gold clusters around xA and xT adopt triangular geometries, whereas irregular structures are obtained in the case of gold clusters complexed around xC and xG. The lengths of the bonds between atoms in the x-bases increase on gold complexation. The aromatic character of the x-bases also increases on gold complexation except for the five-member rings. A significant charge transfer from the x-base to gold atoms is seen in these complexes. Second-order interactions are observed in addition to direct covalent bonds between gold atoms and x-bases.

Quantum Chemical Studies of Structures and Binding in Noncanonical RNA Base pairs: The Trans Watson-Crick:Watson-Crick Family

Abstract The trans Watson-Crick/Watson-Crick family of base pairs represent a geometric class that play important structural and possible functional roles in the ribosome, tRNA, and other functional RNA molecules. They nucleate base triplets and quartets, participate as loop closing terminal base pairs in hair pin motifs and are also responsible for several tertiary interactions that enable sequentially distant regions to interact with each other in RNA molecules. Eleven representative examples spanning nine systems belonging to this geometric family of RNA base pairs, having widely different occurrence statistics in the PDB database, were studied at the HF/6–31G (d, p) level using Morokuma decomposition, Atoms in Molecules as well as Natural Bond Orbital methods in the optimized gas phase geometries and in their crystal structure geometries, respectively. The BSSE and deformation energy corrected interaction energy values for the optimized geometries are compared with the corresponding values in the crystal geometries of the base pairs. For non protonated base pairs in their optimized geometry, these values ranged from −8.19 kcal/mol to −21.84 kcal/mol and compared favorably with those of canonical base pairs. The interaction energies of these base pairs, in their respective crystal geometries, were, however, lesser to varying extents and in one case, that of A:A W:W trans, it was actually found to be positive. The variation in RMSD between the two geometries was also large and ranged from 0.32–2.19 Å. Our analysis shows that the hydrogen bonding characteristics and interaction energies obtained, correlated with the nature and type of hydrogen bonds between base pairs; but the occurrence frequencies, interaction energies, and geometric variabilities were conspicuous by the absence of any apparent correlation. Instead, the nature of local interaction energy hyperspace of different base pairs as inferred from the degree of their respective geometric variability could be correlated with the identities of free and bound hydrogen bond donor/acceptor groups present in interacting bases in conjunction with their tertiary and neighboring group interaction potentials in the global context. It also suggests that the concept of isostericity alone may not always determine covariation potentials for base pairs, particularly for those which may be important for RNA dynamics. These considerations are more important than the absolute values of the interaction energies in their respective optimized geometries in rationalizing their occurrences in functional RNAs. They highlight the importance of revising some of the existing DNA based structure analysis approaches and may have significant implications for RNA structure and dynamics, especially in the context of structure prediction algorithms.

On the role of the cis Hoogsteen:sugar-edge family of base pairs in platforms and triplets-quantum chemical insights into RNA structural biology.

Base pairs belonging to the cis Hoogsteen:sugar-edge (H:S) family play important structural roles in folded RNA molecules. Several of these are present in internal loops, where they are involved in interactions leading to planar dinucleotide platforms which stabilize higher order structures such as base triplets and quartets. We report results of analysis of 30 representative examples spanning 16 possible base pair combinations, with several of them showing multimodality of base pairing geometry. The geometries of 23 of these base pairs were modeled directly from coordinates extracted from RNA crystal structures. The other seven were predicted structures which were modeled on the basis of observed isosteric analogues. After appropriate satisfaction of residual valencies, these structures were relaxed using the B3LYP/6-31G(d,p) method and interaction energies were derived at the RIMP2/aug-cc-pVDZ level of theory. The geometries for each of the studied base pairs have been characterized in terms of the number and nature of H-bonds, rmsd values observed on optimization, base pair geometrical parameters, and sugar pucker analysis. In addition to its evaluation, the nature of intermolecular interaction in these complexes was also analyzed using Morokuma decomposition. The gas phase interaction energies range between -5.2 and -20.6 kcal/mol and, in contrast to the H:S trans base pairs, show enhanced relative importance of the electron correlation component, indicative of the greater role of dispersion energy in stabilization of these base pairs. The rich variety of hydrogen bonding pattern, involving the flexible sugar edge, appears to hold the key to several features of structural motifs, such as planarity and propensity to participate in triplets, observed in this family of base pairs. This work explores these aspects by integrating database analysis, and detailed base pairing geometry analysis at the atomistic level, with ab initio computation of interaction energies. The study, involving alternative classification of base pairs and triplets, provides insights into intrinsic properties of these base pairs and their possible structural and functional roles.

Genetic algorithm optimization of laser pulses for molecular quantum state excitation

Conventionally optimal control theory has been used in the theoretical design of laser pulses through the direct variation in the electric field of the laser pulse as a function of time. This often leads to designed laser pulses which contain a broad and seemingly arbitrary frequency structure that varies in time in a manner which may be difficult to realize experimentally. In contrast, the experimental design of laser pulses has used a genetic algorithm (GA) approach, varying only those laser parameters actually available to the experimentalist. We investigate in this paper the possibility of using GA optimization methods in the theoretical design of laser pulses to bring about quantum state transitions in molecules. This allows us to select only a small limited number of parameters to vary and to choose these parameters so that they correspond to those available to the experimentalist. In the paper we apply our methods to the vibrational-rotational excitation of the HF molecule. We choose a small limited number of frequencies and vary only the associated electric field amplitudes and pulse envelopes. We show that laser pulses designed in this way can lead to very high transition probabilities.

Modeling the noncovalent interactions at the metabolite binding site in purine riboswitches

We present gas phase quantum chemical studies on the metabolite binding interactions in two important purine riboswitches, the adenine and guanine riboswitches, at the B3LYP/6-31G(d,p) level of theory. In order to gain insights into the strucutral basis of their discriminative abilities of regulating gene expression, the structural properties and binding energies for the gas phase optimized geometries of the metabolite bound binding pocket are analyzed and compared with their respective crystal geometries. Kitaura-Morokuma analysis has been carried out to calculate and decompose the interaction energy into various components. NBO and AIM analysis has been carried out to understand the strength and nature of binding of the individual aptamer bases with their respective purine metabolites. The Y74 base, U in case of adenine riboswitch and C in case of guanine riboswitch constitutes the only differentiating element between the two binding pockets. As expected, with W:W cis G:C74 interaction contributing more than 50% of the total binding energy, the interaction energy for metabolite binding as calculated for guanine (-46.43 Kcal/mol) is nearly double compared to the corresponding value for that of adenine (-24.73 Kcal/mol) in the crystal context. Variations in the optimized geometries for different models and comparison of relative contribution to metabolite binding involving four conserved bases reveal the possible role of U47:U51 W:H trans pair in the conformational transition of the riboswitch from the metabolite free to metabolite bound state. Our results are also indicative of significant contributions from stacking and magnesium ion interactions toward cooperativity effects in metabolite recognition.

On the role of Hoogsteen:Hoogsteen interactions in RNA: Ab initio investigations of structures and energies

We use a combination of database analysis and quantum chemical studies to investigate the role of cis and trans Hoogsteen:Hoogsteen (H:H) base pairs and associated higher-order structures in RNA. We add three new examples to the list of previously identified base-pair combinations belonging to these families and, in addition to contextual classification and characterization of their structural and energetic features, we compare their interbase interaction energies and propensities toward participation in triplets and quartets. We find that some base pairs, which are nonplanar in their isolated minimum energy geometries, attain planarity and stability upon triplet formation. A:A H:H trans is the most frequent H:H combination in RNA structures. This base pair occurs at many distinct positions in known rRNA structures, where it helps in the interaction of ribosomal domains in the 50S subunit. It is also present as a part of tertiary interaction in tRNA structures. Although quantum chemical studies suggest an intrinsically nonplanar geometry for this base pair in isolated form, it has the tendency to attain planar geometry in RNA crystal structures by forming higher-order tertiary interactions or in the presence of additional base-phosphate interactions. The tendency of this base pair to form such additional interactions may be helpful in bringing together different segments of RNA, thus making it suitable for the role of facilitator for RNA folding. This also explains the high occurrence frequency of this base pair among all H:H interactions.

A Theoretical Study on Interaction of Small Gold Clusters Aun (n = 4, 6, 8) with xDNA Base Pairs

Abstract xDNA constitutes a novel class of size expanded synthetic nucleic acids in which one of the bases of the base pairs is larger than the natural DNA bases. These expanded bases are called x-bases. In this paper, we investigate the hydrogen bonding characteristics and relevant molecular properties of model complexes (xA…T)-Aun, (xT…A)-Aun, (xG…C)-Aun, and (xC…G)-Aun (n = 4, 6, 8) consisting of xDNA base pairs and gold clusters, in order to study the nature of gold-xDNA binding. We offer detailed characterization of their different aspects, viz., structural, electronic and spectroscopic, effect of gold cluster size, aromaticity, and planarity using quantum mechanics based density functional theory (DFT). Significant charge transfer is seen between the gold clusters and x-base pairs. Gold complexation is found to affect the interbase hydrogen bonding in these complexes. In addition to anchor bonds, X-H…Au type of hydrogen bonding interactions are also found to contribute to the gold-base pair binding in these complexes.

Design of an infrared laser pulse to control the multiphoton dissociation of the Fe-CO bond in CO-heme compounds

Optimal control theory is used to design a laser pulse for the multiphoton dissociation of the Fe-CO bond in the CO-heme compounds. The study uses a hexacoordinated iron-porphyrin-imidazole-CO complex in its ground electronic state as a model for CO liganded to the heme group. The potential energy and dipole moment surfaces for the interaction of the CO ligand with the heme group are calculated using density functional theory. Optimal control theory, combined with a time-dependent quantum dynamical treatment of the laser-molecule interaction, is then used to design a laser pulse capable of efficiently dissociating the CO-heme complex model. The genetic algorithm method is used within the mathematical framework of optimal control theory to perform the optimization process. This method provides good control over the parameters of the laser pulse, allowing optimized pulses with simple time and frequency structures to be designed. The dependence of photodissociation yield on the choice of initial vibrational state and of initial laser field parameters is also investigated. The current work uses a reduced dimensionality model in which only the Fe-C and C-O stretching coordinates are explicitly taken into account in the time-dependent quantum dynamical calculations. The limitations arising from this are discussed in Sec. IV.

Design of Optimal Laser Fields to Control Vibrational Excitations in Carboxy-myoglobin

Optimal control theory is applied to obtain infrared laser pulses for the selective vibrational excitation of a two mathematical dimensional model of carboxy-myoglobin. Density functional theory is used to obtain the potential energy and dipole moment surfaces of the active site model. The Conjugate gradient method is employed to optimize the cost functional and to obtain the optimized laser pulses. Optimized laser fields are found which give virtually 100% excitation probability to preselected vibrational levels.

Design of laser pulses for selective vibrational excitation of the N6-H bond of adenine and adenine-thymine base pair using optimal control theory

Time dependent quantum dynamics and optimal control theory are used for selective vibrational excitation of the N6-H (amino N-H) bond in free adenine and in the adenine-thymine (A-T) base pair. For the N6-H bond in free adenine we have used a one dimensional model while for the hydrogen bond, N6-H(A)...O4(T), present in the A-T base pair, a two mathematical dimensional model is employed. The conjugate gradient method is used for the optimization of the field dependent cost functional. Optimal laser fields are obtained for selective population transfer in both the model systems, which give virtually 100% excitation probability to preselected vibrational levels. The effect of the optimized laser field on the other hydrogen bond, N1(A)...H-N3(T), present in A-T base pair is also investigated.