Refine
Document Type
- Article (5)
Language
- English (5) (remove)
Has Fulltext
- yes (5)
Is part of the Bibliography
- no (5) (remove)
Keywords
- Computational biophysics (1)
- Coxsackievirus B3 (1)
- Energy transfer (1)
- Infrared spectroscopy (1)
- Proteins (1)
- RNA hairpin (1)
- RNA jepit rambut (1)
- ion treatment (1)
- md simulations (1)
- perlakuan ion (1)
Vibrational energy transfer (VET) is essential for protein function. It is responsible for efficient energy dissipation in reaction sites, and has been linked to pathways of allosteric communication. While it is understood that VET occurs via backbone as well as via non-covalent contacts, little is known about the competition of these two transport channels, which determines the VET pathways. To tackle this problem, we equipped the β-hairpin fold of a tryptophan zipper with pairs of non-canonical amino acids, one serving as a VET injector and one as a VET sensor in a femtosecond pump probe experiment. Accompanying extensive non-equilibrium molecular dynamics simulations combined with a master equation analysis unravel the VET pathways. Our joint experimental/computational endeavor reveals the efficiency of backbone vs. contact transport, showing that even if cutting short backbone stretches of only 3 to 4 amino acids in a protein, hydrogen bonds are the dominant VET pathway.
Molecular dynamics has been employed to study the effect of ion treatment on the stability of 14-nucleotide RNA hairpin of Coxsackievirus B3. Three AMBER force fields were used: AMBER94, AMBER98, and AMBER99, which showed no significant structural difference of the hairpin. Thereafter, we applied two different long-range electrostatic treatments that were reaction field and PME methods, and calculated the distribution of ions around the hairpin. Although the structural stabilities of the MD simulations using both methods were similar in 0.14 M Na+, ion environment around the hairpin was notably different. In particular, structural stabilition of the hairpin with increasing ion concentration and with ion Mg2+ cannot be accommodated by simulations using reaction field method. Furthermore, the MD simulations using PME method suggested the strong similarity in structural and dynamical properties of the hairpin with 0.14 M Na+, 0.50 M Na+, 1,03 M Na+, and 0.08 M Mg2+ concentrations. However, the simulations revealed different ion occupations of Na+ and Mg2+.
Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA's secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent ({approx}1 µs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.
In a combined NMR/MD study, the temperature-dependent changes in the conformation of two members of the RNA YNMG-tetraloop motif (cUUCGg and uCACGg) have been investigated at temperatures of 298, 317 and 325 K. The two members have considerable different thermal stability and biological functions. In order to address these differences, the combined NMR/MD study was performed. The large temperature range represents a challenge for both, NMR relaxation analysis (consistent choice of effective bond length and CSA parameter) and all-atom MD simulation with explicit solvent (necessity to rescale the temperature). A convincing agreement of experiment and theory is found. Employing a principle component analysis of the MD trajectories, the conformational distribution of both hairpins at various temperatures is investigated. The ground state conformation and dynamics of the two tetraloops are indeed found to be very similar. Furthermore, both systems are initially destabilized by a loss of the stacking interactions between the first and the third nucleobase in the loop region. While the global fold is still preserved, this initiation of unfolding is already observed at 317 K for the uCACGg hairpin but at a significantly higher temperature for the cUUCGg hairpin.
To facilitate the measurement of intramolecular distances in solvated RNA systems, a combination of spin-labeling, electron paramagnetic resonance (EPR), and molecular dynamics (MD) simulation is presented. The fairly rigid spin label 2,2,5,5-tetramethyl-pyrrolin-1-yloxyl-3-acetylene (TPA) was base and site specifically introduced into RNA through a Sonogashira palladium catalyzed crosscoupling on column. For this purpose 5-iodouridine, 5-iodo-cytidine and 2-iodo-adenosine phosphoramidites were synthesized and incorporated into RNA-sequences. Application of the recently developed ACE (R) chemistry presented the main advantage to limit the reduction of the nitroxide to an amine during the oligonucleotide automated synthesis and thus to increase substantially the reliability of the synthesis and the yield of labeled oligonucleotides. 4-Pulse Electron Double Resonance (PELDOR) was then successfully used to measure the intramolecular spin–spin distances in six doubly labeled RNA-duplexes. Comparison of these results with our previous work on DNA showed that A- and B-Form can be differentiated. Using an all-atom force field with explicit solvent, MD simulations gave results in good agreement with the measured distances and indicated that the RNA A-Form was conserved despite a local destabilization effect of the nitroxide label. The applicability of the method to more complex biological systems is discussed.