Polypyridine-Based Helical Amide Foldamer Channels: Rapid Transport of Water and Protons with High Ion Rejection.

Synthetic strategies that enable rapid construction of covalent organic nanotubes with an angstrom-scale tubular pore remain scarcely reported. Reported here is a remarkably simple and mild one-pot polymerization protocol, employing POCl3 as the polymerization agent. This protocol efficiently generates polypyridine amide foldamer-based covalent organic nanotubes with a 2.8 nm length at a yield of 50 %. Trapping single-file water chains in the 2.8 Štubular cavity, rich in hydrogen-bond donors and acceptors, these tubular polypyridine ensembles rapidly and selectively transport water at a rate of 1.6×109  H2 O⋅S-1 ⋅channel-1 and protons at a speed as fast as gramicidin A, with a high rejection of ions.

OnionNet: a Multiple-Layer Intermolecular-Contact-Based Convolutional Neural Network for Protein-Ligand Binding Affinity Prediction.

Computational drug discovery provides an efficient tool for helping large-scale lead molecule screening. One of the major tasks of lead discovery is identifying molecules with promising binding affinities toward a target, a protein in general. The accuracies of current scoring functions that are used to predict the binding affinity are not satisfactory enough. Thus, machine learning or deep learning based methods have been developed recently to improve the scoring functions. In this study, a deep convolutional neural network model (called OnionNet) is introduced; its features are based on rotation-free element-pair-specific contacts between ligands and protein atoms, and the contacts are further grouped into different distance ranges to cover both the local and nonlocal interaction information between the ligand and the protein. The prediction power of the model is evaluated and compared with other scoring functions using the comparative assessment of scoring functions (CASF-2013) benchmark and the v2016 core set of the PDBbind database. The robustness of the model is further explored by predicting the binding affinities of the complexes generated from docking simulations instead of experimentally determined PDB structures.

Distinct Dynamic and Conformational Features of Human STING in Response to 2'3'-cGAMP and c-di-GMP.

The human stimulator of interferon genes protein (hSTING) can bind cyclic dinucleotides (CDNs) to activate the production of type I interferons and inflammatory cytokines. These CDNs can be either bacterial secondary messengers, 3'3'-CDNs, or endogenous 2'3'-cGAMP. cGAMP, with a unique 2'-5' bond, is the most potent activator of hSTING among all CDNs. However, current understanding of the molecular principles underlying the unique ability of 2'3'-cGAMP to potently activate hSTINGs other than 3'3'-CDNs remains incomplete. In this work, molecular dynamics simulations were used to provide an atomistic picture of the binding of 2'3'-cGAMP and one 3'3'-CDN (c-di-GMP) to hSTING. The results suggest that hSTING binds more strongly to 2'3'-cGAMP than to c-di-GMP, which prefers to bind with a more open and flexible state of hSTING. Finally, a potential "dock-lock-anchor" mechanism is proposed for the activation of hSTING upon the binding of a potent ligand. It is believed that deep insights into understanding the binding of hSTING with 3'3'-CDNs and the endogenous 2'3'-cGAMP would help to establish the principles underlying powerful 2'3'-cGAMP signaling and the nature of hSTING activation, as well as related drug design.

Exploring the structure characteristics and major channels of cytochrome P450 2A6, 2A13, and 2E1 with pilocarpine.

The majority of cytochromes P450 play a critical role in metabolism of endogenous and exogenous substrates, some of its products are carcinogens. Therefore, inhibition of P450 enzymes activity can promote the detoxification and elimination of chemical carcinogens. In this study, molecular dynamics (MD) simulations and adaptive steered molecular dynamics (ASMD) simulations were performed to explore the structure features and channel dynamics of three P450 isoforms 2A6, 2A13, and 2E1 bound with the common inhibitor pilocarpine. The binding free energy results combined with the PMF calculations give a reasonable ranking of binding affinity, which are consistent with the experimental data. Our results uncover how a sequence divergence of different CYP2 enzymes causes individual variations in major channel selections. On the basis of channel bottleneck and energy decomposition analysis, we propose a gating mechanism of their respective major channels in three enzymes, which may be attributed to a reversal of Phe209 in CYP2A6/2A13, as well as the rotation of Phe116 and Phe298 in CYP2E1. The hydrophobic residues not only make strong hydrophobic interactions with inhibitor, but also act as gatekeeper to regulate the opening of channel. The present study provides important insights into the structure-function relationships of three cytochrome P450s and the molecular basis for development of potent inhibitors.

Studying the recognition mechanism of TcaR and ssDNA using molecular dynamic simulations.

The transcription regulator teicoplanin-associate locus regulator (TcaR) plays a vital role in interfering with ssDNA replication and resisting ssDNA phage invasion. Although recent studies demonstrated that TcaR had strong interaction with ssDNA, the dynamics and interaction mechanism of dimeric TcaR bound to ssDNA have not been rationalized at the atomic level. In our study, MD simulations combined with MM-GB/SA calculations were employed to study recognition mechanism between TcaR and ssDNA. The results illuminate that electrostatic interaction is the main driving force for the binding process. We put forward that six anchoring residues (Arg70, Arg71, Ser188, Gln191, Arg221 and Arg222) play a vital role in stabilizing the ssDNA by forming strong hydrogen bond and salt bridge interactions. TcaR undergoes the asymmetric conformational changes at the wHTH domain upon binding to ssDNA. This may be attributed to the changing of electrostatic potential, enhanced contacts and salt bridge interaction. The present study provides new insights into the recognition mechanism of TcaR bound to ssDNA, which could contribute to understanding about the multiple TcaR functions in staphylococci enrich our understanding of MarR family.

Investigation of ligand selectivity in CYP3A7 by molecular dynamics simulations.

Cytochrome P450 (CYP) 3A7 plays a crucial role in the biotransformation of the metabolized endogenous and exogenous steroids. To compare the metabolic capabilities of CYP3A7-ligands complexes, three endogenous ligands were selected, namely dehydroepiandrosterone (DHEA), estrone, and estradiol. In this study, a three-dimensional model of CYP3A7 was constructed by homology modeling using the crystal structure of CYP3A4 as the template and refined by molecular dynamics simulation (MD). The docking method was adopted, combined with MD simulation and the molecular mechanics generalized born surface area method, to probe the ligand selectivity of CYP3A7. These results demonstrate that DHEA has the highest binding affinity, and the results of the binding free energy were in accordance with the experimental conclusion that estrone is better than estradiol. Moreover, several key residues responsible for substrate specificity were identified on the enzyme. Arg372 may be the most important residue due to the low interaction energies and the existence of hydrogen bond with DHEA throughout simulation. In addition, a cluster of Phe residues provides a hydrophobic environment to stabilize ligands. This study provides insights into the structural features of CYP3A7, which could contribute to further understanding of related protein structures and dynamics.

Structural features and dynamic investigations of the membrane-bound cytochrome P450 17A1.

Cytochrome P450 (CYP) 17A1 is a dual-function monooxygenase with a critical role in the synthesis of many human steroid hormones. The enzyme is an important target for treatment of breast and prostate cancers that proliferate in response to estrogens and androgens. Despite the crystallographic structures available for CYP17A1, no membrane-bound structural features of this enzyme at atomic level are available. Accumulating evidence has indicated that the interactions between bounded CYPs and membrane could contribute to the recruitment of lipophilic substrates. To this end, we have investigated the effects on structural characteristics in the presence of the membrane for CYP17A1. The MD simulation results demonstrate a spontaneous insertion process of the enzyme to the lipid. Two predominant modes of CYP17A1 in the membrane are captured, characterized by the depths of insertion and orientations of the enzyme to the membrane surface. The measured heme tilt angles show good consistence with experimental data, thereby verifying the validity of the structural models. Moreover, conformational changes induced by the membrane might have impact on the accessibility of the active site to lipophilic substrates. The dynamics of internal aromatic gate formed by Trp220 and Phe224 are suggested to regulate tunnel opening motions. The knowledge of the membrane binding characteristics could guide future experimental and computational works on membrane-bound CYPs so that various investigations of CYPs in their natural, lipid environment rather than in artificially solubilized forms may be achieved.