Abstract

The effect of the linker length (–CH2–) connecting the diene and the dienophile in macrocyclic trienes on their transannular Diels–Alder (TADA) reactivity has been investigated using density functional theory (DFT) calculations. The relationship between the conformational properties of these reactants and their reaction energy barriers was examined and a quantitative relationship has been obtained. The transition state energy barriers were found to increase with an increase in the linker length, which is in contrast to the expected trend. The conformational preferences of the triene reactants were studied using replica exchange molecular dynamics (REMD) simulations. These calculations reveal that longer linkers lead to a decreased probability of the occurrence of conformations with the diene and dienophile parts of the system vicinal to each other, and thus lower reactivity of systems with long linkers. Excellent correlations between the transition state energies and the cumulative probabilities corresponding to a short distance between the reactive sites, along with the ability of the diene to exist in the s-cisoid form, were observed.

Abstract

Archaeal chromatin proteins share molecular and functional similarities with both bacterial and eukaryotic chromatin proteins. These proteins play an important role in functionally organizing the genomic DNA into a compact nucleoid. Cren7 and Sul7 are two crenarchaeal nucleoidassociated proteins, which are structurally homologous, but not conserved at the sequence level. Cocrystal structures have shown that these two proteins induce a sharp bend on binding to DNA. In this study, we have investigated the architectural properties of these proteins using atomic force microscopy, molecular dynamics simulations and magnetic tweezers. We demonstrate that Cren7 and Sul7 both compact DNA molecules to a similar extent. Using a theoretical model, we quantify the number of individual proteins bound to the DNA as a function of protein concentration and show that forces up to 3.5 pN do not affect this binding. Moreover, we investigate the flexibility of the bending angle induced by Cren7 and Sul7 and show that the protein–DNA complexes differ in flexibility from analogous bacterial and eukaryotic DNA-bending proteins.

Abstract

Vpu is an 81-residue protein encoded by human immunodeficiency virus type I (HIV-1) that facilitates viral release from host cells. The protein has a cytoplasmic domain and a helical transmembrane (TM) domain, of which the latter oligomerizes to form cation-specific ion channels. The number of TM domains that constitute the channel is still unclear, with experimental studies indicating the existence of a variety of oligomeric states. In this study, we have examined the possibility of Vpu to exist in tetra-, penta-, and hexameric states using comprehensive molecular dynamics (MD) simulations. By modeling the TM domain as an ideal α-helix, we carried out replica-exchange MD simulations in an implicit membrane environment for obtaining suitable starting structures, which were then subjected to extensive MD simulations in a fully hydrated lipid bilayer environment. The results show that the pentameric form is the most stable oligomeric state (the tetramer and hexamer models lose their initial channel-like structure), with helices in the pentamer being held together by strong van der Waals interactions. Hydrogen bonds between lipid headgroups and basic/hydrophilic residues on the protein are stronger in the pentamer than in the tetramer or the hexamer, indicating that these interactions might play a role in adhering the pentamer to the membrane. Free energy calculations using umbrella sampling technique have been performed to examine the potassium ion permeation through the pentameric pore. The results show a high free energy barrier corresponding to ion transport indicating weak ion channel activity in agreement with previous experimental biophysical studies.

Abstract

Chemical modifications on the backbone of the oligonucleotides have been shown to improve several factors such as increased cellular uptake, altered duplex conformation and stability. Modified oligonucleotides have also been shown to be potential antisense agents. Some of the common modifications include 2'-O-methyl, 2'-O,4'-C-methylene--D-ribofuranosyl (LNA) and phosphoroamidate. In general, the modifications at ribose sugar of the oligonucleotide attribute special features to the duplexes than other modifications. Molecular dynamics (MD) simulations are valuable tools to study their dynamics, conformational preferences, stabilities etc. Here, we have employed different simulation techniques such as classical MD, replica exchange MD to understand the impact of modification on the duplex properties. The simulations revealed several important aspects related to duplex stability and relation with their base content. Our results also suggest that the differential changes observed in nucleic acids depends on the position, nature of the modified site. This study also provides a molecular level picture for the origin of structural changes arise due to different types of modifications which is very important in designing the most fruitful antisense oligonucleotides.

Abstract

SAM-III riboswitch present in the 5’-UTR of mRNA that translates to SAM synthetase in certain bacteria is involved in regulating the biosynthesis of S-adenosylmethionine (SAM) by acting as an on-off switch for the translation process. We have used molecular dynamics (MD) and steered MD simulations to understand the ligand recognition and the switching mechanisms. Computations on model systems indicated that the riboswitch does not recognize SAM over a structurally related analog, S-adenosylhomocysteine (SAH). Nonetheless, the RNA uses nonspecific interactions to stabilize SAM in its binding pocket compared to SAH. Calculations were done in the absence of the ligand to examine the conformational changes responsible for the transition between on and off states. While the on- state is capable of sampling a large conformational space, formation of duplexlike structure comprising the adenine part of the ligand and sequentially distant nucleotides seems to be crucial for the on- to off- state transition.

Abstract

t has been proposed that combining chemical and thermal denaturation along with mechanical unfolding enables better understanding of the thermodynamics of protein folding, and associated pathways. We have used molecular dynamics (MD) simulations, and umbrella sampling technique to study the effect of thermal and chemical perturbations on the unfolding of Trp-cage miniprotein (TC5b). Free energy profiles for the unfolding of the minicage protein were obtained at two different temperatures, and in presence and absence of urea. The distance between the penultimate terminal residues was used as the reaction coordinate for calculating the free energies. The unfolding involves opening of the Trp-cage, and separation of the secondary structure elements followed by the unfolding of the α-helix. Several energetic and geometric analyses including inter- and intramolecular interactions, and radial distribution functions and their relevance to the folding pathways will be discussed. Similar study on the titin I27 protein will also be presented.

Abstract

: Density functional B3LYP method was used to investigate the preference of intra and intermolecular cyclizations of linear tripeptides containing tetrahydrofuran amino acids. Two distinct model pathways were conceived for the cyclization reaction and all possible transition states and intermediates were located. Analysis of the energetics indicate intermolecular cyclization being favored by both thermodynamic and kinetic control. Geometric and NBO analyses were performed to explain the trends obtained along both the reaction pathways. Conceptual density functional theory based reactive indices also show that reaction pathways leading to intermolecular cyclization of the tripeptides are relatively more facile compared to intramolecular cyclization.

Abstract

Understanding protein folding thermodynamics and kinetics is a central issue in molecular biology. Molecular dynamics (MD) simulations are being extensively used for this purpose. However, direct comparison between simulations and experiments requires both, an accurate description of the system and extensive sampling of the conformational space. Variants of MD have been developed to overcome sampling issues. Here, we use umbrella sampling technique to quantitatively compute the differences in thermodynamic quantities of the unfolding process of Trp-cage miniprotein TC5b (PDB ID: 1L2Y) due to changes in temperature as well as the presence/absence of chemical denaturants. Urea at a concentration of 8M has been used as the chemical denaturant. WHAM was used to obtain the unbiased free energy profiles. The free energy profile can be categorised into three parts – the opening of the Trp-cage, separation of secondary structure elements and distortion of alpha-helix. An increase in temperature or presence of urea effects only the opening of the Trp-cage while the separation of the secondary structure elements and distortion of alpha-helix is independent of temperature and urea. Contour plots were constructed to substantiate the results obtained from free energy profiles.

Abstract

Ribonucleases act as biological counter weights by playing a major role in gene expression. Ribonuclease H is one of them which performs the endonucleotic cleaveage reaction of the 3'-O-P linkage of RNA strand of DNA/RNA hybrids. These enzymes show specificity towards DNA/RNA hybrids but not pure structures. One of the reason for this differentiation could be mixed A/Bconformation of hybrid compared to pure structures and there might be other reasons. There is a debate concerning the substrate recognition mechanism whether this is because of induced fit type mechanism or conformational rearrangement of monomers in order to come close or may be non bonded interactions which will put the enzyme and substrate together in order to perform the hydrolysis reaction. While the role of metal ions in the reaction has been recognised, the mechanism is not obvious. For differentating the structural changes, we performed long molecular dynamics (MD) simulations on ribonuclease H enzyme of B.halodurans with its substrate and their monomers. The calculations show that the large conformational changes in both the monomers favours the reaction. Since the protein hydrophobic core is very stable, these changes will take place at surface level. The principal component analysis also supports the large conformational changes in protein and hybrid. Apart form these changes, the protein binding increases the plasticity of the hybrids. Overall the calculations show that structural deformations in both hybrid and enzyme brings them close to each other. Free energy calculations have been performed to study the proteinhybrid interactions

Abstract

Transannular Diels Alder reaction (TADA), a combination of both intramolecular and Diels-Alder reactions, is a powerful synthetic organic reaction mechanism used in synthesizing polycyclic compounds with high degree of chemo-, regio- and stereospecificity. TADA reactions occur in (x+y+2)-membered triene macrocycles, which contain both the diene and the dienophile moieties, to form A.B.C[x+6+y] type tricyclic compounds. The TADA reactivity of a series of 14-membered macrocyclic rings have been studied using the density functional B3LYP level of theory. The reactants and the transition states are capable of existing in different conformations. The key conformers of the transition states were studied at the B3LYP level. We have done force field parameters for the trienes, and performed detailed replica exchange molecular dynamics simulations to identify all possible isomers of the reactants. We established a linear correlation between the activation energies and the extent of sampling of the reactants in certain conformational states. In this presentation, we also will present the effect of the lengths of the linkers (3 to 5 methylene units) connecting the diene and dienophile units on the TADA reactivities.