MAVIS: merging, annotation, validation, and illustration of structural variants.

Summary: Reliably identifying genomic rearrangements and interpreting their impact is a key step in understanding their role in human cancers and inherited genetic diseases. Many short read algorithmic approaches exist but all have appreciable false negative rates. A common approach is to evaluate the union of multiple tools increasing sensitivity, followed by filtering to retain specificity. Here we describe an application framework for the rapid generation of structural variant consensus, unique in its ability to visualize the genetic impact and context as well as process both genome and transcriptome data. Availability and implementation: http://mavis.bcgsc.ca. Supplementary information: Supplementary data are available at Bioinformatics online.

A platform for oncogenomic reporting and interpretation.

Manual interpretation of variants remains rate limiting in precision oncology. The increasing scale and complexity of molecular data generated from comprehensive sequencing of cancer samples requires advanced interpretative platforms as precision oncology expands beyond individual patients to entire populations. To address this unmet need, we introduce a Platform for Oncogenomic Reporting and Interpretation (PORI), comprising an analytic framework that facilitates the interpretation and reporting of somatic variants in cancer. PORI integrates reporting and graph knowledge base tools combined with support for manual curation at the reporting stage. PORI represents an open-source platform alternative to commercial reporting solutions suitable for comprehensive genomic data sets in precision oncology. We demonstrate the utility of PORI by matching 9,961 pan-cancer genome atlas tumours to the graph knowledge base, calculating therapeutically informative alterations, and making available reports describing select individual samples.

Investigation of the binding of a carbohydrate-mimetic peptide to its complementary anticarbohydrate antibody by STD-NMR spectroscopy and molecular-dynamics simulations.

Saturation transfer difference (STD)-NMR spectroscopy was used to probe experimentally the bioactive solution conformation of the carbohydrate mimic MDWNMHAA 1 of the O-polysaccharide of Shigella flexneri Y when bound to its complementary antibody, mAb SYA/J6. Molecular dynamics simulations using the ZymeCAD™ Molecular Dynamics platform were also undertaken to give a more accurate picture of the conformational flexibility and the possibilities for bound ligand conformations. The ligand topology, or the dynamic epitope, was mapped with the CORCEMA-ST (COmplete Relaxation and Conformational Exchange Matrix Analysis of Saturation Transfer) program that calculates a total matrix analysis of relaxation and exchange effects to generate predicted STD-NMR intensities from simulation. The comparison of these predicted STD enhancements with experimental data was used to select a representative binding mode. A protocol that employed theoretical STD effects calculated at snapshots during the entire course of a molecular dynamics (MD) trajectory of the peptide bound to the Fv portion of the antibody, and not the averaged atomic positions of receptor-ligand complexes, was also examined. In addition, the R factor was calculated on the basis of STD (fit) to avoid T1 bias, and an effective R factor, R(eff), was defined such that if the calculated STD (fit) for proton k was within error of the experimental STD (fit) for proton k, then that calculated STD (fit) for proton k was not included in the calculation of the R factor. This protocol was effective in deriving the antibody-bound solution conformation of the peptide which also differed from the bound conformation determined by X-ray crystallography; however, several discrepancies between experimental and calculated STD (fit) values were observed. The bound conformation was therefore further refined with a simulated annealing refinement protocol known as STD-NMR intensity-restrained CORCEMA optimization (SICO) to give a more accurate representation of the bound peptide epitope. Further optimization was required in this case, but a satisfactory correlation between experimental and calculated STD values was obtained. Attempts were also made to obtain STD enhancements with a synthetic pentasaccharide hapten, corresponding to the O-polysaccharide, while bound to the antibody. However, unfavorable kinetics of binding in this system prevented sufficient STD build-up, which, in turn, hindered a rigorous analysis via full STD build-up curves.

WaterLOGSY NMR experiments in conjunction with molecular-dynamics simulations identify immobilized water molecules that bridge peptide mimic MDWNMHAA to anticarbohydrate antibody SYA/J6.

X-ray crystallographic data of the carbohydrate mimic MDWNMHAA when bound to an anti-Shigella flexneri Y mAb SYA/J6 indicate the immobilization of water molecules, that is, the presence of "bound" waters, in the active site. Water Ligand Observed via Gradient Spectroscopy (WaterLOGSY) was used in conjunction with saturation transfer difference (STD)-NMR spectroscopy to probe the existence of immobilized water molecules in the complex of MDWNMHAA 1 bound to mAb SYA/J6. Molecular dynamics simulations using the ZymeCAD Molecular Dynamics platform were then used to specify the likely locations of these water molecules. Of note, those waters involved in providing complementarity between the peptide and mAb SYA/J6 remained throughout the course of the simulation. Together, the experimental and computational protocols have been used to identify the bound water molecules present in the antibody-peptide complex.

Substrate directs enzyme dynamics by bridging distal sites: UDP-galactopyranose mutase.

UDP-Galactopyranose mutase (UGM) is a flavoenzyme that catalyzes interconversion of UDP-galactopyranose (UDP-Galp) and UDP-galactofuranose (UDP-Galf); its activity depends on FAD redox state. The enzyme is vital to many pathogens, not native to mammals, and is an important drug target. We have probed binding of substrate, UDP-Galp, and UDP to wild type and W160A UGM from K. pneumoniae, and propose that substrate directs recognition loop dynamics by bridging distal FAD and W160 sites; W160 interacts with uracil of the substrate and is functionally essential. Enhanced Trp fluorescence upon substrate binding to UGM indicates conformational changes remote from the binding site because the fluorescence is unchanged upon binding to W70F/W290F UGM where W160 is the sole Trp. MD simulations map these changes to recognition loop closure to coordinate substrate. This requires galactose-FAD interactions as Trp fluorescence is unchanged upon substrate binding to oxidized UGM, or binding of UDP to either form of the enzyme, and MD show heightened recognition loop mobility in complexes with UDP. Consistent with substrate-directed loop closure, UDP binds 10-fold more tightly to oxidized UGM, yet substrate binds tighter to reduced UGM. This requires the W160-U interaction because redox-switched binding affinity of substrate reverses in the W160A mutant where it only binds when oxidized. Without the anchoring W160-U interaction, an alternative binding mode for UDP is detected, and STD-NMR experiments show simultaneous binding of UDP-Galp and UDP to different subsites in oxidized W160A UGM: Substrate no longer directs recognition loop dynamics to coordinate tight binding to the reduced enzyme.

Investigation of binding of UDP-Galf and UDP-[3-F]Galf to UDP-galactopyranose mutase by STD-NMR spectroscopy, molecular dynamics, and CORCEMA-ST calculations.

UDP-galactopyranose mutase (UGM) is the key enzyme involved in the biosynthesis of Galf. UDP-Galp and UDP-Galf are two natural substrates of UGM. A protocol that combines the use of STD-NMR spectroscopy, molecular modeling, and CORCEMA-ST calculations was applied to the investigation of the binding of UDP-Galf and its C3-fluorinated analogue to UGM from Klebsiella pneumoniae. UDP-Galf and UDP-[3-F]Galf were bound to UGM in a manner similar to that of UDP-Galp. The interconversions of UDP-Galf and UDP-[3-F]Galf to their galactopyranose counterparts were catalyzed by the reduced (active) UGM with different catalytic efficiencies, as observed by NMR spectroscopy. The binding affinities of UDP-Galf and UDP-[3-F]Galf were also compared with those of UDP-Galp and UDP by competition STD-NMR experiments. When UGM was in the oxidized (inactive) state, the binding affinities of UDP-Galf, UDP-Galp, and UDP-[3-F]Galf were of similar magnitudes and were lower than that of UDP. However, when UGM was in the reduced state, UDP-Galp had higher binding affinity compared with UDP. Molecular dynamics (MD) simulations indicated that the "open" mobile loop in UGM "closes" upon binding of the substrates. Combined MD simulations and STD-NMR experiments were used to create models of UGM with UDP-Galf and UDP-[3-F]Galf as bound ligands. Calculated values of saturation-transfer effects with CORCEMA-ST (complete relaxation and conformational exchange matrix analysis of saturation transfer) were compared to the experimental STD effects and permitted differentiation between two main conformational families of the bound ligands. Taken together, these results are used to rationalize the different rates of catalytic turnover of UDP-Galf and UDP-[3-F]Galf and shed light on the mechanism of action of UGM.

Can PISEMA experiments be used to extract structural parameters for mobile beta-barrels?

The effect of mobility on 15N chemical shift/15N-(1)H dipolar coupling (PISEMA) solid state NMR experiments applied to macroscopically oriented beta-barrels is assessed using molecular dynamics simulation data of the NalP autotransporter domain embedded in a DMPC bilayer. In agreement with previous findings for alpha-helices, the fast librational motion of the peptide planes is found to have a considerable effect on the calculated PISEMA spectra. In addition, the dependence of the chemical shift anisotropy (CSA) and dipolar coupling parameters on the calculated spectra is evaluated specifically for the beta-barrel case. It is found that the precise choice of the value of the CSA parameters sigma11, sigma22 and sigma33 has only a minor effect, whereas the choice of the CSA parameter theta shifts the position of the peaks by up to 20 ppm and changes the overall shape of the spectrum significantly. As was found for alpha-helices, the choice of the NH bond distance has a large effect on the dipolar coupling constant used for the calculations. Overall, it is found that the alternating beta-strands in the barrel occupy distinct regions of the PISEMA spectra, forming patterns which may prove useful in peak assignment.