Abstract

Riboswitches are RNA-based genetic control elements that function via a conformational transition mechanism when a specific target molecule binds to its binding pocket. To facilitate an atomic detail interpretation of experimental investigations on the role of the adenine ligand on the conformational properties and kinetics of folding of the add adenine riboswitch, we performed molecular dynamics simulations in both the presence and the absence of the ligand. In the absence of ligand, structural deviations were observed in the J23 junction and the P1 stem. Destabilization of the P1 stem in the absence of ligand involves the loss of direct stabilizing interactions with the ligand, with additional contributions from the J23 junction region. The J23 junction of the riboswitch is found to be more flexible, and the tertiary contacts among the junction regions are altered in the absence of the adenine ligand; results suggest that the adenine ligand associates and dissociates from the riboswitch in the vicinity of J23. Good agreement was obtained with the experimental data with the results indicating dynamic behavior of the adenine ligand on the nanosecond time scale to be associated with the dynamic behavior of hydrogen bonding with the riboswitch. Results also predict that direct interactions of the adenine ligand with U74 of the riboswitch are not essential for stable binding although it is crucial for its recognition. The possibility of methodological artifacts and force-field inaccuracies impacting the present observations was checked by additional molecular dynamics simulations in the presence of 2,6-diaminopurine and in the crystal environment. Copyright (c) 2010 Elsevier Ltd. All rights reserved.

Abstract

Sac7d belongs to a family of chromosomal proteins, which are crucial for thermal stabilization of DNA at higher growth temperatures. It is capable of binding DNA nonspecifically, and is responsible for the increase in the melting temperature of DNA in the bound form up to 85 °C. Molecular dynamics (MD) simulations were performed at different temperatures on two protein-DNA complexes of Sac7d. Various structural and energetic parameters were calculated to examine the DNA stability and to investigate the conformational changes in DNA and the protein-DNA interactions. Room temperature simulations indicated very good agreement with the experimental structures. The protein structure is nearly unchanged at both 300 and 360 K, and only up to five base pairs of the DNA are stabilized by Sac7d at 360 K. However, the MD simulations on DNA alone systems show that they lose their helical structures at 360 K further supporting the role of Sac7d in stabilizing the oligomers. At higher temperatures (420 and 480 K), DNA undergoes denaturation in the presence and the absence of the protein. The DNA molecules were found to undergo B- to A-form transitions consistent with experimental studies, and the extent of these transitions are examined in detail. The extent of sampling B- and A-form regions was found to show temperature and sequence dependence. Multiple MD simulations yielded similar results validating the proposed model. Interaction energy calculations corresponding to protein-DNA binding indicates major contribution due to DNA backbone, explaining the nonspecific interactions of Sac7d.

Abstract

One of the most efficient chemical transformations in the natural product synthesis of steroid like molecules is the synthetic methodology that combines the pericyclic and intramolecular reactions termed as transannular Diels-Alder reaction (TADA). TADA reactions occur in (x+y+2)-membered triene macrocycles, which contain both the diene and the dienophile, to form A.B.C [x.6.y] type tricyclic compounds. Initially, a validation study was carried out to identify an appropriate method for studying this reaction pathway. Several methods including HF, post-HF methods such as Mpx and CC, and density functional B3LYP level were considered, and we found that B3LYP method used along with a reasonably smaller basis set (6-31G(d)) is adequate enough. We have extensively studied the TADA reactions of six possible cis/trans-isomers of 14-membered macrocycles, and correlated their reactivities with experimental data wherever applicable. We also have studied the effect of ring size on the TADA reactivity by performing calculations on 552, 562, 572, 652, 672, 752, 762 and 772 systems using the B3LYP level employing the 6-31G(d) and the cc-pVTZ basis sets. We have shown that the activation energy and the reaction energy decrease with respect to the increase in the ring size. In addition to the activation energies, the conformational preferences of the reactants were suggestd to be crucial for the reaction to take place, and on its stereochemical outcome. We have studied their conformational properties using replica exchange molecular dynamics simulations. Initially, we obtained force field parameters for the macrocyclic trienes based on MP2 calculations on six model systems. The analysis concentrated on the dihedral angle corresponding to the central bond of the diene moiety, and the distance between the reaction centers. The extent of conformational sampling of each of the macrocycles is compared with the reactivities. The presentation will discuss the TADA reactivity in terms of activation energies corresponding to the transition states of the reactions, and also based on the conformational properties of the reactants.

Abstract

Hyperthermophilic proteins are naturally occurring proteins that are capable of retaining their structure, and performing their normal functions at temperatures above 80º C. Identification of the structural determinants that are crucial for such a remarkable stability is of fundamental importance, and has potential applications in biocatalysis. Molecular dynamics simulations on two such proteins (Sac7d and Sso7d) illustrated the importance of ion pair interactions for achieving thermostability.1 Previous studies have indicated that the formation of ion pair interactions does not contribute to stability due to high desolvation penalty, while others have argued in favor of stabilization. We have performed umbrella sampling simulations on model systems at different temperatures, which showed that with increasing temperature, the stability of the ion pair interactions become highly favorable. Thus, protein folding involving association of an ion pair may lead to stabilization at enhanced temperature close to the growth temperature of the hyperthemophiles. We have also investigated the effect of protein-DNA interactions, and their role in DNA stabilization at these conditions.2 Canonical B- to A-form transitions were observed in the protein-DNA complexes, and the protein-DNA interaction energies are used to explain the sequence nonspecific nature of the binding activity. The crucial nonbonded interactions that facilitate extraordinary thermal stability of proteins and DNAs in thermophiles will be discussed in this presentation.

Abstract

The enlarging HIV/AIDS pandemic and the virus’ ability to develop resistance to anti-retroviral therapy requires the development of new drugs based on novel targets. Viruses interact extensively with the host at the cellular and molecular levels, which is necessary for the expression of viral proteins and replication of viral genomes. At another level, these interactions also help the virus evade host immune responses. Besides the prototypic retroviral proteins (Gag, Pol and Env), HIV-1 also expresses two regulatory proteins (Tat and Rev) and four accessory proteins (Nef, Vif, Vpr and Vpu). The accessory proteins are not required for viral replication in vitro, but are critical for pathogenesis and disease development in the host. We have focused our efforts on two of these – Nef and Vpu, with an aim to understand their interactions with host proteins and the functional consequences of these interactions for HIV-1 infection and pathogenesis. The Nef protein binds to multiple host proteins to optimize the host cell environment and to enable infected cells to evade the host immune response. A major property of Nef is to downregulate surface expression of the CD4, MHC I and MHC II proteins, by increasing the endocytosis of these molecules. We also discovered that Nef binds to the cytoplasmic domains of the B7 family of costimulatory proteins – CD80 and CD86; this results in their removal from the cell surface and sequestration to the Golgi compartment (1-3). We have now developed an in vitro screen based on purified recombinant Nef proteins and a 20-mer synthetic peptide from the cytoplasmic domain of CD80 and CD86. A chemical library comprising of 1064 compounds was screened and ~7% compounds were obtained as hits, which inhibited the Nef-CD80/CD86 interaction. These compounds are now being evaluated in a cellbased assay and an independent dose dependent response assay. The top hits will be used in a combinatorial chemical approach to synthesize a large array of new molecules for further screening using the in vitro and cell-based screens. The HIV-1 Vpu protein is required for the release of new virions from the surface of infected cells; for this the transmembrane domain of Vpu is critical. Using the Topology Data Bank of Transmembrane Proteins (TOPDB) database we extracted 13 alpha-helical sequences that closely resemble the Vpu-TM domain. We have taken the top two hits and have constructed recombinant Vpu proteins in which the native TM domain was replaced by the identified sequences. These recombinant proteins were successfully expressed in HEK293 cells, and are now being tested for their effects on virion release as well as their ability to inhibit Vpu activity. We describe here two paradigms to screen for anti-HIV activity based on proteinprotein interactions. Results will be presented for both of these.

Abstract

Thermal, chemical and force induced unfolding are three primary experimental approaches that have been successfully used to study protein folding and associated pathways. Recent experimental studies have demonstrated the importance of the use of combination of these techniques to study various aspects of protein folding. We have addressed this issue using two different model systems. The effect of temperature on the mechanical unfolding of titin I27 immunoglobulin domain has been investigated using steered molecular dynamics (MD) simulations performed both at constant velocity and constant force. The maximum force of unfolding decreases with respect to the increase in the temperature, but the extension corresponding to this force remains similar. However, we observed that the maximum probable unfolding force steadily decreases with an increase in temperature (see below). Several computations done on each of the temperatures consistently reproduce these effects. The results will be discussed in light of the available experimental data, and molecular level details on the conformational changes that effect such temperature dependence will be presented. In another study, we have done rigorous free energy calculations on the unfolding of Trp-mini cage protein using umbrella sampling MD simulations. Calculations were done with respect to increase in the temperature, and with respect to the concentration of the urea solution. The change in the free energy profiles, and the effect of external perturbations, and the atomistic details of the unfolding will be discussed. Temperature

Abstract

Design of biomimetic systems is an alternative method of finding biologically active molecules. Peptides incorporating tetrahydrofuran amino acids (TAA) adopt secondary folding patterns and thus mimic natural peptides. The preferred conformation of TAA derived linear and cyclized peptides depends on the stereochemistry of the tetrahydrofuran (THF) ring. Control of stereochemistry of THF ring in the structure of the linear tripeptide Boc-TAA-Leu-Val- OMe containing (2R,5S)-cis TAA, (2S,5R)-cis TAA, (2R,5R)-trans TAA and (2S,5S)-trans TAA is investigated using structural parameters and energetics data obtained from density functional theory (DFT) based calculations at B3LYP/6-31G(d,p) level of theory. We found that the conformations associated with -, -turn structures are stabilized by intramolecular hydrogen bonding interactions. Kinetic control of reactions leading to intra and inter molecular cyclization of tripeptide TAA-Gly-Gly containing (2S,5R)-cis TAA and (2S,5S)-trans TAA is studied at the same level of theory in gas phase.[1] We observe that kinetic control favors cyclodimerization instead of intramolecular cyclization.

Abstract

Macrocyclic trienes are important precursors for the synthesis of several biologically relevant molecules such as steroids. These molecules undergo (4+2) cycloaddition to form tricyclic compounds, which has four newly formed stereocenters. The stereochemistry at these four carbon centers has been shown to be largely controlled by the reaction energy barriers, and the conformational preferences of the macrocycle. Molecular dynamics (MD) simulation is a highly efficient approach for studying conformational properties of chemical molecules and biological macromolecules. MP2/6-31G* method was used to obtain potential energy surfaces corresponding to appropriate model compounds and used as the target data. The force field parameters were refined to accurately represent the target data involving all the model compounds chosen. Using the refined molecular mechanics force field, MD simulations were performed on six possible 14-membered symmetric trienes at different temperatures. We have also used replica exchange MD (RexMD) simulations with eight replicas between 300 and 370 K. The molecules were found to sample a larger conformational space in the RexMD simulation compared to the regular MD simulations. The probability distributions corresponding to the dihedral angles, and distances between the reaction centers have been analyzed. The relationship between the conformational sampling of the macrocyclic trienes vis-à-vis their reactivity will be discussed. Additionally, clustering of the structures from the MD trajectories was done to identify all unique conformations using a new approach.

Abstract

Sac7d is a protein present in Sulfolobus acidocaldarius, a hyperthermophilic archaeon, which is stable at extreme conditions. It is a chromosomal protein, which is capable of binding DNA thereby increasing its melting temperature. It was suggested to stabilize central 3 to 5 base pair steps in DNA. It binds to the minor groove of the DNA causing a kink in the helix without any sequential preference. We have used molecular dynamics (MD) simulations to study the binary complexes of Sac7d bound to DNA at 300 and 360 K. The objective is to understand the atomic detailed mechanism of the capability of Sac7d in facilitating DNA stability at higher temperature. As has been shown before, the structure of the protein is highly stable both at 300 and 360 K. The DNA duplex was found to be stable at 300 K and is similar to the X-ray crystal structure. Consistent with experiments, at least the central four base pairs of the DNA were found to be base paired during most of the simulation at 360 K. The dynamic properties of the protein in the protein-DNA complexes remained unchanged compared to when the protein is unbound. Interestingly, B- to A-form transitions occur in both the DNA molecules when bound to Sac7d, which was confirmed by sugar puckering angles, and backbone dihedral angles. Detailed examination of the protein-DNA interactions reveals that the stabilization of the DNA molecules is mostly due to the favorable interaction of the protein with the backbone of the DNA.

Abstract

Base flipping is a process by which one of the bases of the DNA swings out from its helical position to an extrahelical form. Several proteins such as methyltransferases, glycosylases and endonucleases employ this phenomenon in order to gain access to the target base for chemical modification. Even in the absence of protein, base pair opening occurs in nucleic acids in solution. X-ray crystal studies, fluorescence spectroscopy and NMR imino proton exchange experiments have been extensively employed to study the structures and dynamics involved in this fascinating event. We have employed molecular dynamics simulations to study base flipping both in presence and absence of the protein. The free energy pathways obtained using potential of mean force calculations will be presented. The base pair opening dynamics of AT and GC base pairs was investigated by calculating the equilibrium constants corresponding to the equilibrium between the base pair open and closed states, and compared with experimental data. Our calculations showed that the contribution to the overall equilibrium constants of the base pair opening is primarily from the purine bases. The methodological issues such as sufficient sampling and force field parameters will be addressed. In another study, the role of THR250 HhaI methyltransferase in base flipping will be discussed. THR250 is a highly conserved residue in methyltransferases and was proposed to be crucial for the correct placement of flipped cytosine for methylation. Previous computational studies have shown that the protein actively facilitates the base flipping process. Free energy calculations of base flipping in the presence of the wild type and T250G, T250A and T250D mutants of the protein were computed. The free energy profiles, structural changes and the energetic parameters for the base flipping process will be presented.