Power generation from ambient humidity using protein nanowires.

Harvesting energy from the environment offers the promise of clean power for self-sustained systems1,2. Known technologies-such as solar cells, thermoelectric devices and mechanical generators-have specific environmental requirements that restrict where they can be deployed and limit their potential for continuous energy production3-5. The ubiquity of atmospheric moisture offers an alternative. However, existing moisture-based energy-harvesting technologies can produce only intermittent, brief (shorter than 50 seconds) bursts of power in the ambient environment, owing to the lack of a sustained conversion mechanism6-12. Here we show that thin-film devices made from nanometre-scale protein wires harvested from the microbe Geobacter sulfurreducens can generate continuous electric power in the ambient environment. The devices produce a sustained voltage of around 0.5 volts across a 7-micrometre-thick film, with a current density of around 17 microamperes per square centimetre. We find the driving force behind this energy generation to be a self-maintained moisture gradient that forms within the film when the film is exposed to the humidity that is naturally present in air. Connecting several devices linearly scales up the voltage and current to power electronics. Our results demonstrate the feasibility of a continuous energy-harvesting strategy that is less restricted by location or environmental conditions than other sustainable approaches.

Toward Accurate Simulation of Coupling between Protein Secondary Structure and Phase Separation.

Intrinsically disordered proteins (IDPs) frequently mediate phase separation that underlies the formation of a biomolecular condensate. Together with theory and experiment, efficient coarse-grained (CG) simulations have been instrumental in understanding the sequence-specific phase separation of IDPs. However, the widely used Cα-only models are limited in capturing the peptide nature of IDPs, particularly backbone-mediated interactions and effects of secondary structures, in phase separation. Here, we describe a hybrid resolution (HyRes) protein model toward a more accurate description of the backbone and transient secondary structures in phase separation. With an atomistic backbone and coarse-grained side chains, HyRes can semiquantitatively capture the residue helical propensity and overall chain dimension of monomeric IDPs. Using GY-23 as a model system, we show that HyRes is efficient enough for the direct simulation of spontaneous phase separation and, at the same time, appears accurate enough to resolve the effects of single His to Lys mutations. HyRes simulations also successfully predict increased ß-structure formation in the condensate, consistent with available experimental CD data. We further utilize HyRes to study the phase separation of TPD-43, where several disease-related mutants in the conserved region (CR) have been shown to affect residual helicities and modulate the phase separation propensity as measured by the saturation concentration. The simulations successfully recapitulate the effect of these mutants on the helicity and phase separation propensity of TDP-43 CR. Analyses reveal that the balance between backbone and side chain-mediated interactions, but not helicity itself, actually determines phase separation propensity. These results support that HyRes represents an effective protein model for molecular simulation of IDP phase separation and will help to elucidate the coupling between transient secondary structures and phase separation.

Backbone interactions and secondary structures in phase separation of disordered proteins.

Intrinsically disordered proteins (IDPs) are one of the major drivers behind the formation and characteristics of biomolecular condensates. Due to their inherent flexibility, the backbones of IDPs are significantly exposed, rendering them highly influential and susceptible to biomolecular phase separation. In densely packed condensates, exposed backbones have a heightened capacity to interact with neighboring protein chains, which might lead to strong coupling between the secondary structures and phase separation and further modulate the subsequent transitions of the condensates, such as aging and fibrillization. In this mini-review, we provide an overview of backbone-mediated interactions and secondary structures within biomolecular condensates to underscore the importance of protein backbones in phase separation. We further focus on recent advances in experimental techniques and molecular dynamics simulation methods for probing and exploring the roles of backbone interactions and secondary structures in biomolecular phase separation involving IDPs.

Incorporating physics to overcome data scarcity in predictive modeling of protein function: A case study of BK channels.

Machine learning has played transformative roles in numerous chemical and biophysical problems such as protein folding where large amount of data exists. Nonetheless, many important problems remain challenging for data-driven machine learning approaches due to the limitation of data scarcity. One approach to overcome data scarcity is to incorporate physical principles such as through molecular modeling and simulation. Here, we focus on the big potassium (BK) channels that play important roles in cardiovascular and neural systems. Many mutants of BK channel are associated with various neurological and cardiovascular diseases, but the molecular effects are unknown. The voltage gating properties of BK channels have been characterized for 473 site-specific mutations experimentally over the last three decades; yet, these functional data by themselves remain far too sparse to derive a predictive model of BK channel voltage gating. Using physics-based modeling, we quantify the energetic effects of all single mutations on both open and closed states of the channel. Together with dynamic properties derived from atomistic simulations, these physical descriptors allow the training of random forest models that could reproduce unseen experimentally measured shifts in gating voltage, ∆V1/2, with a RMSE ~ 32 mV and correlation coefficient of R ~ 0.7. Importantly, the model appears capable of uncovering nontrivial physical principles underlying the gating of the channel, including a central role of hydrophobic gating. The model was further evaluated using four novel mutations of L235 and V236 on the S5 helix, mutations of which are predicted to have opposing effects on V1/2 and suggest a key role of S5 in mediating voltage sensor-pore coupling. The measured ∆V1/2 agree quantitatively with prediction for all four mutations, with a high correlation of R = 0.92 and RMSE = 18 mV. Therefore, the model can capture nontrivial voltage gating properties in regions where few mutations are known. The success of predictive modeling of BK voltage gating demonstrates the potential of combining physics and statistical learning for overcoming data scarcity in nontrivial protein function prediction.

Assessing risk of recurrent small bowel obstruction after non-operative management in patients with history of intra-abdominal surgery: a population-based comprehensive analysis in Taiwan.

BACKGROUND: Despite a significant 30% ten-year readmission rate for SBO patients, investigations into recurrent risk factors after non-operative management are scarce. The study aims to generate a risk factor scoring system, the 'Small Bowel Obstruction Recurrence Score' (SBORS), predicting 6-month recurrence of small bowel obstruction (SBO) after successful non-surgical management in patients who have history of intra-abdominal surgery. METHODS: We analyzed data from patients aged ≥ 18 with a history of intra-abdominal surgery and diagnosed with SBO (ICD-9 code: 560, 568) and were successful treated non-surgically between 2004 and 2008. Participants were divided into model-derivation (80%) and validation (20%) group. RESULTS: We analyzed 23,901 patients and developed the SBORS based on factors including the length of hospital stay > 4 days, previous operations > once, hemiplegia, extra-abdominal and intra-abdominal malignancy, esophagogastric surgery and intestino-colonic surgery. Scores > 2 indicated higher rates and risks of recurrence within 6 months (12.96% vs. 7.27%, OR 1.898, p < 0.001 in model-derivation group, 12.60% vs. 7.05%, OR 1.901, p < 0.001 in validation group) with a significantly increased risk of mortality and operative events for recurrent episodes. The SBORS model demonstrated good calibration and acceptable discrimination, with an area under curve values of 0.607 and 0.599 for the score generation and validation group, respectively. CONCLUSIONS: We established the effective 'SBORS' to predict 6-month SBO recurrence risk in patients who have history of intra-abdominal surgery and have been successfully managed non-surgically for the initial obstruction event. Those with scores > 2 face higher recurrence rates and operative risks after successful non-surgical management.

Inner pore hydration free energy controls the activation of big potassium channels.

Hydrophobic gating is an emerging mechanism in regulation of protein ion channels where the pore remains physically open but becomes dewetted to block ion permeation. Atomistic molecular dynamics simulations have played a crucial role in understanding hydrophobic gating by providing the molecular details to complement mutagenesis and structural studies. However, existing studies rely on direct simulations and do not quantitatively describe how the sequence and structural changes may control the delicate liquid-vapor equilibrium of confined water in the pore of the channel protein. To address this limitation, we explore two enhanced sampling methods, namely metadynamics and umbrella sampling, to derive free-energy profiles of pore hydration in both the closed and open states of big potassium (BK) channels, which are important in cardiovascular and neural systems. It was found that metadynamics required substantially longer sampling times and struggled to generate stably converged free-energy profiles due to the slow dynamics of cooperative pore water diffusion even in the barrierless limit. Using umbrella sampling, well-converged free-energy profiles can be readily generated for the wild-type BK channels as well as three mutants with pore-lining mutations experimentally known to dramatically perturb the channel gating voltage. The results show that the free energy of pore hydration faithfully reports the gating voltage of the channel, providing further support for hydrophobic gating in BK channels. Free-energy analysis of pore hydration should provide a powerful approach for quantitative studies of how protein sequence, structure, solution conditions, and/or drug binding may modulate hydrophobic gating in ion channels.

Predicting in-hospital mortality risk for perforated peptic ulcer surgery: the PPUMS scoring system and the benefit of laparoscopic surgery: a population-based study.

BACKGROUND: The major treatment for perforated peptic ulcers (PPU) is surgery. It remains unclear which patient may not get benefit from surgery due to comorbidity. This study aimed to generate a scoring system by predicting mortality for patients with PPU who received non-operative management (NOM) and surgical treatment. METHOD: We extracted the admission data of adult (≥ 18 years) patients with PPU disease from the NHIRD database. We randomly divided patients into 80% model derivation and 20% validation cohorts. Multivariate analysis with a logistic regression model was applied to generate the scoring system, PPUMS. We then apply the scoring system to the validation group. RESULT: The PPUMS score ranged from 0 to 8 points, composite with age (< 45: 0 points, 45-65: 1 point, 65-80: 2 points, > 80: 3 points), and five comorbidities (congestive heart failure, severe liver disease, renal disease, history of malignancy, and obesity: 1 point each). The areas under ROC curve were 0.785 and 0.787 in the derivation and validation groups. The in-hospital mortality rates in the derivation group were 0.6% (0 points), 3.4% (1 point), 9.0% (2 points), 19.0% (3 points), 30.2% (4 points), and 45.9% when PPUMS > 4 point. Patients with PPUMS > 4 had a similar in-hospital mortality risk between the surgery group [laparotomy: odds ratio (OR) = 0.729, p = 0.320, laparoscopy: OR = 0.772, p = 0.697] and the non-surgery group. We identified similar results in the validation group. CONCLUSION: PPUMS scoring system effectively predicts in-hospital mortality for perforated peptic ulcer patients. It factors in age and specific comorbidities is highly predictive and well-calibrated with a reliable AUC of 0.785-0.787. Surgery, no matter laparotomy or laparoscope, significantly reduced mortality for scores < = 4. However, patients with a score > 4 did not show this difference, calling for tailored approaches to treatment based on risk assessment. Further prospective validation is suggested.

Amyloid assembly is dominated by misregistered kinetic traps on an unbiased energy landscape.

Atomistic description of protein fibril formation has been elusive due to the complexity and long time scales of the conformational search. Here, we develop a multiscale approach combining numerous atomistic simulations in explicit solvent to construct Markov State Models (MSMs) of fibril growth. The search for the in-register fully bound fibril state is modeled as a random walk on a rugged two-dimensional energy landscape defined by ß-sheet alignment and hydrogen-bonding states, whereas transitions involving states without hydrogen bonds are derived from kinetic clustering. The reversible association/dissociation of an incoming peptide and overall growth kinetics are then computed from MSM simulations. This approach is applied to derive a parameter-free, comprehensive description of fibril elongation of Aß16-22 and how it is modulated by phenylalanine-to-cyclohexylalanine (CHA) mutations. The trajectories show an aggregation mechanism in which the peptide spends most of its time trapped in misregistered ß-sheet states connected by weakly bound states twith short lifetimes. Our results recapitulate the experimental observation that mutants CHA19 and CHA1920 accelerate fibril elongation but have a relatively minor effect on the critical concentration for fibril growth. Importantly, the kinetic consequences of mutations arise from cumulative effects of perturbing the network of productive and nonproductive pathways of fibril growth. This is consistent with the expectation that nonfunctional states will not have evolved efficient folding pathways and, therefore, will require a random search of configuration space. This study highlights the importance of describing the complete energy landscape when studying the elongation mechanism and kinetics of protein fibrils.

Machine Learning Generation of Dynamic Protein Conformational Ensembles.

Machine learning has achieved remarkable success across a broad range of scientific and engineering disciplines, particularly its use for predicting native protein structures from sequence information alone. However, biomolecules are inherently dynamic, and there is a pressing need for accurate predictions of dynamic structural ensembles across multiple functional levels. These problems range from the relatively well-defined task of predicting conformational dynamics around the native state of a protein, which traditional molecular dynamics (MD) simulations are particularly adept at handling, to generating large-scale conformational transitions connecting distinct functional states of structured proteins or numerous marginally stable states within the dynamic ensembles of intrinsically disordered proteins. Machine learning has been increasingly applied to learn low-dimensional representations of protein conformational spaces, which can then be used to drive additional MD sampling or directly generate novel conformations. These methods promise to greatly reduce the computational cost of generating dynamic protein ensembles, compared to traditional MD simulations. In this review, we examine recent progress in machine learning approaches towards generative modeling of dynamic protein ensembles and emphasize the crucial importance of integrating advances in machine learning, structural data, and physical principles to achieve these ambitious goals.

Activation of TMEM16F by inner gate charged mutations and possible lipid/ion permeation mechanisms.

Transmembrane protein 16F (TMEM16F) is a ubiquitously expressed Ca2+-activated phospholipid scramblase that also functions as a largely non-selective ion channel. Though recent structural studies have revealed the closed and intermediate conformations of mammalian TMEM16F (mTMEM16F), the open and conductive state remains elusive. Instead, it has been proposed that an open hydrophilic pathway may not be required for lipid scrambling. We previously identified an inner activation gate, consisting of F518, Y563, and I612, and showed that charged mutations of the inner gate residues led to constitutively active mTMEM16F scrambling. Herein, atomistic simulations show that lysine substitution of F518 and Y563 can indeed lead to spontaneous opening of the permeation pore in the Ca2+-bound state of mTMEM16F. Dilation of the pore exposes hydrophilic patches in the upper pore region, greatly increases the pore hydration level, and enables lipid scrambling. The putative open state of mTMEM16F resembles the active state of fungal scramblases and is a meta-stable state for the wild-type protein in the Ca2+-bound state. Therefore, mTMEM16F may be capable of supporting the canonical in-groove scrambling mechanism in addition to the out-of-groove one. Further analysis reveals that the in-groove phospholipid and ion transduction pathways of mTMEM16F overlap from the intracellular side up to the inner gate but diverge from each other with different exits to the extracellular side of membrane.

Physics-based modeling provides predictive understanding of selectively promiscuous substrate binding by Hsp70 chaperones.

To help cells cope with protein misfolding and aggregation, Hsp70 molecular chaperones selectively bind a variety of sequences ("selective promiscuity"). Statistical analyses from substrate-derived peptide arrays reveal that DnaK, the E. coli Hsp70, binds to sequences containing three to five branched hydrophobic residues, although otherwise the specific amino acids can vary considerably. Several high-resolution structures of the substrate -binding domain (SBD) of DnaK bound to peptides reveal a highly conserved configuration of the bound substrate and further suggest that the substrate-binding cleft consists of five largely independent sites for interaction with five consecutive substrate residues. Importantly, both substrate backbone orientations (N- to C- and C- to N-) allow essentially the same backbone hydrogen-bonding and side-chain interactions with the chaperone. In order to rationalize these observations, we performed atomistic molecular dynamics simulations to sample the interactions of all 20 amino acid side chains in each of the five sites of the chaperone in the context of the conserved substrate backbone configurations. The resulting interaction energetics provide the basis set for deriving a predictive model that we call Paladin (Physics-based model of DnaK-Substrate Binding). Trained using available peptide array data, Paladin can distinguish binders and nonbinders of DnaK with accuracy comparable to existing predictors and further predicts the detailed configuration of the bound sequence. Tested using existing DnaK-peptide structures, Paladin correctly predicted the binding register in 10 out of 13 substrate sequences that bind in the N- to C- orientation, and the binding orientation in 16 out of 22 sequences. The physical basis of the Paladin model provides insight into the origins of how Hsp70s bind substrates with a balance of selectivity and promiscuity. The approach described here can be extended to other Hsp70s where extensive peptide array data is not available.

Toward Accurate Coarse-Grained Simulations of Disordered Proteins and Their Dynamic Interactions.

Intrinsically disordered proteins (IDPs) play crucial roles in cellular regulatory networks and are now recognized to often remain highly dynamic even in specific interactions and assemblies. Accurate description of these dynamic interactions is extremely challenging using atomistic simulations because of the prohibitive computational cost. Efficient coarse-grained approaches could offer an effective solution to overcome this bottleneck if they could provide an accurate description of key local and global properties of IDPs in both unbound and bound states. The recently developed hybrid-resolution (HyRes) protein model has been shown to be capable of providing a semiquantitative description of the secondary structure propensities of IDPs. Here, we show that greatly improved description of global structures and transient interactions can be achieved by introducing a solvent-accessible surface area-based implicit solvent term followed by reoptimization of effective interaction strengths. The new model, termed HyRes II, can semiquantitatively reproduce a wide range of local and global structural properties of a set of IDPs of various lengths and complexities. It can also distinguish the level of compaction between folded proteins and IDPs. In particular, applied to the disordered N-terminal transactivation domain (TAD) of tumor suppressor p53, HyRes II is able to recapitulate various nontrivial structural properties compared to experimental results, some of them to a level of accuracy that is almost comparable to results from atomistic explicit solvent simulations. Furthermore, we demonstrate that HyRes II can be used to simulate the dynamic interactions of TAD with the DNA-binding domain of p53, generating structural ensembles that are highly consistent with existing NMR data. We anticipate that HyRes II will provide an efficient and relatively reliable tool toward accurate coarse-grained simulations of dynamic protein interactions.

Previous unilateral inguinal hernia repair increase risk of new developed inguinal hernia: a nationwide Longitudinal Cohort Study in Asian male adult patients.

INTRODUCTION: Our study aims to identify that patients who received hernia repair previously did have higher risk of occurrence of newly developed inguinal hernia, named as a contralateral inguinal hernia (CIH), than patients who never received inguinal hernia surgery before. MATERIALS AND METHODS: We collected data from the National Health Insurance Research Database (NHIRD) of Taiwan retrospectively. In the study cohort, 64,089 Asian male adults who underwent primary unilateral inguinal hernia repair during 2003-2008 were included using ICD-9 diagnostic and surgical codes. Another 64,089 male adults without hernia repair history were included as control group via propensity score match. RESULTS: The median follow-up period is 93.53 months. After multivariate analysis, the risk of newly developed inguinal hernia in unilateral inguinal hernia (UIH) repair cohort was significantly higher (adjusted HR 6.364, 95% CI 6.012-6.737, P < 0.001) than the control group. In subgroup analysis, patients without mesh repair (adjusted HR 6.706, P < 0.001) and patients with mesh repair (adjusted HR 5.559, P < 0.001) both showed higher risk of developing newly developed inguinal hernia which needs repair. CONCLUSIONS: Asian men with UIH repair history had a higher risk of developing new inguinal hernia at the contralateral site, namely CIH, than the general population. The surgeon should inform the possibility of CIH after initial herniorrhaphy, therefore, monitoring the occurrence of CIH is necessary.

The utility of intraoperative endoscopy to assist novice surgeons in the detection of gastric stenosis during laparoscopic sleeve gastrectomy.

BACKGROUND: Laparoscopic sleeve gastrectomy (LSG) is a commonly performed bariatric surgery. Gastric stenosis and leaks are 2 major complications associated with LSG and revision surgery might be needed. Herein, we report our experience of intraoperative endoscopy (IOE) to evaluate stenosis and leaks during LSG. METHODS: LSG was performed by three surgeons. Patients who underwent LSG and IOE between January 2016 and March 2020 were enrolled and assigned to two groups: group 1 (1st-30th LSG case for each surgeon) and group 2 (> 30th LSG for each surgeon). Patients' anthropometric and biochemical data pre- and post-LSG, as well as IOE findings and follow-up esophagogastroduodenoscopy records were reviewed. RESULTS: In total, 352 patients were enrolled including 90 patients in group 1 and 262 patients in group 2. Three out of 352 patients (0.9%) were found to have stenosis by IOE, which was related to tightly gastropexy stitch or reinforcement stitch, all of which were in group 1. Stenosis was resolved after removal of the stitch during LSG. The incidence of gastric stenosis detected by IOE was 3.3% (3/90) and 0% (0/262) in group 1 and group 2, respectively (P = 0.003). No leakage was found in this study and no patient developed clinical or endoscopic stenosis after LSG. CONCLUSIONS: The existing evidence showed that IOE can help detect gastric stenosis during LSG, especially for novice surgeons, and the stenosis could be resolved during operation.

Accelerating atomistic simulations of proteins using multiscale enhanced sampling with independent tempering.

Efficient sampling of the conformational space is essential for quantitative simulations of proteins. The multiscale enhanced sampling (MSES) method accelerates atomistic sampling by coupling it to a coarse-grained (CG) simulation. Bias from coupling to the CG model is removed using Hamiltonian replica exchange, such that one could benefit simultaneously from the high accuracy of atomistic models and fast dynamics of CG ones. Here, we extend MSES to allow independent control of the effective temperatures of atomistic and CG simulations, by directly scaling the atomistic and CG Hamiltonians. The new algorithm, named MSES with independent tempering (MSES-IT), supports more sophisticated Hamiltonian and temperature replica exchange protocols to further improve the sampling efficiency. Using a small but nontrivial ß-hairpin, we show that setting the effective temperature of CG model in all conditions to its melting temperature maximizes structural transition rates at the CG level and promotes more efficient replica exchange and diffusion in the condition space. As the result, MSES-IT drive faster reversible transitions at the atomic level and leads to significant improvement in generating converged conformational ensembles compared to the original MSES scheme.

A case report of delayed lower intestinal bleeding after organophosphate poisoning.

BACKGROUND: Organophosphate poisoning is a serious issue and it results in significant casualties in developing countries. Since agriculture remains an important and necessary sector of human society and organophosphate are commonly used in agriculture, it is difficult to prevent organophosphate poisoning. Gastrointestinal bleeding is not a common but life threatening symptom of organophosphate poisoning. We report a rare case of gastrointestine bleeding due to organophosphate poisoning. CASE PRESENTATION: A 78-year-old woman presented to our hospital approximately 12 h after ingesting a mouthful of organophosphate and benzodiazepines in a suicide attempt. Six weeks after successful medical treatment for respiratory failure, she developed recurring melena. Colonoscopy and esophagogastroduodenoscopy findings were negative for ulcers or bleeding. Enteroscopy revealed severe circumferential ulcers with luminal narrowing 10 cm proximal to the ileocecal valve. The patient underwent a 100-cm ileum resection after failed medical treatment and recovered uneventfully. The resected terminal ileum demonstrated severe inflammation and a sharp transitional zone between the healthy and injured mucosa approximately 50 cm proximal to the ileocecal valve. Pathological examination revealed an injured mucosa with inflammatory cell infiltration and structural damage. This case highlights a rare event of OP poisoning with late-onset lower gastrointestinal bleeding, which prolonged the patient's recovery course and parenteral alimentation period. CONCLUSION: We report a rare case of a patient with organophosphate poisoning, with late-onset lower GI tract bleeding, which raised clinical awareness regarding the organophosphate poisoning that induce intestinal symptoms.

The long-term risk of cardiovascular events in patients following bariatric surgery compared to a non-surgical population with obesity and the general population: a comprehensive national cohort study.

PURPOSE: This study evaluated the differences in long-term cardiovascular events between obese patients who received bariatric surgery (BS), those who did not, and the general population (GP). METHODS: Between 2003 and 2008, patients with severe obesity, aged 18-55 years, were divided into the non-surgical (NS) and BS groups and were included in this retrospective study. We also extracted data of healthy civilians defined as the GP. The incidence of cardiovascular events, including myocardial infarction, intracranial hemorrhage, epidural hemorrhage, ischemic stroke, and transient ischemic attack, was defined as the primary end point. Patients were followed up either until the end of 2013, upon reaching the primary end point, or death. RESULTS: After propensity score matching, 1436 patients were included in both the BS and NS groups, and 4829 subjects were enrolled as the GP. Of these, 57 (3.9%), 10 (0.6%), and 30 (0.62%) subjects in the NS, BS, and GP, respectively, experienced cardiovascular events. Multivariate analysis revealed that patients with BS had a significantly lower risk of cardiovascular events (HR = 0.168; 95% CI 0.085-0.328; p < 0.001) than those in the NS group, but it was not significantly different in the BS group compared with the GP (HR = 1.202; 95% CI 0.585-2.471; p = 0.617). CONCLUSION: Long-term risk of cardiovascular events was significantly lower among patients who had BS, compared to the NS obese patients. Thus, the long-term cardiovascular risk between those who received BS and the GP had no significant difference, in a retrospective view.

Positional Isomers of a Non-Nucleoside Substrate Differentially Affect Myosin Function.

Molecular motors have evolved to transduce chemical energy from ATP into mechanical work to drive essential cellular processes, from muscle contraction to vesicular transport. Dysfunction of these motors is a root cause of many pathologies necessitating the need for intrinsic control over molecular motor function. Herein, we demonstrate that positional isomerism can be used as a simple and powerful tool to control the molecular motor of muscle, myosin. Using three isomers of a synthetic non-nucleoside triphosphate, we demonstrate that myosin's force- and motion-generating capacity can be dramatically altered at both the ensemble and single-molecule levels. By correlating our experimental results with computation, we show that each isomer exerts intrinsic control by affecting distinct steps in myosin's mechanochemical cycle. Our studies demonstrate that subtle variations in the structure of an abiotic energy source can be used to control the force and motility of myosin without altering myosin's structure.

Accelerating the Generalized Born with Molecular Volume and Solvent Accessible Surface Area Implicit Solvent Model Using Graphics Processing Units.

The generalized Born with molecular volume and solvent accessible surface area (GBMV2/SA) implicit solvent model provides an accurate description of molecular volume and has the potential to accurately describe the conformational equilibria of structured and disordered proteins. However, its broader application has been limited by the computational cost and poor scaling in parallel computing. Here, we report an efficient implementation of both the electrostatic and nonpolar components of GBMV2/SA on graphics processing unit (GPU) within the CHARMM/OpenMM module. The GPU-GBMV2/SA is numerically equivalent to the original CPU-GBMV2/SA. The GPU acceleration offers ~60- to 70-fold speedup on a single NVIDIA TITAN X (Pascal) graphics card for molecular dynamic simulations of both folded and unstructured proteins of various sizes. The current implementation can be further optimized to achieve even greater acceleration with minimal reduction on the numerical accuracy. The successful development of GPU-GBMV2/SA greatly facilitates its application to biomolecular simulations and paves the way for further development of the implicit solvent methodology. © 2019 Wiley Periodicals, Inc.

Hydrophobic dewetting in gating and regulation of transmembrane protein ion channels.

Water is at the heart of almost all biological phenomena, without which no life that we know of would have been possible. It is a misleadingly complex liquid that exists in near coexistence with the vapor phase under ambient conditions. Confinement within a hydrophobic cavity can tip this balance enough to drive a cooperative dewetting transition. For a nanometer-scale pore, the dewetting transition leads to a stable dry state that is physically open but impermeable to ions. This phenomenon is often referred to as hydrophobic gating. Numerous transmembrane protein ion channels have now been observed to utilize hydrophobic gating in their activation and regulation. Here, we review recent theoretical, simulation, and experimental studies that together have started to establish the principles of hydrophobic gating and discuss how channels of various sizes, topologies, and biological functions can utilize these principles to control the thermodynamic properties of water within their interior pores for gating and regulation. Exciting opportunities remain in multiple areas, particularly on direct experimental detection of hydrophobic dewetting in biological channels and on understanding how the cell may control the hydrophobic gating in regulation of ion channels.