In this work we studied the folding process of the hybrid-1 type human telomeric DNA G-quadruplex with solvent and ions explicitly modeled. Enabled by the powerful bias-exchange metadynamics and large-scale conventional molecular dynamic simulations, the free energy landscape of this G-DNA was obtained for the first time and four folding intermediates were identified, including a triplex and a basically formed quadruplex. The simulations also provided atomistic pictures for the structures and cation binding patterns of the intermediates. The results showed that the structure formation and cation binding are cooperative and mutually supporting each other. The syn/anti reorientation dynamics of the intermediates was also investigated. It was found that the nucleotides usually take correct syn/anti configurations when they form native and stable hydrogen bonds with the others, while fluctuating between two configurations when they do not. Misfolded intermediates with wrong syn/anti configurations were observed in the early intermediates but not in the later ones. Based on the simulations, we also discussed the roles of the non-native interactions. Besides, the formation process of the parallel conformation in the first two G-repeats and the associated reversal loop were studied. Based on the above results, we proposed a folding pathway for the hybrid-1 type G-quadruplex with atomistic details, which is new and more complete compared with previous ones. The knowledge gained for this type of G-DNA may provide a general insight for the folding of the other G-quadruplexes.
Cathelicidins are a large family of cationic antimicrobial peptides (AMPs) found in mammals with broad spectrum antimicrobial activity. LL-37 is the sole amphipathic α-helical AMP from human Cathelicidins family. In addition to its bactericidal capability, LL-37 has antiviral, anti-tumor, and immunoregulatory activity. Despite many experimental studies, its molecular mechanism of action is not yet fully understood. Here, we performed three independent molecular dynamics simulations (600 ns or more) of a LL-37 peptide in the presence of 256 lipid bilayers with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) mimicking bacterial and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) mimicking mammalian membranes. We found that LL-37 can be quickly absorbed onto the POPG bilayer without loss of its helical conformation in the core region and with the helix lying in parallel to the bilayer. The POPG bilayer was deformed. In contrast, LL-37 is slower in reaching the POPC surface and loss much of its helical conformation during the interaction with the bilayer. LL-37 only partially entered the POPC bilayer without significant deformation of the membrane. The observed difference for different bilayers is largely due to the fact that LL-37 is positively charged, POPG is negatively charged, and POPC is neutral. Our simulation results demonstrated the initial stage of disruption of the bacterial membrane by LL-37 in atomic details. Comparison to experimental results on LL-37 and simulation studies in other systems was made.
To understand the role of ribose G-quartets and how they affect the properties of G-quadruplex structures, we studied three systems in which one, two, three, or four deoxyribose G-quartets were substituted with ribose G-quartets. These systems were a parallel DNA intramolecular G-quadruplex, d(TTGGGTGGGTTGGGTGGGTT), and two tetramolecular G-quadruplexes, d(TGGGT) and d(TGGGGT). Thermal denaturation experiments revealed that ribose G-quartets have position-dependent and cumulative effects on G-quadruplex stability. An unexpected destabilization was observed when rG quartets were presented at the 5'-end of the G stack. This observation challenges the general belief that RNA residues stabilize G-quadruplexes. Furthermore, in contrast to past proposals, hydration is not the main factor determining the stability of our RNA/DNA chimeric G-quadruplexes. Interestingly, the presence of rG residues in a central G-quartet facilitated the formation of additional tetramolecular G-quadruplex topologies showing positive circular dichroism signals at 295 nm. 2D NMR analysis of the tetramolecular TGgGGT (lowercase letter indicates ribose) indicates that Gs in the 5'-most G-quartet adopt the syn conformation. These analyses highlight several new aspects of the role of ribose G-quartets on G-quadruplex structure and stability, and demonstrate that the positions of ribose residues are critical for tuning G-quadruplex properties.
Long non-coding RNAs (lncRNAs) play important functional roles in many diverse biological processes. However, not all expressed lncRNAs are functional. Thus, it is necessary to manually collect all experimentally validated functional lncRNAs (EVlncRNA) with their sequences, structures, and functions annotated in a central database. The first release of such a database (EVLncRNAs) was made using the literature prior to 1 May 2016. Since then (till 15 May 2020), 19 245 articles related to lncRNAs have been published. In EVLncRNAs 2.0, these articles were manually examined for a major expansion of the data collected. Specifically, the number of annotated EVlncRNAs, associated diseases, lncRNA-disease associations, and interaction records were increased by 260%, 320%, 484% and 537%, respectively. Moreover, the database has added several new categories: 8 lncRNA structures, 33 exosomal lncRNAs, 188 circular RNAs, and 1079 drug-resistant, chemoresistant, and stress-resistant lncRNAs. All records have checked against known retraction and fake articles. This release also comes with a highly interactive visual interaction network that facilitates users to track the underlying relations among lncRNAs, miRNAs, proteins, genes and other functional elements. Furthermore, it provides links to four new bioinformatics tools with improved data browsing and searching functionality. EVLncRNAs 2.0 is freely available at https://www.sdklab-biophysics-dzu.net/EVLncRNAs2/.
Structure-based models or Gō-like models, which are built from one or multiple particular experimental structures, have been successfully applied to the folding of proteins and RNAs. Recently, a variant termed the hybrid atomistic model advances the description of backbone and side chain interactions of traditional structure-based models, by borrowing the description of local interactions from classical force fields. In this study, we assessed the validity of this model in the folding problem of human telomeric DNA G-quadruplex, where local dihedral terms play important roles. A two-state model was developed and a set of molecular dynamics simulations was conducted to study the folding dynamics of sequence Htel24, which was experimentally validated to adopt two different (3 + 1) hybrid G-quadruplex topologies in K solution. Consistent with the experimental observations, the hybrid-1 conformation was found to be more stable and the hybrid-2 conformation was kinetically more favored. The simulations revealed that the hybrid-2 conformation folded in a higher cooperative manner, which may be the reason why it was kinetically more accessible. Moreover, by building a Markov state model, a two-quartet G-quadruplex state and a misfolded state were identified as competing states to complicate the folding process of Htel24. Besides, the simulations also showed that the transition between hybrid-1 and hybrid-2 conformations may proceed an ensemble of hairpin structures. The hybrid atomistic structure-based model reproduced the kinetic partitioning folding dynamics of Htel24 between two different folds, and thus can be used to study the complex folding processes of other G-quadruplex structures.
How RNA sequences fold to specific tertiary structures is one of the key problems for understanding their dynamics and functions. Here, we study the folding process of an H-type RNA pseudoknot by performing a large-scale all-atom MD simulation and bias-exchange metadynamics. The folding free energy landscapes are obtained and several folding intermediates are identified. It is suggested that the folding occurs via multiple mechanisms, including a step-wise mechanism starting either from the first helix or the second, and a cooperative mechanism with both helices forming simultaneously. Despite of the multiple mechanism nature, the ensemble folding kinetics estimated from a Markov state model is single-exponential. It is also found that the correlation between folding and binding of metal ions is significant, and the bound ions mediate long-range interactions in the intermediate structures. Non-native interactions are found to be dominant in the unfolded state and also present in some intermediates, possibly hinder the folding process of the RNA.
In this work we develop an approach for predicting the tertiary structures of RNA fragments by combining an RNA nucleobase discrete state (RNAnbds) model, a sequential Monte Carlo method, and a statistical potential. The RNAnbds model is designed for optimizing the configuration of nucleobases with respect to their preceding ones along the sequence and their spatial neighbors, in contrast to previous works that focus on RNA backbones. The tests of our approach with the fragments taken from a small RNA pseudoknot and a 23S ribosome RNA show that for short fragments (<10 nucleotides), the root mean square deviations (RMSDs) between the predicted and the experimental ones are generally smaller than 3 Å; for slightly longer fragments (10-15 nucleotides), most RMSDs are smaller than 4 Å. The comparison of our method with another physics-based predictor with a testing set containing nine loops shows that ours is superior in both accuracy and efficiency. Our approach is useful in facilitating RNA three-dimensional structure prediction as well as loop modeling. It also holds the promise of providing insight into the structural ensembles of RNA loops.
Thermodynamics of the permeation of amino acids from water to lipid bilayers is an important first step for understanding the mechanism of cell-permeating peptides and the thermodynamics of membrane protein structure and stability. In this work, we employed bias-exchange metadynamics simulations to simulate the membrane permeation of all 20 amino acids from water to the center of a dipalmitoylphosphatidylcholine (DPPC) membrane (consists of 256 lipids) by using both directional and torsion angles for conformational sampling. The overall accuracy for the free energy profiles obtained is supported by significant correlation coefficients (correlation coefficient at 0.5–0.6) between our results and previous experimental or computational studies. The free energy profiles indicated that (1) polar amino acids have larger free energy barriers than nonpolar amino acids; (2) negatively charged amino acids are the most difficult to enter into the membrane; and (3) conformational transitions for many amino acids during membrane crossing is the key for reduced free energy barriers. These results represent the first set of simulated free energy profiles of membrane crossing for all 20 amino acids.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
334 Leonard St
Brooklyn, NY 11211
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.