Simulation-predicted and -explained inheritance model of pathogenicity confirmed by transgenic mice models.

Variants in the gap junction beta-2 (GJB2) gene are the most common cause of hereditary hearing impairment. However, how GJB2 variants lead to local physicochemical and structural changes in the hexameric ion channels of connexin 26 (Cx26), resulting in hearing impairment, remains elusive. In this study, using molecular dynamics (MD) simulations, we showed that detached inner-wall N-terminal "plugs" aggregated to reduce the channel ion flow in a highly prevalent V37I variant in humans. To examine the predictive ability of the computational platform, an artificial mutant, V37M, of which the effect was previously unknown in hearing loss, was created. Microsecond simulations showed that homo-hexameric V37M Cx26 hemichannels had an abnormal affinity between the inner edge and N-termini to block the narrower side of the cone-shaped Cx26, while the most stable hetero-hexameric channels did not. From the perspective of the conformational energetics of WT and variant Cx26 hexamers, we propose that unaffected carriers could result from a conformational predominance of the WT and pore-shrinkage-incapable hetero-hexamers, while mice with homozygous variants can only harbor an unstable and dysfunctional N-termini-blocking V37M homo-hexamer. Consistent with these predictions, homozygous V37M transgenic mice exhibited apparent hearing loss, but not their heterozygous counterparts, indicating a recessive inheritance mode. Reduced channel conductivity was found in Gjb2V37M/V37M outer sulcus and Claudius cells but not in Gjb2WT/WT cells. We view that the current computational platform could serve as an assessment tool for the pathogenesis and inheritance of GJB2-related hearing impairments and other diseases caused by connexin dysfunction.

Full-Privacy Secured Search Engine Empowered by Efficient Genome-Mapping Algorithms.

Since the 90s, keyword-based search engines have been the only option for people to locate relevant web content through a simple query comprising one to a few keywords. These engines, whether free or paid, retained users' search queries and preferences, often to deliver targeted ads. Additionally, user-uploaded articles for plagiarism detection can further be stored as part of service providers' expanding databases for profit. Essentially, users could not search without exposing their queries to these providers. We present a new solution here: a method for searching the internet using a full article as a query without disclosing the content. Our Sapiens Aperio Veritas Engine (S.A.V.E.) uses an encoding scheme and an FM-index search, borrowed from next-generation human genome sequencing. Each word in a user's query is transformed into one of 12 "amino acids" to create a pseudo-biological sequence (PBS) on the user's device. Plagiarism checks are done by users submitting their locally created PBSs to our cloud service. This detects identical content in our database, which includes all English and Chinese Wikipedia articles and Open Access journals up to April 2021. PBSs, longer than 12 "amino acids", show accurate results with less than 0.8% false positives. Performance-wise, S.A.V.E. runs at a similar genome-mapping speed as Bowtie and is >5 orders faster than BLAST. With both standard and private modes, S.A.V.E. offers a revolutionary, privacy-first search and plagiarism check system. We believe this sets an exciting precedent for future search engines prioritizing user confidentiality. S.A.V.E. can be accessed at https://dyn.life.nthu.edu.tw/SAVE/.

Helical structure motifs made searchable for functional peptide design.

The systematic design of functional peptides has technological and therapeutic applications. However, there is a need for pattern-based search engines that help locate desired functional motifs in primary sequences regardless of their evolutionary conservation. Existing databases such as The Protein Secondary Structure database (PSS) no longer serves the community, while the Dictionary of Protein Secondary Structure (DSSP) annotates the secondary structures when tertiary structures of proteins are provided. Here, we extract 1.7 million helices from the PDB and compile them into a database (Therapeutic Peptide Design database; TP-DB) that allows queries of compounded patterns to facilitate the identification of sequence motifs of helical structures. We show how TP-DB helps us identify a known purification-tag-specific antibody that can be repurposed into a diagnostic kit for Helicobacter pylori. We also show how the database can be used to design a new antimicrobial peptide that shows better Candida albicans clearance and lower hemolysis than its template homologs. Finally, we demonstrate how TP-DB can suggest point mutations in helical peptide blockers to prevent a targeted tumorigenic protein-protein interaction. TP-DB is made available at http://dyn.life.nthu.edu.tw/design/ .

An Effective and Safe Enkephalin Analog for Antinociception.

Opioids account for 69,000 overdose deaths per annum worldwide and cause serious side effects. Safer analgesics are urgently needed. The endogenous opioid peptide Leu-Enkephalin (Leu-ENK) is ineffective when introduced peripherally due to poor stability and limited membrane permeability. We developed a focused library of Leu-ENK analogs containing small hydrophobic modifications. N-pivaloyl analog KK-103 showed the highest binding affinity to the delta opioid receptor (68% relative to Leu-ENK) and an extended plasma half-life of 37 h. In the murine hot-plate model, subcutaneous KK-103 showed 10-fold improved anticonception (142%MPE·h) compared to Leu-ENK (14%MPE·h). In the formalin model, KK-103 reduced the licking and biting time to ~50% relative to the vehicle group. KK-103 was shown to act through the opioid receptors in the central nervous system. In contrast to morphine, KK-103 was longer-lasting and did not induce breathing depression, physical dependence, and tolerance, showing potential as a safe and effective analgesic.

An Efficient Timer and Sizer of Biomacromolecular Motions.

Life ticks as fast as how proteins move. Computationally expensive molecular dynamics simulation has been the only theoretical tool to gauge the time and sizes of these motions, though barely to their slowest ends. Here, we convert a computationally cheap elastic network model (ENM) into a molecular timer and sizer to gauge the slowest functional motions of structured biomolecules. Quasi-harmonic analysis, fluctuation profile matching, and the Wiener-Khintchine theorem are used to define the "time periods," t, for anharmonic principal components (PCs), which are validated by nuclear magnetic resonance (NMR) order parameters. The PCs with their respective "time periods" are mapped to the eigenvalues (λENM) of the corresponding ENM modes. Thus, the power laws t(ns) = 56.1λENM-1.6 and σ2(Å2) = 32.7λENM-3.0 can be established allowing the characterization of the timescales of NMR-resolved conformers, crystallographic anisotropic displacement parameters, and important ribosomal motions, as well as motional sizes of the latter.

Resolution-exchanged structural modeling and simulations jointly unravel that subunit rolling underlies the mechanism of programmed ribosomal frameshifting.

MOTIVATION: Programmed ribosomal frameshifting (PRF) is widely used by viruses and bacteria to produce different proteins from a single mRNA template. How steric hindrance of a PRF-stimulatory mRNA structure transiently modifies the conformational dynamics of the ribosome, and thereby allows tRNA slippage, remains elusive. RESULTS: Here, we leverage linear response theories and resolution-exchanged simulations to construct a structural/dynamics model that connects and rationalizes existing structural, single-molecule and mutagenesis data by resolution-exchanged structural modelling and simulations. Our combined theoretical techniques provide a temporal and spatial description of PRF with unprecedented mechanistic details. We discover that ribosomal unfolding of the PRF-stimulating pseudoknot exerts resistant forces on the mRNA entrance of the ribosome, and thereby drives 30S subunit rolling. Such motion distorts tRNAs, leads to tRNA slippage, and in turn serves as a delicate control of cis-element's unwinding forces over PRF. AVAILABILITY AND IMPLEMENTATION: All the simulation scripts and computational implementations of our methods/analyses (including linear response theory) are included in the bioStructureM suite, provided through GitHub at https://github.com/Yuan-Yu/bioStructureM. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DynOmics: dynamics of structural proteome and beyond.

DynOmics (dynomics.pitt.edu) is a portal developed to leverage rapidly growing structural proteomics data by efficiently and accurately evaluating the dynamics of structurally resolved systems, from individual molecules to large complexes and assemblies, in the context of their physiological environment. At the core of the portal is a newly developed server, ENM 1.0, which permits users to efficiently generate information on the collective dynamics of any structure in PDB format, user-uploaded or database-retrieved. ENM 1.0 integrates two widely used elastic network models (ENMs)-the Gaussian Network Model (GNM) and the Anisotropic Network Model (ANM), extended to take account of molecular environment. It enables users to assess potentially functional sites, signal transduction or allosteric communication mechanisms, and protein-protein and protein-DNA interaction poses, in addition to delivering ensembles of accessible conformers reconstructed at atomic details based on the global modes of motions predicted by the ANM. The 'environment' is defined in a flexible manner, from lipid bilayer and crystal contacts, to substrate or ligands bound to a protein, or surrounding subunits in a multimeric structure or assembly. User-friendly interactive features permit users to easily visualize how the environment alter the intrinsic dynamics of the query systems. ENM 1.0 can be accessed at http://enm.pitt.edu/ or http://dyn.life.nthu.edu.tw/oENM/.

Structural Insights into Substrate Recognition by Clostridium difficile Sortase.

Sortases function as cysteine transpeptidases that catalyze the covalent attachment of virulence-associated surface proteins into the cell wall peptidoglycan in Gram-positive bacteria. The substrate proteins targeted by sortase enzymes have a cell wall sorting signal (CWSS) located at the C-terminus. Up to date, it is still not well understood how sortases with structural resemblance among different classes and diverse species of bacteria achieve substrate specificity. In this study, we focus on elucidating the molecular basis for specific recognition of peptide substrate PPKTG by Clostridium difficile sortase B (Cd-SrtB). Combining structural studies, biochemical assays and molecular dynamics simulations, we have constructed a computational model of Cd-SrtBΔN26-PPKTG complex and have validated the model by site-directed mutagensis studies and fluorescence resonance energy transfer (FRET)-based assay. Furthermore, we have revealed that the fourth amino acid in the N-terminal direction from cleavage site of PPKTG forms specific interaction with Cd-SrtB and plays an essential role in configuring the peptide to allow more efficient substrate-specific cleavage by Cd-SrtB.

Structure and function of chicken interleukin-1 beta mutants: uncoupling of receptor binding and in vivo biological activity.

Receptor-binding and subsequent signal-activation of interleukin-1 beta (IL-1ß) are essential to immune and proinflammatory responses. We mutated 12 residues to identify sites important for biological activity and/or receptor binding. Four of these mutants with mutations in loop 9 (T117A, E118K, E118A, E118R) displayed significantly reduced biological activity. Neither T117A nor E118K mutants substantially affected receptor binding, whereas both mutants lack the IL-1ß signaling in vitro but can antagonize wild-type (WT) IL-1ß. Crystal structures of T117A, E118A, and E118K revealed that the secondary structure or surface charge of loop 9 is dramatically altered compared with that of wild-type chicken IL-1ß. Molecular dynamics simulations of IL-1ß bound to its receptor (IL-1RI) and receptor accessory protein (IL-1RAcP) revealed that loop 9 lies in a pocket that is formed at the IL-1RI/IL-1RAcP interface. This pocket is also observed in the human ternary structure. The conformations of above mutants in loop 9 may disrupt structural packing and therefore the stability in a chicken IL-1ß/IL-1RI/IL-1RAcP signaling complex. We identify the hot spots in IL-1ß that are essential to immune responses and elucidate a mechanism by which IL-1ß activity can be inhibited. These findings should aid in the development of new therapeutics that neutralize IL-1 activity.

iGNM 2.0: the Gaussian network model database for biomolecular structural dynamics.

Gaussian network model (GNM) is a simple yet powerful model for investigating the dynamics of proteins and their complexes. GNM analysis became a broadly used method for assessing the conformational dynamics of biomolecular structures with the development of a user-friendly interface and database, iGNM, in 2005. We present here an updated version, iGNM 2.0 http://gnmdb.csb.pitt.edu/, which covers more than 95% of the structures currently available in the Protein Data Bank (PDB). Advanced search and visualization capabilities, both 2D and 3D, permit users to retrieve information on inter-residue and inter-domain cross-correlations, cooperative modes of motion, the location of hinge sites and energy localization spots. The ability of iGNM 2.0 to provide structural dynamics data on the large majority of PDB structures and, in particular, on their biological assemblies makes it a useful resource for establishing the bridge between structure, dynamics and function.

Comparative study on two methods of back propagation network test in TCM syndrome typing of depression / 中国中西医结合杂志

<p><b>OBJECTIVE</b>To compare the effect of the two methods of back propagation network (BPN) test on TCM syndrome typing of depression.</p><p><b>METHODS</b>Test was carried out by two methods as following: (1) Cross train-test method: 1731 patients with depression typed to 5 syndrome types were randomly divided into 2 groups, and they were trained and tested in turn; (2) Round-Robin method: Test was conducted in an altered cycle mode, that is, in a cycle, one out of the 1731 patients were selected to be tested, while the others were trained, the next cycle started when the test on the selected patient was finished and another one for test was selected. In this way, one cycle after the other, until all patients had been tested.</p><p><b>RESULT</b>The total training sensitivity of the two methods was 97.9% and 98.2% respectively, and the total testing sensitivity was 72.7% and 74.2% respectively.</p><p><b>CONCLUSION</b>(1) The five TCM syndrome types of depression could be well differentiated by BPN, which is valuable for TCM syndrome typing in certain extent; (2) The sensitivity of Round-Robin method is slightly higher than that of Cross train-test method, but in comparison between them no remarkable significance was shown.</p>