Search

Seriously ill patients infected with severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) and hospitalized in intensive care units (ICUs) are commonly given a combination of drugs, a process known as multi‐drug treatment. After extracting data on drug–drug interactions with clinical relevance from available online platforms, we hypothesize that an overall interaction map can be generated for all drugs administered. Furthermore, by combining this approach with simulations of cellular biochemical pathways, we may be able to explain the general clinical outcome. Finally, we postulate that by applying this strategy retrospectively to a cohort of patients hospitalized in ICU, a prediction of the timing of developing acute kidney injury (AKI) could be made. Whether or not this approach can be extended to other diseases is uncertain. Still, we believe it represents a valuable pharmacological insight to help improve clinical outcomes for severely ill patients. [ABSTRACT FROM AUTHOR]

Inside cells, various biological systems work cooperatively for homeostasis and self-replication. These systems do not work independently as they compete for shared elements like ATP and NADH. However, it has been believed that such competition is not a problem in codependent biological systems such as the energy-supplying glycolysis and the energy-consuming translation system. In this study, we biochemically reconstituted the coupling system of glycolysis and translation using purified elements and found that the competition for ATP between glycolysis and protein synthesis interferes with their coupling. Both experiments and simulations revealed that this interference is derived from a metabolic tug-of-war between glycolysis and translation based on their reaction rates, which changes the threshold of the initial substrate concentration for the success coupling. By the metabolic tug-of-war, translation energized by strong glycolysis is facilitated by an exogenous ATPase, which normally inhibits translation. These findings provide chemical insights into the mechanism of competition among biological systems in living cells and provide a framework for the construction of synthetic metabolism in vitro.

Guanylate binding proteins (GBPs) belong to the dynamin-related superfamily and exhibit various functions in the fight against infections. The functions of the human guanylate binding protein 1 (hGBP1) are tightly coupled to GTP hydrolysis and dimerization. Despite known crystal structures of the hGBP1 monomer and GTPase domain dimer, little is known about the dynamics of hGBP1. To gain a mechanistic understanding of hGBP1, we performed sub-millisecond multi-resolution molecular dynamics simulations of both the hGBP1 monomer and dimer. We found that hGBP1 is a highly flexible protein that undergoes a hinge motion similar to the movements observed for other dynamin-like proteins. Another large-scale motion was observed for the C-terminal helix α13, providing a molecular view for the α13–α13 distances previously reported for the hGBP1 dimer. Most of the loops of the GTPase domain were found to be flexible, revealing why GTP binding is needed for hGBP1 dimerization to occur. [ABSTRACT FROM AUTHOR]

The ability of DENV2 to display different morphologies (hence different antigenic properties) complicates vaccine and therapeutics development. Previous studies showed most strains of laboratory adapted DENV2 particles changed from smooth to “bumpy” surfaced morphology when the temperature is switched from 29°C at 37°C. Here we identified five envelope (E) protein residues different between two alternative passage history DENV2 NGC strains exhibiting smooth or bumpy surface morphologies. Several mutations performed on the smooth DENV2 infectious clone destabilized the surface, as observed by cryoEM. Molecular dynamics simulations demonstrated how chemically subtle substitution at various positions destabilized dimeric interactions between E proteins. In contrast, three out of four DENV2 clinical isolates showed a smooth surface morphology at 37°C, and only at high fever temperature (40°C) did they become “bumpy”. These results imply vaccines should contain particles representing both morphologies. For prophylactic and therapeutic treatments, this study also informs on which types of antibodies should be used at different stages of an infection, i.e., those that bind to monomeric E proteins on the bumpy surface or across multiple E proteins on the smooth surfaced virus. [ABSTRACT FROM AUTHOR]

Choline geranate (also described as Choline And GEranic acid, or CAGE) has been developed as a novel biocompatible antiseptic material capable of penetrating skin and aiding the transdermal delivery of co-administered antibiotics. The antibacterial properties of CAGE were analyzed against 24 and 72 hour old biofilms of 11 clinically isolated ESKAPE pathogens (defined as Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumanii, Pseudomonas aeruginosa, and Enterobacter sp, respectively), including multidrug resistant (MDR) isolates. CAGE was observed to eradicate in vitro biofilms at concentrations as low as 3.56 mM (0.156% v:v) in as little as 2 hours, which represents both an improved potency and rate of biofilm eradication relative to that reported for most common standard-of-care topical antiseptics in current use. In vitro time-kill studies on 24 hour old Staphylococcus aureus biofilms indicate that CAGE exerts its antibacterial effect upon contact and a 0.1% v:v solution reduced biofilm viability by over three orders of magnitude (a 3log10 reduction) in 15 minutes. Furthermore, disruption of the protective layer of exopolymeric substances in mature biofilms of Staphylococcus aureus by CAGE (0.1% v:v) was observed in 120 minutes. Insight into the mechanism of action of CAGE was provided with molecular modeling studies alongside in vitro antibiofilm assays. The geranate ion and geranic acid components of CAGE are predicted to act in concert to integrate into bacterial membranes, affect membrane thinning and perturb membrane homeostasis. Taken together, our results show that CAGE demonstrates all properties required of an effective topical antiseptic and the data also provides insight into how its observed antibiofilm properties may manifest. [ABSTRACT FROM AUTHOR]

The voltage-gated potassium channel KCNQ1 (KV7.1) assembles with the KCNE1 accessory protein to generate the slow delayed rectifier current, IKS, which is critical for membrane repolarization as part of the cardiac action potential. Loss-of-function (LOF) mutations in KCNQ1 are the most common cause of congenital long QT syndrome (LQTS), type 1 LQTS, an inherited genetic predisposition to cardiac arrhythmia and sudden cardiac death. A detailed structural understanding of KCNQ1 is needed to elucidate the molecular basis for KCNQ1 LOF in disease and to enable structure-guided design of new anti-arrhythmic drugs. In this work, advanced structural models of human KCNQ1 in the resting/closed and activated/open states were developed by Rosetta homology modeling guided by newly available experimentally-based templates: X. leavis KCNQ1 and various resting voltage sensor structures. Using molecular dynamics (MD) simulations, the capacity of the models to describe experimentally established channel properties including state-dependent voltage sensor gating charge interactions and pore conformations, PIP2 binding sites, and voltage sensor–pore domain interactions were validated. Rosetta energy calculations were applied to assess the utility of each model in interpreting mutation-evoked KCNQ1 dysfunction by predicting the change in protein thermodynamic stability for 50 experimentally characterized KCNQ1 variants with mutations located in the voltage-sensing domain. Energetic destabilization was successfully predicted for folding-defective KCNQ1 LOF mutants whereas wild type-like mutants exhibited no significant energetic frustrations, which supports growing evidence that mutation-induced protein destabilization is an especially common cause of KCNQ1 dysfunction. The new KCNQ1 Rosetta models provide helpful tools in the study of the structural basis for KCNQ1 function and can be used to generate hypotheses to explain KCNQ1 dysfunction. [ABSTRACT FROM AUTHOR]

Cellular DNA is daily exposed to several damaging agents causing a plethora of DNA lesions. As a first aid to restore DNA integrity, several enzymes got specialized in damage recognition and lesion removal during the process called base excision repair (BER). A large number of DNA damage types and several different readers of nucleic acids lesions during BER pathway as well as two sub-pathways were considered in the definition of a model using the Petri net framework. The intuitive graphical representation in combination with precise mathematical analysis methods are the strong advantages of the Petri net-based representation of biological processes and make Petri nets a promising approach for modeling and analysis of human BER. The reported results provide new information that will aid efforts to characterize in silico knockouts as well as help to predict the sensitivity of the cell with inactivated repair proteins to different types of DNA damage. The results can also help in identifying the by-passing pathways that may lead to lack of pronounced phenotypes associated with mutations in some of the proteins. This knowledge is very useful when DNA damage-inducing drugs are introduced for cancer therapy, and lack of DNA repair is desirable for tumor cell death. [ABSTRACT FROM AUTHOR]

T cells use their T-cell receptors (TCRs) to scan other cells for antigenic peptides presented by MHC molecules (pMHC). If a TCR encounters a pMHC, it can trigger a signalling pathway that could lead to the activation of the T cell and the initiation of an immune response. It is currently not clear how the binding of pMHC to the TCR initiates signalling within the T cell. One hypothesis is that conformational changes in the TCR lead to further downstream signalling. Here we investigate four different TCRs in their free state as well as in their pMHC bound state using large scale molecular simulations totalling 26 000 ns. We find that the dynamical features within TCRs differ significantly between unbound TCR and TCR/pMHC simulations. However, apart from expected results such as reduced solvent accessibility and flexibility of the interface residues, these features are not conserved among different TCR types. The presence of a pMHC alone is not sufficient to cause cross-TCR-conserved dynamical features within a TCR. Our results argue against models of TCR triggering involving conserved allosteric conformational changes. [ABSTRACT FROM AUTHOR]

Thirteen tubulin protofilaments, made of αβ-tubulin heterodimers, interact laterally to produce cytoskeletal microtubules. Microtubules exhibit the striking property of dynamic instability, manifested in their intermittent growth and shrinkage at both ends. This behavior is key to many cellular processes, such as cell division, migration, maintenance of cell shape, etc. Although assembly and disassembly of microtubules is known to be linked to hydrolysis of a guanosine triphosphate molecule in the pocket of β-tubulin, detailed mechanistic understanding of corresponding conformational changes is still lacking. Here we take advantage of the recent generation of in-microtubule structures of tubulin to examine the properties of protofilaments, which serve as important microtubule assembly and disassembly intermediates. We find that initially straight tubulin protofilaments, relax to similar non-radially curved and slightly twisted conformations. Our analysis further suggests that guanosine triphosphate hydrolysis primarily affects the flexibility and conformation of the inter-dimer interface, without a strong impact on the shape or flexibility of αβ-heterodimer. Inter-dimer interfaces are significantly more flexible compared to intra-dimer interfaces. We argue that such a difference in flexibility could be key for distinct stability of the plus and minus microtubule ends. The higher flexibility of the inter-dimer interface may have implications for development of pulling force by curving tubulin protofilaments during microtubule disassembly, a process of major importance for chromosome motions in mitosis. [ABSTRACT FROM AUTHOR]

African swine fever virus (ASFV) infection is fatal in domesticated pigs, with a mortality rate approaching 100%. This may result in economic losses and threats to food security. Currently, there are no approved vaccines or antiviral therapies for ASFV. Therefore, in this study, we evaluated congocidine congeners and a tris-benzimidazole as potential inhibitors of ASFV transcription using an in silico approach. We applied redocking of congocidine and docking of its congeners and a tris-benzimidazole to a receptor containing B-DNA with AT-motifs as a target to mimic conserved ASFV late gene promoters. Subsequently, the binding scores of DNA-ligand docked complexes were evaluated and their binding affinity was estimated. Molecular dynamics (MD) simulation was then used to assess ligand behavior within the minor groove. From our results, it is evident the less toxic congocidine congeners and tris-benzimidazole could dock to AT-rich regions significantly. Additionally, the predicted binding affinities had suitable values comparable to other experimentally determined minor groove binders, MD simulation of the docked DNA-ligand complexes and subsequent molecular trajectory visualization further showed that the ligands remained embedded in the minor groove during the time course of simulation, indicating that these ligands may have potential applications in abrogating ASFV transcription. [ABSTRACT FROM AUTHOR]

Electric fields can be a powerful tool to interact with enzymes or proteins, with an intriguing perspective to allow protein manipulation. Recently, researchers have focused the interest on intracellular enzyme modifications triggered by the application of nanosecond pulsed electric fields. These findings were also supported by theoretical predictions from molecular dynamics simulations focussing on significant variations in protein secondary structures. In this work, a theoretical study utilizing molecular dynamics simulations is proposed to explore effects of electric fields of high intensity and very short nanosecond duration applied to the superoxide dismutase (Cu/Zn-SOD or SOD-1), an important enzyme involved in the cellular antioxidant defence mechanism. The effects of 100-nanosecond pulsed electric fields, with intensities ranging from 108 to 7x108 V/m, on a single SOD1 enzyme are presented. We demonstrated that the lowest intensity of 108 V/m, although not inducing structural changes, can produce electrostatic modifications on the reaction centre of the enzyme, as apparent from the dipolar response and the electric field distribution of the protein active site. Electric pulses above 5x108 V/m produced a fast transition between the folded and a partially denatured state, as inferred by the secondary structures analysis. Finally, for the highest field intensity used (7x108 V/m), a not reversible transition toward an unfolded state was observed. [ABSTRACT FROM AUTHOR]

Starting from 1972, colchicine is known as the most useful drug for prevention of familial Mediterranean fever attacks. However, some patients do not respond to colchicine treatment, even taken in high doses. Despite the fact, that different hypotheses have been proposed, the molecular mechanisms of colchicine resistance are not completely clear. It is generally known, that colchicine binds β-tubulin and inhibits microtubules polymerization. The β-tubulin gene has SNPs, which lead to amino acid substitutions, and some of them are located in colchicine binding site (CBS). We have assumed, that this SNPs can affect tubulin-colchicine interaction and might be the reason for colchicine resistance. With this in mind, we modeled 7 amino acid substitutions in CBS, performed molecular dynamics simulations of tubulin-colchicine complex and calculated binding energies for every amino acid substitution. Thus, our study shows, that two amino acid substitutions in the β-tubulin, namely A248T and M257V, reduce binding energy for approximately 2-fold. Based on this, we assume, that these amino acid substitutions could be the reason for colchicine resistance. Thus, our study gives a new insight into colchicine resistance mechanism and provides information for designing colchicine alternatives, that could be effective for colchicine resistant patients. [ABSTRACT FROM AUTHOR]

HIV-1 Gag is a large multidomain poly-protein with flexible unstructured linkers connecting its globular subdomains. It is compact when in solution but assumes an extended conformation when assembled within the immature HIV-1 virion. Here, we use molecular dynamics (MD) simulations to quantitatively characterize the intra-domain interactions of HIV-1 Gag. We find that the matrix (MA) domain and the C-terminal subdomain CActd of the CA capsid domain can form a bound state. The bound state, which is held together primarily by interactions between complementary charged and polar residues, stabilizes the compact state of HIV-1 Gag. We calculate the depth of the attractive free energy potential between the MA/ CActd sites and find it to be about three times larger than the dimerization interaction between the CActd domains. Sequence analysis shows high conservation within the newly-found intra-Gag MA/CActd binding site, as well as its spatial proximity to other well known elements of Gag –such as CActd’s SP1 helix region, its inositol hexaphosphate (IP6) binding site and major homology region (MHR), as well as the MA trimerization site. Our results point to a high, but yet undetermined, functional significance of the intra-Gag binding site. Recent biophysical experiments that address the binding specificity of Gag are interpreted in the context of the MA/CActd bound state, suggesting an important role in selective packaging of genomic RNA by Gag. [ABSTRACT FROM AUTHOR]

The formation of a membrane-enveloped virus starts with the assembly of a curved layer of capsid proteins lining the interior of the plasma membrane (PM) of the host cell. This layer develops into a spherical shell (capsid) enveloped by a lipid-rich membrane. In many cases, the budding process stalls prior to the release of the virus. Recently, Brownian dynamics simulations of a coarse-grained model system reproduced protracted pausing and stalling, which suggests that the origin of pausing/stalling is to be found in the physics of the budding process. Here, we propose that the pausing/stalling observed in the simulations can be understood as a purely kinetic phenomenon associated with the neck geometry. A geometrical potential energy barrier develops during the budding that must be overcome by capsid proteins diffusing along the membrane prior to incorporation into the capsid. The barrier is generated by a conflict between the positive Gauss curvature of the assembling capsid and the negative Gauss curvature of the neck region. A continuum theory description is proposed and is compared with the Brownian simulations of the budding of enveloped viruses. [ABSTRACT FROM AUTHOR]

Leishmaniasis caused by obligate intracellular parasites of genus Leishmania is one of the most neglected tropical diseases threatening 350 million people worldwide. Protein kinases have drawn much attention as potential drug targets due to their important role in various cellular processes. In Leishmania sp. mitogen-activated protein kinase 4 is essential for the parasite survival because of its involvement in various regulatory, apoptotic and developmental pathways. The current study reveals the identification of natural inhibitors of L. donovani mitogen-activated protein kinase-4 (LdMPK4). We have performed in silico docking of 110 natural inhibitors of Leishmania parasite that have been reported earlier and identified two compounds Genistein (GEN) and Chrysin (CHY). The homology model of LdMPK4 was developed, followed by binding affinity studies, and pharmacokinetic properties of the inhibitors were calculated by maintaining ATP as a standard molecule. The modelled structure was deposited in the protein model database with PMDB ID: PM0080988. Molecular dynamic simulation of the enzyme-inhibitor complex along with the free energy calculations over 50 ns showed that GEN and CHY are more stable in their binding. These two molecules, GEN and CHY, can be considered as lead molecules for targeting LdMPK4 enzyme and could emerge as potential LdMPK4 inhibitors. [ABSTRACT FROM AUTHOR]

Human genome contains a group of more than a dozen similar genes with diverse biological functions including antiviral, antibacterial and angiogenesis activities. The characterized gene products of this group show significant sequence similarity and a common structural fold associated with binding and cleavage of ribonucleic acid (RNA) substrates. Therefore, these proteins have been categorized as members of human pancreatic-type ribonucleases (hRNases). hRNases differ in cell/tissue localization and display distinct substrate binding preferences and a wide range of ribonucleolytic catalytic efficiencies. Limited information is available about structural and dynamical properties that influence this diversity among these homologous RNases. Here, we use computer simulations to characterize substrate interactions, electrostatics and dynamical properties of hRNases 1–7 associated with binding to two nucleotide substrates (ACAC and AUAU). Results indicate that even with complete conservation of active-site catalytic triad associated with ribonucleolytic activity, these enzymes show significant differences in substrate interactions. Detailed characterization suggests that in addition to binding site electrostatic and van der Waals interactions, dynamics of distal regions may also play a role in binding. Another key insight is that a small difference in temperature of 300 K (used in experimental studies) and 310 K (physiological temperature) shows significant changes in enzyme-substrate interactions. [ABSTRACT FROM AUTHOR]

RAS is the founding member of a superfamily of GTPases and regulates signaling pathways involved in cellular growth control. While recent studies have shown that the activation state of RAS can be controlled by lysine ubiquitylation and acetylation, the existence of lysine methylation of the RAS superfamily GTPases remains unexplored. In contrast to acetylation, methylation does not alter the side chain charge and it has been challenging to deduce its impact on protein structure by conventional amino acid substitutions. Herein, we investigate lysine methylation on RAS and RAS-related GTPases. We developed GoMADScan (Go language-based Modification Associated Database Scanner), a new user-friendly application that scans and extracts posttranslationally modified peptides from databases. The GoMADScan search on PhosphoSitePlus databases identified methylation of conserved lysine residues in the core GTPase domain of RAS superfamily GTPases, including residues corresponding to RAS Lys-5, Lys-16, and Lys-117. To follow up on these observations, we immunoprecipitated endogenous RAS from HEK293T cells, conducted mass spectrometric analysis and found that RAS residues, Lys-5 and Lys-147, undergo dimethylation and monomethylation, respectively. Since mutations of Lys-5 have been found in cancers and RASopathies, we set up molecular dynamics (MD) simulations to assess the putative impact of Lys-5 dimethylation on RAS structure. Results from our MD analyses predict that dimethylation of Lys-5 does not significantly alter RAS conformation, suggesting that Lys-5 methylation may alter existing protein interactions or create a docking site to foster new interactions. Taken together, our findings uncover the existence of lysine methylation as a novel posttranslational modification associated with RAS and the RAS superfamily GTPases, and putative impact of Lys-5 dimethylation on RAS structure. [ABSTRACT FROM AUTHOR]

Mammals have two types of photoreceptors, rods and cones. While rods are exceptionally sensitive and mediate vision at very low illumination levels, cones operate in daylight and are responsible for the bulk of visual perception in most diurnal animals, including humans. Yet the mechanisms of phototransduction in cones is understudied, largely due to unavailability of pure cone outer segment (COS) preparations. Here we present a novel mathematical model of cone phototransduction that explicitly takes into account complex cone geometry and its multiple physical scales, faithfully reproduces features of the cone response, and is orders of magnitude more efficient than the standard 3D diffusion model. This is accomplished through the mathematical techniques of homogenization and concentrated capacity. The homogenized model is then computationally implemented by finite element method. This homogenized model permits one to analyze the effects of COS geometry on visual transduction and lends itself to performing large numbers of numerical trials, as required for parameter analysis and the stochasticity of rod and cone signal transduction. Agreement between the nonhomogenized, (i.e., standard 3D), and homogenized diffusion models is reported along with their simulation times and memory costs. Virtual expression of rod biochemistry on cone morphology is also presented for understanding some of the characteristic differences between rods and cones. These simulations evidence that 3D cone morphology and ion channel localization contribute to biphasic flash response, i.e undershoot. The 3D nonhomogenized and homogenized models are contrasted with more traditional and coarser well-stirred and 1D longitudinal diffusion models. The latter are single-scale and do not explicitly account for the multi-scale geometry of the COS, unlike the 3D homogenized model. We show that simpler models exaggerate the magnitude of the current suppression, yield accelerated time to peak, and do not predict the local concentration of cGMP at the ionic channels. [ABSTRACT FROM AUTHOR]

Recent evidence has raised in discussion the possibility that cannabidiol can act as a negative allosteric modulator of the cannabinoid type 1 receptor. Here we have used computational methods to study the modulation exerted by cannabidiol on the effects of delta-9-tetrahydrocannabinol in the cannabinoid receptor type 1 and the possibility of direct receptor blockade. We propose a putative allosteric binding site that is located in the N-terminal region of receptor, partially overlapping the orthosteric binding site. Molecular dynamics simulations reveled a coordinated movement involving the outward rotation of helixes 1 and 2 and subsequent expansion of the orthosteric binding site upon cannabidiol binding. Finally, changes in the cytoplasmic region and high helix 8 mobility were related to impaired receptor internalization. Together, these results offer a possible explanation to how cannabidiol can directly modulate effects of delta-9-tetrahydrocannabinol on the cannabinoid receptor type 1. [ABSTRACT FROM AUTHOR]

Nucleotides ratcheted through the biomolecular pores of nanopore sequencers generate raw picoamperage currents, which are segmented into step-current level signals representing the nucleotide sequence. These ‘squiggles’ are a noisy, distorted representation of the underlying true stepped current levels due to experimental and algorithmic factors. We were interested in developing a simulation model to support a white-box approach to identify common distortions, rather than relying on commonly used black box neural network techniques for basecalling nanopore signals. Dynamic time warped-space averaging (DTWA) techniques can generate a consensus from multiple noisy signals without introducing key feature distortions that occur with standard averaging. As a preprocessing tool, DTWA could provide cleaner and more accurate current signals for direct RNA or DNA analysis tools. However, DTWA approaches need modification to take advantage of the a-priori knowledge regarding a common, underlying gold-standard RNA / DNA sequence. Using experimental data, we derive a simulation model to provide known squiggle distortion signals to assist in validating the performance of analysis tools such as DTWA. Simulation models were evaluated by comparing mocked and experimental squiggle characteristics from one Enolase mRNA squiggle group produced by an Oxford MinION nanopore sequencer, and cross-validated using other Enolase, Sequin R1_71_1 and Sequin R2_55_3 mRNA studies. New techniques identified high inserted but low deleted base rates, generating consistent x1.7 squiggle event to base called ratios. Similar probability density and cumulative distribution functions, PDF and CDF, were found across all studies. Experimental PDFs were not the normal distributions expected if squiggle distortion arose from segmentation algorithm artefacts, or through individual nucleotides randomly interacting with individual nanopores. Matching experimental and mocked CDFs required the assumption that there are unique features associated with individual raw-current data streams. Z-normalized signal-to-noise ratios suggest intrinsic sensor limitations being responsible for half the gold standard and noisy squiggle DTW differences. [ABSTRACT FROM AUTHOR]

Studies on the Bin-Amphiphysin-Rvs (BAR) domain have advanced a fundamental understanding of how proteins deform membrane. We previously showed that a BAR domain in tandem with a Pleckstrin Homology (PH domain) underlies the assembly of ACAP1 (Arfgap with Coil-coil, Ankryin repeat, and PH domain I) into an unusual lattice structure that also uncovers a new paradigm for how a BAR protein deforms membrane. Here, we initially pursued computation-based refinement of the ACAP1 lattice to identify its critical protein contacts. Simulation studies then revealed how ACAP1, which dimerizes into a symmetrical structure in solution, is recruited asymmetrically to the membrane through dynamic behavior. We also pursued electron microscopy (EM)-based structural studies, which shed further insight into the dynamic nature of the ACAP1 lattice assembly. As ACAP1 is an unconventional BAR protein, our findings broaden the understanding of the mechanistic spectrum by which proteins assemble into higher-ordered structures to achieve membrane deformation. [ABSTRACT FROM AUTHOR]

Carrier proteins are four-helix bundles that covalently hold metabolites and secondary metabolites, such as fatty acids, polyketides and non-ribosomal peptides. These proteins mediate the production of many pharmaceutically important compounds including antibiotics and anticancer agents. Acyl carrier proteins (ACPs) can be found as part of a multi-domain polypeptide (Type I ACPs), or as part of a multiprotein complex (Type II). Here, the main focus is on ACP2 and ACP3, domains from the type I trans-AT polyketide synthase MmpA, which is a core component of the biosynthetic pathway of the antibiotic mupirocin. During molecular dynamics simulations of their apo, holo and acyl forms ACP2 and ACP3 both form a substrate-binding surface-groove. The substrates bound to this surface-groove have polar groups on their acyl chain exposed and forming hydrogen bonds with the solvent. Bulky hydrophobic residues in the GXDS motif common to all ACPs, and similar residues on helix III, appear to prohibit the formation of a deep tunnel in type I ACPs and type II ACPs from polyketide synthases. In contrast, the equivalent positions in ACPs from type II fatty acid synthases, which do form a deep solvent-excluded substrate-binding tunnel, have the small residue alanine. During simulation, ACP3 with mutations I61A L36A W44L forms a deep tunnel that can fully bury a saturated substrate in the core of the ACP, in contrast to the surface groove of the wild type ACP3. Similarly, in the ACP from E. coli fatty acid synthase, a type II ACP, mutations can change ligand binding from being inside a deep tunnel to being in a surface groove, thus demonstrating how changing a few residues can modify the possibilities for ligand binding. [ABSTRACT FROM AUTHOR]

Prions are unusual protein assemblies that propagate their conformationally-encoded information in absence of nucleic acids. The first prion identified, the scrapie isoform (PrPSc) of the cellular prion protein (PrPC), caused epidemic and epizootic episodes []. Most aggregates of other misfolding-prone proteins are amyloids, often arranged in a Parallel-In-Register-β-Sheet (PIRIBS) [] or β-solenoid conformations []. Similar folding models have also been proposed for PrPSc, although none of these have been confirmed experimentally. Recent cryo-electron microscopy (cryo-EM) and X-ray fiber-diffraction studies provided evidence that PrPSc is structured as a 4-rung β-solenoid (4RβS) [, ]. Here, we combined different experimental data and computational techniques to build the first physically-plausible, atomic resolution model of mouse PrPSc, based on the 4RβS architecture. The stability of this new PrPSc model, as assessed by Molecular Dynamics (MD) simulations, was found to be comparable to that of the prion forming domain of Het-s, a naturally-occurring β-solenoid. Importantly, the 4RβS arrangement allowed the first simulation of the sequence of events underlying PrPC conversion into PrPSc. This study provides the most updated, experimentally-driven and physically-coherent model of PrPSc, together with an unprecedented reconstruction of the mechanism underlying the self-catalytic propagation of prions. [ABSTRACT FROM AUTHOR]

The human serotonin transporter hSERT facilitates the reuptake of its endogenous substrate serotonin from the synaptic cleft into presynaptic neurons after signaling. Reuptake regulates the availability of this neurotransmitter and therefore hSERT plays an important role in balancing human mood conditions. In 2016, the first 3D structures of this membrane transporter were reported in an inhibitor-bound, outward-open conformation. These structures revealed valuable information about interactions of hSERT with antidepressant drugs. Nevertheless, the question remains how serotonin facilitates the specific conformational changes that open and close pathways from the synapse and to the cytoplasm as required for transport. Here, we present a serotonin-bound homology model of hSERT in an outward-occluded state, a key intermediate in the physiological cycle, in which the interactions with the substrate are likely to be optimal. Our approach uses two template structures and includes careful refinement and comprehensive computational validation. According to microsecond-long molecular dynamics simulations, this model exhibits interactions between the gating residues in the extracellular pathway, and these interactions differ from those in an outward-open conformation of hSERT bound to serotonin. Moreover, we predict several features of this state by monitoring the intracellular gating residues, the extent of hydration, and, most importantly, protein-ligand interactions in the central binding site. The results illustrate common and distinct characteristics of these two transporter states and provide a starting point for future investigations of the transport mechanism in hSERT. [ABSTRACT FROM AUTHOR]

γ-Aminobutyric acid (GABA) is the major inhibitory neurotransmitter in the vertebrate and invertebrate nervous system. GABAA receptors are activated by GABA and their agonists, and modulated by a wide variety of recognized drugs, including barbiturates, anesthetics, and benzodiazepines. The phenols propofol, thymol, chlorothymol, carvacrol and eugenol act as positive allosteric modulators on GABAA-R receptor. These GABAergic phenols interact with the lipid membrane, therefore, their anesthetic activity could be the combined result of their specific activity (with receptor proteins) as well as nonspecific interactions (with surrounding lipid molecules) modulating the supramolecular organization of the receptor environment. Therefore, we aimed to contribute to a description of the molecular events that occur at the membrane level as part of the mechanism of general anesthesia, using a molecular dynamic simulation approach. Equilibrium molecular dynamics simulations indicate that the presence of GABAergic phenols in a DPPC bilayer orders lipid acyl chains for carbons near the interface and their effect is not significant at the bilayer center. Phenols interacts with the polar interface of phospholipid bilayer, particularly forming hydrogen bonds with the glycerol and phosphate group. Also, potential of mean force calculations using umbrella sampling show that propofol partition is mainly enthalpic driven at the polar region and entropic driven at the hydrocarbon chains. Finally, potential of mean force indicates that propofol partition into a gel DPPC phase is not favorable. Our in silico results were positively contrasted with previous experimental data. [ABSTRACT FROM AUTHOR]

Profilin 1 (PFN1) protein plays key roles in neuronal growth and differentiation, membrane trafficking, and regulation of the actin cytoskeleton. Four natural variants of PFN1 were described as related to ALS, the most common adult-onset motor neuron disorder. However, the pathological mechanism of PFN1 in ALS is not yet completely understood. The goal of this work is to thoroughly analyze the effects of the ALS-related mutations on PFN1 structure and function using computational simulations. Here, PhD-SNP, PMUT, PolyPhen-2, SIFT, SNAP, SNPS&GO, SAAP, nsSNPAnalyzer, SNPeffect4.0 and I-Mutant2.0 were used to predict the functional and stability effects of PFN1 mutations. ConSurf was used for the evolutionary conservation analysis, and GROMACS was used to perform the MD simulations. The mutations C71G, M114T, and G118V, but not E117G, were predicted as deleterious by most of the functional prediction algorithms that were used. The stability prediction indicated that the ALS-related mutations could destabilize PFN1. The ConSurf analysis indicated that the mutation C71G, M114T, E117G, and G118V occur in highly conserved positions. The MD results indicated that the studied mutations could affect the PFN1 flexibility at the actin and PLP-binding domains, and consequently, their intermolecular interactions. It may be therefore related to the functional impairment of PFN1 upon C71G, M114T, E117G and G118V mutations, and their involvement in ALS development. We also developed a database, SNPMOL (), containing the results presented on this paper for biologists and clinicians to exploit PFN1 and its natural variants. [ABSTRACT FROM AUTHOR]

The cell cycle is the fundamental process of cell populations, it is regulated by environmental cues and by intracellular checkpoints. Cell cycle variability in clonal cell population is caused by stochastic processes such as random partitioning of cellular components to progeny cells at division and random interactions among biomolecules in cells. One of the important biological questions is how the dynamics at the cell cycle scale, which is related to family dependencies between the cell and its descendants, affects cell population behavior in the long-run. We address this question using a “mechanistic” model, built based on observations of single cells over several cell generations, and then extrapolated in time. We used cell pedigree observations of NIH 3T3 cells including FUCCI markers, to determine patterns of inheritance of cell-cycle phase durations and single-cell protein dynamics. Based on that information we developed a hybrid mathematical model, involving bifurcating autoregression to describe stochasticity of partitioning and inheritance of cell-cycle-phase times, and an ordinary differential equation system to capture single-cell protein dynamics. Long-term simulations, concordant with in vitro experiments, demonstrated the model reproduced the main features of our data and had homeostatic properties. Moreover, heterogeneity of cell cycle may have important consequences during population development. We discovered an effect similar to genetic drift, amplified by family relationships among cells. In consequence, the progeny of a single cell with a short cell cycle time had a high probability of eventually dominating the population, due to the heritability of cell-cycle phases. Patterns of epigenetic heritability in proliferating cells are important for understanding long-term trends of cell populations which are either required to provide the influx of maturing cells (such as hematopoietic stem cells) or which started proliferating uncontrollably (such as cancer cells). [ABSTRACT FROM AUTHOR]

Cadherins are homophilic cell-cell adhesion molecules whose aberrant expression has often been shown to correlate with different stages of tumor progression. In this work, we investigate the interaction of two peptidomimetic ligands with the extracellular portion of human E-cadherin using a combination of NMR and computational techniques. Both ligands have been previously developed as mimics of the tetrapeptide sequence Asp1-Trp2-Val3-Ile4 of the cadherin adhesion arm, and have been shown to inhibit E-cadherin-mediated adhesion in epithelial ovarian cancer cells with millimolar potency. To sample a set of possible interactions of these ligands with the E-cadherin extracellular portion, STD-NMR experiments in the presence of two slightly different constructs, the wild type E-cadherin-EC1-EC2 fragment and the truncated E-cadherin-(Val3)-EC1-EC2 fragment, were carried out at three temperatures. Depending on the protein construct, a different binding epitope of the ligand and also a different temperature effect on STD signals were observed, both suggesting an involvement of the Asp1-Trp2 protein sequence among all the possible binding events. To interpret the experimental results at the atomic level and to probe the role of the cadherin adhesion arm in the dynamic interaction with the peptidomimetic ligand, a computational protocol based on docking calculations and molecular dynamics simulations was applied. In agreement with NMR data, the simulations at different temperatures unveil high variability/dynamism in ligand-cadherin binding, thus explaining the differences in ligand binding epitopes. In particular, the modulation of the signals seems to be dependent on the protein flexibility, especially at the level of the adhesive arm, which appears to participate in the interaction with the ligand. Overall, these results will help the design of novel cadherin inhibitors that might prevent the swap dimer formation by targeting both the Trp2 binding pocket and the adhesive arm residues. [ABSTRACT FROM AUTHOR]

The calcium calmodulin-dependent protein kinase II (CaMKII) is a dodecameric holoenzyme important for encoding memory. Its activation, triggered by binding of calcium-calmodulin, persists autonomously after calmodulin dissociation. One (receiver) kinase captures and subsequently phosphorylates the regulatory domain peptide of a donor kinase forming a chained dimer as the first stage of autonomous activation. Protein dynamics simulations examined the conformational changes triggered by dimer formation and phosphorylation, aimed to provide a molecular rationale for human mutations that result in learning disabilities. Ensembles generated from X-ray crystal structures were characterized by network centrality and community analysis. Mutual information related collective motions to local fragment dynamics encoded with a structural alphabet. Implicit solvent tCONCOORD conformational ensembles revealed the dynamic architecture of Inactive kinase domains was co-opted in the activated dimer but the network hub shifted from the nucleotide binding cleft to the captured peptide. Explicit solvent molecular dynamics (MD) showed nucleotide and substrate binding determinants formed coupled nodes in long-range signal relays between regulatory peptides in the dimer. Strain in the extended captured peptide was balanced by reduced flexibility of the receiver kinase C-lobe core. The relays were organized around a hydrophobic patch between the captured peptide and a key binding helix. The human mutations aligned along the relays. Thus, these mutations could disrupt the allosteric network alternatively, or in addition, to altered binding affinities. Non-binding protein sectors distant from the binding sites mediated the allosteric signalling; providing possible targets for inhibitor design. Phosphorylation of the peptide modulated the dielectric of its binding pocket to strengthen the patch, non-binding sectors, domain interface and temporal correlations between parallel relays. These results provide the molecular details underlying the reported positive kinase cooperativity to enrich the discussion on how autonomous activation by phosphorylation leads to long-term behavioural effects. [ABSTRACT FROM AUTHOR]

Single-molecule techniques for protein sequencing are making headway towards single-cell proteomics and are projected to propel our understanding of cellular biology and disease. Yet, single cell proteomics presents a substantial unmet challenge due to the unavailability of protein amplification techniques, and the vast dynamic-range of protein expression in cells. Here, we describe and computationally investigate the feasibility of a novel approach for single-protein identification using tri-color fluorescence and plasmonic-nanopore devices. Comprehensive computer simulations of denatured protein translocation processes through the nanopores show that the tri-color fluorescence time-traces retain sufficient information to permit pattern-recognition algorithms to correctly identify the vast majority of proteins in the human proteome. Importantly, even when taking into account realistic experimental conditions, which restrict the spatial and temporal resolutions as well as the labeling efficiency, and add substantial noise, a deep-learning protein classifier achieves 97% whole-proteome accuracies. Applying our approach for protein datasets of clinical relevancy, such as the plasma proteome or cytokine panels, we obtain ~98% correct protein identification. This study suggests the feasibility of a method for accurate and high-throughput protein identification, which is highly versatile and applicable. [ABSTRACT FROM AUTHOR]

G protein-coupled receptors (GPCRs) control cellular signaling and responses. Many of these GPCRs are modulated by cholesterol and polyunsaturated fatty acids (PUFAs) which have been shown to co-exist with saturated lipids in ordered membrane domains. However, the lipid compositions of such domains extracted from the brain cortex tissue of individuals suffering from GPCR-associated neurological disorders show drastically lowered levels of PUFAs. Here, using free energy techniques and multiscale simulations of numerous membrane proteins, we show that the presence of the PUFA DHA helps helical multi-pass proteins such as GPCRs partition into ordered membrane domains. The mechanism is based on hybrid lipids, whose PUFA chains coat the rough protein surface, while the saturated chains face the raft environment, thus minimizing perturbations therein. Our findings suggest that the reduction of GPCR partitioning to their native ordered environments due to PUFA depletion might affect the function of these receptors in numerous neurodegenerative diseases, where the membrane PUFA levels in the brain are decreased. We hope that this work inspires experimental studies on the connection between membrane PUFA levels and GPCR signaling. [ABSTRACT FROM AUTHOR]

Author summary Improvements in technology often drive scientific discovery. Therefore, research requires sustained investment in the latest equipment and training for the researchers who are going to use it. Prioritising and administering infrastructure investment is challenging because future needs are difficult to predict. In the past, highly computationally demanding research was associated primarily with particle physics and astronomy experiments. However, as biology becomes more quantitative and bioscientists generate more and more data, their computational requirements may ultimately exceed those of physical scientists. Computation has always been central to bioinformatics, but now imaging experiments have rapidly growing data processing and storage requirements. There is also an urgent need for new modelling and simulation tools to provide insight and understanding of these biophysical experiments. Bioscience communities must work together to provide the software and skills training needed in their areas. Research-active institutions need to recognise that computation is now vital in many more areas of discovery and create an environment where it can be embraced. The public must also become aware of both the power and limitations of computing, particularly with respect to their health and personal data. [ABSTRACT FROM AUTHOR]

Fluorescence lifetime imaging microscopy (FLIM) can assess cell’s metabolism through the fluorescence of the co-enzymes NADH and FAD, which exhibit a double-exponential decay, with components related to free and protein-bound conditions. In vivo real time clinical imaging applications demand fast acquisition. As photodamage limits excitation power, this is best achieved using wide-field techniques, like time-gated FLIM, and algorithms that require few images to calculate the decay parameters. The rapid lifetime determination (RLD) algorithm requires only four images to analyze a double-exponential decay. Using computational simulations, we evaluated the accuracy and precision of RLD when measuring endogenous fluorescence lifetimes and metabolic free to protein-bound ratios, for total counts per pixel (TC) lower than 104. The simulations were based on a time-gated FLIM instrument, accounting for its instrument response function, gain and noise. While the optimal acquisition setting depends on the values being measured, the accuracy of the free to protein-bound ratio α2/α1 is stable for low gains and gate separations larger than 1000 ps, while its precision is almost constant for gate separations between 1500 and 2500 ps. For the gate separations and free to protein-bound ratios considered, the accuracy error can be as high as 30% and the precision error can reach 60%. Precision errors lower than 10% cannot be obtained. The best performance occurs for low camera gains and gate separations near 1800 ps. When considering the narrow physiological ranges for the free to protein-bound ratio, the precision errors can be confined to an interval between 10% and 20%. RLD is a valid option when for real time FLIM. The simulations and methodology presented here can be applied to any time-gated FLIM instrument and are useful to obtain the accuracy and precision limits for RLD in the demanding conditions of TC lower than 104. [ABSTRACT FROM AUTHOR]

Gastrotropin, the intracellular carrier of bile salts in the small intestine, binds two ligand molecules simultaneously in its internal cavity. The molecular rearrangements required for ligand entry are not yet fully clear. To improve our understanding of the binding process we combined molecular dynamics simulations with previously published structural and dynamic NMR parameters. The resulting ensembles reveal two distinct modes of barrel opening with one corresponding to the transition between the apo and holo states, whereas the other affecting different protein regions in both ligation states. Comparison of the calculated structures with NMR-derived parameters reporting on slow conformational exchange processes suggests that the protein undergoes partial unfolding along a path related to the second mode of the identified barrel opening motion. [ABSTRACT FROM AUTHOR]

The development of high-throughput biotechnologies allows the collection of omics data to study the biological mechanisms underlying complex diseases at different levels, such as genomics, epigenomics, and transcriptomics. However, each technology is designed to collect a specific type of omics data. Thus, the association between a disease and one type of omics data is usually tested individually, but this strategy is suboptimal. To better articulate biological processes and increase the consistency of variant identification, omics data from various platforms need to be integrated. In this report, we introduce an approach that uses a modified Fisher’s method (denoted as Omnibus-Fisher) to combine separate p-values of association testing for a trait and SNPs, DNA methylation markers, and RNA sequencing, calculated by kernel machine regression into an overall gene-level p-value to account for correlation between omics data. To consider all possible disease models, we extend Omnibus-Fisher to an optimal test by using perturbations. In our simulations, a usual Fisher’s method has inflated type I error rates when directly applied to correlated omics data. In contrast, Omnibus-Fisher preserves the expected type I error rates. Moreover, Omnibus-Fisher has increased power compared to its optimal version when the true disease model involves all types of omics data. On the other hand, the optimal Omnibus-Fisher is more powerful than its regular version when only one type of data is causal. Finally, we illustrate our proposed method by analyzing whole-genome genotyping, DNA methylation data, and RNA sequencing data from a study of childhood asthma in Puerto Ricans. [ABSTRACT FROM AUTHOR]

The qPCR method provides an inexpensive, rapid method for estimating relative telomere length across a set of biological samples. Like all laboratory methods, it involves some degree of measurement error. The estimation of relative telomere length is done subjecting the actual measurements made (the Cq values for telomere and a control gene) to non-linear transformations and combining them into a ratio (the TS ratio). Here, we use computer simulations, supported by mathematical analysis, to explore how errors in measurement affect qPCR estimates of relative telomere length, both in cross-sectional and longitudinal data. We show that errors introduced at the level of Cq values are magnified when the TS ratio is calculated. If the errors at the Cq level are normally distributed and independent of true telomere length, those in the TS ratio are positively skewed and proportional to true telomere length. The repeatability of the TS ratio declines sharply with increasing error in measurement of the Cq values for telomere and/or control gene. In simulated longitudinal data, measurement error alone can produce a pattern of low correlation between successive measures of relative telomere length, coupled with a strong negative dependency of the rate of change on initial relative telomere length. Our results illustrate the importance of reducing measurement error: a small increase in error in Cq values can have large consequences for the power and interpretability of qPCR estimates of relative telomere length. The findings also illustrate the importance of characterising the measurement error in each dataset—coefficients of variation are generally unhelpful, and researchers should report standard deviations of Cq values and/or repeatabilities of TS ratios—and allowing for the known effects of measurement error when interpreting patterns of TS ratio change over time. [ABSTRACT FROM AUTHOR]

Angiotensin-I converting enzyme (ACE) is a two-domain dipeptidylcarboxypeptidase involved in regulating blood pressure via the kallikrein-kininand renin-angiotensin-aldosterone complex. Therefore, ACE is a key drug target for the treatment of cardiovascular system diseases. At present many works are focus on searching for new inhibitory peptides of ACE to control the blood pressure. In order to exploit the interactions between ACE and its inhibitors, molecular dynamics simulations were used. The results showed that (a) the secondary structures of the three inhibitor-protein complexes did not change significantly; (b) root-mean-square deviation (RMSD), radius of gyration (Rg), and solvent-accessible surface area (SASA) values of Leu-Ile-Val-Thr (LIVT)-ACE complexes were significantly higher than that of other systems; (c) the backbone movement of LIVT was vigorous in Asp300-Val350, compared with that in Tyr-Leu-Val-Pro-His (YLVPH) and Tyr-Leu-Val-Arg(YLVR), as shown by the center-of-mass distance; and (d) the backbone movement of Asp300-Val350 may contribute to the interaction between ACE and its inhibitors. Our theoretical results will be helpful to further the design of specific inhibitors of ACE. [ABSTRACT FROM AUTHOR]

The actin family of cytoskeletal proteins is essential to the physiology of virtually all archaea, bacteria, and eukaryotes. While X-ray crystallography and electron microscopy have revealed structural homologies among actin-family proteins, these techniques cannot probe molecular-scale conformational dynamics. Here, we use all-atom molecular dynamic simulations to reveal conserved dynamical behaviors in four prokaryotic actin homologs: MreB, FtsA, ParM, and crenactin. We demonstrate that the majority of the conformational dynamics of prokaryotic actins can be explained by treating the four subdomains as rigid bodies. MreB, ParM, and FtsA monomers exhibited nucleotide-dependent dihedral and opening angles, while crenactin monomer dynamics were nucleotide-independent. We further show that the opening angle of ParM is sensitive to a specific interaction between subdomains. Steered molecular dynamics simulations of MreB, FtsA, and crenactin dimers revealed that changes in subunit dihedral angle lead to intersubunit bending or twist, suggesting a conserved mechanism for regulating filament structure. Taken together, our results provide molecular-scale insights into the nucleotide and polymerization dependencies of the structure of prokaryotic actins, suggesting mechanisms for how these structural features are linked to their diverse functions. [ABSTRACT FROM AUTHOR]

Lamellar and hexagonal lipid structures are of particular importance in the biological processes such as membrane fusion and budding. Atomistic simulations of formation of these phases and transitions between them are computationally prohibitive, hence development of coarse-grained models is an important part of the methodological development in this area. Here we apply systematic bottom-up coarse-graining to model different phase structures formed by 1,2-dioleoylphosphatidylethanolamine (DOPE) lipid molecules. We started from atomistic simulations of DOPE lipids in water carried out at two different water/lipid molar ratio corresponding to the lamellar Lα and inverted hexagonal HII structures at low and high lipid concentrations respectively. The atomistic trajectories were mapped to coarse-grained trajectories, in which each lipid was represented by 14 coarse-grained sites. Then the inverse Monte Carlo method was used to compute the effective coarse-grained potentials which for the coarse-grain model reproduce the same structural properties as the atomistic simulations. The potentials derived from the low concentration atomistic simulation were only able to form a bilayer structure, while both Lα and HII lipid phases were formed in simulations with potentials obtained at high concentration. The typical atomistic configurations of lipids at high concentration combine fragments of both lamellar and non-lamellar structures, that is reflected in the extracted coarse-grained potentials which become transferable and can form a wide range of structures including the inverted hexagonal, bilayer, tubule, vesicle and micellar structures. [ABSTRACT FROM AUTHOR]

Combining structural proteomics experimental data with computational methods is a powerful tool for protein structure prediction. Here, we apply a recently-developed approach for de novo protein structure determination based on the incorporation of short-distance crosslinking data as constraints in discrete molecular dynamics simulations (CL-DMD) for the determination of conformational ensemble of the intrinsically disordered protein α-synuclein in the solution. The predicted structures were in agreement with hydrogen-deuterium exchange, circular dichroism, surface modification, and long-distance crosslinking data. We found that α-synuclein is present in solution as an ensemble of rather compact globular conformations with distinct topology and inter-residue contacts, which is well-represented by movements of the large loops and formation of few transient secondary structure elements. Non-amyloid component and C-terminal regions were consistently found to contain β-structure elements and hairpins. [ABSTRACT FROM AUTHOR]

We propose and theoretically study an approach to massively parallel single molecule peptide sequencing, based on single molecule measurement of the kinetics of probe binding (Havranek, et al., 2013) to the N-termini of immobilized peptides. Unlike previous proposals, this method is robust to both weak and non-specific probe-target affinities, which we demonstrate by applying the method to a range of randomized affinity matrices consisting of relatively low-quality binders. This suggests a novel principle for proteomic measurement whereby highly non-optimized sets of low-affinity binders could be applicable for protein sequencing, thus shifting the burden of amino acid identification from biomolecular design to readout. Measurement of probe occupancy times, or of time-averaged fluorescence, should allow high-accuracy determination of N-terminal amino acid identity for realistic probe sets. The time-averaged fluorescence method scales well to weakly-binding probes with dissociation constants of tens or hundreds of micromolar, and bypasses photobleaching limitations associated with other fluorescence-based approaches to protein sequencing. We argue that this method could lead to an approach with single amino acid resolution and the ability to distinguish many canonical and modified amino acids, even using highly non-optimized probe sets. This readout method should expand the design space for single molecule peptide sequencing by removing constraints on the properties of the fluorescent binding probes. [ABSTRACT FROM AUTHOR]

Novel chemical controls are needed that selectively target human, animal, and plant parasitic nematodes with reduced adverse effects on the host or the environment. We hypothesize that the phosphodiesterase (PDE) enzyme family represents a potential target for development of novel nematicides and anthelmintics. To test this, we identified six PDE families present in the nematode phylum that are orthologous to six of the eleven human PDE families. We characterized the binding interactions of family-selective PDE inhibitors with human and C. elegans PDE4 in conjunction with molecular dynamics (MD) simulations to evaluate differences in binding interactions of these inhibitors within the PDE4 catalytic domain. We observed that roflumilast (human PDE4-selective inhibitor) and zardaverine (selective for human PDE3 and PDE4) were 159- and 77-fold less potent, respectively, in inhibiting C. elegans PDE4. The pan-specific PDE inhibitor isobutyl methyl xanthine (IBMX) had similar affinity for nematode and human PDE4. Of 32 residues within 5 Å of the ligand binding site, five revealed significant differences in non-bonded interaction energies (van der Waals and electrostatic interaction energies) that could account for the differential binding affinities of roflumilast and zardaverine. One site (Phe506 in the human PDE4D3 amino acid sequence corresponding to Tyr253 in C. elegans PDE4) is predicted to alter the binding conformation of roflumilast and zardaverine (but not IBMX) into a less energetically favorable state for the nematode enzyme. The pharmacological differences in sensitivity to PDE4 inhibitors in conjunction with differences in the amino acids comprising the inhibitor binding sites of human and C. elegans PDE4 catalytic domains together support the feasibility of designing the next generation of anthelmintics/nematicides that could selectively bind to nematode PDEs. [ABSTRACT FROM AUTHOR]

ChIP-seq (Chromatin Immunoprecipitation followed by sequencing) is a high-throughput technique to identify genomic regions that are bound in vivo by a particular protein, e.g., a transcription factor (TF). Biological factors, such as chromatin state, indirect and cooperative binding, as well as experimental factors, such as antibody quality, cross-linking, and PCR biases, are known to affect the outcome of ChIP-seq experiments. However, the relative impact of these factors on inferences made from ChIP-seq data is not entirely clear. Here, via a detailed ChIP-seq simulation pipeline, ChIPulate, we assess the impact of various biological and experimental sources of variation on several outcomes of a ChIP-seq experiment, viz., the recoverability of the TF binding motif, accuracy of TF-DNA binding detection, the sensitivity of inferred TF-DNA binding strength, and number of replicates needed to confidently infer binding strength. We find that the TF motif can be recovered despite poor and non-uniform extraction and PCR amplification efficiencies. The recovery of the motif is, however, affected to a larger extent by the fraction of sites that are either cooperatively or indirectly bound. Importantly, our simulations reveal that the number of ChIP-seq replicates needed to accurately measure in vivo occupancy at high-affinity sites is larger than the recommended community standards. Our results establish statistical limits on the accuracy of inferences of protein-DNA binding from ChIP-seq and suggest that increasing the mean extraction efficiency, rather than amplification efficiency, would better improve sensitivity. The source code and instructions for running ChIPulate can be found at . [ABSTRACT FROM AUTHOR]

The fundamental aim of structural analyses in biophysics is to reveal a mutual relation between a molecule's dynamic structure and its physiological function. Small-angle X-ray scattering (SAXS) is an experimental technique for structural characterization of macromolecules in solution and enables time-resolved analysis of conformational changes under physiological conditions. As such experiments measure spatially averaged low-resolution scattering intensities only, the sparse information obtained is not sufficient to uniquely reconstruct a three-dimensional atomistic model. Here, we integrate the information from SAXS into molecular dynamics simulations using computationally efficient native structure-based models. Dynamically fitting an initial structure towards a scattering intensity, such simulations produce atomistic models in agreement with the target data. In this way, SAXS data can be rapidly interpreted while retaining physico-chemical knowledge and sampling power of the underlying force field. We demonstrate our method's performance using the example of three protein systems. Simulations are faster than full molecular dynamics approaches by more than two orders of magnitude and consistently achieve comparable accuracy. Computational demands are reduced sufficiently to run the simulations on commodity desktop computers instead of high-performance computing systems. These results underline that scattering-guided structure-based simulations provide a suitable framework for rapid early-stage refinement of structures towards SAXS data with particular focus on minimal computational resources and time. [ABSTRACT FROM AUTHOR]