Two-dimensional ZnO/BlueP van der Waals heterostructure used for visible-light driven water splitting: A first-principles study.

Solar driven water splitting for hydrogen generation has been considered as an important method for collecting clean energy. Herein, based on first-principles calculations, we propose that ZnO/BlueP van der Waals heterostructure can realize overall water splitting reaction for hydrogen generation. Strikingly, the band-gap of 1.83 eV is appropriate, and band alignments straddle the water redox potentials, ensuring the occurrence of hydrogenevolutionreaction and oxygen evolution reaction. Charge density distribution and carrier mobility exhibit significant charge separation and transfer. Visible-light response is improved compared with those of the isolated monolayers. Moreover, hydrogenevolutionreaction is actually realized on the ZnO layer, while oxygen evolution reaction is implemented on the BlueP layer. Through the investigation of the adsorption and dissociation reactions of H2O, we observe that two neighboring H*s prefer to combine to form H2 by overcoming a lowered energy barrier of 0.75 eV. Strain effect indicates that the lateral compressive strain of -4% to 0% and the vertical tensile strain of 0% to +6% can effectively tune band-gap and band alignments. The results indicate that ZnO/BlueP vdW heterostructure is probable highly efficient photoelectric material used for visible-light driven water splitting for hydrogen generation.

The Antibacterial Activity of Quinazoline and Quinazolinone Hybrids.

Bacterial infections cause substantial morbidity and mortality across the world and pose serious threats to humankind. Drug resistance, especially multidrug resistance resulting from different defensive mechanisms in bacteria, is the leading cause of the failure of chemotherapy, making it an urgent need to develop more effective antibacterials. Quinazoline and quinazolinone frameworks have received considerable attention due to their diversified therapeutic potential. In particular, quinazoline/quinazolinone hybrids can exert antibacterial activity through various mechanisms and are useful scaffolds for the discovery of novel antibacterials. This review principally emphasizes the antibacterial potential, structure-activity relationships (SARs), and mechanism of action of quinazoline and quinazolinone hybrids, covering articles published between 2017 and 2021.

Indole Alkaloids, Synthetic Dimers and Hybrids with Potential In Vivo Anticancer Activity.

Indole, a heterocyclic organic compound, is one of the most promising heterocycles found in natural and synthetic sources since its derivatives possess fascinating structural diversity and various therapeutic properties. Indole alkaloids, synthetic dimers and hybrids could act on diverse targets in cancer cells, and consequently, possess potential antiproliferative effects on various cancers both in vitro and in vivo. Vinblastine, midostaurin, and anlotinib as the representative of indole alkaloids, synthetic dimers and hybrids respectively, have already been clinically applied to treat many types of cancers, demonstrating indole alkaloids, synthetic dimers and hybrids are useful scaffolds for the development of novel anticancer agents. Covering articles published between 2010 and 2020, this review emphasizes the recent development of indole alkaloids, synthetic dimers and hybrids with potential in vivo therapeutic application for cancers.

EVLncRNAs 2.0: an updated database of manually curated functional long non-coding RNAs validated by low-throughput experiments.

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/.

Insights into the Kinetic Partitioning Folding Dynamics of the Human Telomeric G-Quadruplex from Molecular Simulations and Machine Learning.

The human telomeric DNA G-quadruplex follows a kinetic partitioning folding mechanism. The underlying folding landscape potentially has many minima separated by high free-energy barriers. However, using current theoretical models to characterize this complex folding landscape has remained a challenging problem. In this study, by developing a hybrid atomistic structure-based model that merges structural information on the hybrid-1, hybrid-2, and chair-type G-quadruplex topologies, we investigated a kinetic partitioning folding process of human telomeric DNA involving three native folds. The model was validated as it reproduced the experimental observation that the hybrid-1 conformation is the major fold and the hybrid-2 conformation is kinetically more accessible. A three-step mechanism was revealed for the formation of the hybrid-1 conformation, while a two-step mechanism was demonstrated for the formation of hybrid-2 and chair-type conformations. Likewise, a class of state in which structures adopted inappropriate combinations of syn/anti guanine nucleotides was found to greatly slow down the folding process. In addition, by employing the XGBoost machine learning algorithm, three interatom distances and six dihedral angles were identified as essential internal coordinates to represent the low-dimensional folding landscape. The strategy of coupling the multibasin model and the machine learning algorithm may be useful to investigate the conformational dynamics of other multistate biomolecules.

Indole/isatin-containing hybrids as potential antibacterial agents.

The emergence and worldwide spread of drug-resistant bacteria have already posed a serious threat to human life, creating the urgent need to develop potent and novel antibacterial drug candidates with high efficacy. Indole and isatin (indole-2,3-dione) present a wide structural and mechanistic diversity, so their derivatives possess various pharmacological properties and occupy a salient place in the development of new drugs. Indole/isatin-containing hybrids, which demonstrate a promising activity against a panel of clinically important Gram-positive and Gram-negative bacteria, are privileged scaffolds for the discovery of novel antibacterial candidates. This review, covering articles published between January 2015 and May 2020, focuses on the development and structure-activity relationship (SAR) of indole/isatin-containing hybrids with potential application for fighting bacterial infections, to facilitate further rational design of novel drug candidates.

Interplay of hydrophobic and hydrophilic interactions in sequence-dependent cell penetration of spontaneous membrane-translocating peptides revealed by bias-exchange metadynamics simulations.

Spontaneous Membrane Translocating Peptides (SMTPs) can translocate silently across the bilayer and, thus, have the best potential to improve the delivery of therapeutic molecules to cells without toxicity. However, how their translocation mechanisms are affected by a specific peptide sequence remains poorly understood. Here, bias-exchange metadynamics simulations were employed to investigate the translocation mechanisms of five SMTPs with the same composition of amino acids (LLRLR, LRLLR, LLLRR, RLLLR, and LRLRL). Simulation results yield sequence-dependent free energy barrier using the FESs along the z-directional distance. An in-depth analysis of sequence-dependent interactions in different regions of the bilayers indicates that the free-energy barrier height of a specific sequence is resulted from the accessibility balance of isolated or clustered hydrophobic residues (L) and hydrophilic residues (R) that leads to different levels of resistance for moving of a peptide into the hydrophobic center of the membrane. At the maximal of the free-energy barrier, all peptides have a conformation parallel to the membrane surface with the barrier height determined by their affinity to the hydrophobic region. The appropriate bilayer perturbation and GDM+ pairing are beneficial for peptide translocation. These results provide an improved microscopic understanding of how peptide sequence influences the translocation efficiency and mechanism.

Self-derived structure-disrupting peptides targeting methionine aminopeptidase in pathogenic bacteria: a new strategy to generate antimicrobial peptides.

Bacterial infection is one of the leading causes of death in young, elderly, and immune-compromised patients. The rapid spread of multi-drug-resistant (MDR) bacteria is a global health emergency and there is a lack of new drugs to control MDR pathogens. We describe a heretofore-unexplored discovery pathway for novel antibiotics that is based on self-targeting, structure-disrupting peptides. We show that a helical peptide, KFF- EcH3, derived from the Escherichia coli methionine aminopeptidase can disrupt secondary and tertiary structure of this essential enzyme, thereby killing the bacterium (including MDR strains). Significantly, no detectable resistance developed against this peptide. Based on a computational analysis, our study predicted that peptide KFF- EcH3 has the strongest interaction with the structural core of the methionine aminopeptidase. We further used our approach to identify peptide KFF- NgH1 to target the same enzyme from Neisseria gonorrhoeae. This peptide inhibited bacterial growth and was able to treat a gonococcal infection in a human cervical epithelial cell model. These findings present an exciting new paradigm in antibiotic discovery using self-derived peptides that can be developed to target the structures of any essential bacterial proteins.-Zhan, J., Jia, H., Semchenko, E. A., Bian, Y., Zhou, A. M., Li, Z., Yang, Y., Wang, J., Sarkar, S., Totsika, M., Blanchard, H., Jen, F. E.-C., Ye, Q., Haselhorst, T., Jennings, M. P., Seib, K. L., Zhou, Y. Self-derived structure-disrupting peptides targeting methionine aminopeptidase in pathogenic bacteria: a new strategy to generate antimicrobial peptides.

Atomistic Analysis of ToxN and ToxI Complex Unbinding Mechanism.

ToxIN is a triangular structure formed by three protein toxins (ToxNs) and three specific noncoding RNA antitoxins (ToxIs). To respond to stimuli, ToxI is preferentially degraded, releasing the ToxN. Thus, the dynamic character is essential in the normal function interactions between ToxN and ToxI. Here, equilibrated molecular dynamics (MD) simulations were performed to study the stability of ToxN and ToxI. The results indicate that ToxI adjusts the conformation of 3' and 5' termini to bind to ToxN. Steered molecular dynamics (SMD) simulations combined with the recently developed thermodynamic integration in 3nD (TI3nD) method were carried out to investigate ToxN unbinding from the ToxIN complex. The potentials of mean force (PMFs) and atomistic pictures suggest the unbinding mechanism as follows: (1) dissociation of the 5' terminus from ToxN, (2) missing the interactions involved in the 3' terminus of ToxI without three nucleotides (G31, A32, and A33), (3) starting to unfold for ToxI, (4) leaving the binding package of ToxN for three nucleotides of ToxI, (5) unfolding of ToxI. This work provides information on the structure-function relationship at the atomistic level, which is helpful for designing new potent antibacterial drugs in the future.

Exploration of the folding dynamics of human telomeric G-quadruplex with a hybrid atomistic structure-based model.

Structure-based models or Go-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.

Molecular Dynamics Simulations of Human Antimicrobial Peptide LL-37 in Model POPC and POPG Lipid Bilayers.

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.

Fast-Folding Pathways of the Thrombin-Binding Aptamer G-Quadruplex Revealed by a Markov State Model.

G-quadruplex structures participate in many important cellular processes. For a better understanding of their functions, knowledge of the mechanism by which they fold into the functional native structures is necessary. In this work, we studied the folding process of the thrombin-binding aptamer G-quadruplex. Enabled by a computational paradigm that couples an advanced sampling method and a Markov state model, four folding intermediates were identified, including an antiparallel G-hairpin, two G-triplex structures, and a double-hairpin conformation. Likewise, a misfolded structure with a nonnative distribution of syn/anti guanines was also observed. Based on these states, a transition path analysis revealed three fast-folding pathways, along which the thrombin-binding aptamer would fold to the native state directly, with no evidence of potential nonnative competing conformations. The results also showed that the TGT-loop plays an important role in the folding process. The findings of this research may provide general insight about the folding of other G-quadruplex structures.

Bias-Exchange Metadynamics Simulation of Membrane Permeation of 20 Amino Acids.

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.

Unexpected Position-Dependent Effects of Ribose G-Quartets in G-Quadruplexes.

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.

Free energy landscape and multiple folding pathways of an H-type RNA pseudoknot.

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.

Atomistic picture for the folding pathway of a hybrid-1 type human telomeric DNA G-quadruplex.

In this work we studied the folding process of the hybrid-1 type human telomeric DNA G-quadruplex with solvent and K(+) 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.

RNA fragment modeling with a nucleobase discrete-state model.

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.