Source author record

Daniel J. Cole

Daniel J. Cole appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.chem-ph Biological Physics cond-mat.str-el Biomolecules math.DG quant-ph Quantitative Methods cond-mat.other physics.comp-ph

Catalog footprint

What is connected

6works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2023arXiv

Riemannian Geometry and Molecular Similarity II: Kähler Quantization

Shape-similarity between molecules is a tool used by chemists for virtual screening, with the goal of reducing the cost and duration of drug discovery campaigns. This paper reports an entirely novel shape descriptor as an alternative to the previously described RGMolSA descriptors \cite{cole2022riemannian}, derived from the theory of Riemannian geometry and Kähler quantization (KQMolSA). The treatment of a molecule as a series of intersecting spheres allows us to obtain the explicit \textit{Riemannian metric} which captures the geometry of the surface, which can in turn be used to calculate a Hermitian matrix $\mathbb{M}$ as a directly comparable surface representation. The potential utility of this method is demonstrated using a series of PDE5 inhibitors considered to have similar shape. The method shows promise in its capability to handle different conformers, and compares well to existing shape similarity methods. The code and data used to produce the results are available at: \url{https://github.com/RPirie96/KQMolSA}.

preprint2022arXiv

Riemannian Geometry and Molecular Surfaces I: Spectrum of the Laplacian

Ligand-based virtual screening aims to reduce the cost and duration of drug discovery campaigns. Shape similarity can be used to screen large databases, with the goal of predicting potential new hits by comparing to molecules with known favourable properties. This paper presents the theory underpinning RGMolSA, a new alignment-free and mesh-free surface-based molecular shape descriptor derived from the mathematical theory of Riemannian geometry. The treatment of a molecule as a series of intersecting spheres allows the description of its surface geometry using the Riemannian metric, obtained by considering the spectrum of the Laplacian. This gives a simple vector descriptor constructed of the weighted surface area and eight non-zero eigenvalues, which capture the surface shape. We demonstrate the potential of our method by considering a series of PDE5 inhibitors that are known to have similar shape as an initial test case. RGMolSA displays promise when compared to existing shape descriptors and in its capability to handle different molecular conformers. The code and data used to produce the results are available via GitHub: https://github.com/RPirie96/RGMolSA.

preprint2020arXiv

ONETEP + TOSCAM: uniting dynamical mean field theory and linear-scaling density functional theory

We introduce the unification of dynamical mean field theory (DMFT) and linear-scaling density functional theory (DFT), as recently implemented in ONETEP, a linear-scaling DFT package, and TOSCAM, a DMFT toolbox. This code can account for strongly correlated electronic behavior while simultaneously including the effects of the environment, making it ideally suited for studying complex and heterogeneous systems containing transition metals and lanthanides, such as metalloproteins. We systematically introduce the necessary formalism, which must account for the non-orthogonal basis set used by ONETEP. In order to demonstrate the capabilities of this code, we apply it to carbon monoxide-ligated iron porphyrin and explore the distinctly quantum-mechanical character of the iron $3d$ electrons during the process of photodissociation.

preprint2014arXiv

Constrained geometric dynamics of the Fenna-Matthews-Olson complex: The role of correlated motion in reducing uncertainty in excitation energy transfer

The Fenna Mathews Olson (FMO) complex of green sulphur bacteria is an example of a photosynthetic pigment protein complex, in which the electronic properties of the pigments are modified by the protein environment to promote efficient excitonic energy transfer from antenna complexes to the reaction centres. Many of the electronic properties of the FMO complex can be extracted from knowledge of the static crystal structure. However, the recent observation and analysis of long lasting quantum dynamics in the FMO complex point to protein dynamics as a key factor in protecting and generating quantum coherence under laboratory conditions. While fast inter and intra molecular vibrations have been investigated extensively, the slow dynamics which effectively determine the optical inhomogeneous broadening of experimental ensembles has received less attention. Our study employs constrained geometric dynamics to study the flexibility in the protein network by efficiently generating the accessible conformational states from the published crystal structure. Statistical and principle component analysis reveal highly correlated low frequency motions between functionally relevant elements, including strong correlations between pigments that are excitonically coupled. Our analysis reveals a hierarchy of structural interactions which enforce these correlated motions, from the level of monomer monomer interfaces right down to the alpha helices, beta sheets and pigments. In addition to inducing strong spatial correlations across the conformational ensemble, we find that the overall rigidity of the FMO complex is exceptionally high. We suggest that these observations support the idea of highly correlated inhomogeneous disorder of the electronic excited states, which is further supported by the remarkably low variance of the excitonic couplings of the conformational ensemble.

preprint2014arXiv

Renormalization of myoglobin-ligand binding energetics by quantum many-body effects

We carry out a first-principles atomistic study of the electronic mechanisms of ligand binding and discrimination in the myoglobin protein. Electronic correlation effects are taken into account using one of the most advanced methods currently available, namely a linear-scaling density functional theory (DFT) approach wherein the treatment of localized iron 3d electrons is further refined using dynamical mean-field theory (DMFT). This combination of methods explicitly accounts for dynamical and multi-reference quantum physics, such as valence and spin fluctuations, of the 3d electrons, whilst treating a significant proportion of the protein (more than 1000 atoms) with density functional theory. The computed electronic structure of the myoglobin complexes and the nature of the Fe-O2 bonding are validated against experimental spectroscopic observables. We elucidate and solve a long standing problem related to the quantum-mechanical description of the respiration process, namely that DFT calculations predict a strong imbalance between O2 and CO binding, favoring the latter to an unphysically large extent. We show that the explicit inclusion of many body-effects induced by the Hund's coupling mechanism results in the correct prediction of similar binding energies for oxy- and carbonmonoxymyoglobin.

preprint2013arXiv

Ligand Discrimination in Myoglobin from Linear-Scaling DFT+U

Myoglobin modulates the binding of diatomic molecules to its heme group via hydrogen-bonding and steric interactions with neighboring residues, and is an important benchmark for computational studies of biomolecules. We have performed calculations on the heme binding site and a significant proportion of the protein environment (more than 1000 atoms) using linear-scaling density functional theory and the DFT+U method to correct for self-interaction errors associated with localized 3d states. We confirm both the hydrogen-bonding nature of the discrimination effect (3.6 kcal/mol) and assumptions that the relative strain energy stored in the protein is low (less than 1 kcal/mol). Our calculations significantly widen the scope for tackling problems in drug design and enzymology, especially in cases where electron localization, allostery or long-ranged polarization influence ligand binding and reaction.