Q1291H-CFTR molecular dynamics simulations and ex vivo theratyping in nasal epithelial models and clinical response to elexacaftor/tezacaftor/ivacaftor in a Q1291H/F508del patient.

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.

S945L-CFTR molecular dynamics, functional characterization and tezacaftor/ivacaftor efficacy in vivo and in vitro in matched pediatric patient-derived cell models.

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.

Structure and dynamics of the von Willebrand Factor C6 domain.

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.

Molecular Dynamics and Theratyping in Airway and Gut Organoids Reveal R352Q-CFTR Conductance Defect.

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.

Molecular dynamics and functional characterization of I37R-CFTR lasso mutation provide insights into channel gating activity.

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.

Gain-of-Function Variant p.Pro2555Arg of von Willebrand Factor Increases Aggregate Size through Altering Stem Dynamics.

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.

Structure-based screening of binding affinities via small-angle X-ray scattering.

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.

Integrative Structural Biology of Protein-RNA Complexes.

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.

Combined Small-Angle X-ray and Neutron Scattering Restraints in Molecular Dynamics Simulations.

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.

Structure, dynamics and roX2-lncRNA binding of tandem double-stranded RNA binding domains dsRBD1,2 of Drosophila helicase Maleless.

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.

Structure and dynamics of the platelet integrin-binding C4 domain of von Willebrand factor.

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.

The role of small-angle scattering in structure-based screening applications.

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.

A General Small-Angle X-ray Scattering-Based Screening Protocol Validated for Protein-RNA Interactions.

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.

Ab Initio Prediction of NMR Spin Relaxation Parameters from Molecular Dynamics Simulations.

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.

Computing the Rotational Diffusion of Biomolecules via Molecular Dynamics Simulation and Quaternion Orientations.

Rotational diffusion (Drot) is a fundamental property of biomolecules that contains information about molecular dimensions and solute-solvent interactions. While ab initio Drot prediction can be achieved by explicit all-atom molecular dynamics simulations, this is hindered by both computational expense and limitations in water models. We propose coarse-grained force fields as a complementary solution, and show that the MARTINI force field with elastic networks is sufficient to compute Drot in >10 proteins spanning 5-157 kDa. We also adopt a quaternion-based approach that computes Drot orientation directly from autocorrelations of best-fit rotations as used in, e.g., RMSD algorithms. Over 2 µs trajectories, isotropic MARTINI+EN tumbling replicates experimental values to within 10-20%, with convergence analyses suggesting a minimum sampling of >50 × τtheor to achieve sufficient precision. Transient fluctuations in anisotropic tumbling cause decreased precision in predictions of axisymmetric anisotropy and rhombicity, the latter of which cannot be precisely evaluated within 2000 × τtheor for GB3. Thus, we encourage reporting of axial decompositions Dx, Dy, Dz to ease comparability between experiment and simulation. Where protein disorder is absent, we observe close replication of MARTINI+EN Drot orientations versus CHARMM22*/TIP3p and experimental data. This work anticipates the ab initio prediction of NMR-relaxation by combining coarse-grained global motions with all-atom local motions.

Disentangling polydispersity in the PCNA-p15PAF complex, a disordered, transient and multivalent macromolecular assembly.

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.

Length of Hospital Stay Prediction at the Admission Stage for Cardiology Patients Using Artificial Neural Network.

For hospitals' admission management, the ability to predict length of stay (LOS) as early as in the preadmission stage might be helpful to monitor the quality of inpatient care. This study is to develop artificial neural network (ANN) models to predict LOS for inpatients with one of the three primary diagnoses: coronary atherosclerosis (CAS), heart failure (HF), and acute myocardial infarction (AMI) in a cardiovascular unit in a Christian hospital in Taipei, Taiwan. A total of 2,377 cardiology patients discharged between October 1, 2010, and December 31, 2011, were analyzed. Using ANN or linear regression model was able to predict correctly for 88.07% to 89.95% CAS patients at the predischarge stage and for 88.31% to 91.53% at the preadmission stage. For AMI or HF patients, the accuracy ranged from 64.12% to 66.78% at the predischarge stage and 63.69% to 67.47% at the preadmission stage when a tolerance of 2 days was allowed.

Structural Properties of Protein-Detergent Complexes from SAXS and MD Simulations.

In experimental studies of solubilized membrane proteins, the detergent corona influences the protein behavior and the resulting measurement. Thus, combinations of experimental techniques with atomistic modeling have been used to resolve corona structural parameters and distributions. Here, we used small-angle X-ray scattering (SAXS) data and molecular dynamics simulations to study a model protein-detergent complex (PDC) consisting of aquaporin-0 and dodecyl-ß-maltoside molecules (ßDDM). The corona morphology of single snapshots was found to be rough, but it is smooth and compacted in 100-ns-scale ensemble averages. Individual snapshots therefore were unable to accurately represent the ensemble information as captured by experimental SAXS. Mimicking of annular lipids by detergent was also observed. SAXS prediction using different published methods was used to identify optimal ßDDM numbers. Explicit-solvent methods predicted best agreement using 290-ßDDM PDCs, but implicit-solvent methods gave unclear predictions due to overcompensation by free solvation-layer density parameters. Thus, ensemble-based approaches and physically motivated constraints will help to extract structural information from SAXS data.

Interpretation of solution x-ray scattering by explicit-solvent molecular dynamics.

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.

Validating solution ensembles from molecular dynamics simulation by wide-angle X-ray scattering data.

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.