OpenAWSEM with Open3SPN2: A fast,flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.     

Structure of a Ca2 +/CaM:Kv7.4 (KCNQ4) B-Helix Complex Provides Insight into M Current Modulation

Calmodulin (CaM) is an important regulator of Kv7.x (KCNQx) voltage-gated potassium channels. Channels from this family produce neuronal M currents and cardiac and auditory I_KS currents and harbor mutations that cause arrhythmias, epilepsy, and deafness. Despite extensive functional characterization, biochemical and structural details of the interaction between CaM and the channel have remained elusive. Here, we show that both apo-CaM and Ca^2 +/CaM bind to the C-terminal tail of the neuronal channel Kv7.4 (KCNQ4), which is involved in both hearing and mechanosensation. Interactions between apo-CaM and the Kv7.4 tail involve two C-terminal tail segments, known as the A and B segments, whereas the interaction between Ca^2 +/CaM and the Kv7.4 C-terminal tail requires only the B segment. Biochemical studies show that the calcium dependence of the CaM:B segment interaction is conserved in all Kv7 subtypes. X-ray crystallographic determination of the structure of the Ca^2 +/CaM:Kv7.4 B segment complex shows that Ca^2 +/CaM wraps around the Kv7.4 B segment, which forms an α-helix, in an antiparallel orientation that embodies a variation of the classic 1-14 Ca^2 +/CaM interaction motif. Taken together with the context of prior studies, our data suggest a model for modulation of neuronal Kv7 channels involving a calcium-dependent conformational switch from an apo-CaM form that bridges the A and B segments to a Ca^2 +/CaM form bound to the B-helix. The structure presented here also provides a context for a number of disease-causing mutations and for further dissection of the mechanisms by which CaM controls Kv7 function.     

Understanding Confounding Effects in Linguistic Coordination: An Information-Theoretic Approach

We suggest an information-theoretic approach for measuring stylistic coordination in dialogues. The proposed measure has a simple predictive interpretation and can account for various confounding factors through proper conditioning. We revisit some of the previous studies that reported strong signatures of stylistic accommodation, and find that a significant part of the observed coordination can be attributed to a simple confounding effect—length coordination. Specifically, longer utterances tend to be followed by longer responses, which gives rise to spurious correlations in the other stylistic features. We propose a test to distinguish correlations in length due to contextual factors (topic of conversation, user verbosity, etc.) and turn-by-turn coordination. We also suggest a test to identify whether stylistic coordination persists even after accounting for length coordination and contextual factors.     

Synergistic Impact of Solvent and Polymer Additives on the Film Formation of Small Molecule Blend Films for Bulk Heterojunction Solar Cells

The addition of polystyrene (PS), a typical insulator, is empirically shown to increase the power conversion efficiencies (PCEs) of a solution‐deposited bulk heterojunction (BHJ) molecular blend film used in solar cell fabrication: p‐DTS(FBTTh₂)₂/PC₇₁BM. The performance is further improved by small quantities of diiodooctane (DIO), an established solvent additive. In this study, how the addition of PS and DIO affects the film formation of this bulk heterojunction blend film are probed via in situ monitoring of absorbance, thickness, and crystallinity. PS and DIO additives are shown to promote donor crystallite formation on different time scales and through different mechanisms. PS‐containing films retain chlorobenzene solvent, extending evaporation time and promoting phase separation earlier in the casting process. This extended time is insufficient to attain the morphology for optimal PCE results before the film sets. Here is where the presence of DIO comes into play: its low vapor pressure further extends the time scale of film evolution and allows for crystalline rearrangement of the donor phase long after casting, ultimately leading to the best BHJ organization.     

16.1% Efficient Hysteresis‐Free Mesostructured Perovskite Solar Cells Based on Synergistically Improved ZnO Nanorod Arrays

Significant efficiency improvements are reported in mesoscopic perovskite solar cells based on the development of a low‐temperature solution‐processed ZnO nanorod (NR) array exhibiting higher NR aspect ratio, enhanced electron density, and substantially reduced work function than conventional ZnO NRs. These features synergistically result in hysteresis‐free, scan‐independent, and stabilized devices with an efficiency of 16.1%. Electron‐rich, nitrogen‐doped ZnO (N:ZnO) NR‐based electron transporting materials (ETMs) with enhanced electron mobility produced using ammonium acetate show consistently higher efficiencies by one to three power points than undoped ZnO NRs. Additionally, the preferential electrostatic interaction between the ­nonpolar facets of N:ZnO and the conjugated polyelectrolyte polyethylenimine (PEI) has been relied on to promote the hydrothermal growth of high aspect ratio NR arrays and substantially improve the infiltration of the perovskite light absorber into the ETM. Using the same interactions, a conformal PEI coating on the electron‐rich high aspect ratio N:ZnO NR arrays is ­successfully applied, resulting in a favorable work function shift and altogether leading to the significant boost in efficiency from <10% up to >16%. These results largely surpass the state‐of‐the‐art PCE of ZnO‐based perovskite solar cells and highlight the benefits of synergistically combining mesoscale control with doping and surface modification.     

Depositional processes of ribbon carbonates in Middle Cambrian of Iran (Deh-Sufiyan Formation,Central Alborz)

A variety of ribbon carbonates of the Deh-Sufiyan Formation (Middle Cambrian) in Central Alborz Range of northern Iran are studied to provide facies characterization and paleoenvironmental interpretation of ribbon carbonates on shallow-marine carbonate platforms. Seven types of ribbon carbonates are divided based mainly on sedimentary structures, ichnofossils, and bed geometry, which represent deposition during different phases of storm-induced processes. The different features of the storm deposits in ribbon carbonates such as hummocky and swaley cross-stratification, planar lamination, and combined-flow-ripple cross-stratification were formed by combined flows. Identification and interpretation of ichnological signatures and the spatial arrangement of succession of sedimentary structures are used to further refine sedimentary interpretations of parameters such as wave energy, substrate properties, variability in sedimentation rates, and proximality-distality trends of a wave-dominated marine ramp sequence. Successions from individual storm events reflect deposition during increasing combined oscillatory and unidirectional flow succeeded by the waning stages. The study provides depositional processes and models of various ribbon carbonates that may be useful for facies interpretation of ribbon rocks elsewhere.     

Sequence-Level Analysis of the Major European Huntington Disease Haplotype

Huntington disease (HD) reflects the dominant consequences of a CAG-repeat expansion in HTT. Analysis of common SNP-based haplotypes has revealed that most European HD subjects have distinguishable HTT haplotypes on their normal and disease chromosomes and that ∼50% of the latter share the same major HD haplotype. We reasoned that sequence-level investigation of this founder haplotype could provide significant insights into the history of HD and valuable information for gene-targeting approaches. Consequently, we performed whole-genome sequencing of HD and control subjects from four independent families in whom the major European HD haplotype segregates with the disease. Analysis of the full-sequence-based HTT haplotype indicated that these four families share a common ancestor sufficiently distant to have permitted the accumulation of family-specific variants. Confirmation of new CAG-expansion mutations on this haplotype suggests that unlike most founders of human disease, the common ancestor of HD-affected families with the major haplotype most likely did not have HD. Further, availability of the full sequence data validated the use of SNP imputation to predict the optimal variants for capturing heterozygosity in personalized allele-specific gene-silencing approaches. As few as ten SNPs are capable of revealing heterozygosity in more than 97% of European HD subjects. Extension of allele-specific silencing strategies to the few remaining homozygous individuals is likely to be achievable through additional known SNPs and discovery of private variants by complete sequencing of HTT. These data suggest that the current development of gene-based targeting for HD could be extended to personalized allele-specific approaches in essentially all HD individuals of European ancestry.