Search

RNA interference, particularly siRNA induced gene silencing is becoming an important avenue of modern therapeutics. The siRNA is delivered to the cells as short double helical RNA which becomes single stranded for forming the RISC complex. Significant experimental evidence is available for most of the steps except the process of the separation of the two strands. We have attempted to understand the pathway for double stranded siRNA (dsRNA) to single stranded (ssRNA) molecules using steered molecular dynamics simulations. As the process is completely unexplored we have applied force from all possible directions restraining all possible residues to convert dsRNA to ssRNA. We found pulling one strand along helical axis direction restraining far end of the other strand demands excessive force for ssRNA formation. Pulling a central residue of one strand, in a direction perpendicular to the helix axis, while keeping the base paired residue fixed requires intermediate force for strand separation. Moreover we found that in this process the force requirement is quite high for the first bubble formation (nucleation energy) and the bubble propagation energies are quite small. We hypothesize this could be the mode of action adopted by the proteins in the cells., Comment: 25 pages, 7 figures

Cytosolic delivery of functional siRNA remains the major challenge to develop siRNA-based therapeutics. Designing clinically safe and effective siRNA transporter to deliver functional siRNA across the plasma and endosomal membrane remains a key hurdle. With the aim of improving endosomal release, we have designed cyclic and linear peptide-based transporters having an Arg-DHis-Arg template. Computational studies show that the Arg-DHis-Arg template is also stabilized by the Arg-His side-chain hydrogen bonding interaction at physiological pH, which dissociates at lower pH. The overall atomistic interactions were examined by molecular dynamics simulations, which indicate that the extent of peptide_siRNA assembly formation depends greatly on physicochemical properties of the peptides. Our designed peptides having the Arg-DHis-Arg template and two lipidic moieties facilitate high yield of intracellular delivery of siRNA. Additionally, unsaturated lipid, linoleic acid moieties were introduced to promote fusogenicity and facilitate endosomal release and cytosolic delivery. Interestingly, such protease-resistant peptides provide serum stability to siRNA and exhibit high efficacy of erk1 and erk2 gene silencing in the triple negative breast cancer (TNBC) cell line. The peptide having two linoleyl moieties demonstrated comparable efficacy with commercial transfection reagent HiPerFect, as evidenced by the erk1 and erk2 gene knockdown experiment. Additionally, our study shows that ERK1/2 silencing siRNA and doxorubicin-loaded gramicidin-mediated combination therapy is more effective than siRNA-mediated gene silencing-based monotherapy for TNBC treatment.

Promoter binding specificity is one of the important characteristics of transcription by Mycobacterium tuberculosis (Mtb) sigma (σ) factors, which remains unexplored due to limited structural evidence. Our previous study on the structural features of Mtb-SigH, consisting of three alpha helices, and its interaction with core RNA polymerase has been extended herein to determine the little known DNA sequence recognition pattern involving its cognate promoters. Herein, high resolution X-ray crystallographic structures of the protein-DNA complexes were inspected to determine the tentative DNA-binding helix of the σ factor. The binding interface in the available crystal structures is found to be populated mainly with specific residues such as Arg, Asn, Lys, Gln, and Ser. We uncovered the helix 3 of Mtb-SigH containing most of these amino acids, which ranged from Arg 64 to Arg 75, forming the predicted active site. The complex of Mtb-SigH:DNA is modelled with 20 promoter sequences. The binding affinity is predicted by scoring these protein-DNA complexes through proximity and interaction parameters obtained by molecular dynamics simulations. The promoters are ranked considering hydrogen bonding, energy of interaction, buried surface area, and distance between centers of masses in interaction with the protein. The ranking is validated through in vitro transcription assays. The trends of these selected promoter interactions have shown variations parallel to the experimental evaluation, emphasizing the success of the active site determination along with screening of the promoter strength. The promoter interaction of Mtb-SigH can be highly beneficial for understanding the regulation of gene expression of a pathogen and also extends a solid platform to predict promoters for other bacterial σ factors.

The first hydration shell of a protein exhibits heterogeneous behavior owing to several attributes, majorly local polarity and structural flexibility as revealed by solvation dynamics of secondary structural elements. We attempt to recognize the change in complex water counteraction generated due to substantial alteration in flexibility during protein complex formation. The investigation is carried out with the signaling lymphocytic activation molecule (SLAM) family of receptors, expressed by an array of immune cells, and interacting with SLAM-associated protein (SAP), composed of one SH2 domain. All atom molecular dynamics simulations are employed to the aqueous solutions of free SAP and SLAM-peptide bound SAP. We observed that water dynamics around different secondary structural elements became highly affected as well as nicely correlated with the SLAM-peptide induced change in structural rigidity obtained by thermodynamic quantification. A few instances of contradictory dynamic features of water to the change in structural flexibility are explained by means of occluded polar residues by the peptide. For βD, EFloop, and BGloop, both structural flexibility and solvent accessibility of the residues confirm the obvious contribution. Most importantly, we have quantified enhanced restriction in water dynamics around the second Fyn-binding site of the SAP due to SAP-SLAM complexation, even prior to the presence of Fyn. This observation leads to a novel argument that SLAM induced more restricted water molecules could offer more water entropic contribution during the subsequent Fyn binding and provide enhanced stability to the SAP-Fyn complex in the signaling cascade. Finally, SLAM induced water counteraction around the second binding site of the SAP sheds light on the allosteric property of the SAP, which becomes an integral part of the underlying signal transduction mechanism.

Emergence of thousands of crystal structures of noncoding RNA molecules indicates its structural and functional diversity. RNA function is based upon a large variety of structural elements which are specifically assembled in the folded molecules. Along with the canonical Watson-Crick base pairs, different orientations of the bases to form hydrogen-bonded non-canonical base pairs have also been observed in the available RNA structures. Frequencies of occurrences of different non-canonical base pairs in RNA indicate their important role to maintain overall structure and functions of RNA. There are several reports on geometry and energetic stabilities of these non-canonical base pairs. However, their stacking geometry and stacking stability with the neighboring base pairs are not well studied. Among the different non-canonical base pairs, the G:U wobble base pair (G:U W:WC) is most frequently observed in the RNA double helices. Using quantum chemical method and available experimental data set we have studied the stacking geometry of G:U W:WC base pair containing dinucleotide sequences in roll-slide parameters hyperspace for different values of twist. This study indicates that the G:U W:WC base pair can stack well with the canonical base pairs giving rise to large interaction energy. The overall preferred stacking geometry in terms of roll, twist and slide for the eleven possible dinucleotide sequences is seen to be quite dependent on their sequences. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 328–338, 2015.

Designing biologically inspired nanoscale molecular assembly with desired functionality is a challenging endeavour. Here we report the designing of fibrin-inspired nanostructured peptide based sealants which facilitate remarkably fast entrapping of blood corpuscles (~28 seconds) in contrast to fibrin (~56 seconds). Our engineered sealants are stabilized by lysine-aspartate ionic interactions and also by Nε(γ-glutamyl) lysine isopeptide bond mediated covalent interaction. Each sealant is formed by two peptides having complementary charges to promote lysine-aspartate ionic interactions and designed isopeptide bond mediated interactions. Computational analysis reveals the isopeptide bond mediated energetically favourable peptide assemblies in sealants 1–3. Our designed sealants 2 and 3 mimic fibrin-mediated clot formation mechanism in presence of transglutaminase enzyme and blood corpuscles. These fibrin-inspired peptides assemble to form sealants having superior hemostatic activities than fibrin. Designed sealants feature mechanical properties, biocompatibility, biodegradability and high adhesive strength. Such nature-inspired robust sealants might be potentially translated into clinics for facilitating efficient blood clotting to handle traumatic coagulopathy and impaired blood clotting.

The p53 protein activation protects the organism from propagation of cells with damaged DNA having oncogenic mutations. In normal cells, activity of p53 is controlled by interaction with MDM2. The well understood p53-MDM2 interaction facilitates design of ligands that could potentially disrupt or prevent the complexation owing to its emergence as an important objective for cancer therapy. However, thermodynamic quantification of the p53-peptide induced structural changes of the MDM2-protein remains an area to be explored. This study attempts to understand the conformational free energy and entropy costs due to this complex formation from the histograms of dihedral angles generated from molecular dynamics simulations. Residue-specific quantification illustrates that, hydrophobic residues of the protein contribute maximum to the conformational thermodynamic changes. Thermodynamic quantification of structural changes of the protein unfold the fact that, p53 binding provides a source of inter-element cooperativity among the protein secondary structural elements, where the highest affected structural elements (α2 and α4) found at the binding site of the protein affects faraway structural elements (β1 and Loop1) of the protein. The communication perhaps involves water mediated hydrogen bonded network formation. Further, we infer that in inhibitory F19A mutation of P53, though Phe19 is important in the recognition process, it has less prominent contribution in the stability of the complex. Collectively, this study provides vivid microscopic understanding of the interaction within the protein complex along with exploring mutation sites, which will contribute further to engineer the protein function and binding affinity.

The signalling lymphocytic activation molecule (SLAM) family of receptors, expressed by an array of immune cells, associate with SLAM-associated protein (SAP)-related molecules, composed of single SH2 domain architecture. SAP activates Src-family kinase Fyn after SLAM ligation, resulting in a SLAM-SAP-Fyn complex, where, SAP binds the Fyn SH3 domain that does not involve canonical SH3 or SH2 interactions. This demands insight into this SAP mediated signalling cascade. Thermodynamics of the conformational changes are extracted from the histograms of dihedral angles obtained from the all-atom molecular dynamics simulations of this structurally well characterized SAP-SLAM complex. The results incorporate the binding induced thermodynamic changes of individual amino acid as well as the secondary structural elements of the protein and the solvent. Stabilization of the peptide partially comes through a strong hydrogen bonding network with the protein, while hydrophobic interactions also play a significant role where the peptide inserts itself into a hydrophobic cavity of the protein. SLAM binding widens SAP's second binding site for Fyn, which is the next step in the signal transduction cascade. The higher stabilization and less fluctuation of specific residues of SAP in the Fyn binding site, induced by SAP-SLAM complexation, emerge as the key structural elements to trigger the recognition of SAP by the SH3 domain of Fyn. The thermodynamic quantification of the protein due to complexation not only throws deeper understanding in the established mode of SAP-SLAM interaction but also assists in the recognition of the relevant residues of the protein responsible for alterations in its activity.

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by B97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 134-147, 2015.

Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence-directed features, such as small roll values for the purine-pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near-zero slide to be most favorable for the purine-pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar-phosphate backbone with C3-endo sugars and this demands C1-C1 distance of about 5.4 angstrom along the chains. Consideration of an energy penalty term for deviation of C1-C1 distance from the mean value, to the recent DFT-D functionals, specifically B97X-D appears to predict reliable energy contour for AU/AU step. Such distance-based penalty improves energy contours for the other purine-pyrimidine sequences also. (c) 2013 Wiley Periodicals, Inc. Biopolymers 101: 107-120, 2014.

A large amount of experimental evidence is available on the effect of magnesium ions on the structure and stability of DNA double helix. Less is known, however, on how these ions affect the stability and dynamics of the molecule. The static time average pictures from X-ray structures or the quantum chemical energy minimized structures lack understanding of the dynamic DNA–ion interaction. The present work addresses these questions by molecular dynamics simulation studies on two DNA duplexes and their interaction with magnesium ions. Results show typical B-DNA character with occasional excursions to deviated states. We detected expected stability of the duplexes in terms of backbone conformations and base pair parameter by the CHARMM-27 force field. Ion environment analysis shows that Mg2+ retains the coordination sphere throughout the simulation with a preference for major groove over minor. An extensive analysis of the influence of the Mg2+ ion shows no evidence of the popular predictions of groove width narrowing by dipositive metal ion. The major groove atoms show higher occupancy and residence time compared to minor groove for magnesium, where no such distinction is found for the charge neutralizing Na+ ions. The determining factor of Mg2+ ion’s choice in DNA binding site evolves as the steric hindrance faced by the bulky hexahydrated cation where wider major groove gets the preference. We have shown that in case of binding of Mg2+ to DNA non electrostatic contributions play a major role. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:5

Genomic DNA of higher organisms exists as extremely long polymers, while in bacteria and other lower organisms it is circular with no terminal base pairs. Temperature-induced melting of the DNA double helix by localized strand separation has been unattainable by molecular dynamic simulations due to more rapid fraying of the terminal base pairs in oligomeric DNA. However, local-sequence-dependent unfolding of the DNA double helix is extremely important for understanding various biochemical phenomena, and can be addressed by simulating a model polymeric DNA duplex. Here, we present simulations of polymeric B-DNA of sequence d(CGCGCGCGAATTCGCGCGCG)2 at elevated temperatures, along with its equivalent oligomeric constructs for comparison. Initiation of temperature-induced DNA melting was observed with higher fluctuations of the central d(AATT) region only in the model polymer. The polymeric construct shows a definite melting start site at the weaker A/T stretch, which propagates slowly through the CG rich regions. The melting is reflected in the hydrogen bond breaking, i.e. basepair opening, and by disruption of stacking interaction between successive basepairs. Melting at higher temperature of the oligomer, however, was only through terminal fraying, as also reported earlier.

Quantitative-structure-toxicity-relationship (QSTR) models are developed for predicting the toxicity (pIGC50) of 252 aliphatic compounds on Tetrahymena pyriformis. The single parameter models with a simple molecular descriptor, the number of atoms in the molecule, provide reasonable results. Better QSTR models with two parameters result when global electrophilicity is used as the second descriptor. In order to tackle both charge-and frontier-controlled reactions the importance of the local electro (nucleo) philicities and atomic charges is also analysed.

DNA within the living cells experiences a diverse range of temperature, ranging from freezing condition to hot spring water. How the structure, the mechanical properties of DNA, and the solvation dynamics around DNA changes with the temperature is important to understand the functionality of DNA under those acute temperature conditions. In that notion, we have carried out molecular dynamics simulations of a DNA oligomer, containing TATA-box sequence for three different temperatures (250 K, 300 K and 350 K). We observed that the structure of the DNA, in terms of backbone torsion angles, sugar pucker, base pair parameters, and base pair step parameters, did not show any unusual properties within the studied range of temperatures, but significant structural alteration was noticed between BI and BII forms at higher temperature. As expected, the flexibility of the DNA, in terms of the torsional rigidity and the bending rigidity is highly temperature dependent, confirming that flexibility increases with increase in temperature. Additionally, the groove widths of the studied DNA showed temperature sensitivity, specifically, the major groove width decreases and the minor groove width increases, respectively, with the increase in temperature. We observed that at higher temperature, water around both the major and the minor groove of the DNA is less structured. However, the water dynamics around the minor groove of the DNA is more restricted as compared to the water around the major groove throughout the studied range of temperatures, without any anomalous behavior.

Deformation of DNA takes place quite often due to binding of small molecules or proteins with DNA. Such deformation is significant due to minor groove binding and, besides electrostatic interactions, other non-covalent interactions may also play an important role in generating such deformation. TATA-box binding protein (TBP) binds to the minor groove of DNA at the TATA box sequence, producing a large-scale deformation in DNA and initiating transcription. In order to observe the interactions of protein residues with DNA in the minor groove that produce the deformation in the DNA structure, we carried out molecular dynamics simulations of the TBP-DNA system. The results reveal consistent partial intercalation of two Phe residues, distorting stacking interactions at two dinucleotide step sites. We carried out calculations based on dispersion-corrected density functional theory to understand the source of such stabilization. We observed favorable interaction energies between the Phe residues and the base pairs with which they interact. We suggest that salt-bridge interactions between the phosphate groups and Lys or Arg residues, along with the intercalation of Phe residues between two base pair stacks, stabilize the kinked and opened-up DNA conformation.

Understanding unwinding and melting of double helical DNA is very important to characterize role of DNA in replication, transcription, translation etc. Sequence dependent melting thermodynamics is used extensively for detecting promoter regions but melting studies are generally done for short oligonucleotides. This study reports several molecular dynamics (MD) simulations of homopolymeric poly(dA).poly(dT) as regular oligonucleotide fragments as well as its corresponding polymeric constructs with water and charge-neutralizing counterions at different temperatures ranging from 300 to 400 K. We have eliminated the end-effect or terminal peeling propensity by employing MD simulation of DNA oligonucleotides in such a manner that gives rise to properties of polymeric DNA of infinite length. The dynamic properties such as basepairing and stacking geometry, groove width, backbone conformational parameters, bending, distribution of counter ions and number of hydrogen bonds of oligomeric and polymeric constructs of poly(dA).poly(dT) have been analyzed. The oligomer shows terminal fraying or peeling effect at temperatures above 340 K. The polymer shows partial melting at elevated temperatures although complete denaturations of basepairs do not take place. The analysis of cross strand hydrogen bonds shows that the number of N–H···O hydrogen bonds increases with increase in temperature while C–H···O hydrogen bond frequencies decrease with temperature. Restructuring of counterions in the minor groove with temperature appear as initiation of melting in duplex structures.

The three-dimensional structure of DNA contains various sequence-dependent structural information, which control many cellular processes in life, such as replication, transcription, DNA repair, etc. For the above functions, DNA double helices need to unwind or melt locally, which is different from terminal melting, as often seen in molecular dynamics (MD) simulations or even in many DNA crystal structures. We have carried out detailed MD simulations of DNA double helices of regular oligonucleotide fragments as well as in polymeric constructs with water and charge-neutralizing counter-ions at several different temperatures. We wanted to eliminate the end-effect or terminal melting propensity by employing MD simulation of DNA oligonucleotides in such a manner that gives rise to properties of polymeric DNA of infinite length. The polymeric construct is expected to allow us to see local melting at elevated temperatures. Comparative structural analysis of oligonucleotides and its corresponding virtual polymer at various temperatures ranging from 300 K to 400 K is discussed. The general behaviour, such as volume expansion coefficients of both the simulations show high similarity, indicating polymeric construct, does not give many artificial constraints. Local melting of a polymer, even at elevated temperature, may need a high nucleation energy that was not available in the short (7 ns) simulations. We expected to observe such nucleation followed by cooperative melting of the polymers in longer MD runs. Such simulations of different polymeric sequences would facilitate us to predict probable melting origins in a polymeric DNA.

Pigeon liver fatty acid synthetase was inactivated irreversibly by 2, 4, 6–trinitrobenzenesulphonic acid (TNBS). Biphasic inactivation of the enzyme was observed with the inhibitor. NADPH provided protection to the enzyme against inactivation by TNBS and the extent of protection increased with NADPH concentration indicating that the essential lysine residues are present at the NADPH binding site. The stoichiometric results with TNBS showed that 4 mol of lysine residues are modified per mole of fatty acid synthetase upon complete inactivation. The rapid reaction of two amino groups per enzyme molecule led to the loss of 60% of the enzyme activity. These approaches suggested that two lysine residues present at the active site are essential for the enzymatic activity of fatty acid synthetase.

Pigeon liver fatty acid synthetase (FAS) was rapidly inactivated by pyridoxal 5'-phosphate (PLP). Assays of the partial activities of the PLP-treated synthetase showed that only the enoyl-CoA reductase was decreased significantly. The inactivation of both the overall activity and enoyl-CoA reductase activity of FAS by PLP could be reversed by dialysis or dilution but not by reduction with sodium borohydride. Malonyl-CoA and acetyl-CoA did not protect the enzyme, whereas NADPH provided 68% protection against PLP-inactivation indicating that PLP modified lysine residues present at or near the co-enzyme binding site. PLP-treated enzyme after reduction with sodium borohydride exhibited fluorescence with a maximum at 397 nm (irradiation at 325 nm). Stoichiometric analysis showed that modification of four lysine residues per enzyme molecule resulted in complete inactivation of the overall and enoyl-CoA reductase activities of FAS. NADPH prevented the inactivation by protecting two of these lysine residues from modification, suggesting the presence of two essential lysine residues per enzyme molecule. These results are consistent with the hypothesis that each subunit of the enzyme contains an enoyl-CoA reductase domain in which a lysine residue, at or near the active site, interacts with NADPH.

Pigeon liver fatty acid synthetase was inactivated by stoichiometric concentrations of diethylpyrocarbonate (DEP). The second order rate constant for the loss of synthetase activity was similar to the value for enoyl-CoA reductase indicating that ethoxyformylation destroys the ability of the enzyme to reduce the unsaturated acyl intermediate, without significant effect on β-ketoacyl reductase activity. NADPH provided protection to the enzyme against inactivation by DEP indicating that essential histidine residues are present at the active site. DEP-modified enzyme showed a characteristic absorption maxima at 240 nm confirming the formation of ethoxyformic histidine. The reversal of inactivation by hydroxylamine and a pKa value of 7.0 obtained from the pH-rate profile for inactivation again confirmed the specificity of DEP for hisidine. Stoichiometric results showed that two moles of histidine residue per mole of enzyme are essential for the activity of FAS.