Regulation of the protein-conducting channel by a bound ribosome

Toward Understanding the Mechanism of Client-Selective Small Molecule Inhibitors of the Sec61 Translocon

The Sec61 translocon mediates the translocation of numerous, newly synthesized precursor proteins into the lumen of the endoplasmic reticulum or their integration into its membrane. Recently, structural biology revealed conformations of idle or substrate-engaged Sec61, and likewise its interactions with the accessory membrane proteins Sec62, Sec63, and TRAP, respectively. Several natural and synthetic small molecules have been shown to block Sec61-mediated protein translocation. Since this is a key step in protein biogenesis, broad inhibition is generally cytotoxic, which may be problematic for a putative drug target. Interestingly, several compounds exhibit client-selective modes of action, such that only translocation of certain precursor proteins was affected. Here, we discuss recent advances of structural biology, molecular modelling, and molecular screening that aim to use Sec61 as feasible drug target.

Structure of the human dopamine transporter and mechanisms of inhibition

The neurotransmitter dopamine has central roles in mood, appetite, arousal and movement1. Despite its importance in brain physiology and function, and as a target for illicit and therapeutic drugs, the human dopamine transporter (hDAT) and mechanisms by which it is inhibited by small molecules and Zn2+ are without a high-resolution structural context. Here we determine the structure of hDAT in a tripartite complex with the competitive inhibitor and cocaine analogue, (-)-2-β-carbomethoxy-3-β-(4-fluorophenyl)tropane2 (β-CFT), the non-competitive inhibitor MRS72923 and Zn2+ (ref. 4). We show how β-CFT occupies the central site, approximately halfway across the membrane, stabilizing the transporter in an outward-open conformation. MRS7292 binds to a structurally uncharacterized allosteric site, adjacent to the extracellular vestibule, sequestered underneath the extracellular loop 4 (EL4) and adjacent to transmembrane helix 1b (TM1b), acting as a wedge, precluding movement of TM1b and closure of the extracellular gate. A Zn2+ ion further stabilizes the outward-facing conformation by coupling EL4 to EL2, TM7 and TM8, thus providing specific insights into how Zn2+ restrains the movement of EL4 relative to EL2 and inhibits transport activity.

GOLEM: Automated and Robust Cryo-EM-Guided Ligand Docking with Explicit Water Molecules

A detailed understanding of ligand-protein interaction is essential for developing rational drug-design strategies. In recent years, technological advances in cryo-electron microscopy (cryo-EM) brought a new era to the structural determination of biological macromolecules and assemblies at high resolution, marking cryo-EM as a promising tool for studying ligand-protein interactions. However, even in high-resolution cryo-EM results, the densities for the bound small-molecule ligands are often of lower quality due to their relatively dynamic and flexible nature, frustrating their accurate coordinate assignment. To address the challenge of ligand modeling in cryo-EM maps, here we report the development of GOLEM (Genetic Optimization of Ligands in Experimental Maps), an automated and robust ligand docking method that predicts a ligand's pose and conformation in cryo-EM maps. GOLEM employs a Lamarckian genetic algorithm to perform a hybrid global/local search for exploring the ligand's conformational, orientational, and positional space. As an important feature, GOLEM explicitly considers water molecules and places them at optimal positions and orientations. GOLEM takes into account both molecular energetics and the correlation with the cryo-EM maps in its scoring function to optimally place the ligand. We have validated GOLEM against multiple cryo-EM structures with a wide range of map resolutions and ligand types, returning ligand poses in excellent agreement with the densities. As a VMD plugin, GOLEM is free of charge and accessible to the community. With these features, GOLEM will provide a valuable tool for ligand modeling in cryo-EM efforts toward drug discovery.

Transcriptomics of Lactobacillus paracasei: metabolism patterns and cellular responses under high-density culture conditions

Lactobacillus paracasei has significant potential for development and application in the environmental field, particularly in addressing malodor pollution. This study aims to investigate the cellular response of L. paracasei B1 under high-density culture conditions. The selected strain has previously shown effective deodorizing and bacteriostatic abilities. Transcriptomics techniques are employed to dissect the nutrient metabolism pattern of L. paracasei B1 and its response mechanism under environmental stress. The study characterizes the functions of key differentially expressed genes during growth before and after optimizing the culture conditions. The optimization of fermentation culture conditions provides a suitable growth environment for L. paracasei B1, inducing an enhancement of its phosphotransferase system for sugar source uptake and maintaining high levels of glycolysis and pyruvate metabolism. Consequently, the strain is able to grow and multiply rapidly. Under acid stress conditions, glycolysis and pyruvate metabolism are inhibited, and L. paracasei B1 generates additional energy through aerobic respiration to meet the energy demand. The two-component system and quorum sensing play roles in the response and regulation of L. paracasei B1 to adverse environments. The strain mitigates oxygen stress damage through glutathione metabolism, cysteine and methionine metabolism, base excision repair, and purine and pyrimidine metabolism. Additionally, the strain enhances lysine synthesis, the alanine, aspartate, and glutamate metabolic pathways, and relies on the ABC transport system to accumulate amino acid-compatible solutes to counteract acid stress and osmotic stress during pH regulation. These findings establish a theoretical basis for the further development and application of L. paracasei B1 for its productive properties.

Structures and membrane interactions of native serotonin transporter in complexes with psychostimulants

The serotonin transporter (SERT) is a member of the SLC6 neurotransmitter transporter family that mediates serotonin reuptake at presynaptic nerve terminals. SERT is the target of both therapeutic antidepressant drugs and psychostimulant substances such as cocaine and methamphetamines, which are small molecules that perturb normal serotonergic transmission by interfering with serotonin transport. Despite decades of studies, important functional aspects of SERT such as the oligomerization state of native SERT and its interactions with potential proteins remain unresolved. Here, we develop methods to isolate SERT from porcine brain (pSERT) using a mild, nonionic detergent, utilize fluorescence-detection size-exclusion chromatography to investigate its oligomerization state and interactions with other proteins, and employ single-particle cryo-electron microscopy to elucidate the structures of pSERT in complexes with methamphetamine or cocaine, providing structural insights into psychostimulant recognition and accompanying pSERT conformations. Methamphetamine and cocaine both bind to the central site, stabilizing the transporter in an outward open conformation. We also identify densities attributable to multiple cholesterol or cholesteryl hemisuccinate (CHS) molecules, as well as to a detergent molecule bound to the pSERT allosteric site. Under our conditions of isolation, we find that pSERT is best described as a monomeric entity, isolated without interacting proteins, and is ensconced by multiple cholesterol or CHS molecules.

Emerging Trends of Computational Chemistry and Molecular Modeling in Froth Flotation: A Review

Froth flotation is the most versatile process in mineral beneficiation, extensively used to concentrate a wide range of minerals. This process comprises mixtures of more or less liberated minerals, water, air, and various chemical reagents, involving a series of intermingled multiphase physical and chemical phenomena in the aqueous environment. Today's main challenge facing the froth flotation process is to gain atomic-level insights into the properties of its inherent phenomena governing the process performance. While it is often challenging to determine these phenomena via trial-and-error experimentations, molecular modeling approaches not only elicit a deeper understanding of froth flotation but can also assist experimental studies in saving time and budget. Thanks to the rapid development of computer science and advances in high-performance computing (HPC) infrastructures, theoretical/computational chemistry has now matured enough to successfully and gainfully apply to tackle the challenges of complex systems. In mineral processing, however, advanced applications of computational chemistry are increasingly gaining ground and demonstrating merit in addressing these challenges. Accordingly, this contribution aims to encourage mineral scientists, especially those interested in rational reagent design, to become familiarized with the necessary concepts of molecular modeling and to apply similar strategies when studying and tailoring properties at the molecular level. This review also strives to deliver the state-of-the-art integration and application of molecular modeling in froth flotation studies to assist either active researchers in this field to disclose new directions for future research or newcomers to the field to initiate innovative works.

Differential dynamics and direct interaction of bound ligands with lipids in multidrug transporter ABCG2

ABCG2 is an ATP-binding cassette (ABC) transporter that extrudes a wide range of xenobiotics and drugs from the cell and contributes to multidrug resistance in cancer cells. Following our recent structural characterization of topotecan-bound ABCG2, here, we present cryo-EM structures of ABCG2 under turnover conditions in complex with a special modulator and slow substrate, tariquidar, in nanodiscs. The structures reveal that similar to topotecan, tariquidar induces two distinct ABCG2 conformations under turnover conditions (turnover-1 and turnover-2). μs-scale molecular dynamics simulations of drug-bound and apo ABCG2 in native-like lipid bilayers, in both topotecan- and tariquidar-bound states, characterize the ligand size as a major determinant of its binding stability. The simulations highlight direct lipid-drug interactions for the smaller topotecan, which exhibits a highly dynamic binding mode. In contrast, the larger tariquidar occupies most of the available volume in the binding pocket, thus leaving little space for lipids to enter the cavity and interact with it. Similarly, when simulating ABCG2 in the apo inward-open state, we also observe spontaneous penetration of phospholipids into the binding cavity. The captured phospholipid diffusion pathway into ABCG2 offers a putative general path to recruit any hydrophobic/amphiphilic substrates directly from the membrane. Our simulations also reveal that ABCG2 rejects cholesterol as a substrate, which is omnipresent in plasma membranes that contain ABCG2. At the same time, cholesterol is found to prohibit the penetration of phospholipids into ABCG2. These molecular findings have direct functional ramifications on ABCG2's function as a transporter.

Small Proteins in Archaea, a Mainly Unexplored World

In recent years, increasing numbers of small proteins have moved into the focus of science. Small proteins have been identified and characterized in all three domains of life, but the majority remains functionally uncharacterized, lack secondary structure, and exhibit limited evolutionary conservation. While quite a few have already been described for bacteria and eukaryotic organisms, the amount of known and functionally analyzed archaeal small proteins is still very limited. In this review, we compile the current state of research, show strategies for systematic approaches for global identification of small archaeal proteins, and address selected functionally characterized examples. Besides, we document exemplarily for one archaeon the tool development and optimization to identify small proteins using genome-wide approaches.

Structures of the TMC-1 complex illuminate mechanosensory transduction

The initial step in the sensory transduction pathway underpinning hearing and balance in mammals involves the conversion of force into the gating of a mechanosensory transduction channel¹. Despite the profound socioeconomic impacts of hearing disorders and the fundamental biological significance of understanding mechanosensory transduction, the composition, structure and mechanism of the mechanosensory transduction complex have remained poorly characterized. Here we report the single-particle cryo-electron microscopy structure of the native transmembrane channel-like protein 1 (TMC-1) mechanosensory transduction complex isolated from Caenorhabditis elegans. The two-fold symmetric complex is composed of two copies each of the pore-forming TMC-1 subunit, the calcium-binding protein CALM-1 and the transmembrane inner ear protein TMIE. CALM-1 makes extensive contacts with the cytoplasmic face of the TMC-1 subunits, whereas the single-pass TMIE subunits reside on the periphery of the complex, poised like the handles of an accordion. A subset of complexes additionally includes a single arrestin-like protein, arrestin domain protein (ARRD-6), bound to a CALM-1 subunit. Single-particle reconstructions and molecular dynamics simulations show how the mechanosensory transduction complex deforms the membrane bilayer and suggest crucial roles for lipid-protein interactions in the mechanism by which mechanical force is transduced to ion channel gating.

Targeting cryptic-orthosteric site of PD-L1 for inhibitor identification using structure-guided approach

Approved mAbs that block the protein-protein interaction (PPI) interface of the PD-1/PD-L1 immune checkpoint axis have led to significant improvements in cancer treatment. Despite having drawbacks of mAbs only few a compounds are reported till date against this axis. Inhibiting PPIs using small molecules has emerged as a significant therapeutic opportunity, demanding for the identification of drug-like molecules at an accelerated pace under the hit-to-lead campaigns. Due to the PD-L1's cross-talk with PD-1/CD80 and its overexpression on cancer cells, as well as the availability of its crystal structures with small molecules, it is an enticing therapeutic target for structure-assisted small molecule design. Furthermore, the selection of chemical databases enriched with focused designing for PPI interfaces is crucial. Therefore, in this study we have utilized the Asinex signature library for structure-assisted virtual screening to find the potential PD-L1 inhibitors by targeting the cryptic PD-L1 interface, followed by induced fit docking for pose refinements in the pocket. The obtained hits were then subjected to interaction fingerprinting and ligand-based drug-likeness investigations in order to evaluate and analyze their drug-like qualities (ADME). Twelve compounds qualified for molecular dynamics simulations, followed by thermodynamic calculations for evaluation of their stability and energetics inside the pocket. Two novel compounds with different chemical moieties have been identified that are consistent throughout the simulation, mimicking the interactions and binding energies with BMS-1166. These compounds appear as potential therapeutic candidates to be explored experimentally, thereby paving the way for the development of novel leads as immunomodulators.

Folding and Insertion of Transmembrane Helices at the ER

In eukaryotic cells, the endoplasmic reticulum (ER) is the entry point for newly synthesized proteins that are subsequently distributed to organelles of the endomembrane system. Some of these proteins are completely translocated into the lumen of the ER while others integrate stretches of amino acids into the greasy 30 Å wide interior of the ER membrane bilayer. It is generally accepted that to exist in this non-aqueous environment the majority of membrane integrated amino acids are primarily non-polar/hydrophobic and adopt an α-helical conformation. These stretches are typically around 20 amino acids long and are known as transmembrane (TM) helices. In this review, we will consider how transmembrane helices achieve membrane integration. We will address questions such as: Where do the stretches of amino acids fold into a helical conformation? What is/are the route/routes that these stretches take from synthesis at the ribosome to integration through the ER translocon? How do these stretches ‘know’ to integrate and in which orientation? How do marginally hydrophobic stretches of amino acids integrate and survive as transmembrane helices?

All-atom molecular dynamics study of hepatitis B virus containing pregenome RNA in solution

Immature hepatitis B virus (HBV) captures nucleotides in its capsid for reverse transcription. The nucleotides and nucleotide analog drugs, which are triphosphorylated and negatively charged in the cell, approach the capsid via diffusion and are absorbed into it. In this study, we performed a long-time molecular dynamics calculation of the entire HBV capsid containing pregenome RNA to investigate the interactions between the capsid and negatively charged substances. Electric field analysis demonstrated that negatively charged substances can approach the HBV capsid by thermal motion, avoiding spikes. The substances then migrate all over the floor of the HBV capsid. Finally, they find pores through which they can pass through the HBV capsid shell. Free energy profiles were calculated along these pores for small ions to understand their permeability through the pores. Anions (Cl^-) showed higher free energy barriers than cations (Na⁺ and K⁺) through all pores, and the permeation rate of Cl^- was eight times slower than that of K⁺ or Na⁺. Furthermore, the ions were more stable in the capsid than in the bulk water. Thus, the HBV capsid exerts ion selectivity for uptake and provides an environment for ions, such as nucleotides and nucleotide analog drugs, to be stabilized within the capsid.

Molecular Modeling of Signal Peptide Recognition by Eukaryotic Sec Complexes

Here, we review recent molecular modelling and simulation studies of the Sec translocon, the primary component/channel of protein translocation into the endoplasmic reticulum (ER) and bacterial periplasm, respectively. Our focus is placed on the eukaryotic Sec61, but we also mention modelling studies on prokaryotic SecY since both systems operate in related ways. Cryo-EM structures are now available for different conformational states of the Sec61 complex, ranging from the idle or closed state over an inhibited state with the inhibitor mycolactone bound near the lateral gate, up to a translocating state with bound substrate peptide in the translocation pore. For all these states, computational studies have addressed the conformational dynamics of the translocon with respect to the pore ring, the plug region, and the lateral gate. Also, molecular simulations are addressing mechanistic issues of insertion into the ER membrane vs. translocation into the ER, how signal-peptides are recognised at all in the translocation pore, and how accessory proteins affect the Sec61 conformation in the co- and post-translational pathways.

Molecular mechanism of prestin electromotive signal amplification

Hearing involves two fundamental processes: mechano-electrical transduction and signal amplification. Despite decades of studies, the molecular bases for both remain elusive. Here, we show how prestin, the electromotive molecule of outer hair cells (OHCs) that senses both voltage and membrane tension, mediates signal amplification by coupling conformational changes to alterations in membrane surface area. Cryoelectron microscopy (cryo-EM) structures of human prestin bound with chloride or salicylate at a common "anion site" adopt contracted or expanded states, respectively. Prestin is ensconced within a perimeter of well-ordered lipids, through which it induces dramatic deformation in the membrane and couples protein conformational changes to the bulk membrane. Together with computational studies, we illustrate how the anion site is allosterically coupled to changes in the transmembrane domain cross-sectional area and the surrounding membrane. These studies provide insight into OHC electromotility by providing a structure-based mechanism of the membrane motor prestin.

TRAM1 protein may support ER protein import by modulating the phospholipid bilayer near the lateral gate of the Sec61-channel

In mammalian cells, one-third of all polypeptides is transported into or through the ER-membrane via the Sec61-channel. While the Sec61-complex facilitates the transport of all polypeptides with amino-terminal signal peptides (SP) or SP-equivalent transmembrane helices (TMH), the translocating chain-associated membrane protein (now termed TRAM1) was proposed to support transport of a subset of precursors. To identify possible determinants of TRAM1 substrate specificity, we systematically identified TRAM1-dependent precursors by analyzing cellular protein abundance changes upon TRAM1 depletion in HeLa cells using quantitative label-free proteomics. In contrast to previous analysis after TRAP depletion, SP and TMH analysis of TRAM1 clients did not reveal any distinguishing features that could explain its putative substrate specificity. To further address the TRAM1 mechanism, live-cell calcium imaging was carried out after TRAM1 depletion in HeLa cells. In additional contrast to previous analysis after TRAP depletion, TRAM1 depletion did not affect calcium leakage from the ER. Thus, TRAM1 does not appear to act as SP- or TMH-receptor on the ER-membrane’s cytosolic face and does not appear to affect the open probability of the Sec61-channel. It may rather play a supportive role in protein transport, such as making the phospholipid bilayer conducive for accepting SP and TMH in the vicinity of the lateral gate of the Sec61-channel.

Abbreviations: ER, endoplasmic reticulum; OST, oligosaccharyltransferase; RAMP, ribosome-associated membrane protein; SP, signal peptide; SR, SRP-receptor; SRP, signal recognition particle; TMH, signal peptide-equivalent transmembrane helix; TRAM, translocating chain-associated membrane protein; TRAP, translocon-associated protein.

The structure of the TOM core complex in the mitochondrial outer membrane

In the past three decades, significant advances have been made in providing the biochemical background of TOM (translocase of the outer mitochondrial membrane)-mediated protein translocation into mitochondria. In the light of recent cryoelectron microscopy-derived structures of TOM isolated from Neurospora crassa and Saccharomyces cerevisiae, the interpretation of biochemical and biophysical studies of TOM-mediated protein transport into mitochondria now rests on a solid basis. In this review, we compare the subnanometer structure of N. crassa TOM core complex with that of yeast. Both structures reveal remarkably well-conserved symmetrical dimers of 10 membrane protein subunits. The structural data also validate predictions of weakly stable regions in the transmembrane β-barrel domains of the protein-conducting subunit Tom40, which signal the existence of β-strands located in interfaces of protein-protein interactions.

How Does the Ribosome Fold the Proteome?

Folding of polypeptides begins during their synthesis on ribosomes. This process has evolved as a means for the cell to maintain proteostasis, by mitigating the risk of protein misfolding and aggregation. The capacity to now depict this cellular feat at increasingly higher resolution is providing insight into the mechanistic determinants that promote successful folding. Emerging from these studies is the intimate interplay between protein translation and folding, and within this the ribosome particle is the key player. Its unique structural properties provide a specialized scaffold against which nascent polypeptides can begin to form structure in a highly coordinated, co-translational manner. Here, we examine how, as a macromolecular machine, the ribosome modulates the intrinsic dynamic properties of emerging nascent polypeptide chains and guides them toward their biologically active structures.

How does Sec63 affect the conformation of Sec61 in yeast?

The Sec complex catalyzes the translocation of proteins of the secretory pathway into the endoplasmic reticulum and the integration of membrane proteins into the endoplasmic reticulum membrane. Some substrate peptides require the presence and involvement of accessory proteins such as Sec63. Recently, a structure of the Sec complex from Saccharomyces cerevisiae, consisting of the Sec61 channel and the Sec62, Sec63, Sec71 and Sec72 proteins was determined by cryo-electron microscopy (cryo-EM). Here, we show by co-precipitation that the Sec61 channel subunit Sbh1 is not required for formation of stable Sec63-Sec61 contacts. Molecular dynamics simulations started from the cryo-EM conformation of Sec61 bound to Sec63 and of unbound Sec61 revealed how Sec63 affects the conformation of Sec61 lateral gate, plug, pore region and pore ring diameter via three intermolecular contact regions. Molecular docking of SRP-dependent vs. SRP-independent signal peptide chains into the Sec61 channel showed that the pore regions affected by presence/absence of Sec63 play a crucial role in positioning the signal anchors of SRP-dependent substrates nearby the lateral gate.

Incorporation of sensing modalities into de novo designed fluorescence-activating proteins

Through the efforts of many groups, a wide range of fluorescent protein reporters and sensors based on green fluorescent protein and its relatives have been engineered in recent years. Here we explore the incorporation of sensing modalities into de novo designed fluorescence-activating proteins, called mini-fluorescence-activating proteins (mFAPs), that bind and stabilize the fluorescent cis-planar state of the fluorogenic compound DFHBI. We show through further design that the fluorescence intensity and specificity of mFAPs for different chromophores can be tuned, and the fluorescence made sensitive to pH and Ca²⁺ for real-time fluorescence reporting. Bipartite split mFAPs enable real-time monitoring of protein-protein association and (unlike widely used split GFP reporter systems) are fully reversible, allowing direct readout of association and dissociation events. The relative ease with which sensing modalities can be incorporated and advantages in smaller size and photostability make de novo designed fluorescence-activating proteins attractive candidates for optical sensor engineering.

Interplay among Structural Stability, Plasticity, and Energetics Determined by Conformational Attuning of Flexible Loops in PD-1

The dynamics and plasticity of the PD-1/PD-L1 axis are the bottlenecks for the discovery of small-molecule antagonists to perturb this interaction interface significantly. Understanding the process of this protein-protein interaction (PPI) is of fundamental biological interest in structure-based drug designing. Food and Drug Administration (FDA)-approved anti-PD-1 monoclonal antibodies (mAbs) are the first-in-class with distinct binding modes to access this axis clinically; however, their mechanistic aspects remain elusive. Here, we have unveiled the interactive interfaces with PD-L1 and mAbs to investigate the native plasticity of PD-1 at global (structural and dynamical) and local (residue side-chain orientations) levels. We found that the structural stability and coordinated C_α movements are increased in the presence of PD-1's binding partners. The rigorous analysis of these PPIs using computational biophysical approaches revealed PD-1's intrinsic plasticity, its concerted loops' movement (BC, FG, and CC'), distal side-chain motions, and the thermodynamic landscape, which are perturbed remarkably from its unbound to bound states. Based on intra-/inter-residues' contact networks and energetics, the hot-spots have been identified that were found to be essential to arrest the dynamical motions of PD-1 significantly for the rational design of therapeutic agents by mimicking the mAbs mechanism.

NGPINT: a next-generation protein-protein interaction software

Mapping protein-protein interactions at a proteome scale is critical to understanding how cellular signaling networks respond to stimuli. Since eukaryotic genomes encode thousands of proteins, testing their interactions one-by-one is a challenging prospect. High-throughput yeast-two hybrid (Y2H) assays that employ next-generation sequencing to interrogate complementary DNA (cDNA) libraries represent an alternative approach that optimizes scale, cost and effort. We present NGPINT, a robust and scalable software to identify all putative interactors of a protein using Y2H in batch culture. NGPINT combines diverse tools to align sequence reads to target genomes, reconstruct prey fragments and compute gene enrichment under reporter selection. Central to this pipeline is the identification of fusion reads containing sequences derived from both the Y2H expression plasmid and the cDNA of interest. To reduce false positives, these fusion reads are evaluated as to whether the cDNA fragment forms an in-frame translational fusion with the Y2H transcription factor. NGPINT successfully recognized 95% of interactions in simulated test runs. As proof of concept, NGPINT was tested using published data sets and it recognized all validated interactions. NGPINT can process interaction data from any biosystem with an available genome or transcriptome reference, thus facilitating the discovery of protein-protein interactions in model and non-model organisms.

Conformational changes in the nucleotide-binding domains of P-glycoprotein induced by ATP hydrolysis

P-glycoprotein (Pgp) is a member of the ABC transporter superfamily with high physiological importance. Pgp nucleotide-binding domains (NBDs) drive the transport cycle through ATP binding and hydrolysis. We use molecular dynamics simulations to investigate the ATP hydrolysis-induced conformational changes in NBDs. Five systems, including all possible ATP/ADP combinations in the NBDs and the APO system, are simulated. ATP/ADP exchange induces conformational changes mostly within the conserved signature motif of the NBDs, resulting in relative orientational changes in the NBDs. Nucleotide removal leads to additional orientational changes in the NBDs, allowing their dissociation. Furthermore, we capture putative hydrolysis-competent configurations in which the conserved glutamate in the Walker-B motif acts as a catalytic base capturing a water molecule likely initiating ATP hydrolysis.

NeoR, a near-infrared absorbing rhodopsin

The Rhizoclosmatium globosum genome encodes three rhodopsin-guanylyl cyclases (RGCs), which are predicted to facilitate visual orientation of the fungal zoospores. Here, we show that RGC1 and RGC2 function as light-activated cyclases only upon heterodimerization with RGC3 (NeoR). RGC1/2 utilize conventional green or blue-light-sensitive rhodopsins (λmax = 550 and 480 nm, respectively), with short-lived signaling states, responsible for light-activation of the enzyme. The bistable NeoR is photoswitchable between a near-infrared-sensitive (NIR, λmax = 690 nm) highly fluorescent state (QF = 0.2) and a UV-sensitive non-fluorescent state, thereby modulating the activity by NIR pre-illumination. No other rhodopsin has been reported so far to be functional as a heterooligomer, or as having such a long wavelength absorption or high fluorescence yield. Site-specific mutagenesis and hybrid quantum mechanics/molecular mechanics simulations support the idea that the unusual photochemical properties result from the rigidity of the retinal chromophore and a unique counterion triad composed of two glutamic and one aspartic acids. These findings substantially expand our understanding of the natural potential and limitations of spectral tuning in rhodopsin photoreceptors.

Di-aryl guanidinium derivatives: Towards improved α2-Adrenergic affinity and antagonist activity

Compounds with excellent receptor engagement displaying α₂-AR antagonist activity are useful not only for therapeutic purposes (e.g. antidepressants), but also to help in the crystallization of this particular GPCR. Therefore, based on our broad experience in the topic, we have prepared eighteen di-aryl (phenyl and/or pyridin-2-yl) mono- or di-substituted guanidines and 2-aminoimidazolines. The in vitro α₂-AR binding affinity experiments in human brain tissue showed the advantage of a 2-aminoimidazolinium cation, a di-arylmethylene core, a conformationally locked pyridin-2-yl-guanidine and a di-substituted guanidinium to achieve good α₂-AR engagement. After different in vitro [³⁵S]GTPγS binding experiments in human prefrontal cortex tissue, it was possible to identify that compounds 7a, 7b and 7c were α₂-AR partial agonist, whereas 8h was a potent α₂-AR antagonist. Docking and MD studies with a model of α_2A-AR and two crystal structures suggest that antagonism is achieved by compounds carrying a di-substituted guanidine which substituent occupy a pocket adjacent to TM5 without engaging S200^5.42 or S204^5.46, and a mono-substituted cationic group, which favorably interacts with E94^2.65.

Protein Transport Across the Bacterial Plasma Membrane by the Sec Pathway

More than a third of all bacterial polypeptides, comprising the 'exportome', are transported to extracytoplasmic locations. Most of the exportome is targeted and inserts into ('membranome') or crosses ('secretome') the plasma membrane. The membranome and secretome use distinct targeting signals and factors, and driving forces, but both use the ubiquitous and essential Sec translocase and its SecYEG protein-conducting channel. Membranome export is co-translational and uses highly hydrophobic N-terminal signal anchor sequences recognized by the signal recognition particle on the ribosome, that also targets C-tail anchor sequences. Translating ribosomes drive movement of these polypeptides through the lateral gate of SecY into the inner membrane. On the other hand, secretome export is post-translational and carries two types of targeting signals: cleavable N-terminal signal peptides and multiple short hydrophobic targeting signals in their mature domains. Secretome proteins remain translocation competent due to occupying loosely folded to completely non-folded states during targeting. This is accomplished mainly by the intrinsic properties of mature domains and assisted by signal peptides and/or chaperones. Secretome proteins bind to the dimeric SecA subunit of the translocase. SecA converts from a dimeric preprotein receptor to a monomeric ATPase motor and drives vectorial crossing of chains through SecY aided by the proton motive force. Signal peptides are removed by signal peptidases and translocated chains fold or follow subsequent trafficking.

Improving the Efficiency of Ligand-Binding Protein Design with Molecular Dynamics Simulations

Custom-designed ligand-binding proteins represent a promising class of macromolecules with exciting applications toward the design of new enzymes or the engineering of antibodies and small-molecule recruited proteins for therapeutic interventions. However, several challenges remain in designing a protein sequence such that the binding site organization results in high affinity interaction with a bound ligand. Here, we study the dynamics of explicitly solvated designed proteins through all-atom molecular dynamics (MD) simulations to gain insight into the causes that lead to the low affinity or instability of most of these designs, despite the prediction of their success by the computational design methodology. Simulations ranging from 500 to 1000 ns per replicate were conducted on 37 designed protein variants encompassing two distinct folds and a range of ligand affinities, resulting in more than 180 μs of combined sampling. The simulations provide retrospective insights into the properties affecting ligand affinity that can prove useful in guiding further steps of design optimization. Features indicate that entropic components are particularly important for affinity, which are not easily incorporated in the empirical models often used in design protocols. Additionally, we demonstrate that the application of machine learning approaches built upon the output from the simulations can help discriminate between successful and failed binders, such that MD could act as a screening step in protein design, resulting in a more efficient process.

Large-scale simulations of nucleoprotein complexes: ribosomes, nucleosomes, chromatin, chromosomes and CRISPR

Recent advances in biotechnology such as Hi-C, CRISPR/Cas9 and ribosome display have placed nucleoprotein complexes at center stage. Understanding the structural dynamics of these complexes aids in optimizing protocols and interpreting data for these new technologies. The integration of simulation and experiment has helped advance mechanistic understanding of these systems. Coarse-grained simulations, reduced-description models, and explicit solvent molecular dynamics simulations yield useful complementary perspectives on nucleoprotein complex structural dynamics. When combined with Hi-C, cryo-EM, and single molecule measurements, these simulations integrate disparate forms of experimental data into a coherent mechanism.

Serotonin transporter-ibogaine complexes illuminate mechanisms of inhibition and transport

The serotonin transporter (SERT) regulates neurotransmitter homeostasis through the sodium- and chloride-dependent recycling of serotonin into presynaptic neurons^1-3. Major depression and anxiety disorders are treated using selective serotonin reuptake inhibitors-small molecules that competitively block substrate binding and thereby prolong neurotransmitter action^2,4. The dopamine and noradrenaline transporters, together with SERT, are members of the neurotransmitter sodium symporter (NSS) family. The transport activities of NSSs can be inhibited or modulated by cocaine and amphetamines^2,3, and genetic variants of NSSs are associated with several neuropsychiatric disorders including attention deficit hyperactivity disorder, autism and bipolar disorder^2,5. Studies of bacterial NSS homologues-including LeuT-have shown how their transmembrane helices (TMs) undergo conformational changes during the transport cycle, exposing a central binding site to either side of the membrane^1,6-12. However, the conformational changes associated with transport in NSSs remain unknown. To elucidate structure-based mechanisms for transport in SERT we investigated its complexes with ibogaine, a hallucinogenic natural product with psychoactive and anti-addictive properties^13,14. Notably, ibogaine is a non-competitive inhibitor of transport but displays competitive binding towards selective serotonin reuptake inhibitors^15,16. Here we report cryo-electron microscopy structures of SERT-ibogaine complexes captured in outward-open, occluded and inward-open conformations. Ibogaine binds to the central binding site, and closure of the extracellular gate largely involves movements of TMs 1b and 6a. Opening of the intracellular gate involves a hinge-like movement of TM1a and the partial unwinding of TM5, which together create a permeation pathway that enables substrate and ion diffusion to the cytoplasm. These structures define the structural rearrangements that occur from the outward-open to inward-open conformations, and provide insight into the mechanism of neurotransmitter transport and ibogaine inhibition.

Acceleration of cryo-EM Flexible Fitting for Large Biomolecular Systems by Efficient Space Partitioning

Flexible fitting is a powerful technique to build the 3D structures of biomolecules from cryoelectron microscopy (cryo-EM) density maps. One popular method is a cross-correlation coefficient-based approach, where the molecular dynamics (MD) simulation is carried out with the biasing potential that includes the cross-correlation coefficient between the experimental and simulated density maps. Here, we propose efficient parallelization schemes for the calculation of the cross-correlation coefficient to accelerate flexible fitting. Our schemes are tested for small, medium, and large biomolecules using CPU and hybrid CPU + GPU architectures. The scheme for the atomic decomposition MD is suitable for small proteins such as Ca²⁺-ATPase with the all-atom Go model, while that for the domain decomposition MD is better for larger systems such as ribosome with the all-atom Go or the all-atom explicit solvent models. Our methods allow flexible fitting for various biomolecules with reasonable computational cost. This approach also connects high-resolution structure refinements with investigation of protein structure-function relationship.

Mechanism by which water and protein electrostatic interactions control proton transfer at the active site of channelrhodopsin

Channelrhodopsins are light-sensitive ion channels whose reaction cycles involve conformation-coupled transfer of protons. Understanding how channelrhodopsins work is important for applications in optogenetics, where light activation of these proteins triggers changes in the transmembrane potential across excitable membranes. A fundamental open question is how the protein environment ensures that unproductive proton transfer from the retinal Schiff base to the nearby carboxylate counterion is avoided in the resting state of the channel. To address this question, we performed combined quantum mechanical/molecular mechanical proton transfer calculations with explicit treatment of the surrounding lipid membrane. The free energy profiles computed for proton transfer to the counterion, either via a direct jump or mediated by a water molecule, demonstrate that, when retinal is all-trans, water and protein electrostatic interactions largely favour the protonated retinal Schiff base state. We identified a conserved lysine group as an essential structural element for the proton transfer energetics in channelrhodopsins.

An Update on Sec61 Channel Functions, Mechanisms, and Related Diseases

The membrane of the endoplasmic reticulum (ER) of nucleated human cells harbors the protein translocon, which facilitates membrane integration or translocation of almost every newly synthesized polypeptide targeted to organelles of the endo- and exocytotic pathway. The translocon comprises the polypeptide-conducting Sec61 channel and several additional proteins and complexes that are permanently or transiently associated with the heterotrimeric Sec61 complex. This ensemble of proteins facilitates ER targeting of precursor polypeptides, modification of precursor polypeptides in transit through the Sec61 complex, and Sec61 channel gating, i.e., dynamic regulation of the pore forming subunit to mediate precursor transport and calcium efflux. Recently, cryoelectron tomography of translocons in native ER membrane vesicles, derived from human cell lines or patient fibroblasts, and even intact cells has given unprecedented insights into the architecture and dynamics of the native translocon and the Sec61 channel. These structural data are discussed in light of different Sec61 channel activities including ribosome receptor function, membrane insertion, and translocation of newly synthesized polypeptides as well as the putative physiological roles of the Sec61 channel as a passive ER calcium leak channel. Furthermore, the structural insights into the Sec61 channel are incorporated into an overview and update on Sec61 channel-related diseases—the Sec61 channelopathies—and novel therapeutic concepts for their treatment.

The Sec System: Protein Export in Escherichia coli

In Escherichia coli, proteins found in the periplasm or the outer membrane are exported from the cytoplasm by the general secretory, Sec, system before they acquire stably folded structure. This dynamic process involves intricate interactions among cytoplasmic and membrane proteins, both peripheral and integral, as well as lipids. In vivo, both ATP hydrolysis and proton motive force are required. Here, we review the Sec system from the inception of the field through early 2016, including biochemical, genetic, and structural data.

Accessing gap-junction channel structure-function relationships through molecular modeling and simulations

Background

Gap junction channels (GJCs) are massive protein channels connecting the cytoplasm of adjacent cells. These channels allow intercellular transfer of molecules up to ~1 kDa, including water, ions and other metabolites. Unveiling structure-function relationships coded into the molecular architecture of these channels is necessary to gain insight on their vast biological function including electrical synapse, inflammation, development and tissular homeostasis. From early works, computational methods have been critical to analyze and interpret experimental observations. Upon the availability of crystallographic structures, molecular modeling and simulations have become a valuable tool to assess structure-function relationships in GJCs. Modeling different connexin isoforms, simulating the transport process, and exploring molecular variants, have provided new hypotheses and out-of-the-box approaches to the study of these important channels.

Methods

Here, we review foundational structural studies and recent developments on GJCs using molecular modeling and simulation techniques, highlighting the methods and the cross-talk with experimental evidence.

Results and discussion

By comparing results obtained by molecular modeling and simulations techniques with structural and functional information obtained from both recent literature and structural databases, we provide a critical assesment of structure-function relationships that can be obtained from the junction between theoretical and experimental evidence.

Use of restrained molecular dynamics to predict the conformations of phosphorylated receiver domains in two-component signaling systems

Two‐component signaling (TCS) is the primary means by which bacteria, as well as certain plants and fungi, respond to external stimuli. Signal transduction involves stimulus‐dependent autophosphorylation of a sensor histidine kinase and phosphoryl transfer to the receiver domain of a downstream response regulator. Phosphorylation acts as an allosteric switch, inducing structural and functional changes in the pathway's components. Due to their transient nature, phosphorylated receiver domains are challenging to characterize structurally. In this work, we provide a methodology for simulating receiver domain phosphorylation to predict conformations that are nearly identical to experimental structures. Using restrained molecular dynamics, phosphorylated conformations of receiver domains can be reliably sampled on nanosecond timescales. These simulations also provide data on conformational dynamics that can be used to identify regions of functional significance related to phosphorylation. We first validated this approach on several well‐characterized receiver domains and then used it to compare the upstream and downstream components of the fungal Sln1 phosphorelay. Our results demonstrate that this technique provides structural insight, obtained in the absence of crystallographic or NMR information, regarding phosphorylation‐induced conformational changes in receiver domains that regulate the output of their associated signaling pathway. To our knowledge, this is the first time such a protocol has been described that can be broadly applied to TCS proteins for predictive purposes. Proteins 2016; 85:155–176. © 2016 Wiley Periodicals, Inc.

Electron tunneling in proteins program

We developed a unique integrated software package (called Electron Tunneling in Proteins Program or ETP) which provides an environment with different capabilities such as tunneling current calculation, semi-empirical quantum mechanical calculation, and molecular modeling simulation for calculation and analysis of electron transfer reactions in proteins. ETP program is developed as a cross-platform client-server program in which all the different calculations are conducted at the server side while only the client terminal displays the resulting calculation outputs in the different supported representations. ETP program is integrated with a set of well-known computational software packages including Gaussian, BALLVIEW, Dowser, pKip, and APBS. In addition, ETP program supports various visualization methods for the tunneling calculation results that assist in a more comprehensive understanding of the tunneling process. © 2016 Wiley Periodicals, Inc.

Advances in the molecular dynamics flexible fitting method for cryo-EM modeling

Molecular Dynamics Flexible Fitting (MDFF) is an established technique for fitting all-atom structures of molecules into corresponding cryo-electron microscopy (cryo-EM) densities. The practical application of MDFF is simple but requires a user to be aware of and take measures against a variety of possible challenges presented by each individual case. Some of these challenges arise from the complexity of a molecular structure or the limited quality of available structural models and densities to be interpreted, while others stem from the intricacies of MDFF itself. The current article serves as an overview of the strategies that have been developed since MDFF's inception to overcome common challenges and successfully perform MDFF simulations.

CHARMM Force-Fields with Modified Polyphosphate Parameters Allow Stable Simulation of the ATP-Bound Structure of Ca(2+)-ATPase

Adenosine triphosphate (ATP) is an indispensable energy source in cells. In a wide variety of biological phenomena like glycolysis, muscle contraction/relaxation, and active ion transport, chemical energy released from ATP hydrolysis is converted to mechanical forces to bring about large-scale conformational changes in proteins. Investigation of structure-function relationships in these proteins by molecular dynamics (MD) simulations requires modeling of ATP in solution and ATP bound to proteins with accurate force-field parameters. In this study, we derived new force-field parameters for the triphosphate moiety of ATP based on the high-precision quantum calculations of methyl triphosphate. We tested our new parameters on membrane-embedded sarcoplasmic reticulum Ca(2+)-ATPase and four soluble proteins. The ATP-bound structure of Ca(2+)-ATPase remains stable during MD simulations, contrary to the outcome in shorter simulations using original parameters. Similar results were obtained with the four ATP-bound soluble proteins. The new force-field parameters were also tested by investigating the range of conformations sampled during replica-exchange MD simulations of ATP in explicit water. Modified parameters allowed a much wider range of conformational sampling compared with the bias toward extended forms with original parameters. A diverse range of structures agrees with the broad distribution of ATP conformations in proteins deposited in the Protein Data Bank. These simulations suggest that the modified parameters will be useful in studies of ATP in solution and of the many ATP-utilizing proteins.

Multi-scaled normal mode analysis method for dynamics simulation of protein-membrane complexes: A case study of potassium channel gating motion correlations

Membrane proteins play critically important roles in many cellular activities such as ions and small molecule transportation, signal recognition, and transduction. In order to fulfill their functions, these proteins must be placed in different membrane environments and a variety of protein-lipid interactions may affect the behavior of these proteins. One of the key effects of protein-lipid interactions is their ability to change the dynamics status of membrane proteins, thus adjusting their functions. Here, we present a multi-scaled normal mode analysis (mNMA) method to study the dynamics perturbation to the membrane proteins imposed by lipid bi-layer membrane fluctuations. In mNMA, channel proteins are simulated at all-atom level while the membrane is described with a coarse-grained model. mNMA calculations clearly show that channel gating motion can tightly couple with a variety of membrane deformations, including bending and twisting. We then examined bi-channel systems where two channels were separated with different distances. From mNMA calculations, we observed both positive and negative gating correlations between two neighboring channels, and the correlation has a maximum as the channel center-to-center distance is close to 2.5 times of their diameter. This distance is larger than recently found maximum attraction distance between two proteins embedded in membrane which is 1.5 times of the protein size, indicating that membrane fluctuation might impose collective motions among proteins within a larger area. The hybrid resolution feature in mNMA provides atomic dynamics information for key components in the system without costing much computer resource. We expect it to be a conventional simulation tool for ordinary laboratories to study the dynamics of very complicated biological assemblies. The source code is available upon request to the authors.

Variation in the ribosome interacting loop of the Sec61α from Giardia lamblia

The interaction between the ribosome and the endoplasmic reticulum-located Sec61 protein translocon is mediated through an arginine residue of Sec61α, which is conserved in all prokaryotic and eukaryotic orthologues characterized to date. Using in silico approaches we report that instead of arginine, this ribosome-interaction function is most likely discharged by a lysine residue in the protist Giardia lamblia. This functional substitution of the R with a K in GlSec61α may have taken place to accommodate a G-rich rRNA.

Investigation of the conformational dynamics of the apo A2A adenosine receptor

The activation/deactivation processes for G-protein coupled receptors (GPCRs) have been computationally studied for several different classes, including rhodopsin, the β2 adrenergic receptor, and the M2 muscarinic receptor. Despite determined cocrystal structures of the adenosine A2A receptor (A2A AR) in complex with antagonists, agonists and an antibody, the deactivation process of this GPCR is not completely understood. In this study, we investigate the convergence of two apo simulations, one starting with an agonist-bound conformation (PDB: 3QAK)(14) and the other starting with an antagonist-bound conformation (PDB: 3EML)(11) . Despite the two simulations not completely converging, we were able to identify distinct intermediate steps of the deactivation process characterized by the movement of Y288(7.53) in the NPxxY motif. We find that Y288(7.53) contributes to the process by forming hydrogen bonds to residues in transmembrane helices 2 and 7 and losing these interactions upon full deactivation. Y197(5.58) also plays a role in the process by forming a hydrogen bond only once the side chain moves from the lipid interface to the middle of the helical bundle.

GENESIS: a hybrid-parallel and multi-scale molecular dynamics simulator with enhanced sampling algorithms for biomolecular and cellular simulations

GENESIS (Generalized-Ensemble Simulation System) is a new software package for molecular dynamics (MD) simulations of macromolecules. It has two MD simulators, called ATDYN and SPDYN. ATDYN is parallelized based on an atomic decomposition algorithm for the simulations of all-atom force-field models as well as coarse-grained Go-like models. SPDYN is highly parallelized based on a domain decomposition scheme, allowing large-scale MD simulations on supercomputers. Hybrid schemes combining OpenMP and MPI are used in both simulators to target modern multicore computer architectures. Key advantages of GENESIS are (1) the highly parallel performance of SPDYN for very large biological systems consisting of more than one million atoms and (2) the availability of various REMD algorithms (T-REMD, REUS, multi-dimensional REMD for both all-atom and Go-like models under the NVT, NPT, NPAT, and NPγT ensembles). The former is achieved by a combination of the midpoint cell method and the efficient three-dimensional Fast Fourier Transform algorithm, where the domain decomposition space is shared in real-space and reciprocal-space calculations. Other features in SPDYN, such as avoiding concurrent memory access, reducing communication times, and usage of parallel input/output files, also contribute to the performance. We show the REMD simulation results of a mixed (POPC/DMPC) lipid bilayer as a real application using GENESIS. GENESIS is released as free software under the GPLv2 licence and can be easily modified for the development of new algorithms and molecular models. WIREs Comput Mol Sci 2015, 5:310–323. doi: 10.1002/wcms.1220

Lateral opening of the bacterial translocon on ribosome binding and signal peptide insertion

Proteins are co-translationally inserted into the bacterial plasma membrane via the SecYEG translocon by lateral release of hydrophobic transmembrane segments into the phospholipid bilayer. The trigger for lateral opening of the translocon is not known. Here we monitor lateral opening by photo-induced electron transfer (PET) between two fluorophores attached to the two SecY helices at the rim of the gate. In the resting translocon, the fluorescence is quenched, consistent with a closed conformation. Ribosome binding to the translocon diminishes PET quenching, indicating opening of the gate. The effect is larger with ribosomes exposing hydrophobic transmembrane segments and vanishes at low temperature. We propose a temperature-dependent dynamic equilibrium between closed and open conformations of the translocon that is shifted towards partially and fully open by ribosome binding and insertion of a hydrophobic peptide, respectively. The combined effects of ribosome and peptide binding allow for co-translational membrane insertion of successive transmembrane segments.

Integral membrane proteins laterally partition from the SecYEG translocon into the phospholipid bilayer. Here, the authors use photo-induced electron transfer to show that ribosome binding induces the opening of the lateral gate, and demonstrate that lateral opening does not happen at low temperature.

Flexibility in the insulin receptor ectodomain enables docking of insulin in crystallographic conformation observed in a hormone-bound microreceptor

Insulin binding to the insulin receptor (IR) is the first key step in initiating downstream signaling cascades for glucose homeostasis in higher organisms. The molecular details of insulin recognition by IR are not yet completely understood, but a picture of hormone/receptor interactions at one of the epitopes (Site 1) is beginning to emerge from recent structural evidence. However, insulin-bound structures of truncated IR suggest that crystallographic conformation of insulin cannot be accommodated in the full IR ectodomain due to steric overlap of insulin with the first two type III fibronectin domains (F1 and F2), which are contributed to the insulin binding-pocket by the second subunit in the IR homodimer. A conformational change in the F1-F2 pair has thus been suggested. In this work, we present an all-atom structural model of complex of insulin and the IR ectodomain, where no structural overlap of insulin with the receptor domains (F1 and F2) is observed. This structural model was arrived at by flexibly fitting parts of our earlier insulin/IR all-atom model into the simulated density maps of crystallized constructs combined with conformational sampling from apo-IR solution conformations. Importantly, our experimentally-consistent model helps rationalize yet unresolved Site.

xMDFF: molecular dynamics flexible fitting of low-resolution X-ray structures

X-ray crystallography remains the most dominant method for solving atomic structures. However, for relatively large systems, the availability of only medium-to-low-resolution diffraction data often limits the determination of all-atom details. A new molecular dynamics flexible fitting (MDFF)-based approach, xMDFF, for determining structures from such low-resolution crystallographic data is reported. xMDFF employs a real-space refinement scheme that flexibly fits atomic models into an iteratively updating electron-density map. It addresses significant large-scale deformations of the initial model to fit the low-resolution density, as tested with synthetic low-resolution maps of D-ribose-binding protein. xMDFF has been successfully applied to re-refine six low-resolution protein structures of varying sizes that had already been submitted to the Protein Data Bank. Finally, via systematic refinement of a series of data from 3.6 to 7 Å resolution, xMDFF refinements together with electrophysiology experiments were used to validate the first all-atom structure of the voltage-sensing protein Ci-VSP.

GPU-accelerated analysis and visualization of large structures solved by molecular dynamics flexible fitting

Hybrid structure fitting methods combine data from cryo-electron microscopy and X-ray crystallography with molecular dynamics simulations for the determination of all-atom structures of large biomolecular complexes. Evaluating the quality-of-fit obtained from hybrid fitting is computationally demanding, particularly in the context of a multiplicity of structural conformations that must be evaluated. Existing tools for quality-of-fit analysis and visualization have previously targeted small structures and are too slow to be used interactively for large biomolecular complexes of particular interest today such as viruses or for long molecular dynamics trajectories as they arise in protein folding. We present new data-parallel and GPU-accelerated algorithms for rapid interactive computation of quality-of-fit metrics linking all-atom structures and molecular dynamics trajectories to experimentally-determined density maps obtained from cryo-electron microscopy or X-ray crystallography. We evaluate the performance and accuracy of the new quality-of-fit analysis algorithms vis-à-vis existing tools, examine algorithm performance on GPU-accelerated desktop workstations and supercomputers, and describe new visualization techniques for results of hybrid structure fitting methods.