Search

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 (threefold) the maximum steady-state actin-activated ATPase activity (k cat) and decreases (eightfold) the actin concentration at which ATPase is one-half maximal (K ATPase). We also found a twofold to fourfold increase in the actin-activated power stroke and phosphate release rate constants at 30 μM actin, which overall enhanced the duty ratio threefold. 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.

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.

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 (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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Researchers developed a method called fast Fourier transform (FFT)-based modeling of atomistic protein–protein interactions applied to cross second virial coefficient B23(FMAPB23) to characterize protein kinase A (PKA) subunits as macromolecular regulators of PKA RIα liquid-liquid phase separation. The cross second virial coefficient B23 was calculated for Mg2+-bound PKAcat [denoted as PKAcat(Mg2+)] and cAMP-bound RIα [denoted as RIa(cAMP)]. Five structures for the first molecule were taken from PDB: 6NO7 chains E and G, 3O7L chain B, and 4NTT chains A and B, while two structures for the second molecule were taken from 4MX3 chains A and B. The results from the total of ten runs were used to calculate a mean value and a standard error of the mean.

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. Płoskoń 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.

We develop a hybrid approach that combines the Monte Carlo (MC) method, a variational implicit-solvent model (VISM), and a binary level-set method for the simulation of biomolecular binding in an aqueous solvent. The solvation free energy for the biomolecular complex is estimated by minimizing the VISM free-energy functional of all possible solute-solvent interfaces that are used as dielectric boundaries. This functional consists of the solute volumetric, solute-solvent interfacial, solute-solvent van der Waals interaction, and electrostatic free energy. A technique of shifting the dielectric boundary is used to accurately predict the electrostatic part of the solvation free energy. Minimizing such a functional in each MC move is made possible by our new and fast binary level-set method. This method is based on the approximation of surface area by the convolution of an indicator function with a compactly supported kernel and is implemented by simple flips of numerical grid cells locally around the solute-solvent interface. We apply our approach to the p53-MDM2 system for which the two molecules are approximated by rigid bodies. Our efficient approach captures some of the poses before the final bound state. All-atom molecular dynamics simulations with most of such poses quickly reach the final bound state. Our work is a new step toward realistic simulations of biomolecular interactions. With further improvement of coarse graining and MC sampling, and combined with other models, our hybrid approach can be used to study the free-energy landscape and kinetic pathways of ligand binding to proteins.

The atomic structure of the complete myosin tail within thick filaments isolated from Lethocerus indicus flight muscle is described and compared to crystal structures of recombinant, human cardiac myosin tail segments. Overall, the agreement is good with three exceptions: the proximal S2, in which the filament has heads attached but the crystal structure doesn't, and skip regions 2 and 4. At the head-tail junction, the tail α-helices are asymmetrically structured encompassing well-defined unfolding of 12 residues for one myosin tail, ∼4 residues of the other, and different degrees of α-helix unwinding for both tail α-helices, thereby providing an atomic resolution description of coiled-coil "uncoiling" at the head-tail junction. Asymmetry is observed in the nonhelical C termini; one C-terminal segment is intercalated between ribbons of myosin tails, the other apparently terminating at Skip 4 of another myosin tail. Between skip residues, crystal and filament structures agree well. Skips 1 and 3 also agree well and show the expected α-helix unwinding and coiled-coil untwisting in response to skip residue insertion. Skips 2 and 4 are different. Skip 2 is accommodated in an unusual manner through an increase in α-helix radius and corresponding reduction in rise/residue. Skip 4 remains helical in one chain, with the other chain unfolded, apparently influenced by the acidic myosin C terminus. The atomic model may shed some light on thick filament mechanosensing and is a step in understanding the complex roles that thick filaments of all species undergo during muscle contraction.

Influenza neuraminidase is an important drug target. Glycans are present on neuraminidase and are generally considered to inhibit antibody binding via their glycan shield. In this work, we studied the effect of glycans on the binding kinetics of antiviral drugs to the influenza neuraminidase. We created all-atom in silico systems of influenza neuraminidase with experimentally derived glycoprofiles consisting of four systems with different glycan conformations and one system without glycans. Using Brownian dynamics simulations, we observe a two- to eightfold decrease in the rate of ligand binding to the primary binding site of neuraminidase due to the presence of glycans. These glycans are capable of covering much of the surface area of neuraminidase, and the ligand binding inhibition is derived from glycans sterically occluding the primary binding site on a neighboring monomer. Our work also indicates that drugs preferentially bind to the primary binding site (i.e., the active site) over the secondary binding site, and we propose a binding mechanism illustrating this. These results help illuminate the complex interplay between glycans and ligand binding on the influenza membrane protein neuraminidase.

2'-deoxy-ATP (dATP) is a naturally occurring small molecule that has shown promise as a therapeutic because it significantly increases cardiac myocyte force development even at low dATP/ATP ratios. To investigate mechanisms by which dATP alters myosin crossbridge dynamics, we used Brownian dynamics simulations to calculate association rates between actin and ADP- or dADP-bound myosin. These rates were then directly incorporated in a mechanistic Monte Carlo Markov Chain model of cooperative sarcomere contraction. A unique combination of increased powerstroke and detachment rates was required to match experimental steady-state and kinetic data for dATP force production in rat cardiac myocytes when the myosin attachment rate in the model was constrained by the results of a Brownian dynamics simulation. Nearest-neighbor cooperativity was seen to contribute to, but not fully explain, the steep relationship between dATP/ATP ratio and steady-state force-development observed at lower dATP concentrations. Dynamic twitch simulations performed using measured calcium transients as inputs showed that the effects of dATP on the crossbridge alone were not sufficient to explain experimentally observed enhancement of relaxation kinetics by dATP treatment. Hence, dATP may also affect calcium handling even at low concentrations. By enabling the effects of dATP on sarcomere mechanics to be predicted, this multi-scale modeling framework may elucidate the molecular mechanisms by which dATP can have therapeutic effects on cardiac contractile dysfunction.

Fatty acid synthases (FASs) and polyketide synthases (PKSs) iteratively elongate and often reduce two-carbon ketide units in de novo fatty acid and polyketide biosynthesis. Cycles of chain extensions in FAS and PKS are initiated by an acyltransferase (AT), which loads monomer units onto acyl carrier proteins (ACPs), small, flexible proteins that shuttle covalently linked intermediates between catalytic partners. Formation of productive ACP-AT interactions is required for catalysis and specificity within primary and secondary FAS and PKS pathways. Here, we use the Escherichia coli FAS AT, FabD, and its cognate ACP, AcpP, to interrogate type II FAS ACP-AT interactions. We utilize a covalent crosslinking probe to trap transient interactions between AcpP and FabD to elucidate the X-ray crystal structure of a type II ACP-AT complex. Our structural data are supported using a combination of mutational, crosslinking, and kinetic analyses, and long-timescale molecular dynamics (MD) simulations. Together, these complementary approaches reveal key catalytic features of FAS ACP-AT interactions. These mechanistic inferences suggest that AcpP adopts multiple, productive conformations at the AT binding interface, allowing the complex to sustain high transacylation rates. Furthermore, MD simulations support rigid body subdomain motions within the FabD structure that may play a key role in AT activity and substrate selectivity.

Carbon-carbon bond forming reactions are essential transformations in natural product biosynthesis. During de novo fatty acid and polyketide biosynthesis, β-ketoacyl-acyl carrier protein (ACP) synthases (KS), catalyze this process via a decarboxylative Claisen-like condensation reaction. KSs must recognize multiple chemically distinct ACPs and choreograph a ping-pong mechanism, often in an iterative fashion. Here, we report crystal structures of substrate mimetic bearing ACPs in complex with the elongating KSs from Escherichia coli, FabF and FabB, in order to better understand the stereochemical features governing substrate discrimination by KSs. Complemented by molecular dynamics (MD) simulations and mutagenesis studies, these structures reveal conformational states accessed during KS catalysis. These data taken together support a gating mechanism that regulates acyl-ACP binding and substrate delivery to the KS active site. Two active site loops undergo large conformational excursions during this dynamic gating mechanism and are likely evolutionarily conserved features in elongating KSs.

Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the finite element method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open.

Advances in imaging methods such as electron microscopy, tomography, and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution of partial differential equations. In this work, we outline the steps for a naïve user to approach the Geometry-preserving Adaptive MeshER software version 2, a mesh generation code written in C++ designed to convert structural data sets to realistic geometric meshes while preserving the underlying shapes. We present two example cases: 1) mesh generation at the subcellular scale as informed by electron tomography and 2) meshing a protein with a structure from x-ray crystallography. We further demonstrate that the meshes generated by the Geometry-preserving Adaptive MeshER software are suitable for use with numerical methods. Together, this collection of libraries and tools simplifies the process of constructing realistic geometric meshes from structural biology data.

We use Brownian-Langevin dynamics principles to derive a coarse-graining multiscale myofilament model that can describe the thin-filament activation process during contraction. The model links atomistic molecular simulations of protein-protein interactions in the thin-filament regulatory unit to sarcomere-level activation dynamics. We first calculate the molecular interaction energy between tropomyosin and actin surface using Brownian dynamics simulations. This energy profile is then generalized to account for the observed tropomyosin transitions between its regulatory stable states. The generalized energy landscape then served as a basis for developing a filament-scale model using Langevin dynamics. This integrated analysis, spanning molecular to thin-filament scales, is capable of tracking the events of the tropomyosin conformational changes as it moves over the actin surface. The tropomyosin coil with flexible overlap regions between adjacent tropomyosins is represented in the model as a system of coupled stochastic ordinary differential equations. The proposed multiscale approach provides a more detailed molecular connection between tropomyosin dynamics, the trompomyosin-actin interaction-energy landscape, and the generated force by the sarcomere.

We introduce CASC: a new, modern, and header-only C++ library which provides a data structure to represent arbitrary dimension abstract simplicial complexes (ASC) with user-defined classes stored directly on the simplices at each dimension. This is accomplished by using the latest C++ language features including variadic template parameters introduced in C++11 and automatic function return type deduction from C++14. Effectively CASC decouples the representation of the topology from the interactions of user data. We present the innovations and design principles of the data structure and related algorithms. This includes a metadata aware decimation algorithm which is general for collapsing simplices of any dimension. We also present an example application of this library to represent an orientable surface mesh.

Acyl carrier proteins (ACP)s transport intermediates through many primary and secondary metabolic pathways. Studying the effect of substrate identity on ACP structure has been hindered by the lability of the thioester bond that attaches acyl substrates to the 4'-phosphopantetheine cofactor of ACP. Here we show that an acyl acyl-carrier protein synthetase (AasS) can be used in real time to shift the hydrolysis equilibrium toward favoring acyl-ACP during solution NMR spectroscopy. Only 0.005 molar equivalents of AasS enables 1 week of stability to palmitoyl-AcpP from Escherichia coli. 2D NMR spectra enabled with this method revealed that the tethered palmitic acid perturbs nearly every secondary structural region of AcpP. This technique will allow previously unachievable structural studies of unstable acyl-ACP species, contributing to the understanding of these complex biosynthetic pathways.

At the center of many complex biosynthetic pathways, the acyl carrier protein (ACP) shuttles substrates to appropriate enzymatic partners to produce fatty acids and polyketides. Carrier proteins covalently tether their cargo via a thioester linkage to a phosphopantetheine cofactor. Due to the labile nature of this linkage, chemoenzymatic methods have been developed that involve replacement of the thioester with a more stable amide or ester bond. We explored the importance of the thioester bond to the structure of the carrier protein by using solution NMR spectroscopy and molecular dynamics simulations. Remarkably, the replacement of sulfur with other heteroatoms results in significant structural changes, thus suggesting more rigorous selections of isosteric substitutes is needed.

Noncoding RNA (ncRNA) has a key role in regulating gene expression, mediating fundamental processes and diseases via a variety of yet unknown mechanisms. Here, we review recent applications of conventional and enhanced Molecular Dynamics (MD) simulations methods to address the mechanistic function of large biomolecular systems that are tightly involved in the ncRNA function and that are of key importance in life sciences. This compendium focuses of three biomolecular systems, namely the CRISPR-Cas9 genome editing machinery, group II intron ribozyme and the ribonucleoprotein complex of the spliceosome, which edit and process ncRNA. We show how the application of a novel accelerated MD simulations method has been key in disclosing the conformational transitions underlying RNA binding in the CRISPR-Cas9 complex, suggesting a mechanism for RNA recruitment and clarifying the conformational changes required for attaining genome editing. As well, we discuss the use of mixed quantum-classical MD simulations in deciphering the catalytic mechanism of RNA splicing as operated by group II intron ribozyme, one of the largest ncRNA structures crystallized so far. Finally, we debate the future challenges and opportunities in the field, discussing the recent application of MD simulations for unraveling the functional biophysics of the spliceosome, a multi-mega Dalton complex of proteins and small nuclear RNAs that performs RNA splicing in humans. This showcase of applications highlights the current talent of MD simulations to dissect atomic-level details of complex biomolecular systems instrumental for the design of finely engineered genome editing machines. As well, this review aims at inspiring future investigations of several other ncRNA regulatory systems, such as micro and small interfering RNAs, which achieve their function and specificity using RNA-based recognition and targeting strategies.

Fatty acid biosynthesis in α- and γ-proteobacteria requires two functionally distinct dehydratases, FabA and FabZ. Here, mechanistic cross-linking facilitates the structural characterization of a stable hexameric complex of six Escherichia coli FabZ dehydratase subunits with six AcpP acyl carrier proteins. The crystal structure sheds light on the divergent substrate selectivity of FabA and FabZ by revealing distinct architectures of the binding pocket. Molecular dynamics simulations demonstrate differential biasing of substrate orientations and conformations within the active sites of FabA and FabZ such that FabZ is preorganized to catalyze only dehydration, while FabA is primed for both dehydration and isomerization.

To maintain the lipid asymmetry of the cell envelope in Gram-negative bacteria, the MlaC protein serves as a lipid transfer factor and delivers phospholipids from the outer to the inner membrane. A strategy of antibiotic discovery is to design a proper compound that can tightly bind to the MlaC protein and inhibit the MlaC function. In this study, we performed virtual screening on multiple MlaC structures obtained from molecular dynamics simulations to identify potential MlaC binders. Our results suggested that clorobiocin is a compound that could bind to the MlaC protein. Through the comparison of the bound geometry between clorobiocin and novobiocin, we pointed out that the methyl-pyrrole group of the noviose sugar in clorobiocin forms hydrophobic interactions with amino acids in the phospholipid binding pocket, which allows the compound to bind deep in the active site. This also explains why clorobiocin shows a tighter binding affinity than novobiocin. Our study highlights a practical path of antibiotic development against Gram-negative bacteria.

Allosteric regulation plays an important role in many biological processes, such as signal transduction, transcriptional regulation, and metabolism. Allostery is rooted in the fundamental physical properties of macromolecular systems, but its underlying mechanisms are still poorly understood. A collection of contributions to a recent interdisciplinary CECAM (Center Européen de Calcul Atomique et Moléculaire) workshop is used here to provide an overview of the progress and remaining limitations in the understanding of the mechanistic foundations of allostery gained from computational and experimental analyses of real protein systems and model systems. The main conceptual frameworks instrumental in driving the field are discussed. We illustrate the role of these frameworks in illuminating molecular mechanisms and explaining cellular processes, and describe some of their promising practical applications in engineering molecular sensors and informing drug design efforts.

Since the discovery of the DNA double helix, the main molecular repository of genetic information, scientists have been struggling to find ways to efficiently manipulate genes. The ability to mark, modify, or regulate specific sequences of DNA in a controlled fashion is of key importance because of the ways that gene editing could be used to improve human life. For example, genetic therapies are being developed to permanently cure cancer and other life-threatening diseases.

Beta-secretase 1 (BACE-1) is an aspartyl protease implicated in the overproduction of β-amyloid fibrils responsible for Alzheimer disease. The process of β-amyloid genesis is known to be pH dependent, with an activity peak between solution pH of 3.5 and 5.5. We have studied the pH-dependent dynamics of BACE-1 to better understand the pH dependent mechanism. We have implemented support for graphics processor unit (GPU) accelerated constant pH molecular dynamics within the AMBER molecular dynamics software package and employed this to determine the relative population of different aspartyl dyad protonation states in the pH range of greatest β-amyloid production, followed by conventional molecular dynamics to explore the differences among the various aspartyl dyad protonation states. We observed a difference in dynamics between double-protonated, mono-protonated, and double-deprotonated states over the known pH range of higher activity. These differences include Tyr 71-aspartyl dyad proximity and active water lifetime. This work indicates that Tyr 71 stabilizes catalytic water in the aspartyl dyad active site, enabling BACE-1 activity.

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. [ABSTRACT FROM AUTHOR]

Abstract Objective Recent advances in electron microscopy have, for the first time, enabled imaging of single cells in 3D at a nanometer length scale resolution. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. However, this will require a system for going from EM images to 3D volume meshes which can be used in finite element simulations. Methods In this paper, we develop an end-to-end pipeline for this task by adapting and extending computer graphics mesh processing and smoothing algorithms. Our workflow makes use of our recently rewritten mesh processing software, GAMer 2 , which implements several mesh conditioning algorithms and serves as a platform to connect different pipeline steps. Results We apply this pipeline to a series of electron micrographs of dendrite morphology explored at three different length scales and show that the resultant meshes are suitable for finite element simulations. Conclusion Our pipeline, which consists of free and open-source community driven tools, is a step towards routine physical simulations of biological processes in realistic geometries. Significance We posit that a new frontier at the intersection of computational technologies and single cell biology is now open. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery.

Nonribosomal peptide synthetases (NRPSs) are responsible for the synthesis of a variety of bioactive natural products with clinical and economic significance. Interestingly, these large multimodular enzyme machineries incorporate nonproteinogenic d-amino acids through the use of auxiliary epimerization domains, converting l-amino acids into d-amino acids that impart into the resulting natural products unique bioactivity and resistance to proteases. Due to the large and complex nature of NRPSs, several questions remain unanswered about the mechanism of the catalytic domain reactions. We have investigated the use of mechanism-based crosslinkers to probe the mechanism of an epimerization domain in gramicidin S biosynthesis. In addition, MD simulations were performed, showcasing the possible roles of catalytic residues within the epimerization domain.

The atomic-level mechanisms that coordinate ligand release from protein pockets are only known for a handful of proteins. Here, we report results from accelerated molecular dynamics simulations for benzene dissociation from the buried cavity of the T4 lysozyme Leu99Ala mutant (L99A). In these simulations, benzene is released through a previously characterized, sparsely populated room-temperature excited state of the mutant, explaining the coincidence for experimentally measured benzene off rate and apo protein slow-timescale NMR relaxation rates between ground and excited states. The path observed for benzene egress is a multistep ligand migration from the buried cavity to ultimate release through an opening between the F/G-, H-, and I-helices and requires a number of cooperative multiresidue and secondary-structure rearrangements within the C-terminal domain of L99A. These rearrangements are identical to those observed along the ground state to excited state transitions characterized by molecular dynamic simulations run on the Anton supercomputer. Analyses of the molecular properties of the residues lining the egress path suggest that protein surface electrostatic potential may play a role in the release mechanism. Simulations of wild-type T4 lysozyme also reveal that benzene-egress-associated dynamics in the L99A mutant are potentially exaggerations of the substrate-processivity-related dynamics of the wild type.