Mathematical modelling of DNA
Course will be taught online. Lectures (video, audio, notes) and exercises will be uploaded week by week, usually on Mondays. An interactive zoom session will be held on each Friday (15.15 to 17.00) to answer questions based on previous week (or weeks) lecture video and exercise session(s).
Raushan Singh, Rahul Sharma
1st and 2nd year courses in math or physics, (or with teacher's permission).
Helpful although not required
This course is designed to be an introduction, within the particular context of DNA, to the interplay between analysis, computation and experiment that makes up the process called mathematical modelling. In addition to students whose primary interest is in DNA, the syllabus is intended for students wishing an introduction to the modelling process in general, and the course will describe a number of widely encountered mathematical and computational techniques.
The course will be a detailed introduction to the cgDNA and cgDNA+ sequence-dependent coarse grain models of DNA, including both how to use them to predict various biologically pertinent sequence-dependent expectations with an associated Monte Carlo code, and all the extensive underlying applied mathematics necessary to estimate their parameter sets from a library of Molecular Dynamics simulations. The cgDNA models are a research tool that have their own web page . The course will work through the details of publications described on that page, specifically, , ,  and  below.
The course has five chapters.
0) Introduction to DNA and a brief overview of its coarse grain models.
1) The sequence-dependent cgDNA models.
2) Monte Carlo methods for sampling cgDNA model equilibrium distributions and application to DNA persistence lengths.
3) Parameter estimation for the cgDNA models from Molecular Dynamics time series.
4) Equality constrained nonlinear optimisation with application to computing cgDNA equilibria.
Note: Login id to access lecture notes is "moddna" and password will be communicated through email.
Description of the basic structure of DNA, and multiscaling (or coarse graining) approaches. The need for a tertiary structure model of DNA, i.e. a sequence-dependent coarse grain model. Overview of the cgDNA coarse grain model to predict a Gaussian PDF for the configuration distribution of a DNA fragment of given sequence.
Here are the lecture notes and video of this week. Here you have supplementary material for this first lecture., eg some pre-recordded videos describing some Molecular Dynamics (or MD) very fine grain simulations of DNA.
Definiton of Inputs (configuration coordinate w, sequence S and parameter set P) and Output (a Gaussian pdf, with its ground state and stiffness matrix, on the configuration space coordinates w) of a cgDNA model. Motivation of why it is reasonable to expect an equilibrium configuration space pdf to exist, and why it is reasonable to estimate coarse grain parameter sets P with a few intensive fine grain atomisitc Molecular Dynamics simulations of dsDNA plus solvent bath.
Here are the lecture notes and video and audio of this week. There is no friday video this week because there was no specific questions just a general back ground discussion. Note that in the lecture of this week the symbol "N" is mistakenly used to denote two different things namely the number of base-pairs in the sequence S, and also the dimension of the configuration space coordinate vector w (and the groundstate \mu and size of the stifness matrix K). Should have been two different symbols. The length of the sequence is certainly related to the dimension of the configuration space coordinate vector (see future lectures) but they are never the same.
Two lectures preparing for the specification of our dsDNA coordinate system using chains of rigid bodies. Lecture 1, the matrix groups SO(3) and SE(3). Lecture 2, interpretation of SO(3) and SE(3) as rigid body configurations and relative rigid body displacements between pairs of rigid body configurations.
Here are the lecture notes and video and audio of Lecture 1.
Here are the lecture notes and video and audio of Lecture 2.
|Week 4||Symmetric coordinates for relative SE(3) displacements between a pair of rigid bodies, and the importance of introducing a mid-point frame. Definition of mid frame involves (principal) square root of rotation matrix, but Euclidean average of origins. (more detail of choices in Exercise Session 4). With components of relative translation expressed in mid-frame, and components of Cayley vector of relative expression expressed in any of the three frames R, R^+, R^-, the transformation on the coordinates corresponding to reversing roles of + and - bodies is u <-> -u, v<-> -v. A related transformation for the Crick-Watson strand symmetry will be used in our coordinate system for double stranded DNA, but it is different because for the moment note that there is no account taken of the additional feature that the Crick and Watson embedding rules for frame into base are different. See Qu 2 in Series 4. Overview of material on Cayley transforms and Cayley vectors already introduced in Exercise sets. Geometrical view of the relations between axis-angle, Euler-parameter/quaternion and Cayley vector parametrisations of SO(3). Here are the lecture notes and video and audio of this week. In the following lectures and exercises we will make use of various matrix factorisations. A brief summary of the results we will use are provided in this PDF . Most or indeed all of the factorisations should be familiar to you.|
|Week 5||Generalisation of symmetrized relative SE(3) coordinates between a pair of rigid bodies to the cgDNA+ model internal coordinates and the associated tree structure for a quadruple chain of rigid bodies, with intra and inter and 5'-phosphate coordinates. Watson or reading strand, and the re-embedding of frames on the Crick strand to avoid rotations through angles close to \pi. Definition of base-pair and junction frames as mid-frames. cgDNA+ model intra and inter configuration coordinates: translations expressed in mid-frames (base-pair frame between two base frames for intras, junction frames between two base-pair frames for inters) and Cayley vectors of relative rotations for both intra and inter relative rotations (with matrix multiplication on the right). 5'-Phosphate Cayley vectors and translations in associated base frame (before re-embedding for the Crick bases). Here are the lecture notes and video and audio of this week.|
|Week 6||Further (and final) details on cgDNA+ coordinate vectors w. Non dimensionalisation and scaling of cgDNA+ coordinates. Transformation under Crick-Watson change of reading strand for cgDNA+ coordinates (more detailed treatment in exercise session). Short introduction to using the cgDNAweb server. Here are the lecture notes and video and audio of this week. Here are the lecture notes and video and audio of short introduction to cgDNAweb+.|
|Week 7||Definitions and assumptions underlying the cgDNA and cgDNA+ rigid base and double nucleotide chain coarse grain model free energies and their associated Gaussian PDFs: a) (five) nearest-neighbour base/ nucleotide interactions, plus b) dimer sequence-dependence of parameter set blocks. Leads to a Gaussian model where the stiffness matrix has a banded structure with overlapping blocks. Description of the Crick-Watson symmetry properties of predicted sequence dependent ground-states and stiffness matrices, and the (sufficient) Crick-Watson symmetry properties of the elements of the parameter set. The latter properties reduce the total number of independent blocks in the parameter set to 16 independent 5'XY cases and 10 independent XY dimer steps, 6/12 in non-palindromic steps and 4 palindromic sets that have to satisfy self-symmetry conditions. Here are the lecture notes (of two lecture) of this week. Here are the video and audio of Lecture 1. Here are the video and audio of Lecture 2.|
|Week 8||Start of Chapter 2: What can be done with the cgDNA model? Brief discussion of i) probabilities and looping experiments, and longer discussion of ii) expectations, specifically correlations along a polymer chain. Numerical approximations of both from an ensemble of configurations generated by an appropriate (direct sampling for our multivariate Gaussian pdf as opposed to Metropolis or Markov chain sampling that are necessary for more complicated pdfs) Monte Carlo code e.g. cgDNApmc, counting hits and misses for i), and averaging over an ensemble as a simple quadrature rule for ii). First mention of expectations leading to persistence lengths. Importance for efficiency of MC of bandedness of the stiffness matrix and use of the associated banded Cholesky factorisation to diagonalise. Here are the lecture notes and video and audio of this week.|
|Week 9||End of Chapter 2: More detailed discussion on expectations leading to tangent-tangent and Flory persistence lengths. Correlations of relative frame rotations and translations along a chain using homogeneous coordinates in SE(3) and the associated matrix multiplication. Simplifications when junction statistics are independent (the I.D. case), and when the chain is uniform (the I.I.D. case). Exponential decay of frame rotation correlations as the index difference grows, and convergence of the translation block to the Flory persistence vector. Exercises comparison with cgDNAmc data for poly(XY) ie sequences with close to intrinsically straight ground states, and cgDNAmc data for lambda-phage sequences some with significantly bent ground states, where shape factorised semi-log tan-tan plots remain close to linear is made in the exercise session. Here are the lecture notes and video and audio of this week.|
|Week 10||Start of Chapter 3: Parameter Estimation in the cgDNA/+ model. Start estimation of oligomer based mean and centred covariance from MD time series data for the cgDNA coarse grain variables. This week describes the Maximum Likelihood approach to obtain estimates of oligomer-based Gaussian pdfs from an ensemble of configuration snapshots. Next week we consider an alternative Maximum Entropy approach. Cases both with and without imposed banded sparsity pattern in the stiffness matrix. Note: this weeks video lecture and notes are precisely the same as last years course, in which there was no mention of the cgDNA+ model, as opposed to the cgDNA model which is of course mentioned. Re-using the old video works as the material of this week and next is quite generally applicable to estimating parameters in any Gaussian pdf with an overlapping block structure in the stiffness matrix. The only difference between cgDNA and cgDNA+ for the purpose of this weeks lecture is that cgDNA is 18x18 overlapping blocks with 6x6 overlaps, while cgDNA+ is 42x42 overlapping blocks with 18x18 overlaps (and two end blocks different size), and these differences in block size make no difference to the general mathematical approach that is discussed. The 2021 exercises are adapted from the 2020 ones to make explicit mention of both cgDNA+ and cgDNA cases. Here are the lecture notes and video and audio of this week.|
|Week 11||Introduction of entropy for a continuous pdf with respect to an associated measure, and relative entropy (or Kullback-Leibler divergence) between two continuous pdfs with associated measure. Contrast between between discrete probability distributions and continuous probability density functions or pdfs, and differences between the cases of bounded and unbounded domains. Jensen inequality to prove that entropy minimizing (or maximizing depending on sign convention) pdf on a bounded domain is uniform with respect to the measure. Jensen inequality to prove that relative entropy is always non-negative. Jaynes max entropy principle to characterize Boltzmann pdfs as pdfs maximising (or minimising in our sign convention) entropy subject to constrained of values of some moments, and sufficiency of the associated first-order necessary conditions. In the case of constrained first and second moments on R^n, characterisation of Gaussian pdfs as maximum entropy distributions. Here are the lecture notes and video and audio of this week.|
|Week 12||For the single (short) lecture of this week. Topic is general further properties of Kullback-Leibler divergence. Specifically a) KL between two factorisable pdfs is sum of KL between each of the factors, which motivates the scaling of KL per degree of freedom, b) closed form expression for KL between two Gaussians (which makes the use of KL in high dimensional model parameter estimation feasible), and c) the two ways KL can be used in parameter estimation. Here are the lecture notes and video and audio for lecture of this week.|
|Week 13||The final two lectures in the course (both long, partially compensated by the single short lecture of last week). Lecture 1, designing a training library for the cgDNA+ model, some aspects of Molecular Dynamics simulations, and using palindromy to assess convergence of time series of snapshots. Lecture 2: assessing validity of cgDNA+ model assumptions on palindromic training set data. Using Fisher information (or second variation of KL divergence) in numerics of cgDNA+ parameter set estimation. Constructing positive definite paraneter set blocks using non-uniqueness of optimal parameter sets. Here are the lecture notes and video and audio for lecture 1 of this week. Here are the lecture notes and video and audio for lecture 2 of this week.|
Summary and description of the exercices
This document contains an overview and a description of all the exercises given so far.
-  A DNA Coarse-Grain Rigid Base Model and Parameter Estimation from Molecular Dynamics Simulations , D. Petkevičiūtė Thesis #5520, EPFL, (2012).
-  cgDNA: a software package for the prediction of sequence-dependent coarse-grain free energies of B-form DNA , D. Petkevičiūtė, M. Pasi, O. Gonzalez and J. H. Maddocks Nucleic Acids Research 42, no. 20 (2014), p. e153, (2014) .
-  A sequence-dependent rigid-base model of DNA , O. Gonzalez, D. Petkevičiūtė, and J. H. Maddocks, Journal of Chemical Physics 138, no. 5 (2013), p. 055122 1-28 .
-  Sequence-dependent persistence lengths of DNA , J. S. Mitchell, J. Glowacki, A. E. Grandchamp, R. S. Manning and J. H. Maddocks, Journal of Chemical Theory and Computation, no. 13 (2017), p. 1539-1555 .
-  Absolute versus relative entropy parameter estimation in a coarse-grain model of DNA , O. Gonzalez, M. Pasi, D. Petkevičiūtė, J. Glowacki, J.H. Maddocks, Multiscale Modeling and Simulation 15, no. 3 (2017), p. 1073 - 1107 .
-  A sequence-dependent coarse-grain model of B-DNA with explicit description of bases and phosphate groups parametrised from large scale Molecular Dynamics simulations , A. Patelli Thesis #9552, EPFL, (2019).
-  On the stochastic modeling of rigid body systems with application to polymer dynamics , J. Walter, O. Gonzalez and J.H. Maddocks, SIAM Multiscale Modeling and Simulation 8, no. 3 (2010), p. 1018 – 1053.
References for general books on DNA.
 Understanding DNA, The molecule & how it work C. R. Calladine, H. R. Drew, B. F. Luisi, A. A. Travers, Third Edition, 2004, Academic Press, ISBN 9780121550893 .
Summary: Understanding DNA explains, step by step, why DNA forms specific structures, the form of these structures and how they fundamentally affect the biological processes of transcription and replication.
 Unraveling Dna: The Most Important Molecule Of Life M. D. Frank-Kamenetskii, Revised and Updated Edition, 1997, Perseus Publishing, ISBN 9780201155846.
Summary: A curious blend of history, biographical details to cover the development of molecular biology from the influence of physicists earlier in the century, through the central dogma of molecular biology to discussion of social issues raised by genetic engineering.
 DNA topology A. D. Bates & A. Maxwell, 2005, Oxford University Press, ISBN 9780198506553.
Summary: A clear, concise explanation of the relevance of supercoiling and catenation in the context of biological activity of the DNA molecule.
 DNA structure and Function R. R. Sinden, 1994, Academic Press, ISBN 9780126457506.
Summary: a timely and comprehensive resource, that provides a simple yet comprehensive introduction to nearly all aspects of DNA structure. It also explains current ideas on the biological significance of classic and alternative DNA conformations.