Dilated cardiomyopathy mutation in beta-cardiac myosin enhances actin activation of the power stroke and phosphate release.

Inherited mutations in human beta-cardiac myosin (M2ß) can lead to severe forms of heart failure. The E525K mutation in M2ß is associated with dilated cardiomyopathy (DCM) and was found to stabilize the interacting heads motif (IHM) and autoinhibited super-relaxed (SRX) state in dimeric heavy meromyosin. However, in monomeric M2ß subfragment 1 (S1) we found that E525K enhances (3-fold) the maximum steady-state actin-activated ATPase activity (kcat) and decreases (6-fold) the actin concentration at which ATPase is one-half maximal (KATPase). We also found a 3 to 4-fold increase in the actin-activated power stroke and phosphate release rate constants at 30 µM actin, which overall enhanced the duty ratio 3-fold. Loaded motility assays revealed that the enhanced intrinsic motor activity translates to increased ensemble force in M2ß S1. Glutamate 525, located near the actin binding region in the so-called activation loop, is highly conserved and predicted to form a salt-bridge with another conserved residue (lysine 484) in the relay helix. Enhanced sampling molecular dynamics simulations predict that the charge reversal mutation disrupts the E525-K484 salt-bridge, inducing conformations with a more flexible relay helix and a wide phosphate release tunnel. Our results highlight a highly conserved allosteric pathway associated with actin activation of the power stroke and phosphate release and suggest an important feature of the autoinhibited IHM is to prevent this region of myosin from interacting with actin. The ability of the E525K mutation to stabilize the IHM likely overrides the enhanced intrinsic motor properties, which may be key to triggering DCM pathogenesis.

Glycosylation and Crowded Membrane Effects on Influenza Neuraminidase Stability and Dynamics.

All protein simulations are conducted with varying degrees of simplification, oftentimes with unknown ramifications about how these simplifications affect the interpretability of the results. In this work, we investigated how protein glycosylation and lateral crowding effects modulate an array of properties characterizing the stability and dynamics of influenza neuraminidase. We constructed three systems: (1) glycosylated neuraminidase in a whole virion (i.e., crowded membrane) environment, (2) glycosylated neuraminidase in its own lipid bilayer, and (3) unglycosylated neuraminidase in its own lipid bilayer. We saw that glycans tend to stabilize the protein structure and reduce its conformational flexibility while restricting the solvent movement. Conversely, a crowded membrane environment encouraged exploration of the free energy landscape and a large-scale conformational change, while making the protein structure more compact. Understanding these effects informs what factors one must consider in attempting to recapture the desired level of physical accuracy.

Glycosylation and Crowded Membrane Effects on Influenza Neuraminidase Stability and Dynamics.

All protein simulations are conducted with varying degrees of simplifications, oftentimes with unknown ramifications on how these simplifications affect the interpretability of the results. In this work we investigated how protein glycosylation and lateral crowding effects modulate an array of properties characterizing the stability and dynamics of influenza neuraminidase. We constructed three systems: 1) Glycosylated neuraminidase in a whole virion (i.e. crowded membrane) environment 2) Glycosylated neuraminidase in its own lipid bilayer 3) Unglycosylated neuraminidase in its own lipid bilayer. We saw that glycans tend to stabilize the protein structure and reduce its conformational flexibility while restricting solvent movement. Conversely, a crowded membrane environment encouraged exploration of the free energy landscape and a large scale conformational change while making the protein structure more compact. Understanding these effects informs what factors one must consider while attempting to recapture the desired level of physical accuracy.

Interface Engineering of Carrier-Protein-Dependent Metabolic Pathways.

Carrier-protein-dependent metabolic pathways biosynthesize fatty acids, polyketides, and non-ribosomal peptides, producing metabolites with important pharmaceutical, environmental, and industrial properties. Recent findings demonstrate that these pathways rely on selective communication mechanisms involving protein-protein interactions (PPIs) that guide enzyme reactivity and timing. While rational design of these PPIs could enable pathway design and modification, this goal remains a challenge due to the complex nature of protein interfaces. Computational methods offer an encouraging avenue, though many score functions fail to predict experimental observables, leading to low success rates. Here, we improve upon the Rosetta score function, leveraging experimental data through iterative rounds of computational prediction and mutagenesis, to design a hybrid fatty acid-non-ribosomal peptide initiation pathway. By increasing the weight of the electrostatic score term, the computational protocol proved to be more predictive, requiring fewer rounds of iteration to identify mutants with high in vitro activity. This allowed efficient design of new PPIs between a non-ribosomal peptide synthetase adenylation domain, PltF, and a fatty acid synthase acyl carrier protein, AcpP, as validated by activity and structural studies. This method provides a promising platform for customized pathway design, establishing a standard for carrier-protein-dependent pathway engineering through PPI optimization.

Multiscale computational modeling of the effects of 2'-deoxy-ATP on cardiac muscle calcium handling.

2'-Deoxy-ATP (dATP), a naturally occurring near analog of ATP, is a well-documented myosin activator that has been shown to increase contractile force, improve pump function, and enhance lusitropy in the heart. Calcium transients in cardiomyocytes with elevated levels of dATP show faster calcium decay compared with cardiomyocytes with basal levels of dATP, but the mechanisms behind this are unknown. Here, we design and utilize a multiscale computational modeling framework to test the hypothesis that dATP acts on the sarcoendoplasmic reticulum calcium-ATPase (SERCA) pump to accelerate calcium re-uptake into the sarcoplasmic reticulum during cardiac relaxation. Gaussian accelerated molecular dynamics simulations of human cardiac SERCA2A in the E1 apo, ATP-bound and dATP-bound states showed that dATP forms more stable contacts in the nucleotide binding pocket of SERCA and leads to increased closure of cytosolic domains. These structural changes ultimately lead to changes in calcium binding, which we assessed using Brownian dynamics simulations. We found that dATP increases calcium association rate constants to SERCA and that dATP binds to apo SERCA more rapidly than ATP. Using a compartmental ordinary differential equation model of human cardiomyocyte excitation-contraction coupling, we found that these increased association rate constants contributed to the accelerated rates of calcium transient decay observed experimentally. This study provides clear mechanistic evidence of enhancements in cardiac SERCA2A pump function due to interactions with dATP.

PERIOD phosphorylation leads to feedback inhibition of CK1 activity to control circadian period.

PERIOD (PER) and Casein Kinase 1δ regulate circadian rhythms through a phosphoswitch that controls PER stability and repressive activity in the molecular clock. CK1δ phosphorylation of the familial advanced sleep phase (FASP) serine cluster embedded within the Casein Kinase 1 binding domain (CK1BD) of mammalian PER1/2 inhibits its activity on phosphodegrons to stabilize PER and extend circadian period. Here, we show that the phosphorylated FASP region (pFASP) of PER2 directly interacts with and inhibits CK1δ. Co-crystal structures in conjunction with molecular dynamics simulations reveal how pFASP phosphoserines dock into conserved anion binding sites near the active site of CK1δ. Limiting phosphorylation of the FASP serine cluster reduces product inhibition, decreasing PER2 stability and shortening circadian period in human cells. We found that Drosophila PER also regulates CK1δ via feedback inhibition through the phosphorylated PER-Short domain, revealing a conserved mechanism by which PER phosphorylation near the CK1BD regulates CK1 kinase activity.

Capturing Differences in the Regulation of LRRK2 Dynamics and Conformational States by Small Molecule Kinase Inhibitors.

Mutations in the human leucine rich repeat protein kinase-2 (LRRK2) create risk factors for Parkinson's disease, and pathological functions of LRRK2 are often correlated with aberrant kinase activity. Past research has focused on developing selective LRRK2 kinase inhibitors. In this study, we combined enhanced sampling simulations with HDX-MS to characterize the inhibitor-induced dynamic changes and the allosteric communications within the C-terminal domains of LRRK2, LRRK2RCKW. We find that the binding of MLi-2 (a type I kinase inhibitor) stabilizes a closed kinase conformation and reduces the global dynamics of LRRK2RCKW, leading to a more compact LRRK2RCKW structure. In contrast, the binding of Rebastinib (a type II kinase inhibitor) stabilizes an open kinase conformation, which promotes a more extended LRRK2RCKW structure. By probing the distinct effects of the type I and type II inhibitors, key interdomain interactions are found to regulate the communication between the kinase domain and the GTPase domain. The intermediate states revealed in our simulations facilitate the efforts toward in silico design of allosteric modulators that control LRRK2 conformations and potentially mediate the oligomeric states of LRRK2 and its interactions with other proteins.

Integrating comparative modeling and accelerated simulations reveals conformational and energetic basis of actomyosin force generation.

Muscle contraction is performed by arrays of contractile proteins in the sarcomere. Serious heart diseases, such as cardiomyopathy, can often be results of mutations in myosin and actin. Direct characterization of how small changes in the myosin-actin complex impact its force production remains challenging. Molecular dynamics (MD) simulations, although capable of studying protein structure-function relationships, are limited owing to the slow timescale of the myosin cycle as well as a lack of various intermediate structures for the actomyosin complex. Here, employing comparative modeling and enhanced sampling MD simulations, we show how the human cardiac myosin generates force during the mechanochemical cycle. Initial conformational ensembles for different myosin-actin states are learned from multiple structural templates with Rosetta. This enables us to efficiently sample the energy landscape of the system using Gaussian accelerated MD. Key myosin loop residues, whose substitutions are related to cardiomyopathy, are identified to form stable or metastable interactions with the actin surface. We find that the actin-binding cleft closure is allosterically coupled to the myosin motor core transitions and ATP-hydrolysis product release from the active site. Furthermore, a gate between switch I and switch II is suggested to control phosphate release at the prepowerstroke state. Our approach demonstrates the ability to link sequence and structural information to motor functions.

Stomatal CO2/bicarbonate sensor consists of two interacting protein kinases, Raf-like HT1 and non-kinase-activity requiring MPK12/MPK4.

The continuing rise in the atmospheric carbon dioxide (CO2) concentration causes stomatal closing, thus critically affecting transpirational water loss, photosynthesis, and plant growth. However, the primary CO2 sensor remains unknown. Here, we show that elevated CO2 triggers interaction of the MAP kinases MPK4/MPK12 with the HT1 protein kinase, thus inhibiting HT1 kinase activity. At low CO2, HT1 phosphorylates and activates the downstream negatively regulating CBC1 kinase. Physiologically relevant HT1-mediated phosphorylation sites in CBC1 are identified. In a genetic screen, we identify dominant active HT1 mutants that cause insensitivity to elevated CO2. Dominant HT1 mutants abrogate the CO2/bicarbonate-induced MPK4/12-HT1 interaction and HT1 inhibition, which may be explained by a structural AlphaFold2- and Gaussian-accelerated dynamics-generated model. Unexpectedly, MAP kinase activity is not required for CO2 sensor function and CO2-triggered HT1 inhibition and stomatal closing. The presented findings reveal that MPK4/12 and HT1 together constitute the long-sought primary stomatal CO2/bicarbonate sensor upstream of the CBC1 kinase in plants.

Essential Role of Loop Dynamics in Type II NRPS Biomolecular Recognition.

Non-ribosomal peptides play a critical role in the clinic as therapeutic agents. To access more chemically diverse therapeutics, non-ribosomal peptide synthetases (NRPSs) have been targeted for engineering through combinatorial biosynthesis; however, this has been met with limited success in part due to the lack of proper protein-protein interactions between non-cognate proteins. Herein, we report our use of chemical biology to enable X-ray crystallography, molecular dynamics (MD) simulations, and biochemical studies to elucidate binding specificities between peptidyl carrier proteins (PCPs) and adenylation (A) domains. Specifically, we determined X-ray crystal structures of a type II PCP crosslinked to its cognate A domain, PigG and PigI, and of PigG crosslinked to a non-cognate PigI homologue, PltF. The crosslinked PCP-A domain structures possess large protein-protein interfaces that predominantly feature hydrophobic interactions, with specific electrostatic interactions that orient the substrate for active site delivery. MD simulations of the PCP-A domain complexes and unbound PCP structures provide a dynamical evaluation of the transient interactions formed at PCP-A domain interfaces, which confirm the previously hypothesized role of a PCP loop as a crucial recognition element. Finally, we demonstrate that the interfacial interactions at the PCP loop 1 region can be modified to control PCP binding specificity through gain-of-function mutations. This work suggests that loop conformational preferences and dynamism account for improved shape complementary in the PCP-A domain interactions. Ultimately, these studies show how crystallographic, biochemical, and computational methods can be used to rationally re-engineer NRPSs for non-cognate interactions.

Architecture and self-assembly of the jumbo bacteriophage nuclear shell.

Bacteria encode myriad defences that target the genomes of infecting bacteriophage, including restriction-modification and CRISPR-Cas systems1. In response, one family of large bacteriophages uses a nucleus-like compartment to protect its replicating genomes by excluding host defence factors2-4. However, the principal composition and structure of this compartment remain unknown. Here we find that the bacteriophage nuclear shell assembles primarily from one protein, which we name chimallin (ChmA). Combining cryo-electron tomography of nuclear shells in bacteriophage-infected cells and cryo-electron microscopy of a minimal chimallin compartment in vitro, we show that chimallin self-assembles as a flexible sheet into closed micrometre-scale compartments. The architecture and assembly dynamics of the chimallin shell suggest mechanisms for its nucleation and growth, and its role as a scaffold for phage-encoded factors mediating macromolecular transport, cytoskeletal interactions, and viral maturation.

Investigating Intrinsically Disordered Proteins With Brownian Dynamics.

Intrinsically disordered proteins (IDPs) have recently become systems of great interest due to their involvement in modulating many biological processes and their aggregation being implicated in many diseases. Since IDPs do not have a stable, folded structure, however, they cannot be easily studied with experimental techniques. Hence, conducting a computational study of these systems can be helpful and be complementary with experimental work to elucidate their mechanisms. Thus, we have implemented the coarse-grained force field for proteins (COFFDROP) in Browndye 2.0 to study IDPs using Brownian dynamics (BD) simulations, which are often used to study large-scale motions with longer time scales and diffusion-limited molecular associations. Specifically, we have checked our COFFDROP implementation with eight naturally occurring IDPs and have investigated five (Glu-Lys)25 IDP sequence variants. From measuring the hydrodynamic radii of eight naturally occurring IDPs, we found the ideal scaling factor of 0.786 for non-bonded interactions. We have also measured the entanglement indices (average C α distances to the other chain) between two (Glu-Lys)25 IDP sequence variants, a property related to molecular association. We found that entanglement indices decrease for all possible pairs at excess salt concentration, which is consistent with long-range interactions of these IDP sequence variants getting weaker at increasing salt concentration.

Lipoprotein-associated phospholipase A2: A paradigm for allosteric regulation by membranes.

Lipoprotein-associated phospholipase A2 (Lp-PLA2) associates with low- and high-density lipoproteins in human plasma and specifically hydrolyzes circulating oxidized phospholipids involved in oxidative stress. The association of this enzyme with the lipoprotein's phospholipid monolayer to access its substrate is the most crucial first step in its catalytic cycle. The current study demonstrates unequivocally that a significant movement of a major helical peptide region occurs upon membrane binding, resulting in a large conformational change upon Lp-PLA2 binding to a phospholipid surface. This allosteric regulation of an enzyme's activity by a large membrane-like interface inducing a conformational change in the catalytic site defines a unique dimension of allosterism. The mechanism by which this enzyme associates with phospholipid interfaces to select and extract a single phospholipid substrate molecule and carry out catalysis is key to understanding its physiological functioning. A lipidomics platform was employed to determine the precise substrate specificity of human recombinant Lp-PLA2 and mutants. This study uniquely elucidates the association mechanism of this enzyme with membranes and its resulting conformational change as well as the extraction and binding of specific oxidized and short acyl-chain phospholipid substrates. Deuterium exchange mass spectrometry coupled with molecular dynamics simulations was used to define the precise specificity of the subsite for the oxidized fatty acid at the sn-2 position of the phospholipid backbone. Despite the existence of several crystal structures of this enzyme cocrystallized with inhibitors, little was understood about Lp-PLA2's specificity toward oxidized phospholipids.

Gaussian-Accelerated Molecular Dynamics with the Weighted Ensemble Method: A Hybrid Method Improves Thermodynamic and Kinetic Sampling.

Gaussian-accelerated molecular dynamics (GaMD) is a well-established enhanced sampling method for molecular dynamics simulations that effectively samples the potential energy landscape of the system by adding a boost potential, which smoothens the surface and lowers the energy barriers between states. GaMD is unable to give time-dependent properties such as kinetics directly. On the other hand, the weighted ensemble (WE) method can efficiently sample transitions between states with its many weighted trajectories, which directly yield rates and pathways. However, convergence to equilibrium conditions remains a challenge for the WE method. Hence, we have developed a hybrid method that combines the two methods, wherein GaMD is first used to sample the potential energy landscape of the system and WE is subsequently used to further sample the potential energy landscape and kinetic properties of interest. We show that the hybrid method can sample both thermodynamic and kinetic properties more accurately and quickly compared to using either method alone.

Data for molecular dynamics simulations of Escherichia coli cytochrome bd oxidase with the Amber force field.

Cytochrome bd-type quinol oxidase is an important metalloenzyme that allows many bacteria to survive in low oxygen conditions. Since bd oxidase is found in many prokaryotes but not in eukaryotes, it has emerged as a promising bacterial drug target. Examples of organisms containing bd oxidases include the Mycobacterium tuberculosis (Mtb) bacterium that causes tuberculosis (TB) in humans, the Vibrio cholerae bacterium that causes cholera, the Pseudomonas aeruginosa bacterium that contributes to antibiotic resistance and sepsis, and the Campylobacter jejuni bacterium that causes food poisoning. Escherichia coli (E. coli) is another organism exhibiting the cytochrome bd oxidase. Since it has the highest sequence identity to Mtb (36%) and we are ultimately interested in finding drug targets for TB, we have built parameters for the E. coli bd oxidase (Protein Data Bank ID number: 6RKO) that are compatible with the all-atom Amber ff14SB force field for molecular dynamics (MD) simulations. Specifically, we built parameters for the three heme cofactors present in all species of bacterial cytochrome bd-type oxidases (heme b 558 , heme b 595 , and heme d) along with their axial ligands. This data report includes the parameter and library files that can be used with Amber's LEaP program to generate input files for MD simulations using the Amber software package. We also provide the PDB data files of the initial model both by itself and solvated with TIP3P water molecules and counterions.

A glycan gate controls opening of the SARS-CoV-2 spike protein.

SARS-CoV-2 infection is controlled by the opening of the spike protein receptor binding domain (RBD), which transitions from a glycan-shielded 'down' to an exposed 'up' state to bind the human angiotensin-converting enzyme 2 receptor and infect cells. While snapshots of the 'up' and 'down' states have been obtained by cryo-electron microscopy and cryo-electron tomagraphy, details of the RBD-opening transition evade experimental characterization. Here over 130 µs of weighted ensemble simulations of the fully glycosylated spike ectodomain allow us to characterize more than 300 continuous, kinetically unbiased RBD-opening pathways. Together with ManifoldEM analysis of cryo-electron microscopy data and biolayer interferometry experiments, we reveal a gating role for the N-glycan at position N343, which facilitates RBD opening. Residues D405, R408 and D427 also participate. The atomic-level characterization of the glycosylated spike activation mechanism provided herein represents a landmark study for ensemble pathway simulations and offers a foundation for understanding the fundamental mechanisms of SARS-CoV-2 viral entry and infection.

Elucidation of Cryptic and Allosteric Pockets within the SARS-CoV-2 Main Protease.

The SARS-CoV-2 pandemic has rapidly spread across the globe, posing an urgent health concern. Many quests to computationally identify treatments against the virus rely on in silico small molecule docking to experimentally determined structures of viral proteins. One limit to these approaches is that protein dynamics are often unaccounted for, leading to overlooking transient, druggable conformational states. Using Gaussian accelerated molecular dynamics to enhance sampling of conformational space, we identified cryptic pockets within the SARS-CoV-2 main protease, including some within regions far from the active site. These simulations sampled comparable dynamics and pocket volumes to conventional brute force simulations carried out on two orders of magnitude greater timescales.

Decoding allosteric regulation by the acyl carrier protein.

Enzymes in multistep metabolic pathways utilize an array of regulatory mechanisms to maintain a delicate homeostasis [K. Magnuson, S. Jackowski, C. O. Rock, J. E. Cronan, Jr, Microbiol. Rev. 57, 522-542 (1993)]. Carrier proteins in particular play an essential role in shuttling substrates between appropriate enzymes in metabolic pathways. Although hypothesized [E. Ploskon et al., Chem. Biol. 17, 776-785 (2010)], allosteric regulation of substrate delivery has never before been demonstrated for any acyl carrier protein (ACP)-dependent pathway. Studying these mechanisms has remained challenging due to the transient and dynamic nature of protein-protein interactions, the vast diversity of substrates, and substrate instability [K. Finzel, D. J. Lee, M. D. Burkart, ChemBioChem 16, 528-547 (2015)]. Here we demonstrate a unique communication mechanism between the ACP and partner enzymes using solution NMR spectroscopy and molecular dynamics to elucidate allostery that is dependent on fatty acid chain length. We demonstrate that partner enzymes can allosterically distinguish between chain lengths via protein-protein interactions as structural features of substrate sequestration are translated from within the ACP four-helical bundle to the protein surface, without the need for stochastic chain flipping. These results illuminate details of cargo communication by the ACP that can serve as a foundation for engineering carrier protein-dependent pathways for specific, desired products.