Skip to content

Latest commit

 

History

History
566 lines (489 loc) · 18 KB

nucl_canon_trj_preview.md

File metadata and controls

566 lines (489 loc) · 18 KB
Raw

NcH2Arun_1 - Canonical nucleosome with full tails (PDB ID 5AV9)

Back

<style> .slidecontainer { width: 100%; } 
  .slider {
    -webkit-appearance: none;
    width: 100%;
    height: 10px;
    background: #d3d3d3;
    outline: none;
    opacity: 0.7;
    -webkit-transition: .2s;
    transition: opacity .2s;
  }

  .slider:hover {
    opacity: 1;
  }

  .slider::-webkit-slider-thumb {
    -webkit-appearance: none;
    appearance: none;
    width: 15px;
    height: 15px;
    background: #222222;
    cursor: pointer;
  }

  .slider::-moz-range-thumb {
    width: 25px;
    height: 25px;
    background: #4CAF50;
    cursor: pointer;
  }

</style>
<script src="https://d3js.org/d3.v4.js"></script> <script src="https://unpkg.com/[email protected]/dist/ngl.js"></script> <script src="https://code.jquery.com/jquery-3.5.1.min.js" integrity="sha256-9/aliU8dGd2tb6OSsuzixeV4y/faTqgFtohetphbbj0=" crossorigin="anonymous"></script> <script> var pdb="trj/canon_nucl_for_web.pdb" var xtc="trj/canon_nucl_for_web.xtc" var csvfile="dat/cH2A_dist_unwrap.csv" var trjstep = 10; $(document).ready(function() { window.stage = new NGL.Stage("viewport0", { backgroundColor: "#FFFFFF" }); window.stage.loadFile(pdb).then(function(nucl) { var aspectRatio = 2; var radius = 1.5; var radiusScale = 4; 
  var hyper_scheme = NGL.ColormakerRegistry.addSelectionScheme([
    ["orange", ".CA"],
    ["blue", "_N"],
    ["red", "_O"],
    ["grey", "*"]
  ], "DA");
  var residues = NGL.ColormakerRegistry.addSelectionScheme([
    ["blue", "40 and ARG"],
["lightcyan", "42 and ARG"],
["cyan", "49 and ARG"],
    ["green", "41 and TYR"],
["pink", "39 and HIS"],
    ["white", "*"]
  ], "protein");
  
  var shape = new NGL.Shape( "Axes" );
		shape.addArrow( [ 0, 0, -35 ], [ 0, 0, 35 ], [ 0.04, 0.8, 0.03 ], 2.0 );
  shape.addArrow( [ -58, 0, 0 ], [ 58, 0, 0 ], [ 0.8,0.06 ,0.102  ], 2.0 );
  shape.addArrow( [ 0, -55, 0 ], [ 0, 55, 0 ], [ 0.3, 0.3, 0.3 ], 2.0 );
  window.axes = stage.addComponentFromObject( shape );
  window.axes.addRepresentation( "buffer" );
  window.axes.autoView();
  window.axes.setVisibility(false);
  stage.animationControls.rotate([ 0, 1, 0, 0 ],0);
  stage.setParameters({cameraType: "orthographic"});

  

  window.arg_lys_selection = nucl.addRepresentation('hyperball', {
    "sele": "(ARG or LYS) and not _H",
    color: hyper_scheme,
    radius: 3.5
  });
  window.arg_lys_selection.setVisibility(false);
  
  window.dna_latch_selection = nucl.addRepresentation('hyperball', {
    "sele": "39-49 and (:A or :E) and not _H",
    color: residues,
    radius: 3.5
  });
  window.dna_latch_selection.setVisibility(false);
  
  //H3 39-49 - назвать H3 39-49 DNA latch

    window.replaced_AA_selection = nucl.addRepresentation('hyperball', {
    "sele": "(10 or 16 or 40 or 99 or 123 or 124 or 126 or 128) and (:C or :G) and not _H",
    color: hyper_scheme,
    radius: 3.5
  });
  window.replaced_AA_selection.setVisibility(false);


  nucl.addRepresentation('cartoon', {
    "sele": ":A :E",
    "color": 0x020AED,
    "aspectRatio": aspectRatio,
    'radiusScale': radiusScale,
    'radiusType': 'sstruc',
    "capped": true,
    'subdiv': 10,
    'diffuseInterior': false,
    'useInteriorColor': false
  });
  nucl.addRepresentation('cartoon', {
    "sele": ":B :F",
    "color": "green",
    "aspectRatio": aspectRatio,
    'radiusScale': radiusScale,
    'radiusType': 'sstruc',
    "capped": true,
    'subdiv': 10,
    'diffuseInterior': false,
    'useInteriorColor': false
  });
  nucl.addRepresentation('cartoon', {
    "sele": ":C :G",
    "color": 0xE0F705,
    "aspectRatio": aspectRatio,
    'radiusScale': radiusScale,
    'radiusType': 'sstruc',
    "capped": true,
    'subdiv': 10,
    'diffuseInterior': false,
    'useInteriorColor': false
  });
  nucl.addRepresentation('cartoon', {
    "sele": ":D :H",
    "color": 0xCE0000,
    "aspectRatio": aspectRatio,
    'radiusScale': radiusScale,
    'radiusType': 'sstruc',
    "capped": true,
    'subdiv': 10,
    'diffuseInterior': false,
    'useInteriorColor': false
  });
  nucl.addRepresentation('spacefill', {
    "sele": ":I and -72 and .C1'",
    "color": "steelblue",
    "radius":5,
    'diffuseInterior': false,
    'useInteriorColor': false
  });
  nucl.addRepresentation('spacefill', {
    "sele": ":J and -72 and .C1'",
    "color": "orange",
    "radius":5,
    'diffuseInterior': false,
    'useInteriorColor': false
  });
  window.nucl_cartoon = nucl.addRepresentation('cartoon', {
    "sele": "nucleic",
    "color": 'grey',
    "aspectRatio": aspectRatio,
    "radius": radius,
    "radiusSegments": 1,
    "capped": 0
  });
  window.nucl_base = nucl.addRepresentation('base', {
    "sele": "nucleic",
    "color": 'grey'
  });

  NGL.autoLoad(xtc).then(function(frames) {
    nucl.addTrajectory(frames);
    window.traj = nucl.trajList[0].trajectory;
    var player = new NGL.TrajectoryPlayer(window.traj, {
      start: 0,
      timeout: 1,
      mode: "once",
      interpolateType: "spline",
      step: 1,
      interpolateStep: 5
    });
    window.traj.signals.frameChanged.add(function() {
      var fnum = window.traj.currentFrame;
      $('#myRange')[0].value = fnum;
      $("#frame_counter")[0].innerHTML = (fnum * trjstep).toFixed(2);
      tooltipLine.attr('stroke', 'black')
        .attr('x1', x(fnum))
        .attr('x2', x(fnum))
        .attr('y1', 0)
        .attr('y2', height);
    });



    //window.player.play();
    $('#myRange')[0].setAttribute('max', window.traj.frames.length - 1)

  });
  nucl.autoView();

});


var slider = document.getElementById("myRange");
var output = document.getElementById("frame_counter");
output.innerHTML = slider.value;
window.slider = slider;

slider.oninput = function() {
  //output.innerHTML = (this.value*trjstep).toFixed(2);
  window.traj.player.pause();
  window.traj.setFrame(this.value);
  //tooltipLine.attr('stroke', 'black')
  //  .attr('x1', x(this.value))
  //  .attr('x2', x(this.value))
  //  .attr('y1', 0)
  //  .attr('y2', height);
}


$('input[type=checkbox][name=ref_str_check]').on('change', toggle_reference_structure);
$('input[type=checkbox][name=arg_lys_check]').on('change', toggle_lys_arg_visibility);
$('input[type=checkbox][name=latch_check]').on('change', toggle_latch_visibility);
$('input[type=checkbox][name=highlight_DA_check]').on('change', toggle_DA_highlight);
$('input[type=checkbox][name=ortho_check]').on('change', toggle_orthographic);
$('input[type=checkbox][name=axes_check]').on('change', toggle_axes_visibility);
$('input[type=checkbox][name=replaced_AA_check]').on('change', toggle_replaced_AA_visibility);



function load_reference_structure() {
  window.ref_str_mol = window.stage.loadFile(pdb).then(function(nucl) {
    window.ref_str_mol = nucl;
    window.ref_str_repr_prot = nucl.addRepresentation('cartoon', {
      "sele": 'protein',
      "color": '#29d6d9',
      "aspectRatio": 2,
      'radiusScale': 4.1,
      'radiusType': 'sstruc',
      "aspectRatio": 2,
      "radiusSegments": 1,
      "capped": 0

    });
    window.ref_str_repr_nucl = nucl.addRepresentation('cartoon', {
      "sele": 'nucleic',
      "color": '#29d6d9',
      "aspectRatio": 2,
      "radius": 1.51,
      "radiusSegments": 1,
      "capped": 0
    });
    window.ref_str_repr_base = nucl.addRepresentation('base', {
      "sele": 'nucleic',
      "color": '#29d6d9'
    });
  })
}

function toggle_reference_structure() {
  var state = $(this).is(":checked");
  if (typeof window.ref_str_repr_prot == "undefined") {
    load_reference_structure();

  } else {
    window.ref_str_repr_prot.setVisibility(state);
    window.ref_str_repr_nucl.setVisibility(state);
    window.ref_str_repr_base.setVisibility(state)
  }


}

function toggle_DA_highlight() {
  var state = $(this).is(":checked");
  if (state === false) {
    window.nucl_cartoon.setColor('grey')
    window.nucl_base.setColor('grey')
  } else {
    var scheme = NGL.ColormakerRegistry.addSelectionScheme([
      ["pink", "DA"],
      ["grey", "*"]
    ], "DA");
    window.nucl_cartoon.setColor(scheme)
    window.nucl_base.setColor(scheme)

  }
}

function toggle_lys_arg_visibility() {
  var state = $(this).is(":checked");
  window.arg_lys_selection.setVisibility(state);
}
function toggle_latch_visibility() {
  var state = $(this).is(":checked");
  window.dna_latch_selection.setVisibility(state);
}

function toggle_replaced_AA_visibility() {
  var state = $(this).is(":checked");
  window.replaced_AA_selection.setVisibility(state);
}

function toggle_axes_visibility() {
  var state = $(this).is(":checked");
  window.axes.setVisibility(state);
}

function toggle_orthographic() {
  var state = $(this).is(":checked");
  if ($(this).is(":checked")){
  	window.stage.setParameters({cameraType: "orthographic"})
    } else {
    window.stage.setParameters({cameraType: "perspective"})
    }
}


var margin = {
    top: 10,
    right: 30,
    bottom: 40,
    left: 60
  },
  width = 800 - margin.left - margin.right,
  height = 200 - margin.top - margin.bottom;

// append the svg object to the body of the page
var svg = d3.select("#my_dataviz")
  .append("svg")
  .attr("width", width + margin.left + margin.right)
  .attr("height", height + margin.top + margin.bottom)
  .append("g")
  .attr("transform",
    "translate(" + margin.left + "," + margin.top + ")");
const tooltipLine = svg.append('line');
var x = d3.scaleLinear();
var y = d3.scaleLinear();

 //Read the data
d3.csv(csvfile,

  // When reading the csv, I must format variables:

  // Now I can use this dataset:
  function(data) {
    data.forEach(function(d) {
      d.Frame = d.Frame/10;
    });
    // Add X axis --> it is a date format

    x.domain([0, d3.max(data, function(d) {
        return +d.Frame;
      })])
      .range([0, width])

    svg.append("g")
      .attr("transform", "translate(0," + height + ")")
      .attr("class", "axis")
      .call(d3.axisBottom(x)
        .tickFormat(function(d) {
          return d*10;
        }))

    // Add Y axis

    y.domain([0, Math.max(d3.max(data, function(d) {
        return +d.prox;
      }),
      d3.max(data, function(d) {
        return +d.dist;
      }))
      ])
      .range([height, 0]);
    svg.append("g")
      .call(d3.axisLeft(y));


    // Add the line
    svg.append("path")
      .datum(data)
      .attr("fill", "none")
      .attr("stroke", "steelblue")
      .attr("opacity", 0.4)
      .attr("stroke-width", 2)
      .attr("d", d3.line()
        .x(function(d) {
          return x(d.Frame)
        })
        .y(function(d) {
          return y(d.prox)
        })
      )
    svg.append("path")
      .datum(data)
      .attr("fill", "none")
      .attr("stroke", "orange")
      .attr("opacity", 0.4)
      .attr("stroke-width", 2)
      .attr("d", d3.line()
        .x(function(d) {
          return x(d.Frame)
        })
        .y(function(d) {
          return y(d.dist)
        })
      )
     svg.append("path")
      .datum(data)
      .attr("fill", "none")
      .attr("stroke", "steelblue")
      .attr("stroke-width", 3)
      .attr("d", d3.line()
        .x(function(d) {
          return x(d.Frame)
        })
        .y(function(d) {
          return y(d.prox_filtered)
        })
      )
     svg.append("path")
      .datum(data)
      .attr("fill", "none")
      .attr("stroke", "orange")
      .attr("stroke-width", 3)
      .attr("d", d3.line()
        .x(function(d) {
          return x(d.Frame)
        })
        .y(function(d) {
          return y(d.dist_filtered)
        })
      )
      svg.append("text")
      .attr("class", "x label")
      .attr("text-anchor", "end")
      .attr("x", width-width/2)
      .attr("y", height + 35)
      .text("Time, ns");

      svg.append("text")
      .attr("class", "y label")
      .attr("text-anchor", "end")
      .attr("y", -45)
      .attr("dy", ".75em")
      .attr("transform", "rotate(-90)")
      .text("Unwrapped base pairs");

    tipBox = svg.append('rect')
      .attr('width', width)
      .attr('height', height)
      .attr('opacity', 0)
      .on('mousemove', drawTooltip)
      .on('mouseout', removeTooltip);
    const tooltip = d3.select('#tooltip');


    function removeTooltip() {
      if (tooltip) tooltip.style('display', 'none');
    }

    function drawTooltip() {
      const frame = Math.floor((x.invert(d3.mouse(tipBox.node())[0])));
      window.traj.player.pause();
      window.traj.setFrame(frame);

      tooltipLine.attr('stroke', 'black')
        .attr('x1', x(frame))
        .attr('x2', x(frame))
        .attr('y1', 0)
        .attr('y2', height);
      tooltip.html('Proximal unwrap (smoothed): ' + Math.floor(data[frame].prox_filtered) + 'bp<br>Distal unwrap (smoothed): ' + Math.floor(data[frame].dist_filtered) + 'bp')
        .style('display', 'block')
        .style('left', d3.event.pageX + 20)
        .style('top', d3.event.pageY - 20)
      /*        .selectAll()
              .data(data[frame])
              .append('div')
              .html(d => 1 ); */
    }

 })

});

</script>

H3

H4

H2A

H2B

DNA


Blue and orange balls mark distal and proximal entry/exit bais pairs.

Show starting state 
<input class="form-check-input " type="checkbox" name="ortho_check" value="" id="ortho_check" checked="true">
<label class="form-check-label " for="ortho_check">
  Orthographic
</label>

<input class="form-check-input " type="checkbox" name="axes_check" value="" id="axes_check">
<label class="form-check-label " for="ortho_check">
  Show axes
</label>



<input class="form-check-input " type="checkbox" name="arg_lys_check" value="" id="arg_lys_check">
<label class="form-check-label " for="arg_lys_check">
  Show ARG LYS
</label> 
<!--
<input class="form-check-input " type="checkbox" name="latch_check" value="" id="latch_check">
<label class="form-check-label " for="latch_check">
  Show H3 39-49 DNA latch
</label> -->

<br>
<input class="form-check-input " type="checkbox" name="replaced_AA_check" value="" id="replaced_AA_check">
<label class="form-check-label " for="replaced_AA_check">
  Show H2A.J substitutions
</label>
    
<input class="form-check-input " type="checkbox" name="highlight_DA_check" value="" id="highlight_DA_check">
<label class="form-check-label " for="highlight_DA_check">
  Highlight ADE
</label>


<div id="viewport0" style="height:500px; border: thin solid black"></div>
<div class="slidecontainer">
  <button type="submit" class="btn" name="play_button" data-toggle="button" id='play' onclick='window.traj.player.play();'>Play</button>
  <button type="submit" class="btn" name="play_button" data-toggle="button" id='pause' onclick='window.traj.player.pause();'>Pause</button>
  <input type="range" min="0" max="100" value="0" class="slider" id="myRange">
  <p>Time: <span id="frame_counter"></span> ns</p>

</div>
<h4>Number of detached DNA base pairs from the nucleosome core (excluding DNA linkers)</h4>
<div id='tooltip' style='position:absolute;background-color:lightgray;padding:5px'></div>
<div id="my_dataviz"></div>