A computational study to assess the pathogenicity of single or combinations of missense variants on respiratory Complex I

Variants found in the respiratory complex I (CI) subunit genes encoded by mitochondrial DNA can cause severe genetic diseases. However, it is difficult to establish a priori whether a single or a combination of CI variants may impact oxidative phosphorylation. Here we propose a computational approach based on coarse-grained molecular dynamics simulations. One of the primary CI variants (m.14484T>C/MT-ND6) associated with the Leber hereditary optic neuropathy was used as a test case. This variant was investigated alone or in combination with two additional rare CI variants whose role remains uncertain. We found that the primary variant stiffens CI dynamics in the crucial E-channel region and that one of the other variants, located in the vicinity of the primary one, further worsens the stiffening. This approach may be extended to other variants candidate to exert a pathogenic impact on CI function, or to investigate the interaction of multiple variants. Teaser Molecular dynamics is able to predict the functional impact of variants hitting respiratory complex I mitochondrial genes.

[1]  J. Hirst,et al.  Investigation of hydrated channels and proton pathways in a high-resolution cryo-EM structure of mammalian complex I , 2023, Science advances.

[2]  Ville R. I. Kaila,et al.  Quinone Catalysis Modulates Proton Transfer Reactions in the Membrane Domain of Respiratory Complex I , 2023, Journal of the American Chemical Society.

[3]  C. Sue,et al.  Low disease risk and penetrance in Leber hereditary optic neuropathy , 2022, American journal of human genetics.

[4]  D. Mackey,et al.  Is the disease risk and penetrance in Leber hereditary optic neuropathy actually low? , 2022, American journal of human genetics.

[5]  G. Pierdominici-Sottile,et al.  On the Uses of PCA to Characterise Molecular Dynamics Simulations of Biological Macromolecules: Basics and Tips for an Effective Use. , 2022, Chemphyschem : a European journal of chemical physics and physical chemistry.

[6]  J. Hirst,et al.  Making the leap from structure to mechanism: are the open states of mammalian complex I identified by cryoEM resting states or catalytic intermediates? , 2022, Current opinion in structural biology.

[7]  G. Hummer,et al.  Small ionic radii limit time step in Martini 3 molecular dynamics simulations. , 2022, The Journal of chemical physics.

[8]  A. Padavannil,et al.  The Mysterious Multitude: Structural Perspective on the Accessory Subunits of Respiratory Complex I , 2022, Frontiers in Molecular Biosciences.

[9]  Ville R. I. Kaila,et al.  Deactivation blocks proton pathways in the mitochondrial complex I , 2021, Proceedings of the National Academy of Sciences.

[10]  Vivek Sharma,et al.  High-resolution structure and dynamics of mitochondrial complex I—Insights into the proton pumping mechanism , 2021, bioRxiv.

[11]  Anna L. Duncan,et al.  Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins , 2021, Science Advances.

[12]  Filipe M. Sousa,et al.  Coarse-Grained Parameterization of Nucleotide Cofactors and Metabolites: Protonation Constants, Partition Coefficients, and Model Topologies , 2021, J. Chem. Inf. Model..

[13]  L. Sazanov,et al.  The coupling mechanism of mammalian respiratory complex I , 2020, Science.

[14]  D. Gallagher,et al.  Key role of quinone in the mechanism of respiratory complex I , 2020, Nature Communications.

[15]  F. Simonelli,et al.  A Martini coarse-grained model of the calcein fluorescent dye , 2018, Journal of Physics D: Applied Physics.

[16]  Jumin Lee,et al.  CHARMM‐GUI Martini Maker for modeling and simulation of complex bacterial membranes with lipopolysaccharides , 2017, J. Comput. Chem..

[17]  Ben M. Webb,et al.  Comparative Protein Structure Modeling Using MODELLER , 2016, Current protocols in bioinformatics.

[18]  Helgi I. Ingólfsson,et al.  CHARMM-GUI Martini Maker for Coarse-Grained Simulations with the Martini Force Field. , 2015, Journal of chemical theory and computation.

[19]  Helgi I Ingólfsson,et al.  The power of coarse graining in biomolecular simulations , 2013, Wiley interdisciplinary reviews. Computational molecular science.

[20]  W F Drew Bennett,et al.  Improved Parameters for the Martini Coarse-Grained Protein Force Field. , 2013, Journal of chemical theory and computation.

[21]  M. Zeviani,et al.  Rare Primary Mitochondrial DNA Mutations and Probable Synergistic Variants in Leber’s Hereditary Optic Neuropathy , 2012, PloS one.

[22]  R. E. Kocher.  Principle , 2011 .

[23]  Hyeon Joo,et al.  OPM database and PPM web server: resources for positioning of proteins in membranes , 2011, Nucleic Acids Res..

[24]  Elizabeth J. Denning,et al.  MDAnalysis: A toolkit for the analysis of molecular dynamics simulations , 2011, J. Comput. Chem..

[25]  Rouslan G. Efremov,et al.  The architecture of respiratory complex I , 2010, Nature.

[26]  Zongfu Cao,et al.  The novel G10680A mutation is associated with complete penetrance of the LHON/T14484C family. , 2009, Mitochondrion.

[27]  Taehoon Kim,et al.  CHARMM‐GUI: A web‐based graphical user interface for CHARMM , 2008, J. Comput. Chem..

[28]  R. Larson,et al.  The MARTINI Coarse-Grained Force Field: Extension to Proteins. , 2008, Journal of chemical theory and computation.

[29]  G. Meer,et al.  Membrane lipids: where they are and how they behave , 2008, Nature Reviews Molecular Cell Biology.

[30]  Patrick F Chinnery,et al.  Clinical expression of Leber hereditary optic neuropathy is affected by the mitochondrial DNA-haplogroup background. , 2007, American journal of human genetics.

[31]  M. Parrinello,et al.  Canonical sampling through velocity rescaling. , 2007, The Journal of chemical physics.

[32]  A. Sali,et al.  Statistical potential for assessment and prediction of protein structures , 2006, Protein science : a publication of the Protein Society.

[33]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[34]  P. Nguyen,et al.  Energy landscape of a small peptide revealed by dihedral angle principal component analysis , 2004, Proteins.

[35]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

[36]  V. Carelli,et al.  Mitochondrial dysfunction as a cause of optic neuropathies , 2004, Progress in Retinal and Eye Research.

[37]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[38]  K. Kuczera,et al.  Transitions from alpha to pi helix observed in molecular dynamics simulations of synthetic peptides. , 2000, Biochemistry.

[39]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[40]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[41]  G Vriend,et al.  The essential dynamics of thermolysin: Confirmation of the hinge‐bending motion and comparison of simulations in vacuum and water , 1995, Proteins.

[42]  J. Thornton,et al.  PROCHECK: a program to check the stereochemical quality of protein structures , 1993 .

[43]  D. Mackey,et al.  A variant of Leber hereditary optic neuropathy characterized by recovery of vision and by an unusual mitochondrial genetic etiology. , 1992, American journal of human genetics.

[44]  N. Howell,et al.  Leber hereditary optic neuropathy: involvement of the mitochondrial ND1 gene and evidence for an intragenic suppressor mutation. , 1991, American journal of human genetics.

[45]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[46]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .