Whole-Cell Models and Simulations in Molecular Detail.

Comprehensive data about the composition and structure of cellular components have enabled the construction of quantitative whole-cell models. While kinetic network-type models have been established, it is also becoming possible to build physical, molecular-level models of cellular environments. This review outlines challenges in constructing and simulating such models and discusses near- and long-term opportunities for developing physical whole-cell models that can connect molecular structure with biological function.

Non-Abelian topological order and anyons on a trapped-ion processor.

Non-Abelian topological order is a coveted state of matter with remarkable properties, including quasiparticles that can remember the sequence in which they are exchanged1-4. These anyonic excitations are promising building blocks of fault-tolerant quantum computers5,6. However, despite extensive efforts, non-Abelian topological order and its excitations have remained elusive, unlike the simpler quasiparticles or defects in Abelian topological order. Here we present the realization of non-Abelian topological order in the wavefunction prepared in a quantum processor and demonstrate control of its anyons. Using an adaptive circuit on Quantinuum's H2 trapped-ion quantum processor, we create the ground-state wavefunction of D4 topological order on a kagome lattice of 27 qubits, with fidelity per site exceeding 98.4 per cent. By creating and moving anyons along Borromean rings in spacetime, anyon interferometry detects an intrinsically non-Abelian braiding process. Furthermore, tunnelling non-Abelions around a torus creates all 22 ground states, as well as an excited state with a single anyon-a peculiar feature of non-Abelian topological order. This work illustrates the counterintuitive nature of non-Abelions and enables their study in quantum devices.

Structural basis of transcription: RNA polymerase II substrate binding and metal coordination using a free-electron laser.

Catalysis and translocation of multisubunit DNA-directed RNA polymerases underlie all cellular mRNA synthesis. RNA polymerase II (Pol II) synthesizes eukaryotic pre-mRNAs from a DNA template strand buried in its active site. Structural details of catalysis at near-atomic resolution and precise arrangement of key active site components have been elusive. Here, we present the free-electron laser (FEL) structures of a matched ATP-bound Pol II and the hyperactive Rpb1 T834P bridge helix (BH) mutant at the highest resolution to date. The radiation-damage-free FEL structures reveal the full active site interaction network, including the trigger loop (TL) in the closed conformation, bonafide occupancy of both site A and B Mg2+, and, more importantly, a putative third (site C) Mg2+ analogous to that described for some DNA polymerases but not observed previously for cellular RNA polymerases. Molecular dynamics (MD) simulations of the structures indicate that the third Mg2+ is coordinated and stabilized at its observed position. TL residues provide half of the substrate binding pocket while multiple TL/BH interactions induce conformational changes that could allow translocation upon substrate hydrolysis. Consistent with TL/BH communication, a FEL structure and MD simulations of the T834P mutant reveal rearrangement of some active site interactions supporting potential plasticity in active site function and long-distance effects on both the width of the central channel and TL conformation, likely underlying its increased elongation rate at the expense of fidelity.

Transferable deep generative modeling of intrinsically disordered protein conformations.

Intrinsically disordered proteins have dynamic structures through which they play key biological roles. The elucidation of their conformational ensembles is a challenging problem requiring an integrated use of computational and experimental methods. Molecular simulations are a valuable computational strategy for constructing structural ensembles of disordered proteins but are highly resource-intensive. Recently, machine learning approaches based on deep generative models that learn from simulation data have emerged as an efficient alternative for generating structural ensembles. However, such methods currently suffer from limited transferability when modeling sequences and conformations absent in the training data. Here, we develop a novel generative model that achieves high levels of transferability for intrinsically disordered protein ensembles. The approach, named idpSAM, is a latent diffusion model based on transformer neural networks. It combines an autoencoder to learn a representation of protein geometry and a diffusion model to sample novel conformations in the encoded space. IdpSAM was trained on a large dataset of simulations of disordered protein regions performed with the ABSINTH implicit solvent model. Thanks to the expressiveness of its neural networks and its training stability, idpSAM faithfully captures 3D structural ensembles of test sequences with no similarity in the training set. Our study also demonstrates the potential for generating full conformational ensembles from datasets with limited sampling and underscores the importance of training set size for generalization. We believe that idpSAM represents a significant progress in transferable protein ensemble modeling through machine learning.

Robust Measurements of n-Point Correlation Functions of Driven-Dissipative Quantum Systems on a Digital Quantum Computer.

We propose and demonstrate a unified hierarchical method to measure n-point correlation functions that can be applied to driven, dissipative, or otherwise open or nonequilibrium quantum systems. In this method, the time evolution of the system is repeatedly interrupted by interacting an ancilla qubit with the system through a controlled operation, and measuring the ancilla immediately afterward. We discuss the robustness of this method as compared to other ancilla-based interferometric techniques (such as the Hadamard test), and highlight its advantages for near-term quantum simulations of open quantum systems. We implement the method on a quantum computer in order to measure single-particle Green's functions of a driven-dissipative fermionic system. This Letter shows that dynamical correlation functions for driven-dissipative systems can be robustly measured with near-term quantum computers.

Varying molecular interactions explain aspects of crowder-dependent enzyme function of a viral protease.

Biochemical processes in cells, including enzyme-catalyzed reactions, occur in crowded conditions with various background macromolecules occupying up to 40% of cytoplasm's volume. Viral enzymes in the host cell also encounter such crowded conditions as they often function at the endoplasmic reticulum membranes. We focus on an enzyme encoded by the hepatitis C virus, the NS3/4A protease, which is crucial for viral replication. We have previously found experimentally that synthetic crowders, polyethylene glycol (PEG) and branched polysucrose (Ficoll), differently affect the kinetic parameters of peptide hydrolysis catalyzed by NS3/4A. To gain understanding of the reasons for such behavior, we perform atomistic molecular dynamics simulations of NS3/4A in the presence of either PEG or Ficoll crowders and with and without the peptide substrates. We find that both crowder types make nanosecond long contacts with the protease and slow down its diffusion. However, they also affect the enzyme structural dynamics; crowders induce functionally relevant helical structures in the disordered parts of the protease cofactor, NS4A, with the PEG effect being more pronounced. Overall, PEG interactions with NS3/4A are slightly stronger but Ficoll forms more hydrogen bonds with NS3. The crowders also interact with substrates; we find that the substrate diffusion is reduced much more in the presence of PEG than Ficoll. However, contrary to NS3, the substrate interacts more strongly with Ficoll than with PEG crowders, with the substrate diffusion being similar to crowder diffusion. Importantly, crowders also affect the substrate-enzyme interactions. We observe that both PEG and Ficoll enhance the presence of substrates near the active site, especially near catalytic H57 but Ficoll crowders increase substrate binding more than PEG molecules.

Characterization of RNA polymerase II trigger loop mutations using molecular dynamics simulations and machine learning.

Catalysis and fidelity of multisubunit RNA polymerases rely on a highly conserved active site domain called the trigger loop (TL), which achieves roles in transcription through conformational changes and interaction with NTP substrates. The mutations of TL residues cause distinct effects on catalysis including hypo- and hyperactivity and altered fidelity. We applied molecular dynamics simulation (MD) and machine learning (ML) techniques to characterize TL mutations in the Saccharomyces cerevisiae RNA Polymerase II (Pol II) system. We did so to determine relationships between individual mutations and phenotypes and to associate phenotypes with MD simulated structural alterations. Using fitness values of mutants under various stress conditions, we modeled phenotypes along a spectrum of continual values. We found that ML could predict the phenotypes with 0.68 R2 correlation from amino acid sequences alone. It was more difficult to incorporate MD data to improve predictions from machine learning, presumably because MD data is too noisy and possibly incomplete to directly infer functional phenotypes. However, a variational auto-encoder model based on the MD data allowed the clustering of mutants with different phenotypes based on structural details. Overall, we found that a subset of loss-of-function (LOF) and lethal mutations tended to increase distances of TL residues to the NTP substrate, while another subset of LOF and lethal substitutions tended to confer an increase in distances between TL and bridge helix (BH). In contrast, some of the gain-of-function (GOF) mutants appear to cause disruption of hydrophobic contacts among TL and nearby helices.

Stability and deformation of biomolecular condensates under the action of shear flow.

Biomolecular condensates play a key role in cytoplasmic compartmentalization and cell functioning. Despite extensive research on the physico-chemical, thermodynamic, or crowding aspects of the formation and stabilization of the condensates, one less studied feature is the role of external perturbative fluid flow. In fact, in living cells, shear stress may arise from streaming or active transport processes. Here, we investigate how biomolecular condensates are deformed under different types of shear flows. We first model Couette flow perturbations via two-way coupling between the condensate dynamics and fluid flow by deploying Lattice Boltzmann Molecular Dynamics. We then show that a simplified approach where the shear flow acts as a static perturbation (one-way coupling) reproduces the main features of the condensate deformation and dynamics as a function of the shear rate. With this approach, which can be easily implemented in molecular dynamics simulations, we analyze the behavior of biomolecular condensates described through residue-based coarse-grained models, including intrinsically disordered proteins and protein/RNA mixtures. At lower shear rates, the fluid triggers the deformation of the condensate (spherical to oblated object), while at higher shear rates, it becomes extremely deformed (oblated or elongated object). At very high shear rates, the condensates are fragmented. We also compare how condensates of different sizes and composition respond to shear perturbation, and how their internal structure is altered by external flow. Finally, we consider the Poiseuille flow that realistically models the behavior in microfluidic devices in order to suggest potential experimental designs for investigating fluid perturbations in vitro.

Biosynthesis and trafficking of heme o and heme a: new structural insights and their implications for reaction mechanisms and prenylated heme transfer.

Aerobic respiration is a key energy-producing pathway in many prokaryotes and virtually all eukaryotes. The final step of aerobic respiration is most commonly catalyzed by heme-copper oxidases embedded in the cytoplasmic or mitochondrial membrane. The majority of these terminal oxidases contain a prenylated heme (typically heme a or occasionally heme o) in the active site. In addition, many heme-copper oxidases, including mitochondrial cytochrome c oxidases, possess a second heme a cofactor. Despite the critical role of heme a in the electron transport chain, the details of the mechanism by which heme b, the prototypical cellular heme, is converted to heme o and then to heme a remain poorly understood. Recent structural investigations, however, have helped clarify some elements of heme a biosynthesis. In this review, we discuss the insight gained from these advances. In particular, we present a new structural model of heme o synthase (HOS) based on distance restraints from inferred coevolutionary relationships and refined by molecular dynamics simulations that are in good agreement with the experimentally determined structures of HOS homologs. We also analyze the two structures of heme a synthase (HAS) that have recently been solved by other groups. For both HOS and HAS, we discuss the proposed catalytic mechanisms and highlight how new insights into the heme-binding site locations shed light on previously obtained biochemical data. Finally, we explore the implications of the new structural data in the broader context of heme trafficking in the heme a biosynthetic pathway and heme-copper oxidase assembly.

Multi-state modeling of G-protein coupled receptors at experimental accuracy.

The family of G-protein coupled receptors (GPCRs) is one of the largest protein families in the human genome. GPCRs transduct chemical signals from extracellular to intracellular regions via a conformational switch between active and inactive states upon ligand binding. While experimental structures of GPCRs remain limited, high-accuracy computational predictions are now possible with AlphaFold2. However, AlphaFold2 only predicts one state and is biased toward either the active or inactive conformation depending on the GPCR class. Here, a multi-state prediction protocol is introduced that extends AlphaFold2 to predict either active or inactive states at very high accuracy using state-annotated templated GPCR databases. The predicted models accurately capture the main structural changes upon activation of the GPCR at the atomic level. For most of the benchmarked GPCRs (10 out of 15), models in the active and inactive states were closer to their corresponding activation state structures. Median RMSDs of the transmembrane regions were 1.12 Å and 1.41 Å for the active and inactive state models, respectively. The models were more suitable for protein-ligand docking than the original AlphaFold2 models and template-based models. Finally, our prediction protocol predicted accurate GPCR structures and GPCR-peptide complex structures in GPCR Dock 2021, a blind GPCR-ligand complex modeling competition. We expect that high accuracy GPCR models in both activation states will promote understanding in GPCR activation mechanisms and drug discovery for GPCRs. At the time, the new protocol paves the way towards capturing the dynamics of proteins at high-accuracy via machine-learning methods.

Entanglement from Tensor Networks on a Trapped-Ion Quantum Computer.

The ability to selectively measure, initialize, and reuse qubits during a quantum circuit enables a mapping of the spatial structure of certain tensor-network states onto the dynamics of quantum circuits, thereby achieving dramatic resource savings when simulating quantum systems with limited entanglement. We experimentally demonstrate a significant benefit of this approach to quantum simulation: the entanglement structure of an infinite system-specifically the half-chain entanglement spectrum-is conveniently encoded within a small register of "bond qubits" and can be extracted with relative ease. Using Honeywell's model H0 quantum computer equipped with selective midcircuit measurement and reset, we quantitatively determine the near-critical entanglement entropy of a correlated spin chain directly in the thermodynamic limit and show that its phase transition becomes quickly resolved upon expanding the bond-qubit register.

Short disordered protein segment regulates cross-species transmission of a yeast prion.

Soluble prion proteins contingently encounter foreign prion aggregates, leading to cross-species prion transmission. However, how its efficiency is regulated by structural fluctuation of the host soluble prion protein remains unsolved. In the present study, through the use of two distantly related yeast prion Sup35 proteins, we found that a specific conformation of a short disordered segment governs interspecies prion transmissibility. Using a multidisciplinary approach including high-resolution NMR and molecular dynamics simulation, we identified critical residues within this segment that allow interspecies prion transmission in vitro and in vivo, by locally altering dynamics and conformation of soluble prion proteins. Remarkably, subtle conformational differences caused by a methylene group between asparagine and glutamine sufficed to change the short segment structure and substantially modulate the cross-seeding activity. Thus, our findings uncover how conformational dynamics of the short segment in the host prion protein impacts cross-species prion transmission. More broadly, our study provides mechanistic insights into cross-seeding between heterologous proteins.

Clustering and dynamics of crowded proteins near membranes and their influence on membrane bending.

Atomistic molecular dynamics simulations of concentrated protein solutions in the presence of a phospholipid bilayer are presented to gain insights into the dynamics and interactions at the cytosol-membrane interface. The main finding is that proteins that are not known to specifically interact with membranes are preferentially excluded from the membrane, leaving a depletion zone near the membrane surface. As a consequence, effective protein concentrations increase, leading to increased protein contacts and clustering, whereas protein diffusion becomes faster near the membrane for proteins that do occasionally enter the depletion zone. Since protein-membrane contacts are infrequent and short-lived in this study, the structure of the lipid bilayer remains largely unaffected by the crowded protein solution, but when proteins do contact lipid head groups, small but statistically significant local membrane curvature is induced, on average.

Crowding affects structural dynamics and contributes to membrane association of the NS3/4A complex.

Using molecular dynamics simulations, we describe how crowded environments affect the internal dynamics and diffusion of the hepatitis C virus proteases NS3/4A. This protease plays a key role in viral replication and is successfully used as a target for antiviral treatment. The NS3 enzyme requires a peptide cofactor, called NS4A, with its central part interacting with the NS3 ß-sheet, and flexible, protruding terminal tails that are unstructured in water solution. The simulations describe the enzyme and water molecules at atomistic resolution, whereas crowders are modeled via either all-atom or coarse-grained models to emphasize different aspects of crowding. Crowders reflect the polyethylene glycol (PEG) molecules used in the experiments to mimic the crowded surrounding. A bead-shell model of folded coarse-grained PEG molecules considers mainly the excluded volume effect, whereas all-atom PEG models afford more protein-like crowder interactions. Circular dichroism spectroscopy experiments of the NS4A N-terminal tail show that a helical structure is formed in the presence of PEG crowders. The simulations suggest that crowding may assist in the formation of an NS4A helical fragment, positioned exactly where a transmembrane helix would fold upon the NS4A contact with the membrane. In addition, partially interactive PEGs help the NS4A N-tail to detach from the protease surface, thus enabling the process of helix insertion and potentially helping the virus establish a replication machinery needed to produce new viruses. Results point to an active role of crowding in assisting structural changes in disordered protein fragments that are necessary for their biological function.

Physics-based protein structure refinement in the era of artificial intelligence.

Protein structure refinement is the last step in protein structure prediction pipelines. Physics-based refinement via molecular dynamics (MD) simulations has made significant progress during recent years. During CASP14, we tested a new refinement protocol based on an improved sampling strategy via MD simulations. MD simulations were carried out at an elevated temperature (360 K). An optimized use of biasing restraints and the use of multiple starting models led to enhanced sampling. The new protocol generally improved the model quality. In comparison with our previous protocols, the CASP14 protocol showed clear improvements. Our approach was successful with most initial models, many based on deep learning methods. However, we found that our approach was not able to refine machine-learning models from the AlphaFold2 group, often decreasing already high initial qualities. To better understand the role of refinement given new types of models based on machine-learning, a detailed analysis via MD simulations and Markov state modeling is presented here. We continue to find that MD-based refinement has the potential to improve AI predictions. We also identified several practical issues that make it difficult to realize that potential. Increasingly important is the consideration of inter-domain and oligomeric contacts in simulations; the presence of large kinetic barriers in refinement pathways also continues to present challenges. Finally, we provide a perspective on how physics-based refinement could continue to play a role in the future for improving initial predictions based on machine learning-based methods.

New parallel computing algorithm of molecular dynamics for extremely huge scale biological systems.

In this paper, we address high performance extreme-scale molecular dynamics (MD) algorithm in the GENESIS software to perform cellular-scale molecular dynamics (MD) simulations with more than 100,000 CPU cores. It includes (1) the new algorithm of real-space nonbonded interactions maximizing the performance on ARM CPU architecture, (2) reciprocal-space nonbonded interactions minimizing communicational cost, (3) accurate temperature/pressure evaluations that allows a large time step, and (4) effective parallel file inputs/outputs (I/O) for MD simulations of extremely huge systems. The largest system that contains 1.6 billion atoms was simulated using MD with a performance of 8.30 ns/day on Fugaku supercomputer. It extends the available size and time of MD simulations to answer unresolved questions of biomacromolecules in a living cell.

Experimental accuracy in protein structure refinement via molecular dynamics simulations.

Refinement is the last step in protein structure prediction pipelines to convert approximate homology models to experimental accuracy. Protocols based on molecular dynamics (MD) simulations have shown promise, but current methods are limited to moderate levels of consistent refinement. To explore the energy landscape between homology models and native structures and analyze the challenges of MD-based refinement, eight test cases were studied via extensive simulations followed by Markov state modeling. In all cases, native states were found very close to the experimental structures and at the lowest free energies, but refinement was hindered by a rough energy landscape. Transitions from the homology model to the native states require the crossing of significant kinetic barriers on at least microsecond time scales. A significant energetic driving force toward the native state was lacking until its immediate vicinity, and there was significant sampling of off-pathway states competing for productive refinement. The role of recent force field improvements is discussed and transition paths are analyzed in detail to inform which key transitions have to be overcome to achieve successful refinement.

High-accuracy protein structures by combining machine-learning with physics-based refinement.

Protein structure prediction has long been available as an alternative to experimental structure determination, especially via homology modeling based on templates from related sequences. Recently, models based on distance restraints from coevolutionary analysis via machine learning to have significantly expanded the ability to predict structures for sequences without templates. One such method, AlphaFold, also performs well on sequences where templates are available but without using such information directly. Here we show that combining machine-learning based models from AlphaFold with state-of-the-art physics-based refinement via molecular dynamics simulations further improves predictions to outperform any other prediction method tested during the latest round of CASP. The resulting models have highly accurate global and local structures, including high accuracy at functionally important interface residues, and they are highly suitable as initial models for crystal structure determination via molecular replacement.

The endoplasmic reticulum acetyltransferases ATase1/NAT8B and ATase2/NAT8 are differentially regulated to adjust engagement of the secretory pathway.

Nε-lysine acetylation of nascent glycoproteins within the endoplasmic reticulum (ER) lumen regulates the efficiency of the secretory pathway. The ER acetylation machinery consists of the membrane transporter, acetyl-CoA transporter 1 (AT-1/SLC33A1), and two acetyltransferases, ATase1/NAT8B and ATase2/NAT8. Dysfunctional ER acetylation is associated with severe neurological diseases with duplication of AT-1/SLC33A1 being associated with autism spectrum disorder, intellectual disability, and dysmorphism. Neuron-specific AT-1 over-expression in the mouse alters neuron morphology and function, causing an autism-like phenotype, indicating that ER acetylation plays a key role in neurophysiology. As such, characterizing the molecular mechanisms that regulate the acetylation machinery could reveal critical information about its biology. By using structure-biochemistry approaches, we discovered that ATase1 and ATase2 share enzymatic properties but differ in that ATase1 is post-translationally regulated via acetylation. Furthermore, gene expression studies revealed that the promoters of AT-1, ATase1, and ATase2 contain functional binding sites for the neuron-related transcription factors cAMP response element-binding protein and the immediate-early genes c-FOS and c-JUN, and that ATase1 and ATase2 exhibit additional modes of transcriptional regulation relevant to aging and Alzheimer's disease. In vivo rodent gene expression experiments revealed that Atase2 is specifically induced following activity-dependent events. Finally, over-expression of either ATase1 or ATase2 was sufficient to increase the engagement of the secretory pathway in PC12 cells. Our results indicate important regulatory roles for ATase1 and ATase2 in neuron function with induction of ATase2 expression potentially serving as a critical event that adjusts the efficiency of the secretory pathway for activity-dependent neuronal functions.

Accelerating the Generalized Born with Molecular Volume and Solvent Accessible Surface Area Implicit Solvent Model Using Graphics Processing Units.

The generalized Born with molecular volume and solvent accessible surface area (GBMV2/SA) implicit solvent model provides an accurate description of molecular volume and has the potential to accurately describe the conformational equilibria of structured and disordered proteins. However, its broader application has been limited by the computational cost and poor scaling in parallel computing. Here, we report an efficient implementation of both the electrostatic and nonpolar components of GBMV2/SA on graphics processing unit (GPU) within the CHARMM/OpenMM module. The GPU-GBMV2/SA is numerically equivalent to the original CPU-GBMV2/SA. The GPU acceleration offers ~60- to 70-fold speedup on a single NVIDIA TITAN X (Pascal) graphics card for molecular dynamic simulations of both folded and unstructured proteins of various sizes. The current implementation can be further optimized to achieve even greater acceleration with minimal reduction on the numerical accuracy. The successful development of GPU-GBMV2/SA greatly facilitates its application to biomolecular simulations and paves the way for further development of the implicit solvent methodology. © 2019 Wiley Periodicals, Inc.