CO2 Response Screen in Grass Brachypodium Reveals Shows Key Role of a MAP-kinase in CO2-induced Stomatal Closure.Lopez, B.N.K., P.H.O. Ceciliato, F. Rangel, E. Salem, K. Kernig, Y. Takahashi, K, Chow, L. Zhang, M. Sidhom, C. Seitz, R. Sibout, D.L. Laudencia-Chingcuanco, D. Woods, J.A. McCammon, J. Vogel, J.I. SchroederPlant Physiology, to appear (2024).
Plants respond to increased CO2 concentrations through rapid stomatal closure which can contribute to increased water use efficiency. Grasses display faster stomatal responses than eudicots due to dumbbell-shaped guard cells flanked by subsidiary cells working in opposition. However, forward genetic screening for stomatal CO2 signal transduction mutants in grasses has not been reported. The grass model Brachypodium distachyon is closely related to agronomically important cereal crops, sharing largely collinear genomes. To gain insights into CO2 control mechanisms of stomatal movements in grasses, we developed a forward genetics screen with an EMS-mutagenized Brachypodium distachyon M5 generation population using infrared imaging to identify plants with altered canopy leaf temperature at elevated CO2. Among isolated mutants, a “chill1” mutant exhibited cooler leaf temperatures than wildtype Bd21-3 parent control plants after exposure to increased [CO2]. chill1 plants showed strongly impaired high CO2-induced stomatal closure, despite retaining a robust abscisic acid-induced stomatal closing response. Through bulked segregant whole-genome-sequencing analyses followed by analyses of further backcrossed F4 generation plants and generation and characterization of CRISPR-cas9 mutants, chill1 was mapped to a protein kinase, BdMPK5. The chill1 mutation impaired BdMPK5 protein-mediated CO2/HCO3- sensing in vitro. Furthermore, AlphaFold2-directed structural modeling suggests that the identified BdMPK5-D90N chill1 mutant residue is located at the interface with the HT1 Raf-like kinase. BdMPK5 is a key signaling component involved in CO2-induced stomatal movements, potentially functioning as a component of the CO2 sensor in grasses.
Dilated cardiomyopathy mutation in beta-cardiac myosin enhances actin activation of the power stroke and phosphate release.Bodt, S.M.L., J. Ge, W. Ma, D.V. Rasicci, R. Desetty, J.A. McCammon, C.M. Yengo.PNAS Nexus to appear (2024).
Multiscale modeling shows that 2-deoxy-ATP improves ventricular function by modulating calcium cycling, crossbridge kinetics, and myosin recruitment.Teitgen, A.E., M.T. Hock, K.J. McCabe, M.C. Childers, G.A. Huber, B. Marzbane, J.A. McCammon, D.A. Beard, M. Regnier, A.D. McCulloch.Proc. Natl. Acad. Sci. USA, to appear (2024).
2'-deoxy-ATP (dATP) improves cardiac function by increasing the rate of crossbridge cycling and Ca[Formula: see text] transient decay. However, the mechanisms of these effects and how therapeutic responses to dATP are achieved when dATP is only a small fraction of the total ATP pool remain poorly understood. Here, we used a multiscale computational modeling approach to analyze the mechanisms by which dATP improves ventricular function. We integrated atomistic simulations of prepowerstroke myosin and actomyosin association, filament-scale Markov state modeling of sarcomere mechanics, cell-scale analysis of myocyte Ca[Formula: see text] dynamics and contraction, organ-scale modeling of biventricular mechanoenergetics, and systems level modeling of circulatory dynamics. Molecular and Brownian dynamics simulations showed that dATP increases the actomyosin association rate by 1.9 fold. Markov state models predicted that dATP increases the pool of myosin heads available for crossbridge cycling, increasing steady-state force development at low dATP fractions by 1.3 fold due to mechanosensing and nearest-neighbor cooperativity. This was found to be the dominant mechanism by which small amounts of dATP can improve contractile function at myofilament to organ scales. Together with faster myocyte Ca[Formula: see text] handling, this led to improved ventricular contractility, especially in a failing heart model in which dATP increased ejection fraction by 16% and the energy efficiency of cardiac contraction by 1%. This work represents a complete multiscale model analysis of a small molecule myosin modulator from single molecule to organ system biophysics and elucidates how the molecular mechanisms of dATP may improve cardiovascular function in heart failure with reduced ejection fraction.
The Mechanism of Allosteric Regulation of Calcium-Independent Phospholipase A2 by ATP and Calmodulin Binding to the Ankyrin Domain.Mouchlis. V.D., Y.H. Hsu, D. Hayashi, J. Cao, S. Li, J.A. McCammon, E.A. DennisProc. Natl. Acad. Sci. USA 121, to appear (2024).
Group VIA calcium-independent phospholipase A2 (iPLA2) is a member of the PLA2 superfamily that exhibits calcium-independent activity in contrast to the other two major types, secreted phospholipase A2 (sPLA2) and cytosolic phospholipase A2 (cPLA2), which both require calcium for their enzymatic activity. Adenosine triphosphate (ATP) has been reported to allosterically activate iPLA2, and this has now been verified with a lipidomics-based mixed-micelle assay, but its mechanism of action has been unknown. Hydrogen/deuterium exchange mass spectrometry (HDX-MS) was employed to identify ATP interaction peptide regions located within the ankyrin repeat domain at which ATP interacts. Molecular dynamics simulations revealed the mechanism by which ATP binds to its site and the main residues that interact. Site-directed mutagenesis was used to verify the importance of these residues in the role of ATP in regulating iPLA2 activity. Importantly, calcium was found to abolish the enhancing regulatory function of ATP and to promote the inhibitory activity by calmodulin. Given previous evidence that calcium does not bind directly to iPLA2, its effect appears to be indirect via association with ATP and/or calmodulin. Using HDX-MS, we found that calmodulin interacts with the N terminus peptide region of iPLA2 consisting of residues 20 to 28. These two regulatory iPLA2 sites open the road to the development of potential targets for therapeutic intervention.
NetSci: A Library for High Performance Biomolecular Simulation Network Analysis Computation.Stokely, A., L. Votapka, M. Hock, A. Teitgen, J.A. McCammon, A. McCulloch, R.E. Amaro.J. Chem. Info. Model. 64, 7966-7976 (2024).
We present the NetSci program–an open-source scientific software package designed for estimating mutual information (MI) between data sets using GPU acceleration and a k-nearest-neighbor algorithm. This approach significantly enhances calculation speed, achieving improvements of several orders of magnitude over traditional CPU-based methods, with data set size limits dictated only by available hardware. To validate NetSci, we accurately compute MI for an analytically verifiable two-dimensional Gaussian distribution and replicate the generalized correlation (GC) analysis previously conducted on the B1 domain of protein G. We also apply NetSci to molecular dynamics simulations of the Sarcoendoplasmic Reticulum Calcium-ATPase (SERCA) pump, exploring the allosteric mechanisms and pathways influenced by ATP and 2′-deoxy-ATP (dATP) binding. Our analysis reveals distinct allosteric effects induced by ATP compared to dATP, with predicted information pathways from the bound nucleotide to the calcium-binding domain differing based on the nucleotide involved. NetSci proves to be a valuable tool for estimating MI and GC in various data sets and is particularly effective for analyzing intraprotein communication and information transfer.
Structure and Dynamics in Drug Discovery: A Perspective.Wei, H., J.A. McCammon.Nature Portfolio Journals NPJ Drug Discovery 1, 1 (2024).
The roles of computing in structure-based drug discovery are considered, from early studies based on some of the first experimental structures of enzyme-inhibitor complexes, through the use of advanced molecular dynamics simulations and machine learning methods. This perspective aims to explore the history, current trends, and future directions of these methodologies.