C1 domains are independently folded modules that are responsible for targeting their parent proteins to lipid membranes containing diacylglycerol (DAG), a ubiquitous second messenger. The DAG-binding affinities of C1 domains determine the threshold concentration of DAG required for the propagation of the signaling response and the selectivity of this response among the DAG receptors in the cell. The structural information currently available for C1 domains offers little insight into the molecular basis of their differential DAG-binding affinities. In this work, we characterized the C1B domain of Protein Kinase Cα (C1Bα) and its diagnostic mutant, Y123W, using solution NMR methods and molecular dynamics simulations. The mutation did not perturb the C1Bα structure or sub-nanosecond dynamics of the protein backbone, but resulted in a >100-fold increase of DAG binding affinity and substantial change in μs-timescale conformational dynamics, as quantified by NMR rotating-frame relaxation-dispersion methods. The differences in the conformational exchange behavior between the wild-type and Y123W C1Bα were localized to the hinge regions of ligand-binding loops. Molecular dynamics simulations provided insight into the identity of the exchanging conformers and revealed the significance of a particular residue, Gln128, in modulating the geometry of the ligand-binding site. Taken together with the results of binding studies, our findings suggest that the conformational dynamics and preferential partitioning of the tryptophan sidechain into the water-lipid interface are important factors that modulate the DAG-binding properties of C1 domains.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has motivated a widespread effort to understand its epidemiology and pathogenic mechanisms. Modern high-throughput sequencing technology has led to the deposition of vast numbers of SARS-CoV-2 genome sequences in curated repositories, which have been useful in mapping the spread of the virus around the globe. They also provide a unique opportunity to observe virus evolution in real time. Here, we evaluate two sets of SARS-CoV-2 genomic sequences to identify emerging variants within structured cis-regulatory elements of the SARS-CoV-2 genome. Overall, 20 variants are present at a minor allele frequency of at least 0.5%. Several enhance the stability of Stem Loop 1 in the 5ʹ untranslated region (UTR), including a group of co-occurring variants that extend its length. One appears to modulate the stability of the frameshifting pseudoknot between ORF1a and ORF1b, and another perturbs a bi-ss molecular switch in the 3ʹUTR. Finally, 5 variants destabilize structured elements within the 3ʹUTR hypervariable region, including the S2M (stem loop 2 m) selfish genetic element, raising questions as to the functional relevance of these structures in viral replication. Two of the most abundant variants appear to be caused by RNA editing, suggesting host-viral defense contributes to SARS-CoV-2 genome heterogeneity. Our analysis has implications for the development of therapeutics that target viral cis-regulatory RNA structures or sequences.

Profilin-1 (PFN1) plays important roles in modulating actin dynamics through binding both monomeric actin and proteins enriched with polyproline motifs. Mutations in PFN1 have been linked to the neurodegenerative disease amyotrophic lateral sclerosis (ALS). However, whether ALS-linked mutations affect PFN1 function has remained unclear. To address this question, we employed an unbiased proteomics analysis in mammalian cells to identify proteins that differentially interact with mutant and wild-type (WT) PFN1. These studies uncovered differential binding between two ALS-linked PFN1 variants, G118V and M114T, and select formin proteins. Furthermore, both variants augmented formin-mediated actin assembly relative to PFN1 WT. Molecular dynamics simulations revealed mutation-induced changes in the internal dynamic couplings within an alpha helix of PFN1 that directly contacts both actin and polyproline, as well as structural fluctuations within the actin- and polyproline-binding regions of PFN1. These data indicate that ALS-PFN1 variants have the potential for heightened flexibility in the context of the ternary actin–PFN1–polyproline complex during actin assembly. Conversely, PFN1 C71G was more severely destabilized than the other PFN1 variants, resulting in reduced protein expression in both transfected and ALS patient lymphoblast cell lines. Moreover, this variant exhibited loss-of-function phenotypes in the context of actin assembly. Perturbations in actin dynamics and assembly can therefore result from ALS-linked mutations in PFN1. However, ALS-PFN1 variants may dysregulate actin polymerization through different mechanisms that depend upon the solubility and stability of the mutant protein.

Allosteric regulation is an essential function of many proteins that control a variety of different processes such as catalysis, signal transduction, and gene regulation. Structural rearrangements have historically been considered the main means of communication between different parts of a protein. Recent studies have highlighted the importance, however, of changes in protein flexibility as an effective way to mediate allosteric communication across a protein. Scapharca dimeric hemoglobin (HbI) is the simplest possible allosteric system, with cooperative ligand binding between two identical subunits. Thermodynamic equilibrium studies of the binding of oxygen to HbI have shown that cooperativity is an entropically driven effect. The change in entropy of the system observed upon ligand binding may arise from changes in the protein, the ligand, or the water of the system. The goal of this study is to determine the contribution of the change in entropy of the protein backbone to HbI cooperative binding. Molecular dynamics simulations and nuclear magnetic resonance relaxation techniques have revealed that the fast internal motions of HbI contribute to the cooperative binding to carbon monoxide in two ways: (1) by contributing favorably to the free energy of the system and (2) by participating in the cooperative mechanism at the HbI subunit interface. The internal dynamics of the weakly cooperative HbI mutant, F97Y, were also investigated with the same methods. The changes in backbone NH dynamics observed for F97Y HbI upon ligand binding are not as large as for the wild type, in agreement with the reduced cooperativity observed for this mutant. The results of this study indicate that interface flexibility and backbone conformational entropy of HbI participate in and are important for the cooperative mechanism of carbon monoxide binding.

The ability to determine the effect of charge changes on the free energy is necessary for fundamental studies of the electrostatic contribution to protein binding and stability. Currently, calculations of differences in free energy for charge mutations carried out in systems with periodic boundary conditions must include an approximate self-energy correction that can be on the same order of magnitude as the calculated free energy change. Here, a new method for accurately calculating the free energy change associated with any alchemical mutation, regardless of charge, is presented. In this method, paired mutations of opposite charge exactly cancel the self-energy term because of its quadratic charge dependence. Since the self-energy term implicitly cancels within the method, a correction never needs to be applied, and the statistical uncertainty associated is thereby removed. An implementation procedure is described and applied to the free energy of ionic hydration and a charged amino acid mutation.