RESUMO
SUCNR1 is an auto- and paracrine sensor of the metabolic stress signal succinate. Using unsupervised molecular dynamics (MD) simulations (170.400 ns) and mutagenesis across human, mouse, and rat SUCNR1, we characterize how a five-arginine motif around the extracellular pole of TM-VI determines the initial capture of succinate in the extracellular vestibule (ECV) to either stay or move down to the orthosteric site. Metadynamics demonstrate low-energy succinate binding in both sites, with an energy barrier corresponding to an intermediate stage during which succinate, with an associated water cluster, unlocks the hydrogen-bond-stabilized conformationally constrained extracellular loop (ECL)-2b. Importantly, simultaneous binding of two succinate molecules through either a "sequential" or "bypassing" mode is a frequent endpoint. The mono-carboxylate NF-56-EJ40 antagonist enters SUCNR1 between TM-I and -II and does not unlock ECL-2b. It is proposed that occupancy of both high-affinity sites is required for selective activation of SUCNR1 by high local succinate concentrations.
Assuntos
Receptores Acoplados a Proteínas G , Ácido Succínico , Camundongos , Ratos , Animais , Humanos , Ácido Succínico/metabolismo , Receptores Acoplados a Proteínas G/metabolismo , Simulação de Dinâmica Molecular , Succinatos/metabolismo , Estresse FisiológicoRESUMO
Cryptococcus neoformans is a fungal pathogen responsible for cryptococcosis and cryptococcal meningitis. The C. neoformans' capsular polysaccharide and its shed exopolysaccharide function both as key virulence factors and to protect the fungal cell from phagocytosis. Currently, a glycoconjugate of these polysaccharides is being explored as a vaccine to protect against C. neoformans infection. In this study, NOE and J-coupling values from NMR experiments were consistent with a converged structure of the synthetic decasaccharide, GXM10-Ac3, calculated from MD simulations. GXM10-Ac3 was designed as an extension of glucuronoxylomannan (GXM) polysaccharide motif (M2) which is common in the clinically predominant serotype A strains and is recognized by protective forms of GXM-specific monoclonal antibodies. The M2 motif is a hexasaccharide with a three-residue α-mannan backbone, modified by ß-(1â2)-xyloses (Xyl) on the first two mannoses (Man) and a ß-(1â2)-glucuronic acid (GlcA) on the third Man. Combined NMR and MD analyses reveal that GXM10-Ac3 adopts an extended structure, with Xyl/GlcA branches alternating sides along the α-mannan backbone. O-acetyl esters also alternate sides and are grouped in pairs. MD analysis of a twelve M2-repeating unit polymer supports the notion that the GXM10-Ac3 structure is uniformly represented throughout the polysaccharide. This derived GXM model displays high flexibility while maintaining a structural identity, yielding insights to further explore intermolecular interactions between polysaccharides, interactions with anti-GXM mAbs, and the cryptococcal polysaccharide architecture.
Assuntos
Criptococose , Cryptococcus neoformans , Humanos , Mananas , Cryptococcus neoformans/química , Polissacarídeos/química , Criptococose/microbiologia , Espectroscopia de Ressonância Magnética , Anticorpos Monoclonais , Anticorpos AntifúngicosRESUMO
In recent years, self-assembly has emerged as a powerful tool for fabricating functional materials. Since self-assembly is fundamentally determined by the particle interactions in the system, if we can gain full control over these interactions, it would open the door for creating functional materials by design. In this paper, we exploit capillary interactions between colloidal particles at liquid interfaces to create two-dimensional (2D) materials where particle interactions and self-assembly can be fully programmed using particle shape alone. Specifically, we consider colloidal particles which are polygonal plates with homogeneous surface chemistry and undulating edges as this particle geometry gives us precise and independent control over both short-range hard-core repulsions and longer-range capillary interactions. To illustrate the immense potential provided by our system for programming self-assembly, we use minimum energy calculations and Monte Carlo simulations to show that polygonal plates with different in-plane shapes (hexagons, truncated triangles, triangles, squares) and edge undulations of different multipolar order (hexapolar, octopolar, dodecapolar) can be used to create a rich variety of 2D structures, including hexagonal close-packed, honeycomb, Kagome, and quasicrystal lattices. Since the required particle shapes can be readily fabricated experimentally, we can use our colloidal system to control the entire process chain for materials design, from initial design and fabrication of the building blocks, to final assembly of the emergent 2D material.
RESUMO
Solutions of long, flexible polymer molecules are complex fluids that simultaneously exhibit fluid-like and solid-like behavior. When subjected to an external flow, dilute polymer solutions exhibit elastic turbulence-a unique, chaotic flow state absent in Newtonian fluids, like water. Unlike its Newtonian counterpart, elastic turbulence is caused by polymer molecules stretching and aligning in the flow, and can occur at vanishing inertia. While experimental realizations of elastic turbulence are well-documented, there is currently no understanding of its mechanism. Here, we present large-scale direct numerical simulations of elastic turbulence in pressure-driven flows through straight channels. We demonstrate that the transition to elastic turbulence is sub-critical, giving rise to spot-like flow structures that, further away from the transition, eventually spread throughout the domain. We provide evidence that elastic turbulence is organized around unstable coherent states that are localized close to the channel midplane.
RESUMO
Metallic alloys often form phases-known as solid solutions-in which chemical elements are spread out on the same crystal lattice in an almost random manner. The tendency of certain chemical motifs to be more common than others is known as chemical short-range order (SRO), and it has received substantial consideration in alloys with multiple chemical elements present in large concentrations due to their extreme configurational complexity (e.g., high-entropy alloys). SRO renders solid solutions "slightly less random than completely random," which is a physically intuitive picture, but not easily quantifiable due to the sheer number of possible chemical motifs and their subtle spatial distribution on the lattice. Here, we present a multiscale method to predict and quantify the SRO state of an alloy with atomic resolution, incorporating machine learning techniques to bridge the gap between electronic-structure calculations and the characteristic length scale of SRO. The result is an approach capable of predicting SRO length scale in agreement with experimental measurements while comprehensively correlating SRO with fundamental quantities such as local lattice distortions. This work advances the quantitative understanding of solid-solution phases, paving the way for the rigorous incorporation of SRO length scales into predictive mechanical and thermodynamic models.
RESUMO
Clearance of serotonin (5-hydroxytryptamine, 5-HT) from the synaptic cleft after neuronal signaling is mediated by serotonin transporter (SERT), which couples this process to the movement of a Na+ ion down its chemical gradient. After release of 5-HT and Na+ into the cytoplasm, the transporter faces a rate-limiting challenge of resetting its conformation to be primed again for 5-HT and Na+ binding. Early studies of vesicles containing native SERT revealed that K+ gradients can provide an additional driving force, via K+ antiport. Moreover, under appropriate conditions, a H+ ion can replace K+. Intracellular K+ accelerates the resetting step. Structural studies of SERT have identified two binding sites for Na+ ions, but the K+ site remains enigmatic. Here, we show that K+ antiport can drive substrate accumulation into vesicles containing SERT extracted from a heterologous expression system, allowing us to study the residues responsible for K+ binding. To identify candidate binding residues, we examine many cation binding configurations using molecular dynamics simulations, predicting that K+ binds to the so-called Na2 site. Site-directed mutagenesis of residues in this site can eliminate the ability of both K+ and H+ to drive 5-HT accumulation into vesicles and, in patch clamp recordings, prevent the acceleration of turnover rates and the formation of a channel-like state by K+ or H+. In conclusion, the Na2 site plays a pivotal role in orchestrating the sequential binding of Na+ and then K+ (or H+) ions to facilitate 5-HT uptake in SERT.
Assuntos
Simulação de Dinâmica Molecular , Potássio , Proteínas da Membrana Plasmática de Transporte de Serotonina , Sódio , Proteínas da Membrana Plasmática de Transporte de Serotonina/metabolismo , Proteínas da Membrana Plasmática de Transporte de Serotonina/química , Proteínas da Membrana Plasmática de Transporte de Serotonina/genética , Potássio/metabolismo , Sítios de Ligação , Humanos , Sódio/metabolismo , Serotonina/metabolismo , Ligação Proteica , AnimaisRESUMO
The SNAP receptor (SNARE) proteins syntaxin-1, SNAP-25, and synaptobrevin mediate neurotransmitter release by forming tight SNARE complexes that fuse synaptic vesicles with the plasma membranes in microseconds. Membrane fusion is generally explained by the action of proteins on macroscopic membrane properties such as curvature, elastic modulus, and tension, and a widespread model envisions that the SNARE motifs, juxtamembrane linkers, and C-terminal transmembrane regions of synaptobrevin and syntaxin-1 form continuous helices that act mechanically as semirigid rods, squeezing the membranes together as they assemble ("zipper") from the N to the C termini. However, the mechanism underlying fast SNARE-induced membrane fusion remains unknown. We have used all-atom molecular dynamics simulations to investigate this mechanism. Our results need to be interpreted with caution because of the limited number and length of the simulations, but they suggest a model of membrane fusion that has a natural physicochemical basis, emphasizes local molecular events over general membrane properties, and explains extensive experimental data. In this model, the central event that initiates fast (microsecond scale) membrane fusion occurs when the SNARE helices zipper into the juxtamembrane linkers which, together with the adjacent transmembrane regions, promote encounters of acyl chains from both bilayers at the polar interface. The resulting hydrophobic nucleus rapidly expands into stalk-like structures that gradually progress to form a fusion pore, aided by the SNARE transmembrane regions and without clearly discernible intermediates. The propensity of polyunsaturated lipids to participate in encounters that initiate fusion suggests that these lipids may be important for the high speed of neurotransmitter release.
Assuntos
Fusão de Membrana , Proteínas SNARE , Proteínas SNARE/metabolismo , Simulação de Dinâmica Molecular , Proteínas R-SNARE , Sintaxina 1 , Neurotransmissores , LipídeosRESUMO
Clay minerals are implicated in the retention of biomolecules within organic matter in many soil environments. Spectroscopic studies have proposed several mechanisms for biomolecule adsorption on clays. Here, we employ molecular dynamics simulations to investigate these mechanisms in hydrated adsorbate conformations of montmorillonite, a smectite-type clay, with ten biomolecules of varying chemistry and structure, including sugars related to cellulose and hemicellulose, lignin-related phenolic acid, and amino acids with different functional groups. Our molecular modeling captures biomolecule-clay and biomolecule-biomolecule interactions that dictate selectivity and competition in adsorption retention and interlayer nanopore trapping, which we determine experimentally by NMR and X-ray diffraction, respectively. Specific adsorbate structures are important in facilitating the electrostatic attraction and Van der Waals energies underlying the hierarchy in biomolecule adsorption. Stabilized by a network of direct and water-bridged hydrogen bonds, favorable electrostatic interactions drive this hierarchy whereby amino acids with positively charged side chains are preferentially adsorbed on the negatively charged clay surface compared to the sugars and carboxylate-rich aromatics and amino acids. With divalent metal cations, our model adsorbate conformations illustrate hydrated metal cation bridging of carboxylate-containing biomolecules to the clay surface, thus explaining divalent cation-promoted adsorption from our experimental data. Adsorption experiments with a mixture of biomolecules reveal selective inhibition in biomolecule adsorption, which our molecular modeling attributes to electrostatic biomolecule-biomolecule pairing that is more energetically favorable than the biomolecule-clay complex. In sum, our findings highlight chemical and structural features that can inform hypotheses for predicting biomolecule adsorption at water-clay interfaces.
Assuntos
Simulação de Dinâmica Molecular , Água , Argila , Adsorção , Água/química , Eletricidade Estática , Aminoácidos , AçúcaresRESUMO
The estrogen receptor-α (ER) is thought to function only as a homodimer but responds to a variety of environmental, metazoan, and therapeutic estrogens at subsaturating doses, supporting binding mixtures of ligands as well as dimers that are only partially occupied. Here, we present a series of flexible ER ligands that bind to receptor dimers with individual ligand poses favoring distinct receptor conformations-receptor conformational heterodimers-mimicking the binding of two different ligands. Molecular dynamics simulations showed that the pairs of different ligand poses changed the correlated motion across the dimer interface to generate asymmetric communication between the dimer interface, the ligands, and the surface binding sites for epigenetic regulatory proteins. By examining the binding of the same ligand in crystal structures of ER in the agonist vs. antagonist conformers, we also showed that these allosteric signals are bidirectional. The receptor conformer can drive different ligand binding modes to support agonist vs. antagonist activity profiles, a revision of ligand binding theory that has focused on unidirectional signaling from the ligand to the coregulator binding site. We also observed differences in the allosteric signals between ligand and coregulator binding sites in the monomeric vs. dimeric receptor, and when bound by two different ligands, states that are physiologically relevant. Thus, ER conformational heterodimers integrate two different ligand-regulated activity profiles, representing different modes for ligand-dependent regulation of ER activity.
Assuntos
Receptor alfa de Estrogênio , Estrogênios , Simulação de Dinâmica Molecular , Multimerização Proteica , Receptor alfa de Estrogênio/metabolismo , Receptor alfa de Estrogênio/química , Regulação Alostérica , Humanos , Ligantes , Estrogênios/metabolismo , Estrogênios/química , Sítios de Ligação , Ligação Proteica , Conformação ProteicaRESUMO
We unveil the multifractal behavior of Ising spin glasses in their low-temperature phase. Using the Janus II custom-built supercomputer, the spin-glass correlation function is studied locally. Dramatic fluctuations are found when pairs of sites at the same distance are compared. The scaling of these fluctuations, as the spin-glass coherence length grows with time, is characterized through the computation of the singularity spectrum and its corresponding Legendre transform. A comparatively small number of site pairs controls the average correlation that governs the response to a magnetic field. We explain how this scenario of dramatic fluctuations (at length scales smaller than the coherence length) can be reconciled with the smooth, self-averaging behavior that has long been considered to describe spin-glass dynamics.
RESUMO
Atypical Chemokine Receptor 3 (ACKR3) belongs to the G protein-coupled receptor family but it does not signal through G proteins. The structural properties that govern the functional selectivity and the conformational dynamics of ACKR3 activation are poorly understood. Here, we combined hydrogen/deuterium exchange mass spectrometry, site-directed mutagenesis, and molecular dynamics simulations to examine the binding mode and mechanism of action of ACKR3 ligands of different efficacies. Our results show that activation or inhibition of ACKR3 is governed by intracellular conformational changes of helix 6, intracellular loop 2, and helix 7, while the DRY motif becomes protected during both processes. Moreover, we identified the binding sites and the allosteric modulation of ACKR3 upon ß-arrestin 1 binding. In summary, this study highlights the structure-function relationship of small ligands, the binding mode of ß-arrestin 1, the activation dynamics, and the atypical dynamic features in ACKR3 that may contribute to its inability to activate G proteins.
Assuntos
Simulação de Dinâmica Molecular , Ligação Proteica , Receptores CXCR , Humanos , Receptores CXCR/metabolismo , Receptores CXCR/genética , Sítios de Ligação , Conformação Proteica , beta-Arrestina 1/metabolismo , beta-Arrestina 1/genética , Ligantes , Células HEK293 , Mutagênese Sítio-Dirigida , Regulação Alostérica , Relação Estrutura-AtividadeRESUMO
Proton-powered c-ring rotation in mitochondrial ATP synthase is crucial to convert the transmembrane protonmotive force into torque to drive the synthesis of adenosine triphosphate (ATP). Capitalizing on recent cryo-EM structures, we aim at a structural and energetic understanding of how functional directional rotation is achieved. We performed multi-microsecond atomistic simulations to determine the free energy profiles along the c-ring rotation angle before and after the arrival of a new proton. Our results reveal that rotation proceeds by dynamic sliding of the ring over the a-subunit surface, during which interactions with conserved polar residues stabilize distinct intermediates. Ordered water chains line up for a Grotthuss-type proton transfer in one of these intermediates. After proton transfer, a high barrier prevents backward rotation and an overall drop in free energy favors forward rotation, ensuring the directionality of c-ring rotation required for the thermodynamically disfavored ATP synthesis. The essential arginine of the a-subunit stabilizes the rotated configuration through a salt bridge with the c-ring. Overall, we describe a complete mechanism for the rotation step of the ATP synthase rotor, thereby illuminating a process critical to all life at atomic resolution.
Assuntos
ATPases Mitocondriais Próton-Translocadoras , Prótons , ATPases Mitocondriais Próton-Translocadoras/metabolismo , Conformação Proteica , Trifosfato de Adenosina , Rotação , ATPases Translocadoras de Prótons/metabolismoRESUMO
The fat mass and obesity-associated fatso (FTO) protein is a member of the Alkb family of dioxygenases and catalyzes oxidative demethylation of N6-methyladenosine (m6A), N1-methyladenosine (m1A), 3-methylthymine (m3T), and 3-methyluracil (m3U) in single-stranded nucleic acids. It is well established that the catalytic activity of FTO proceeds via two coupled reactions. The first reaction involves decarboxylation of alpha-ketoglutarate (αKG) and formation of an oxyferryl species. In the second reaction, the oxyferryl intermediate oxidizes the methylated nucleic acid to reestablish Fe(II) and the canonical base. However, it remains unclear how binding of the nucleic acid activates the αKG decarboxylation reaction and why FTO demethylates different methyl modifications at different rates. Here, we investigate the interaction of FTO with 5-mer DNA oligos incorporating the m6A, m1A, or m3T modifications using solution NMR, molecular dynamics (MD) simulations, and enzymatic assays. We show that binding of the nucleic acid to FTO activates a two-state conformational equilibrium in the αKG cosubstrate that modulates the O2 accessibility of the Fe(II) catalyst. Notably, the substrates that provide better stabilization to the αKG conformation in which Fe(II) is exposed to O2 are demethylated more efficiently by FTO. These results indicate that i) binding of the methylated nucleic acid is required to expose the catalytic metal to O2 and activate the αKG decarboxylation reaction, and ii) the measured turnover of the demethylation reaction (which is an ensemble average over the entire sample) depends on the ability of the methylated base to favor the Fe(II) state accessible to O2.
Assuntos
Dioxigenase FTO Dependente de alfa-Cetoglutarato , Ferro , Ácidos Cetoglutáricos , Dioxigenase FTO Dependente de alfa-Cetoglutarato/metabolismo , Dioxigenase FTO Dependente de alfa-Cetoglutarato/química , Dioxigenase FTO Dependente de alfa-Cetoglutarato/genética , Ácidos Cetoglutáricos/metabolismo , Ácidos Cetoglutáricos/química , Ferro/metabolismo , Ferro/química , Humanos , Especificidade por Substrato , Adenosina/análogos & derivados , Adenosina/metabolismo , Adenosina/química , Conformação Proteica , Uracila/metabolismo , Uracila/análogos & derivados , Uracila/química , Simulação de Dinâmica Molecular , Timina/análogos & derivadosRESUMO
Mechanical energy, specifically in the form of ultrasound, can induce pressure variations and temperature fluctuations when applied to an aqueous media. These conditions can both positively and negatively affect protein complexes, consequently altering their stability, folding patterns, and self-assembling behavior. Despite much scientific progress, our current understanding of the effects of ultrasound on the self-assembly of amyloidogenic proteins remains limited. In the present study, we demonstrate that when the amplitude of the delivered ultrasonic energy is sufficiently low, it can induce refolding of specific motifs in protein monomers, which is sufficient for primary nucleation; this has been revealed by MD. These ultrasound-induced structural changes are initiated by pressure perturbations and are accelerated by a temperature factor. Furthermore, the prolonged action of low-amplitude ultrasound enables the elongation of amyloid protein nanofibrils directly from natively folded monomeric lysozyme protein, in a controlled manner, until it reaches a critical length. Using solution X-ray scattering, we determined that nanofibrillar assemblies, formed either under the action of sound or from natively fibrillated lysozyme, share identical structural characteristics. Thus, these results provide insights into the effects of ultrasound on fibrillar protein self-assembly and lay the foundation for the potential use of sound energy in protein chemistry.
Assuntos
Amiloide , Muramidase , Amiloide/química , Amiloide/metabolismo , Muramidase/química , Muramidase/metabolismo , Dobramento de Proteína , Temperatura , Ondas Ultrassônicas , Simulação de Dinâmica MolecularRESUMO
Localized deformation and randomly shaped imperfections are salient features of buckling-type instabilities in thin-walled load-bearing structures. However, it is generally agreed that their complex interactions in response to mechanical loading are not yet sufficiently understood, as evidenced by buckling-induced catastrophic failures which continue to today. This study investigates how the intimate coupling between localization mechanisms and geometric imperfections combine to determine the statistics of the pressure required to buckle (the illustrative example of) a hemispherical shell. The geometric imperfections, in the form of a surface, are defined by a random field generated over the nominally hemispherical shell geometry, and the probability distribution of the buckling pressure is computed via stochastic finite element analysis. Monte-Carlo simulations are performed for a wide range of the shell's radius to thickness ratio, as well as the correlation length of the spatial distribution of the imperfection. The results show that over this range, the buckling pressure is captured by the Weibull distribution. In addition, the analyses of the deformation patterns observed during the simulations provide insights into the effects of certain characteristic lengths on the local buckling that triggers global instability. In light of the simulation results, a probabilistic model is developed for the statistics of the buckling load that reveals how the dimensionless radius plays a dual role which remained hidden in previous deterministic analyses. The implications of the present model for reliability-based design of shell structures are discussed.
RESUMO
Proteins play a key role in biological electron transport, but the structure-function relationships governing the electronic properties of peptides are not fully understood. Despite recent progress, understanding the link between peptide conformational flexibility, hierarchical structures, and electron transport pathways has been challenging. Here, we use single-molecule experiments, molecular dynamics (MD) simulations, nonequilibrium Green's function-density functional theory (NEGF-DFT), and unsupervised machine learning to understand the role of secondary structure on electron transport in peptides. Our results reveal a two-state molecular conductance behavior for peptides across several different amino acid sequences. MD simulations and Gaussian mixture modeling are used to show that this two-state molecular conductance behavior arises due to the conformational flexibility of peptide backbones, with a high-conductance state arising due to a more defined secondary structure (beta turn or 310 helices) and a low-conductance state occurring for extended peptide structures. These results highlight the importance of helical conformations on electron transport in peptides. Conformer selection for the peptide structures is rationalized using principal component analysis of intramolecular hydrogen bonding distances along peptide backbones. Molecular conformations from MD simulations are used to model charge transport in NEGF-DFT calculations, and the results are in reasonable qualitative agreement with experiments. Projected density of states calculations and molecular orbital visualizations are further used to understand the role of amino acid side chains on transport. Overall, our results show that secondary structure plays a key role in electron transport in peptides, which provides broad avenues for understanding the electronic properties of proteins.
Assuntos
Simulação de Dinâmica Molecular , Peptídeos , Estrutura Secundária de Proteína , Transporte de Elétrons , Peptídeos/química , Peptídeos/metabolismo , Ligação de HidrogênioRESUMO
Understanding the intracellular dynamics of brain cells entails performing three-dimensional molecular simulations incorporating ultrastructural models that can capture cellular membrane geometries at nanometer scales. While there is an abundance of neuronal morphologies available online, e.g. from NeuroMorpho.Org, converting those fairly abstract point-and-diameter representations into geometrically realistic and simulation-ready, i.e. watertight, manifolds is challenging. Many neuronal mesh reconstruction methods have been proposed; however, their resulting meshes are either biologically unplausible or non-watertight. We present an effective and unconditionally robust method capable of generating geometrically realistic and watertight surface manifolds of spiny cortical neurons from their morphological descriptions. The robustness of our method is assessed based on a mixed dataset of cortical neurons with a wide variety of morphological classes. The implementation is seamlessly extended and applied to synthetic astrocytic morphologies that are also plausibly biological in detail. Resulting meshes are ultimately used to create volumetric meshes with tetrahedral domains to perform scalable in silico reaction-diffusion simulations for revealing cellular structure-function relationships. Availability and implementation: Our method is implemented in NeuroMorphoVis, a neuroscience-specific open source Blender add-on, making it freely accessible for neuroscience researchers.
Assuntos
Simulação por Computador , Neurônios , Neurônios/ultraestrutura , Neurônios/citologia , Modelos Neurológicos , Humanos , Animais , Astrócitos/citologia , Astrócitos/ultraestruturaRESUMO
Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully length messenger RNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in polymerase chain reaction (PCR) amplification, barcode read errors and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.
Assuntos
Genoma , RNA , RNA-Seq , Análise de Sequência de RNA , Simulação por Computador , RNA/genética , Sequenciamento de Nucleotídeos em Larga EscalaRESUMO
Somatic variation is a major type of genetic variation contributing to human diseases including cancer. Of the vast quantities of somatic variants identified, the functional impact of many somatic variants, in particular the missense variants, remains unclear. Lack of the functional information prevents the translation of rich variation data into clinical applications. We previously developed a method named Ramachandran Plot-Molecular Dynamics Simulations (RP-MDS), aiming to predict the function of germline missense variants based on their effects on protein structure stability, and successfully applied to predict the deleteriousness of unclassified germline missense variants in multiple cancer genes. We hypothesized that regardless of their different genetic origins, somatic missense variants and germline missense variants could have similar effects on the stability of their affected protein structure. As such, the RP-MDS method designed for germline missense variants should also be applicable to predict the function of somatic missense variants. In the current study, we tested our hypothesis by using the somatic missense variants in TP53 as a model. Of the 397 somatic missense variants analyzed, RP-MDS predicted that 195 (49.1%) variants were deleterious as they significantly disturbed p53 structure. The results were largely validated by using a p53-p21 promoter-green fluorescent protein (GFP) reporter gene assay. Our study demonstrated that deleterious somatic missense variants can be identified by referring to their effects on protein structural stability.
Assuntos
Mutação de Sentido Incorreto , Estabilidade Proteica , Proteína Supressora de Tumor p53 , Humanos , Proteína Supressora de Tumor p53/genética , Proteína Supressora de Tumor p53/química , Simulação de Dinâmica Molecular , Neoplasias/genética , Conformação ProteicaRESUMO
Cross-linkers play a critical role in capturing protein dynamics in chemical cross-linking mass spectrometry techniques. Various types of cross-linkers with different backbone features are widely used in the study of proteins. However, it is still not clear how the cross-linkers' backbone affect their own structure and their interactions with proteins. In this study, we systematically characterized and compared methylene backbone and polyethylene glycol (PEG) backbone cross-linkers in terms of capturing protein structure and dynamics. The results indicate the cross-linker with PEG backbone have a better ability to capture the inter-domain dynamics of calmodulin, adenylate kinase, maltodextrin binding protein and dual-specificity protein phosphatase. We further conducted quantum chemical calculations and all-atom molecular dynamics simulations to analyze thermodynamic and kinetic properties of PEG backbone and methylene backbone cross-linkers. Solution nuclear magnetic resonance was employed to validate the interaction interface between proteins and cross-linkers. Our findings suggest that the polarity distribution of PEG backbone enhances the accessibility of the cross-linker to the protein surface, facilitating the capture of sites located in dynamic regions. By comprehensively benchmarking with disuccinimidyl suberate (DSS)/bis-sulfosuccinimidyl-suberate(BS3), bis-succinimidyl-(PEG)2 revealed superior advantages in protein dynamic conformation analysis in vitro and in vivo, enabling the capture of a greater number of cross-linking sites and better modeling of protein dynamics. Furthermore, our study provides valuable guidance for the development and application of PEG backbone cross-linkers.