RESUMO
A significant challenge to making targeted cystic fibrosis transmembrane conductance regulator (CFTR) modulator therapies accessible to all individuals with cystic fibrosis (CF) are many mutations in the CFTR gene that can cause CF, most of which remain uncharacterized. Here, we characterized the structural and functional defects of the rare CFTR mutation R352Q, with a potential role contributing to intrapore chloride ion permeation, in patient-derived cell models of the airway and gut. CFTR function in differentiated nasal epithelial cultures and matched intestinal organoids was assessed using an ion transport assay and forskolin-induced swelling assay, respectively. CFTR potentiators (VX-770, GLPG1837, and VX-445) and correctors (VX-809, VX-445, with or without VX-661) were tested. Data from R352Q-CFTR were compared with data of 20 participants with mutations with known impact on CFTR function. R352Q-CFTR has residual CFTR function that was restored to functional CFTR activity by CFTR potentiators but not the corrector. Molecular dynamics simulations of R352Q-CFTR were carried out, which indicated the presence of a chloride conductance defect, with little evidence supporting a gating defect. The combination approach of in vitro patient-derived cell models and in silico molecular dynamics simulations to characterize rare CFTR mutations can improve the specificity and sensitivity of modulator response predictions and aid in their translational use for CF precision medicine.
Assuntos
Regulador de Condutância Transmembrana em Fibrose Cística , Fibrose Cística , Aminofenóis/farmacologia , Cloretos/metabolismo , Fibrose Cística/genética , Regulador de Condutância Transmembrana em Fibrose Cística/genética , Regulador de Condutância Transmembrana em Fibrose Cística/metabolismo , Humanos , Simulação de Dinâmica Molecular , Mutação , Organoides/metabolismoRESUMO
Von Willebrand disease (VWD) is a bleeding disorder with different levels of severity. VWD-associated mutations are located in the von Willebrand factor (VWF) gene, coding for the large multidomain plasma protein VWF with essential roles in hemostasis and thrombosis. On the one hand, a variety of mutations in the C-domains of VWF are associated with increased bleeding upon vascular injury. On the other hand, VWF gain-of-function (GOF) mutations in the C4 domain have recently been identified, which induce an increased risk of myocardial infarction. Mechanistic insights into how these mutations affect the molecular behavior of VWF are scarce and holistic approaches are challenging due to the multidomain and multimeric character of this large protein. Here, we determine the structure and dynamics of the C6 domain and the single nucleotide polymorphism (SNP) variant G2705R in C6 by combining nuclear magnetic resonance spectroscopy, molecular dynamics simulations and aggregometry. Our findings indicate that this mutation mostly destabilizes VWF by leading to a more pronounced hinging between both subdomains of C6. Hemostatic parameters of variant G2705R are close to normal under static conditions, but the missense mutation results in a gain-of-function under flow conditions, due to decreased VWF stem stability. Together with the fact that two C4 variants also exhibit GOF characteristics, our data underline the importance of the VWF stem region in VWF's hemostatic activity and the risk of mutation-associated prothrombotic properties in VWF C-domain variants due to altered stem dynamics.
Assuntos
Fator de von Willebrand , Fator de von Willebrand/genéticaRESUMO
Von Willebrand factor (VWF) is a key player in the regulation of hemostasis by promoting recruitment of platelets to sites of vascular injury. An array of 6 C domains forms the dimeric C-terminal VWF stem. Upon shear force activation, the stem adopts an open conformation allowing the adhesion of VWF to platelets and the vessel wall. To understand the underlying molecular mechanism and associated functional perturbations in disease-related variants, knowledge of high-resolution structures and dynamics of C domains is of paramount interest. Here, we present the solution structure of the VWF C4 domain, which binds to the platelet integrin and is therefore crucial for the VWF function. In the structure, we observed 5 intra- and inter-subdomain disulfide bridges, of which 1 is unique in the C4 domain. The structure further revealed an unusually hinged 2-subdomain arrangement. The hinge is confined to a very short segment around V2547 connecting the 2 subdomains. Together with 2 nearby inter-subdomain disulfide bridges, this hinge induces slow conformational changes and positional alternations of both subdomains with respect to each other. Furthermore, the structure demonstrates that a clinical gain-of-function VWF variant (Y2561) is more likely to have an effect on the arrangement of the C4 domain with neighboring domains rather than impairing platelet integrin binding.
Assuntos
Plaquetas/metabolismo , Integrinas/metabolismo , Fator de von Willebrand/química , Fator de von Willebrand/metabolismo , Sequência de Aminoácidos , Dissulfetos/metabolismo , Humanos , Espectroscopia de Ressonância Magnética , Modelos Moleculares , Ligação Proteica , Domínios Proteicos , Soluções , Relação Estrutura-AtividadeRESUMO
Maleless (MLE) is an evolutionary conserved member of the DExH family of helicases in Drosophila. Besides its function in RNA editing and presumably siRNA processing, MLE is best known for its role in remodelling non-coding roX RNA in the context of X chromosome dosage compensation in male flies. MLE and its human orthologue, DHX9 contain two tandem double-stranded RNA binding domains (dsRBDs) located at the N-terminal region. The two dsRBDs are essential for localization of MLE at the X-territory and it is presumed that this involves binding roX secondary structures. However, for dsRBD1 roX RNA binding has so far not been described. Here, we determined the solution NMR structure of dsRBD1 and dsRBD2 of MLE in tandem and investigated its role in double-stranded RNA (dsRNA) binding. Our NMR and SAXS data show that both dsRBDs act as independent structural modules in solution and are canonical, non-sequence-specific dsRBDs featuring non-canonical KKxAXK RNA binding motifs. NMR titrations combined with filter binding experiments and isothermal titration calorimetry (ITC) document the contribution of dsRBD1 to dsRNA binding in vitro. Curiously, dsRBD1 mutants in which dsRNA binding in vitro is strongly compromised do not affect roX2 RNA binding and MLE localization in cells. These data suggest alternative functions for dsRBD1 in vivo.
Assuntos
Proteínas Cromossômicas não Histona/química , DNA Helicases/química , Proteínas de Drosophila/química , Drosophila melanogaster/química , RNA Longo não Codificante/química , Fatores de Transcrição/química , Sequência de Aminoácidos , Animais , Sítios de Ligação , Proteínas Cromossômicas não Histona/genética , Proteínas Cromossômicas não Histona/metabolismo , Clonagem Molecular , DNA Helicases/genética , DNA Helicases/metabolismo , Mecanismo Genético de Compensação de Dose , Motivo de Ligação ao RNA de Cadeia Dupla , Proteínas de Drosophila/genética , Proteínas de Drosophila/metabolismo , Drosophila melanogaster/genética , Drosophila melanogaster/metabolismo , Escherichia coli/genética , Escherichia coli/metabolismo , Expressão Gênica , Regulação da Expressão Gênica , Vetores Genéticos/química , Vetores Genéticos/metabolismo , Cinética , Masculino , Modelos Moleculares , Conformação de Ácido Nucleico , Ligação Proteica , Conformação Proteica em alfa-Hélice , Conformação Proteica em Folha beta , RNA Longo não Codificante/genética , RNA Longo não Codificante/metabolismo , Proteínas Recombinantes/química , Proteínas Recombinantes/genética , Proteínas Recombinantes/metabolismo , Alinhamento de Sequência , Homologia de Sequência de Aminoácidos , Fatores de Transcrição/genética , Fatores de Transcrição/metabolismoRESUMO
The intrinsically disordered p15PAF regulates DNA replication and repair when interacting with the Proliferating Cell Nuclear Antigen (PCNA) sliding clamp. As many interactions between disordered proteins and globular partners involved in signaling and regulation, the complex between p15PAF and trimeric PCNA is of low affinity, forming a transient complex that is difficult to characterize at a structural level due to its inherent polydispersity. We have determined the structure, conformational fluctuations, and relative population of the five species that coexist in solution by combining small-angle X-ray scattering (SAXS) with molecular modelling. By using explicit ensemble descriptions for the individual species, built using integrative approaches and molecular dynamics (MD) simulations, we collectively interpreted multiple SAXS profiles as population-weighted thermodynamic mixtures. The analysis demonstrates that the N-terminus of p15PAF penetrates the PCNA ring and emerges on the back face. This observation substantiates the role of p15PAF as a drag regulating PCNA processivity during DNA repair. Our study reveals the power of ensemble-based approaches to decode structural, dynamic, and thermodynamic information from SAXS data. This strategy paves the way for deciphering the structural bases of flexible, transient and multivalent macromolecular assemblies involved in pivotal biological processes.
Assuntos
Proteínas de Transporte/química , Antígeno Nuclear de Célula em Proliferação/química , Proteínas de Ligação a DNA , Humanos , Proteínas Intrinsicamente Desordenadas/química , Simulação de Dinâmica Molecular , Complexos Multiproteicos/química , Ressonância Magnética Nuclear Biomolecular , Estrutura Quaternária de Proteína , Proteínas Recombinantes/química , Espalhamento a Baixo Ângulo , Difração de Raios XRESUMO
Small- and wide-angle x-ray scattering (SWAXS) and molecular dynamics (MD) simulations are complementary approaches that probe conformational transitions of biomolecules in solution, even in a time-resolved manner. However, the structural interpretation of the scattering signals is challenging, while MD simulations frequently suffer from incomplete sampling or from a force-field bias. To combine the advantages of both techniques, we present a method that incorporates solution scattering data as a differentiable energetic restraint into explicit-solvent MD simulations, termed SWAXS-driven MD, with the aim to direct the simulation into conformations satisfying the experimental data. Because the calculations fully rely on explicit solvent, no fitting parameters associated with the solvation layer or excluded solvent are required, and the calculations remain valid at wide angles. The complementarity of SWAXS and MD is illustrated using three biological examples, namely a periplasmic binding protein, aspartate carbamoyltransferase, and a nuclear exportin. The examples suggest that SWAXS-driven MD is capable of refining structures against SWAXS data without foreknowledge of possible reaction paths. In turn, the SWAXS data accelerates conformational transitions in MD simulations and reduces the force-field bias.
Assuntos
Aspartato Carbamoiltransferase/química , Carioferinas/química , Simulação de Dinâmica Molecular , Espalhamento a Baixo Ângulo , Difração de Raios X/métodos , Algoritmos , Sequência de Aminoácidos , Dados de Sequência Molecular , Solventes/químicaRESUMO
Wide-angle x-ray scattering (WAXS) experiments of biomolecules in solution have become increasingly popular because of technical advances in light sources and detectors. However, the structural interpretation of WAXS profiles is problematic, partly because accurate calculations of WAXS profiles from structural models have remained challenging. In this work, we present the calculation of WAXS profiles from explicit-solvent molecular dynamics (MD) simulations of five different proteins. Using only a single fitting parameter that accounts for experimental uncertainties because of the buffer subtraction and dark currents, we find excellent agreement to experimental profiles both at small and wide angles. Because explicit solvation eliminates free parameters associated with the solvation layer or the excluded solvent, which would require fitting to experimental data, we minimize the risk of overfitting. We further find that the influence from water models and protein force fields on calculated profiles are insignificant up to q≈15nm(-1). Using a series of simulations that allow increasing flexibility of the proteins, we show that incorporating thermal fluctuations into the calculations significantly improves agreement with experimental data, demonstrating the importance of protein dynamics in the interpretation of WAXS profiles. In addition, free MD simulations up to one microsecond suggest that the calculated profiles are highly sensitive with respect to minor conformational rearrangements of proteins, such as an increased flexibility of a loop or an increase of the radius of gyration by < 1%. The present study suggests that quantitative comparison between MD simulations and experimental WAXS profiles emerges as an accurate tool to validate solution ensembles of biomolecules.
Assuntos
Espectroscopia de Ressonância Magnética/métodos , Simulação de Dinâmica Molecular , Espalhamento a Baixo Ângulo , Difração de Raios X , Aldose-Cetose Isomerases/química , Sequência de Aminoácidos , Citocromos c/química , Dados de Sequência Molecular , Muramidase/química , Ubiquitina/químicaRESUMO
Background: Cystic fibrosis (CF) is caused by a wide spectrum of mutations in the CF transmembrane conductance regulator (CFTR) gene, with some leading to non-classical clinical presentations. We present an integrated in vivo, in silico and in vitro investigation of an individual with CF carrying the rare Q1291H-CFTR allele and the common F508del allele. At age 56 years, the participant had obstructive lung disease and bronchiectasis, qualifying for Elexacaftor/Tezacaftor/Ivacaftor (ETI) CFTR modulator treatment due to their F508del allele. Q1291H CFTR incurs a splicing defect, producing both a normally spliced but mutant mRNA isoform and a misspliced isoform with a premature termination codon, causing nonsense mediated decay. The effectiveness of ETI in restoring Q1291H-CFTR is largely unknown. Methods: We collected clinical endpoint measurements, including forced expiratory volume in 1 s percent predicted (FEV1pp) and body mass index (BMI), and examined medical history. In silico simulations of the Q1291H-CFTR were compared to Q1291R, G551D, and wild-type (WT)-CFTR. We quantified relative Q1291H CFTR mRNA isoform abundance in patient-derived nasal epithelial cells. Differentiated pseudostratified airway epithelial cell models at air liquid interface were created and ETI treatment impact on CFTR was assessed by electrophysiology assays and Western blot. Results: The participant ceased ETI treatment after 3 months due to adverse events and no improvement in FEV1pp or BMI. In silico simulations of Q1291H-CFTR identified impairment of ATP binding similar to known gating mutants Q1291R and G551D-CFTR. Q1291H and F508del mRNA transcripts composed 32.91% and 67.09% of total mRNA respectively, indicating 50.94% of Q1291H mRNA was misspliced and degraded. Mature Q1291H-CFTR protein expression was reduced (3.18% ± 0.60% of WT/WT) and remained unchanged with ETI. Baseline CFTR activity was minimal (3.45 ± 0.25 µA/cm2) and not enhanced with ETI (5.73 ± 0.48 µA/cm2), aligning with the individual's clinical evaluation as a non-responder to ETI. Conclusion: The combination of in silico simulations and in vitro theratyping in patient-derived cell models can effectively assess CFTR modulator efficacy for individuals with non-classical CF manifestations or rare CFTR mutations, guiding personalized treatment strategies and optimizing clinical outcomes.
RESUMO
Cystic Fibrosis (CF) results from over 400 different disease-causing mutations in the CF Transmembrane Conductance Regulator (CFTR) gene. These CFTR mutations lead to numerous defects in CFTR protein function. A novel class of targeted therapies (CFTR modulators) have been developed that can restore defects in CFTR folding and gating. This study aimed to characterize the functional and structural defects of S945L-CFTR and interrogate the efficacy of modulators with two modes of action: gating potentiator [ivacaftor (IVA)] and folding corrector [tezacaftor (TEZ)]. The response to these modulators in vitro in airway differentiated cell models created from a participant with S945L/G542X-CFTR was correlated with in vivo clinical outcomes of that participant at least 12 months pre and post modulator therapy. In this participants' airway cell models, CFTR-mediated chloride transport was assessed via ion transport electrophysiology. Monotherapy with IVA or TEZ increased CFTR activity, albeit not reaching statistical significance. Combination therapy with TEZ/IVA significantly (p = 0.02) increased CFTR activity 1.62-fold above baseline. Assessment of CFTR expression and maturation via western blot validated the presence of mature, fully glycosylated CFTR, which increased 4.1-fold in TEZ/IVA-treated cells. The in vitro S945L-CFTR response to modulator correlated with an improvement in in vivo lung function (ppFEV1) from 77.19 in the 12 months pre TEZ/IVA to 80.79 in the 12 months post TEZ/IVA. The slope of decline in ppFEV1 significantly (p = 0.02) changed in the 24 months post TEZ/IVA, becoming positive. Furthermore, there was a significant improvement in clinical parameters and a fall in sweat chloride from 68 to 28â mmol/L. The mechanism of dysfunction of S945L-CFTR was elucidated by in silico molecular dynamics (MD) simulations. S945L-CFTR caused misfolding of transmembrane helix 8 and disruption of the R domain, a CFTR domain critical to channel gating. This study showed in vitro and in silico that S945L causes both folding and gating defects in CFTR and demonstrated in vitro and in vivo that TEZ/IVA is an efficacious modulator combination to address these defects. As such, we support the utility of patient-derived cell models and MD simulations in predicting and understanding the effect of modulators on CFTR function on an individualized basis.
RESUMO
The multimeric plasma glycoprotein (GP) von Willebrand factor (VWF) is best known for recruiting platelets to sites of injury during primary hemostasis. Generally, mutations in the VWF gene lead to loss of hemostatic activity and thus the bleeding disorder von Willebrand disease. By employing cone and platelet aggregometry and microfluidic assays, we uncovered a platelet GPIIb/IIIa-dependent prothrombotic gain of function (GOF) for variant p.Pro2555Arg, located in the C4 domain, leading to an increase in platelet aggregate size. We performed complementary biophysical and structural investigations using circular dichroism spectra, small-angle X-ray scattering, nuclear magnetic resonance spectroscopy, molecular dynamics simulations on the single C4 domain, and dimeric wild-type and p.Pro2555Arg constructs. C4-p.Pro2555Arg retained the overall structural conformation with minor populations of alternative conformations exhibiting increased hinge flexibility and slow conformational exchange. The dimeric protein becomes disordered and more flexible. Our data suggest that the GOF does not affect the binding affinity of the C4 domain for GPIIb/IIIa. Instead, the increased VWF dimer flexibility enhances temporal accessibility of platelet-binding sites. Using an interdisciplinary approach, we revealed that p.Pro2555Arg is the first VWF variant, which increases platelet aggregate size and shows a shear-dependent function of the VWF stem region, which can become hyperactive through mutations. Prothrombotic GOF variants of VWF are a novel concept of a VWF-associated pathomechanism of thromboembolic events, which is of general interest to vascular health but not yet considered in diagnostics. Thus, awareness should be raised for the risk they pose. Furthermore, our data implicate the C4 domain as a novel antithrombotic drug target.
Assuntos
Mutação com Ganho de Função , Variação Genética , Complexo Glicoproteico GPIIb-IIIa de Plaquetas/metabolismo , Fator de von Willebrand/genética , Mutação com Ganho de Função/genética , Hemostasia , Humanos , Agregação Plaquetária , Domínios Proteicos/genética , Doenças de von Willebrand/sangue , Fator de von Willebrand/metabolismoRESUMO
Characterization of I37R, a mutation located in the lasso motif of the CFTR chloride channel, was conducted by theratyping several CFTR modulators from both potentiator and corrector classes. Intestinal current measurements in rectal biopsies, forskolin-induced swelling (FIS) in intestinal organoids, and short circuit current measurements in organoid-derived monolayers from an individual with I37R/F508del CFTR genotype demonstrated that the I37R-CFTR results in a residual function defect amenable to treatment with potentiators and type III, but not type I, correctors. Molecular dynamics of I37R using an extended model of the phosphorylated, ATP-bound human CFTR identified an altered lasso motif conformation which results in an unfavorable strengthening of the interactions between the lasso motif, the regulatory (R) domain, and the transmembrane domain 2 (TMD2). Structural and functional characterization of the I37R-CFTR mutation increases understanding of CFTR channel regulation and provides a potential pathway to expand drug access to CF patients with ultra-rare genotypes.
RESUMO
Free energy calculations for protein-ligand dissociation have been tested and validated for small ligands (50 atoms or less), but there has been a paucity of studies for larger, peptide-size ligands due to computational limitations. Previously we have studied the energetics of dissociation in a potassium channel-charybdotoxin complex by using umbrella sampling molecular-dynamics simulations, and established the need for carefully chosen coordinates and restraints to maintain the physiological ligand conformation. Here we address the ligand integrity problem further by constructing additional potential of mean forces for dissociation of charybdotoxin using restraints. We show that the large discrepancies in binding free energy arising from simulation artifacts can be avoided by using appropriate restraints on the ligand, which enables determination of the binding free energy within the chemical accuracy. We make several suggestions for optimal choices of harmonic potential parameters and restraints to be used in binding studies of large ligands.
Assuntos
Proteínas de Bactérias/química , Proteínas de Bactérias/metabolismo , Charibdotoxina/química , Charibdotoxina/metabolismo , Canais de Potássio/química , Canais de Potássio/metabolismo , Amidas/química , Fenômenos Biomecânicos , Modelos Moleculares , Ligação Proteica , Conformação Proteica , Propriedades de Superfície , TermodinâmicaRESUMO
Protein-protein and protein-ligand interactions often involve conformational changes or structural rearrangements that can be quantified by solution small-angle X-ray scattering (SAXS). These scattering intensity measurements reveal structural details of the bound complex, the number of species involved and, additionally, the strength of interactions if carried out as a titration. Although a core part of structural biology workflows, SAXS-based titrations are not commonly used in drug discovery contexts. This is because prior knowledge of expected sample requirements, throughput and prediction accuracy is needed to develop reliable ligand screens. This study presents the use of the histidine-binding protein (26â kDa) and other periplasmic binding proteins to benchmark ligand screen performance. Sample concentrations and exposure times were varied across multiple screening trials at four beamlines to investigate the accuracy and precision of affinity prediction. The volatility ratio between titrated scattering curves and a common apo reference is found to most reliably capture the extent of structural and population changes. This obviates the need to explicitly model scattering intensities of bound complexes, which can be strongly ligand-dependent. Where the dissociation constant is within 102 of the protein concentration and the total exposure times exceed 20â s, the titration protocol presented at 0.5â mg ml-1 yields affinities comparable to isothermal titration calorimetry measurements. Estimated throughput ranges between 20 and 100 ligand titrations per day at current synchrotron beamlines, with the limiting step imposed by sample handling and cleaning procedures.
RESUMO
Ribonucleoprotein complexes (RNPs) are central to all processes in the cell. One of the prerequisites to understand how RNPs work is to determine their high-resolution structures. With the recent revolution in cryoelectron microscopy this task has become easier for large RNP machines, such as ribosomes, spliceosomes, and polymerases. However, the transient and highly dynamic nature of many RNPs makes structure determination a challenging task. Thus, an integrative structural and molecular biology approach is required, tackling three key challenges: (1) identification of cognate RNA sequences; (2) collection of structural data by conducting X-ray crystallography, NMR, electron microscopy, small-angle scattering (SAS), and other experiments; and (3) the creation of structural models that integrates all experimental restraints. Given the breadth of expertise required, this review presents an overview of available methods and successful examples with the goal to provide readers with a selection of promising options for structure determination of RNPs.
Assuntos
Ribonucleoproteínas/química , Ribonucleoproteínas/metabolismo , Sequência de Bases , Microscopia Crioeletrônica , Cristalografia por Raios X , Modelos Moleculares , Ribonucleoproteínas/genética , Espalhamento a Baixo ÂnguloRESUMO
Ion channel-toxin complexes are ideal systems for computational studies of protein-ligand interactions, because, in most cases, the channel axis provides a natural reaction coordinate for unbinding of a ligand and a wealth of physiological data is available to check the computational results. We use a recently determined structure of a potassium channel-charybdotoxin complex in molecular dynamics simulations to investigate the mechanism and energetics of unbinding. Pairs of residues on the channel protein and charybdotoxin that are involved in the binding are identified, and their behavior is traced during umbrella-sampling simulations as charybdotoxin is moved away from the binding site. The potential of mean force for the unbinding of charybdotoxin is constructed from the umbrella sampling simulations using the weighted histogram analysis method, and barriers observed are correlated with specific breaking of interactions and influx of water molecules into the binding site. Charybdotoxin is found to undergo conformational changes as a result of the reaction coordinate choice--a nontrivial decision for larger ligands--which we explore in detail, and for which we propose solutions. Agreement between the calculated and the experimental binding energies is obtained once the energetic consequences of these conformational changes are included in the calculations.
Assuntos
Charibdotoxina/metabolismo , Modelos Moleculares , Canais de Potássio/metabolismo , Charibdotoxina/química , Ligantes , Espectroscopia de Ressonância Magnética , Canais de Potássio/química , Ligação Proteica , Conformação Proteica , Termodinâmica , Água/metabolismoRESUMO
Small-angle X-ray and small-angle neutron scattering (SAXS/SANS) provide unique structural information on biomolecules and their complexes in solution. SANS may provide multiple independent data sets by means of contrast variation experiments, that is, by measuring at different D2O concentrations and different perdeuteration conditions of the biomolecular complex. However, even the combined data from multiple SAXS/SANS sets is by far insufficient to define all degrees of freedom of a complex, leading to a significant risk of overfitting when refining biomolecular structures against SAXS/SANS data. Hence, to control against overfitting, the low-information SAXS/SANS data must be complemented by accurate physical models, and, if possible, refined models should be cross-validated against independent data not used during the refinement. We present a method for refining atomic biomolecular structures against multiple sets of SAXS and SANS data using all-atom molecular dynamics simulations. Using the protein citrate synthase and the protein/RNA complex Sxl-Unr-msl2 mRNA as test cases, we demonstrate how multiple SAXS and SANS sets may be used for refinement and cross-validation, thereby excluding overfitting during refinement. For the Sxl-Unr-msl2 complex, we find that perdeuteration of the Unr domain leads to a unique, slightly compacted conformation, whereas other perdeuteration conditions lead to similar solution conformations compared to the nondeuterated state. In line with our previous method for predicting SAXS curves, SANS curves were predicted with explicit-solvent calculations, taking atomic models for both the hydration layer and the excluded solvent into account, thereby avoiding the use of solvent-related fitting parameters and solvent-reduced neutron scattering lengths. We expect the method to be useful for deriving and validating solution structures of biomolecules and soft-matter complexes, and for critically assessing whether multiple SAXS and SANS sets are mutually compatible.
Assuntos
Simulação de Dinâmica Molecular , Difração de Nêutrons , Espalhamento a Baixo Ângulo , Difração de Raios X , Animais , Citrato (si)-Sintase/química , Proteínas de Ligação a DNA/química , Drosophila , Proteínas de Drosophila/química , Difração de Nêutrons/métodos , RNA/química , Proteínas de Ligação a RNA/química , Suínos , Difração de Raios X/métodosRESUMO
Steady-state evoked potentials (SSEPs) elicited by sinusoidal stimuli are predicted from a physiologically-based model, including bielectrode and volume conduction effects. Comparison with visual SSEPs yields constraints on phase and latency of the retinothalamic transfer function that are consistent with experiment. Predictions of phase velocities measured as SSEPs cross the cortex are consistent with low values measured for slow waves in sleep, while resonant behavior induced by corticothalamic loops, especially near the alpha peak, contributes to wide scatter in waking-state phase velocity measurements comparable to effects from volume conduction. The common use of bielectrode derivations to compensate for volume conduction effects is examined and shown to be incomplete, tending to lead to underestimates of phase velocity, especially at low frequencies and near the alpha peak, due to incorrect elimination of true long-wavelength contributions to the SSEP.
Assuntos
Córtex Cerebral/fisiologia , Potenciais Evocados/fisiologia , Modelos Neurológicos , Animais , Eletroencefalografia/métodos , HumanosRESUMO
In many biomolecular interactions, changes in the assembly states and structural conformations of participants can act as a complementary reporter of binding to functional and thermodynamic assays. This structural information is captured by a number of structural biology and biophysical techniques that are viable either as primary screens in small-scale applications or as secondary screens to complement higher throughput methods. In particular, small-angle X-ray scattering (SAXS) reports the average distance distribution between all atoms after orientational averaging. Such information is important when for example investigating conformational changes involved in inhibitory and regulatory mechanisms where binding events do not necessarily cause functional changes. Thus, we summarise here the current and prospective capabilities of SAXS-based screening in the context of other methods that yield structural information. Broad guidelines are also provided to assist readers in preparing screening protocols that are tailored to available X-ray sources.
RESUMO
1H-15N NMR spin relaxation and relaxation dispersion experiments can reveal the time scale and extent of protein motions across the ps-ms range, where the ps-ns dynamics revealed by fundamental quantities R1, R2, and heteronuclear NOE can be well-sampled by molecular dynamics simulations (MD). Although the principles of relaxation prediction from simulations are well-established, numerous NMR-MD comparisons have hitherto focused upon the aspect of order parameters S2 due to common artifacts in the prediction of transient dynamics. We therefore summarize here all necessary components and highlight existing and proposed solutions, such as the inclusion of quantum mechanical zero-point vibrational corrections and separate MD convergence of global and local motions in coarse-grained and all-atom force fields, respectively. For the accuracy of the MD prediction to be tested, two model proteins GB3 and Ubiquitin are used to validate five atomistic force fields against published NMR data supplemented by the coarse-grained force field MARTINI+EN. In Amber and CHARMM-type force fields, quantitative agreement was achieved for structured elements with minimum adjustment of global parameters. Deviations from experiment occur in flexible loops and termini, indicating differences in both the extent and time scale of backbone motions. The lack of systematic patterns and water model dependence suggests that modeling of the local environment limits prediction accuracy. Nevertheless, qualitative accuracy in a 2 µs CHARMM36m Stam2 VHS domain simulation demonstrates the potential of MD-based interpretation in combination with NMR-measured dynamics, increasing the utility of spin relaxation in integrative structural biology.
RESUMO
We present a screening protocol utilizing small-angle X-ray scattering (SAXS) to obtain structural information on biomolecular interactions independent of prior knowledge, so as to complement affinity-based screening and provide leads for further exploration. This protocol categorizes ligand titrations by computing pairwise agreement between curves, and separately estimates affinities by quantifying complex formation as a departure from the linear sum properties of solution SAXS. The protocol is validated by sparse sequence search around the native poly uridine RNA motifs of the two-RRM domain Sex-lethal protein (Sxl). The screening of 35 RNA motifs between 4 to 10 nucleotides reveals a strong variation of resulting complexes, revealed to be preference-switching between 1:1 and 2:2 binding stoichiometries upon addition of structural modeling. Validation of select sequences in isothermal calorimetry and NMR titration retrieves domain-specific roles and function of a guanine anchor. These findings reinforce the suitability of SAXS as a complement in lead identification.