Cyclic Ion Mobility for Hydrogen/Deuterium Exchange-Mass Spectrometry Applications

Hydrogen/deuterium exchange-mass spectrometry (HDX-MS) has emerged as a powerful tool to probe protein dynamics. As a bottom-up technique, HDX-MS provides information at peptide-level resolution, allowing structural localization of dynamic changes. Consequently, the HDX-MS data quality is largely determined by the number of peptides that are identified and monitored after deuteration. Integration of ion mobility (IM) into HDX-MS workflows has been shown to increase the data quality by providing an orthogonal mode of peptide ion separation in the gas phase. This is of critical importance for challenging targets such as integral membrane proteins (IMPs), which often suffer from low sequence coverage or redundancy in HDX-MS analyses. The increasing complexity of samples being investigated by HDX-MS, such as membrane mimetic reconstituted and in vivo IMPs, has generated need for instrumentation with greater resolving power. Recently, Giles et al. developed cyclic ion mobility (cIM), an IM device with racetrack geometry that enables scalable, multipass IM separations. Using one-pass and multipass cIM routines, we use the recently commercialized SELECT SERIES Cyclic IM spectrometer for HDX-MS analyses of four detergent solubilized IMP samples and report its enhanced performance. Furthermore, we develop a novel processing strategy capable of better handling multipass cIM data. Interestingly, use of one-pass and multipass cIM routines produced unique peptide populations, with their combined peptide output being 31 to 222% higher than previous generation SYNAPT G2-Si instrumentation. Thus, we propose a novel HDX-MS workflow with integrated cIM that has the potential to enable the analysis of more complex systems with greater accuracy and speed.

Dynamic coupling of fast channel gating with slow ATP-turnover underpins protein transport through the Sec translocon

The Sec translocon is a highly conserved membrane assembly for polypeptide transport across, or into, lipid bilayers. In bacteria, secretion through the core channel complex-SecYEG in the inner membrane-is powered by the cytosolic ATPase SecA. Here, we use single-molecule fluorescence to interrogate the conformational state of SecYEG throughout the ATP hydrolysis cycle of SecA. We show that the SecYEG channel fluctuations between open and closed states are much faster (~20-fold during translocation) than ATP turnover, and that the nucleotide status of SecA modulates the rates of opening and closure. The SecY variant PrlA4, which exhibits faster transport but unaffected ATPase rates, increases the dwell time in the open state, facilitating pre-protein diffusion through the pore and thereby enhancing translocation efficiency. Thus, rapid SecYEG channel dynamics are allosterically coupled to SecA via modulation of the energy landscape, and play an integral part in protein transport. Loose coupling of ATP-turnover by SecA to the dynamic properties of SecYEG is compatible with a Brownian-rachet mechanism of translocation, rather than strict nucleotide-dependent interconversion between different static states of a power stroke.

A community effort in SARS-CoV-2 drug discovery

The COVID-19 pandemic continues to pose a substantial threat to human lives and is likely to do so for years to come. Despite the availability of vaccines, searching for efficient small-molecule drugs that are widely available, including in low- and middle-income countries, is an ongoing challenge. In this work, we report the results of an open science community effort, the "Billion molecules against COVID-19 challenge", to identify small-molecule inhibitors against SARS-CoV-2 or relevant human receptors. Participating teams used a wide variety of computational methods to screen a minimum of 1 billion virtual molecules against 6 protein targets. Overall, 31 teams participated, and they suggested a total of 639,024 molecules, which were subsequently ranked to find 'consensus compounds'. The organizing team coordinated with various contract research organizations (CROs) and collaborating institutions to synthesize and test 878 compounds for biological activity against proteases (Nsp5, Nsp3, TMPRSS2), nucleocapsid N, RdRP (only the Nsp12 domain), and (alpha) spike protein S. Overall, 27 compounds with weak inhibition/binding were experimentally identified by binding-, cleavage-, and/or viral suppression assays and are presented here. Open science approaches such as the one presented here contribute to the knowledge base of future drug discovery efforts in finding better SARS-CoV-2 treatments.

Plant size, latitude, and phylogeny explain within-population variability in herbivory

Interactions between plants and herbivores are central in most ecosystems, but their strength is highly variable. The amount of variability within a system is thought to influence most aspects of plant-herbivore biology, from ecological stability to plant defense evolution. Our understanding of what influences variability, however, is limited by sparse data. We collected standardized surveys of herbivory for 503 plant species at 790 sites across 116° of latitude. With these data, we show that within-population variability in herbivory increases with latitude, decreases with plant size, and is phylogenetically structured. Differences in the magnitude of variability are thus central to how plant-herbivore biology varies across macroscale gradients. We argue that increased focus on interaction variability will advance understanding of patterns of life on Earth.

A unifying mechanism for protein transport through the core bacterial Sec machinery

Encapsulation and compartmentalization are fundamental to the evolution of cellular life, but they also pose a challenge: how to partition the molecules that perform biological functions-the proteins-across impermeable barriers into sub-cellular organelles, and to the outside. The solution lies in the evolution of specialized machines, translocons, found in every biological membrane, which act both as gate and gatekeeper across and into membrane bilayers. Understanding how these translocons operate at the molecular level has been a long-standing ambition of cell biology, and one that is approaching its denouement; particularly in the case of the ubiquitous Sec system. In this review, we highlight the fruits of recent game-changing technical innovations in structural biology, biophysics and biochemistry to present a largely complete mechanism for the bacterial version of the core Sec machinery. We discuss the merits of our model over alternative proposals and identify the remaining open questions. The template laid out by the study of the Sec system will be of immense value for probing the many other translocons found in diverse biological membranes, towards the ultimate goal of altering or impeding their functions for pharmaceutical or biotechnological purposes.

A molecular evolution algorithm for ligand design in DOCK

As a complement to virtual screening, de novo design of small molecules is an alternative approach for identifying potential drug candidates. Here, we present a new 3D genetic algorithm to evolve molecules through breeding, mutation, fitness pressure, and selection. The method, termed DOCK_GA, builds upon and leverages powerful sampling, scoring, and searching routines previously implemented into DOCK6. Three primary experiments were used during development: Single-molecule evolution evaluated three selection methods (elitism, tournament, and roulette), in four clinically relevant systems, in terms of mutation type and crossover success, chemical properties, ensemble diversity, and fitness convergence, among others. Large scale benchmarking assessed performance across 651 different protein-ligand systems. Ensemble-based evolution demonstrated using multiple inhibitors simultaneously to seed growth in a SARS-CoV-2 target. Key takeaways include: (1) The algorithm is robust as demonstrated by the successful evolution of molecules across a large diverse dataset. (2) Users have flexibility with regards to parent input, selection method, fitness function, and molecular descriptors. (3) The program is straightforward to run and only requires a single executable and input file at run-time. (4) The elitism selection method yields more tightly clustered molecules in terms of 2D/3D similarity, with more favorable fitness, followed by tournament and roulette.

Towards a molecular mechanism underlying mitochondrial protein import through the TOM and TIM23 complexes

Nearly all mitochondrial proteins need to be targeted for import from the cytosol. For the majority, the first port of call is the translocase of the outer membrane (TOM complex), followed by a procession of alternative molecular machines, conducting transport to their final destination. The pre-sequence translocase of the inner membrane (TIM23-complex) imports proteins with cleavable pre-sequences. Progress in understanding these transport mechanisms has been hampered by the poor sensitivity and time resolution of import assays. However, with the development of an assay based on split NanoLuc luciferase, we can now explore this process in greater detail. Here, we apply this new methodology to understand how ∆ψ and ATP hydrolysis, the two main driving forces for import into the matrix, contribute to the transport of pre-sequence-containing precursors (PCPs) with varying properties. Notably, we found that two major rate-limiting steps define PCP import time: passage of PCP across the outer membrane and initiation of inner membrane transport by the pre-sequence - the rates of which are influenced by PCP size and net charge. The apparent distinction between transport through the two membranes (passage through TOM is substantially complete before PCP-TIM engagement) is in contrast with the current view that import occurs through TOM and TIM in a single continuous step. Our results also indicate that PCPs spend very little time in the TIM23 channel - presumably rapid success or failure of import is critical for maintenance of mitochondrial fitness.

Rate-limiting transport of positively charged arginine residues through the Sec-machinery is integral to the mechanism of protein secretion

Transport of proteins across and into membranes is a fundamental biological process with the vast majority being conducted by the ubiquitous Sec machinery. In bacteria, this is usually achieved when the SecY-complex engages the cytosolic ATPase SecA (secretion) or translating ribosomes (insertion). Great strides have been made towards understanding the mechanism of protein translocation. Yet, important questions remain – notably, the nature of the individual steps that constitute transport, and how the proton-motive force (PMF) across the plasma membrane contributes. Here, we apply a recently developed high-resolution protein transport assay to explore these questions. We find that pre-protein transport is limited primarily by the diffusion of arginine residues across the membrane, particularly in the context of bulky hydrophobic sequences. This specific effect of arginine, caused by its positive charge, is mitigated for lysine which can be deprotonated and transported across the membrane in its neutral form. These observations have interesting implications for the mechanism of protein secretion, suggesting a simple mechanism through which the PMF can aid transport by enabling a 'proton ratchet', wherein re-protonation of exiting lysine residues prevents channel re-entry, biasing transport in the outward direction.

Effect of Postactivation Potentiation After Medium vs. High Inertia Eccentric Overload Exercise on Standing Long Jump, Countermovement Jump, and Change of Direction Performance

ABSTRACT

Beato, M, De Keijzer, KL, Leskauskas, Z, Allen, WJ, Dello Iacono, A, and McErlain-Naylor, SA. Effect of postactivation potentiation after medium vs. high inertia eccentric overload exercise on standing long jump, countermovement jump, and change of direction performance. J Strength Cond Res 35(9): 2616-2621, 2021-This study aimed to evaluate the postactivation potentiation (PAP) effects of an eccentric overload (EOL) exercise on vertical and horizontal jumps and change of direction (COD) performance. Twelve healthy physically active male subjects were involved in a crossover study. The subjects performed 3 sets of 6 repetitions of EOL half squats for maximal power using a flywheel ergometer. Postactivation potentiation using an EOL exercise was compared between a medium (M-EOL) vs. high inertia (H-EOL) experimental condition. Long jump (LJ) was recorded at 30 seconds, 3, and 6 minutes after both EOL exercises and compared with baseline values (control). The same procedure was used to assess countermovement jump (CMJ) height and peak power and 5-m COD test (COD-5m). A fully Bayesian statistical approach to provide probabilistic statements was used in this study. Long jump performance reported improvements after M-EOL and H-EOL exercise (Bayes factor [BF10] = 32.7, strong; BF10 = 9.2, moderate), respectively. Countermovement jump height (BF10 = 135.6, extreme; BF10 > 200, extreme), CMJ peak power (BF10 > 200, extreme; BF10 = 56.1, very strong), and COD-5m (BF10 = 55.7, very strong; BF10 = 16.4, strong) reported improvements after M-EOL and H-EOL exercise, respectively. Between analysis did not report meaningful differences in performance between M-EOL and H-EOL exercises. The present outcomes highlight that PAP using an EOL (M-EOL and H-EOL) improves LJ, CMJ height, CMJ peak power, and COD-5m in male athletes. The optimal time window for the PAP effect was found for both EOL conditions from 3 to 6 minutes. However, M-EOL and H-EOL produce similar PAP effect on LJ, CMJ, and COD-5m tasks.

BiasNet: A Model to Predict Ligand Bias Toward GPCR Signaling

Signaling bias is a feature of many G protein-coupled receptor (GPCR) targeting drugs with potential clinical implications. Whether it is therapeutically advantageous for a drug to be G protein biased or β-arrestin biased depends on the context of the signaling pathway. Here, we explored GPCR ligands that exhibit biased signaling to gain insights into scaffolds and pharmacophores that lead to bias. More specifically, we considered BiasDB, a database containing information about GPCR biased ligands, and focused our analysis on ligands which show either a G protein or β-arrestin bias. Five different machine learning models were trained on these ligands using 15 different sets of features. Molecular fragments which were important for training the models were analyzed. Two of these fragments (number of secondary amines and number of aromatic amines) were more prevalent in β-arrestin biased ligands. After training a random forest model on HierS scaffolds, we found five scaffolds, which demonstrated G protein or β-arrestin bias. We also conducted t-SNE clustering, observing correspondence between unsupervised and supervised machine learning methods. To increase the applicability of our work, we developed a web implementation of our models, which can predict bias based on user-provided SMILES, drug names, or PubChem CID. Our web implementation is available at: drugdiscovery.utep.edu/biasnet.

Role of alkylated residues in the tetrapeptide self-assembly-A molecular dynamics study

Designing peptide sequences that self-assemble into well-defined nanostructures can open a new venue for the development of novel drug carriers and molecular contrast agents. Current approaches are often based on a linear block-design of amphiphilic peptides where a hydrophilic peptide chain is terminated by a hydrophobic tail. Here, a new template for a self-assembling tetrapeptide (YXKX, Y = tyrosine, X = alkylated tyrosine, K = lysine) is proposed with two distinct sides relative to the peptide's backbone: alkylated hydrophobic residues on one side and hydrophilic residues on the other side. Using all-atom molecular dynamics simulations, the self-assembly pathway of the tetrapeptide is analyzed for two different concentrations. At both concentrations, tetrapeptides self-assembled into a nanosphere structure. The alkylated tyrosines initialize the self-assembly process via a strong hydrophobic effect and to reduce exposure to the aqueous solvent, they formed a hydrophobic core. The hydrophilic residues occupied the surface of the self-assembled nanosphere. Ordered arrangement of tetrapeptides within the nanosphere with the backbone hydrogen bonding led to a beta sheet formation. Alkyl chain length constrained the size and shape of the nanosphere. This study provides foundation for further exploration of self-assembling structures that are based on peptides with hydrophobic and hydrophilic moieties located on the opposite sides of a peptide backbone.

Inter-membrane association of the Sec and BAM translocons for bacterial outer-membrane biogenesis

The outer-membrane of Gram-negative bacteria is critical for surface adhesion, pathogenicity, antibiotic resistance and survival. The major constituent - hydrophobic β-barrel Outer-Membrane Proteins (OMPs) - are first secreted across the inner-membrane through the Sec-translocon for delivery to periplasmic chaperones, for example SurA, which prevent aggregation. OMPs are then offloaded to the β-Barrel Assembly Machinery (BAM) in the outer-membrane for insertion and folding. We show the Holo-TransLocon (HTL) - an assembly of the protein-channel core-complex SecYEG, the ancillary sub-complex SecDF, and the membrane 'insertase' YidC - contacts BAM through periplasmic domains of SecDF and YidC, ensuring efficient OMP maturation. Furthermore, the proton-motive force (PMF) across the inner-membrane acts at distinct stages of protein secretion: (1) SecA-driven translocation through SecYEG and (2) communication of conformational changes via SecDF across the periplasm to BAM. The latter presumably drives efficient passage of OMPs. These interactions provide insights of inter-membrane organisation and communication, the importance of which is becoming increasingly apparent.

A molecular dual carriageway

In order to enter a cell, an ammonium ion must first dissociate to form an ammonia molecule and a hydrogen ion (a proton), which then pass through the cell membrane separately and recombine inside.

Brain and lesion segmentation in multiple sclerosis using fully convolutional neural networks: A large-scale study

OBJECTIVE

To investigate the performance of deep learning (DL) based on fully convolutional neural network (FCNN) in segmenting brain tissues in a large cohort of multiple sclerosis (MS) patients.

METHODS

We developed a FCNN model to segment brain tissues, including T2-hyperintense MS lesions. The training, validation, and testing of FCNN were based on ~1000 magnetic resonance imaging (MRI) datasets acquired on relapsing-remitting MS patients, as a part of a phase 3 randomized clinical trial. Multimodal MRI data (dual-echo, FLAIR, and T1-weighted images) served as input to the network. Expert validated segmentation was used as the target for training the FCNN. We cross-validated our results using the leave-one-center-out approach.

RESULTS

We observed a high average (95% confidence limits) Dice similarity coefficient for all the segmented tissues: 0.95 (0.92-0.98) for white matter, 0.96 (0.93-0.98) for gray matter, 0.99 (0.98-0.99) for cerebrospinal fluid, and 0.82 (0.63-1.0) for T2 lesions. High correlations between the DL segmented tissue volumes and ground truth were observed (R² > 0.92 for all tissues). The cross validation showed consistent results across the centers for all tissues.

CONCLUSION

The results from this large-scale study suggest that deep FCNN can automatically segment MS brain tissues, including lesions, with high accuracy.

Post-Translational Protein Transport by the Sec Complex

Although it has been studied for 30 years, the mechanism by which secretory proteins are transported post-translationally into the endoplasmic reticulum (ER) has not yet been fully resolved. Recently published structures (Itskanov and Park, Science 2019;363:84-87; Wu, X. et al. Nature 2019;566:136-139) of the heptameric secretory (Sec) complex which mediates post-translational import into the yeast ER shed new light on the process.

Platform for Automated Real-Time High Performance Analytics on Medical Image Data

Biomedical data are quickly growing in volume and in variety, providing clinicians an opportunity for better clinical decision support. Here, we demonstrate a robust platform that uses software automation and high performance computing (HPC) resources to achieve real-time analytics of clinical data, specifically magnetic resonance imaging (MRI) data. We used the Agave application programming interface to facilitate communication, data transfer, and job control between an MRI scanner and an off-site HPC resource. In this use case, Agave executed the graphical pipeline tool GRAphical Pipeline Environment (GRAPE) to perform automated, real-time, quantitative analysis of MRI scans. Same-session image processing will open the door for adaptive scanning and real-time quality control, potentially accelerating the discovery of pathologies and minimizing patient callbacks. We envision this platform can be adapted to other medical instruments, HPC resources, and analytics tools.

ATP-induced asymmetric pre-protein folding as a driver of protein translocation through the Sec machinery

Transport of proteins across membranes is a fundamental process, achieved in every cell by the 'Sec' translocon. In prokaryotes, SecYEG associates with the motor ATPase SecA to carry out translocation for pre-protein secretion. Previously, we proposed a Brownian ratchet model for transport, whereby the free energy of ATP-turnover favours the directional diffusion of the polypeptide (Allen et al., 2016). Here, we show that ATP enhances this process by modulating secondary structure formation within the translocating protein. A combination of molecular simulation with hydrogendeuterium-exchange mass spectrometry and electron paramagnetic resonance spectroscopy reveal an asymmetry across the membrane: ATP-induced conformational changes in the cytosolic cavity promote unfolded pre-protein structure, while the exterior cavity favours its formation. This ability to exploit structure within a pre-protein is an unexplored area of protein transport, which may apply to other protein transporters, such as those of the endoplasmic reticulum and mitochondria.

A molecular dynamics approach towards evaluating osmotic and thermal stress in the extracellular environment

OBJECTIVE

A molecular dynamics approach to understanding fundamental mechanisms of combined thermal and osmotic stress induced by thermochemical ablation (TCA) is presented.

METHODS

Structural models of fibronectin and fibronectin bound to its integrin receptor provide idealized models for the effects of thermal and osmotic stress in the extracellular matrix. Fibronectin binding to integrin is known to facilitate cell survival. The extracellular environment produced by TCA at the lesion boundary was modelled at 37 °C and 43 °C with added sodium chloride (NaCl) concentrations (0, 40, 80, 160, and 320 mM). Atomistic simulations of solvated proteins were performed using the GROMOS96 force field and TIP3P water model. Computational results were compared with the results of viability studies of human hepatocellular carcinoma (HCC) cell lines HepG2 and Hep3B under matching thermal and osmotic experimental conditions.

RESULTS

Cell viability was inversely correlated with hyperthermal and hyperosmotic stresses. Added NaCl concentrations were correlated with a root mean square fluctuation increase of the fibronectin arginylglycylaspartic acid (RGD) binding domain. Computed interaction coefficients demonstrate preferential hydration of the protein model and are correlated with salt-induced strengthening of hydrophobic interactions. Under the combined hyperthermal and hyperosmotic stress conditions (43 °C and 320 mM added NaCl), the free energy change required for fibronectin binding to integrin was less favorable than that for binding under control conditions (37 °C and 0 mM added NaCl).

CONCLUSION

Results quantify multiple measures of structural changes as a function of temperature increase and addition of NaCl to the solution. Correlations between cell viability and stability measures suggest that protein aggregates, non-functional proteins, and less favorable cell attachment conditions have a role in TCA-induced cell stress.

Customizable de novo design strategies for DOCK: Application to HIVgp41 and other therapeutic targets

De novo design can be used to explore vast areas of chemical space in computational lead discovery. As a complement to virtual screening, from-scratch construction of molecules is not limited to compounds in pre-existing vendor catalogs. Here, we present an iterative fragment growth method, integrated into the program DOCK, in which new molecules are built using rules for allowable connections based on known molecules. The method leverages DOCK's advanced scoring and pruning approaches and users can define very specific criteria in terms of properties or features to customize growth toward a particular region of chemical space. The code was validated using three increasingly difficult classes of calculations: (1) Rebuilding known X-ray ligands taken from 663 complexes using only their component parts (focused libraries), (2) construction of new ligands in 57 drug target sites using a library derived from ∼13M drug-like compounds (generic libraries), and (3) application to a challenging protein-protein interface on the viral drug target HIVgp41. The computational testing confirms that the de novo DOCK routines are robust and working as envisioned, and the compelling results highlight the potential utility for designing new molecules against a wide variety of important protein targets. © 2017 Wiley Periodicals, Inc.

Structure-based identification of inhibitors targeting obstruction of the HIVgp41 N-heptad repeat trimer

The viral protein HIVgp41 is an attractive and validated drug target that proceeds through a sequence of conformational changes crucial for membrane fusion, which facilitates viral entry. Prior work has identified inhibitors that interfere with the formation of a required six-helix bundle, composed of trimeric C-heptad (CHR) and N-heptad (NHR) repeat elements, through blocking association of an outer CHR helix or obstructing formation of the inner NHR trimer itself. In this work, we employed similarity-based scoring to identify and experimentally characterize 113 compounds, related to 2 small-molecule inhibitors recently reported by Allen et al. (Bioorg. Med. Chem Lett.2015, 25 2853-59), proposed to act via the NHR trimer obstruction mechanism. The compounds were first tested in an HIV cell-cell fusion assay with the most promising evaluated in a second, more biologically relevant viral entry assay. Of the candidates, compound #11 emerged as the most promising hit (IC₅₀=37.81µM), as a result of exhibiting activity in both assays with low cytotoxicity, as was similarly seen with the known control peptide inhibitor C34. The compound also showed no inhibition of VSV-G pseudotyped HIV entry compared to a control inhibitor suggesting it was specific for HIVgp41. Molecular dynamics simulations showed the predicted DOCK pose of #11 interacts with HIVgp41 in an energetic fashion (per-residue footprints) similar to the four native NHR residues (IQLT) which candidate inhibitors were intended to mimic.

Channel crossing: how are proteins shipped across the bacterial plasma membrane?

The structure of the first protein-conducting channel was determined more than a decade ago. Today, we are still puzzled by the outstanding problem of protein translocation—the dynamic mechanism underlying the consignment of proteins across and into membranes. This review is an attempt to summarize and understand the energy transducing capabilities of protein-translocating machines, with emphasis on bacterial systems: how polypeptides make headway against the lipid bilayer and how the process is coupled to the free energy associated with ATP hydrolysis and the transmembrane protein motive force. In order to explore how cargo is driven across the membrane, the known structures of the protein-translocation machines are set out against the background of the historic literature, and in the light of experiments conducted in their wake. The paper will focus on the bacterial general secretory (Sec) pathway (SecY-complex), and its eukaryotic counterpart (Sec61-complex), which ferry proteins across the membrane in an unfolded state, as well as the unrelated Tat system that assembles bespoke channels for the export of folded proteins.

Unlocking the Bacterial SecY Translocon

The Sec translocon performs protein secretion and membrane protein insertion at the plasma membrane of bacteria and archaea (SecYEG/β), and the endoplasmic reticular membrane of eukaryotes (Sec61). Despite numerous structures of the complex, the mechanism underlying translocation of pre-proteins, driven by the ATPase SecA in bacteria, remains unresolved. Here we present a series of biochemical and computational analyses exploring the consequences of signal sequence binding to SecYEG. The data demonstrate that a signal sequence-induced movement of transmembrane helix 7 unlocks the translocon and that this conformational change is communicated to the cytoplasmic faces of SecY and SecE, involved in SecA binding. Our findings progress the current understanding of the dynamic action of the translocon during the translocation initiation process. The results suggest that the converging effects of the signal sequence and SecA at the cytoplasmic face of SecYEG are decisive for the intercalation and translocation of pre-protein through the SecY channel.

•

Validation of previously observed signal sequence-induced “unlocking” of SecYEG

•

Conformational changes upon SecYEG unlocking are relayed to SecA binding site

•

Unlocking the translocon perturbs the interaction between SecY and SecE

•

Conformational changes distinct between secretion and membrane protein insertion

DOCK 6: Impact of new features and current docking performance

This manuscript presents the latest algorithmic and methodological developments to the structure-based design program DOCK 6.7 focused on an updated internal energy function, new anchor selection control, enhanced minimization options, a footprint similarity scoring function, a symmetry-corrected root-mean-square deviation algorithm, a database filter, and docking forensic tools. An important strategy during development involved use of three orthogonal metrics for assessment and validation: pose reproduction over a large database of 1043 protein-ligand complexes (SB2012 test set), cross-docking to 24 drug-target protein families, and database enrichment using large active and decoy datasets (Directory of Useful Decoys [DUD]-E test set) for five important proteins including HIV protease and IGF-1R. Relative to earlier versions, a key outcome of the work is a significant increase in pose reproduction success in going from DOCK 4.0.2 (51.4%) → 5.4 (65.2%) → 6.7 (73.3%) as a result of significant decreases in failure arising from both sampling 24.1% → 13.6% → 9.1% and scoring 24.4% → 21.1% → 17.5%. Companion cross-docking and enrichment studies with the new version highlight other strengths and remaining areas for improvement, especially for systems containing metal ions. The source code for DOCK 6.7 is available for download and free for academic users at http://dock.compbio.ucsf.edu/.

Structural Basis for Ceramide Recognition and Hydrolysis by Human Neutral Ceramidase

Neutral ceramidase (nCDase) catalyzes conversion of the apoptosis-associated lipid ceramide to sphingosine, the precursor for the proliferative factor sphingosine-1-phosphate. As an enzyme regulating the balance of ceramide and sphingosine-1-phosphate, nCDase is emerging as a therapeutic target for cancer. Here, we present the 2.6-Å crystal structure of human nCDase in complex with phosphate that reveals a striking, 20-Å deep, hydrophobic active site pocket stabilized by a eukaryotic-specific subdomain not present in bacterial ceramidases. Utilizing flexible ligand docking, we predict a likely binding mode for ceramide that superimposes closely with the crystallographically observed transition state analog phosphate. Our results suggest that nCDase uses a new catalytic strategy for Zn(2+)-dependent amidases, and generates ceramide specificity by sterically excluding sphingolipids with bulky headgroups and specifically recognizing the small hydroxyl head group of ceramide. Together, these data provide a foundation to aid drug development and establish common themes for how proteins recognize the bioactive lipid ceramide.

Small molecule inhibitors of HIVgp41 N-heptad repeat trimer formation

Identification of mechanistically novel anti-HIV fusion inhibitors was accomplished using a computer-aided structure-based design approach with the goal of blocking the formation of the N-heptad repeat (NHR) trimer of the viral protein gp41. A virtual screening strategy that included per-residue interaction patterns (footprints) was employed to identify small molecules compatible with putative binding pockets at the internal interface of the NHR helices at the core native viral six-helix bundle. From a screen of ∼2.8 million compounds using the DOCK program, 120 with favorable energetic and footprint overlap characteristics were purchased and experimentally tested leading to two compounds with favorable cell-cell fusion (IC50) and cytotoxicity profiles. Importantly, both hits were identified on the basis of scores containing footprint overlap terms and would not have been identified using the standard DOCK energy function alone. To our knowledge, these compounds represent the first reported small molecules that inhibit viral entry via the proposed NHR-trimer obstruction mechanism.

Steered molecular dynamics identifies critical residues of the Nodamura virus B2 suppressor of RNAi

Nearly all RNA viruses produce double-stranded RNA (dsRNA) during their replication cycles--an important pathogen-associated molecular pattern recognized by the RNA interference (RNAi) pathway in invertebrates and plants. Nodamura virus (NoV) encodes a suppressor of RNA silencing termed B2, which binds to dsRNA and prevents the initiation of RNAi as well as the loading of silencing complexes. Using the published crystal structure of NoV-B2, we performed a series of molecular dynamics (MD) simulations to determine the relative electrostatic and van der Waals contributions of various residues in binding dsRNA, identifying four novel potential interactors: R56, E48, P68 and R69. Additionally, steered MD was used to simulate the binding affinity of NoV-B2 sequences bearing substitutions at positions F49, R56 or R59 to dsRNA, with F49S and R56L/R59L substitutions found to have a significant negative impact on the ability of NoV-B2 to bind dsRNA. NoV RNA1 variants were tested for self-directed replication in both vertebrate (RNAi⁻) and invertebrate (RNAi⁺) cultured cells. Consistent with a role in dsRNA binding, NoV replication in F49C and F49S variant constructs was affected negatively only in RNAi⁺ cells. Thus, we used a combination of MD simulations and experimental mutagenesis to further characterize residues important for NoV-dsRNA interactions.

Implementation of the Hungarian algorithm to account for ligand symmetry and similarity in structure-based design

False negative docking outcomes for highly symmetric molecules are a barrier to the accurate evaluation of docking programs, scoring functions, and protocols. This work describes an implementation of a symmetry-corrected root-mean-square deviation (RMSD) method into the program DOCK based on the Hungarian algorithm for solving the minimum assignment problem, which dynamically assigns atom correspondence in molecules with symmetry. The algorithm adds only a trivial amount of computation time to the RMSD calculations and is shown to increase the reported overall docking success rate by approximately 5% when tested over 1043 receptor–ligand systems. For some families of protein systems the results are even more dramatic, with success rate increases up to 16.7%. Several additional applications of the method are also presented including as a pairwise similarity metric to compare molecules during de novo design, as a scoring function to rank-order virtual screening results, and for the analysis of trajectories from molecular dynamics simulation. The new method, including source code, is available to registered users of DOCK6 (http://dock.compbio.ucsf.edu).

Strategies for lead discovery: application of footprint similarity targeting HIVgp41

A highly-conserved binding pocket on HIVgp41 is an important target for development of anti-viral inhibitors. Holden et al. (Bioorg. Med. Chem. Lett.2012, 22, 3011) recently reported 7 experimentally-verified leads identified through a computational screen to the gp41 pocket in conjunction with a new DOCK scoring method (termed FPS scoring) developed in our laboratory. The method employs molecular footprints based on per-residue van der Waals interactions, electrostatic interactions, or the sum. In this work, we critically examine the gp41 screening results, prioritized using different scoring methods, in terms of two main criteria: (1) ligand pose properties which include footprint and energy score decompositions, MW, number of rotatable bonds, ligand efficiency, formal charge, and volume overlap, and (2) ligand pose stability which includes footprint stability (changes in footprint overlap) and rmsd stability (changes in geometry). Relative to standard DOCK scoring, pose property analyses demonstrate how FPS scoring can be used to identify ligands that mimic a known reference (derived here from the native gp41 substrate), while pose stability analyses demonstrate how FPS scoring can be used to enrich for compounds with greater overall stability during molecular dynamics (MD) simulations. Compellingly, of the 115 compounds tested experimentally, the 7 active compounds, as a group, more closely mimic the footprints made by the reference and show greater MD stability compared to the inactive group. Extensive studies using 116 protein-ligand complexes as controls reveal that ligands in their crystallographic binding pose also maintain higher FPS scores and smaller rmsds than do accompanying decoys, confirming that native poses are indeed 'stable' under the same conditions and that monitoring FPS variability during compound prioritization is likely to be beneficial. Overall, the results suggest the new scoring method will complement current virtual screening approaches for both the identification (FPS-ranking) and prioritization (FPS-stability) of target-compatible molecules in a quantitative and logical way.

Functional convergence of structurally distinct thioesterases from cyanobacteria and plants involved in phylloquinone biosynthesis

The synthesis of phylloquinone (vitamin K1) in photosynthetic organisms requires a thioesterase that hydrolyzes 1,4-dihydroxy-2-naphthoyl-CoA (DHNA-CoA) to release 1,4-dihydroxy-2-naphthoate (DHNA). Cyanobacteria and plants contain distantly related hotdog-fold thioesterases that catalyze this reaction, although the structural basis of these convergent enzymatic activities is unknown. To investigate this, the crystal structures of hotdog-fold DHNA-CoA thioesterases from the cyanobacterium Synechocystis (Slr0204) and the flowering plant Arabidopsis thaliana (AtDHNAT1) were determined. These enzymes form distinct homotetramers and use different active sites to catalyze hydrolysis of DHNA-CoA, similar to the 4-hydroxybenzoyl-CoA (4-HBA-CoA) thioesterases from Pseudomonas and Arthrobacter. Like the 4-HBA-CoA thioesterases, the DHNA-CoA thioesterases contain either an active-site aspartate (Slr0204) or glutamate (AtDHNAT1) that are predicted to be catalytically important. Computational modeling of the substrate-bound forms of both enzymes indicates the residues that are likely to be involved in substrate binding and catalysis. Both enzymes are selective for DHNA-CoA as a substrate, but this selectivity is achieved using divergent predicted binding strategies. The Slr0204 binding pocket is predominantly hydrophobic and closely conforms to DHNA, while that of AtDHNAT1 is more polar and solvent-exposed. Considered in light of the related 4-HBA-CoA thioesterases, these structures indicate that hotdog-fold thioesterases using either an active-site aspartate or glutamate diverged into distinct clades prior to the evolution of strong substrate specificity in these enzymes.

dNTP-dependent conformational transitions in the fingers subdomain of Klentaq1 DNA polymerase: insights into the role of the "nucleotide-binding" state

BACKGROUND

Conformational selection plays a key role in the polymerase cycle.

RESULTS

Klentaq1 exists in conformational equilibrium between three states (open, closed, and “nucleotide-binding”) whose level of occupancy is determined by the bound substrate.

CONCLUSION

The “nucleotide-binding” state plays a pivotal role in the reaction pathway.

SIGNIFICANCE

Direct evidence is provided for the role of a conformationally distinct “nucleotide-binding” state during dNTP incorporation. DNA polymerases are responsible for the accurate replication of DNA. Kinetic, single-molecule, and x-ray studies show that multiple conformational states are important for DNA polymerase fidelity. Using high precision FRET measurements, we show that Klentaq1 (the Klenow fragment of Thermus aquaticus DNA polymerase 1) is in equilibrium between three structurally distinct states. In the absence of nucleotide, the enzyme is mostly open, whereas in the presence of DNA and a correctly base-pairing dNTP, it re-equilibrates to a closed state. In the presence of a dNTP alone, with DNA and an incorrect dNTP, or in elevated MgCl2 concentrations, an intermediate state termed the "nucleotide-binding" state predominates. Photon distribution and hidden Markov modeling revealed fast dynamic and slow conformational processes occurring between all three states in a complex energy landscape suggesting a mechanism in which dNTP delivery is mediated by the nucleotide-binding state. After nucleotide binding, correct dNTPs are transported to the closed state, whereas incorrect dNTPs are delivered to the open state.

Grid-based molecular footprint comparison method for docking and de novo design: application to HIVgp41

Scoring functions are a critically important component of computer-aided screening methods for the identification of lead compounds during early stages of drug discovery. Here, we present a new multigrid implementation of the footprint similarity (FPS) scoring function that was recently developed in our laboratory which has proven useful for identification of compounds which bind to a protein on a per-residue basis in a way that resembles a known reference. The grid-based FPS method is much faster than its Cartesian-space counterpart, which makes it computationally tractable for on-the-fly docking, virtual screening, or de novo design. In this work, we establish that: (i) relatively few grids can be used to accurately approximate Cartesian space footprint similarity, (ii) the method yields improved success over the standard DOCK energy function for pose identification across a large test set of experimental co-crystal structures, for crossdocking, and for database enrichment, and (iii) grid-based FPS scoring can be used to tailor construction of new molecules to have specific properties, as demonstrated in a series of test cases targeting the viral protein HIVgp41. The method is available in the program DOCK6.

Dissection of pilus tip assembly by the FimD usher monomer

Type 1 pili are representative of a class of bacterial surface structures assembled by the conserved chaperone/usher pathway and used by uropathogenic Escherichia coli to attach to bladder cells during infection. The outer membrane assembly platform-the usher-is critical for the formation of pili, catalysing the polymerisation of pilus subunits and enabling the secretion of the nascent pilus. Despite extensive structural characterisation of the usher, a number of questions about its mechanism remain, notably its oligomerisation state, and how it orchestrates the ordered assembly of pilus subunits. We demonstrate here that the FimD usher is able to catalyse in vitro pilus assembly effectively in its monomeric form. Furthermore, by establishing the kinetics of usher-catalysed reactions between various pilus subunits, we establish a complete kinetic model of tip fibrillum assembly, able to account for the order of subunits in native type 1 pili.

Mobility of the SecA 2-helix-finger is not essential for polypeptide translocation via the SecYEG complex

Polypeptide translocation in bacteria, once underway, requires only one copy each of SecA and SecYEG and does not require the mobility of the SecA 2-helix-finger.

The bacterial ATPase SecA and protein channel complex SecYEG form the core of an essential protein translocation machinery. The nature of the conformational changes induced by each stage of the hydrolytic cycle of ATP and how they are coupled to protein translocation are not well understood. The structure of the SecA–SecYEG complex revealed a 2-helix-finger (2HF) of SecA in an ideal position to contact the substrate protein and push it through the membrane. Surprisingly, immobilization of this finger at the edge of the protein channel had no effect on translocation, whereas its imposition inside the channel blocked transport. This analysis resolves the stoichiometry of the active complex, demonstrating that after the initiation process translocation requires only one copy each of SecA and SecYEG. The results also have important implications on the mechanism of energy transduction and the power stroke driving transport. Evidently, the 2HF is not a highly mobile transducing element of polypeptide translocation.

The role of chaperone-subunit usher domain interactions in the mechanism of bacterial pilus biogenesis revealed by ESI-MS

The PapC usher is a β-barrel outer membrane protein essential for assembly and secretion of P pili that are required for adhesion of pathogenic E. coli, which cause the development of pyelonephritis. Multiple protein subunits form the P pilus, the highly specific assembly of which is coordinated by the usher. Despite a wealth of structural knowledge, how the usher catalyzes subunit polymerization and orchestrates a correct and functional order of subunit assembly remain unclear. Here, the ability of the soluble N-terminal (UsherN), C-terminal (UsherC2), and Plug (UsherP) domains of the usher to bind different chaperone-subunit (PapDPapX) complexes is investigated using noncovalent electrospray ionization mass spectrometry. The results reveal that each usher domain is able to bind all six PapDPapX complexes, consistent with an active role of all three usher domains in pilus biogenesis. Using collision induced dissociation, combined with competition binding experiments and dissection of the adhesin subunit, PapG, into separate pilin and adhesin domains, the results reveal why PapG has a uniquely high affinity for the usher, which is consistent with this subunit always being displayed at the pilus tip. In addition, we show how the different soluble usher domains cooperate to coordinate and control efficient pilus assembly at the usher platform. As well as providing new information about the protein-protein interactions that determine pilus biogenesis, the results highlight the power of noncovalent MS to interrogate biological mechanisms, especially in complex mixtures of species.

Pilus biogenesis at the outer membrane of Gram-negative bacterial pathogens

Pili belong to a broad class of bacterial surface structures that play a key role in infection and pathogenicity. The largest and best characterised pilus biogenesis system--the chaperone-usher pathway--is particularly remarkable in its ability to synthesise and display highly organised structures at the outer membrane without any input from endogenous energy sources. The past few years have heralded exciting new developments in our understanding of the structural biology and mechanism of pilus assembly, which are discussed in this review. Such knowledge will be particularly important in the future, as we approach an era of widespread resistance to common antibiotics and require new targets.

Crystal structure of the FimD usher bound to its cognate FimC-FimH substrate

Type 1 pili are the archetypal representative of a widespread class of adhesive multisubunit fibres in Gram-negative bacteria. During pilus assembly, subunits dock as chaperone-bound complexes to an usher, which catalyzes their polymerization and mediates pilus translocation across the outer membrane. We report the crystal structure of the full-length FimD usher bound to the FimC:FimH chaperone:adhesin complex and that of the unbound form of the FimD translocation domain. The FimD:FimC:FimH structure shows FimH inserted inside the FimD 24-stranded β-barrel translocation channel. FimC:FimH is held in place through interactions with the two C-terminal periplasmic domains of FimD, a binding mode confirmed in solution by electron paramagnetic resonance spectroscopy. To accommodate FimH, the usher plug domain is displaced from the barrel lumen to the periplasm, concomitant with a dramatic conformational change in the β-barrel. The N-terminal domain of FimD is observed in an ideal position to catalyse incorporation of a newly recruited chaperone:subunit complex. The FimD:FimC:FimH structure provides unique insights into the pilus subunit incorporation cycle, and captures the first view of a protein transporter in the act of secreting its cognate substrate.

Modeling the relationship between the p53 C-terminal domain and its binding partners using molecular dynamics

Fifty percent of all cancer cases result from mutations of the TP53 gene, which encodes the tumor suppressor p53, and it is hypothesized that the p53-mediated checkpoint pathway is compromised in most of the remaining cases. The p53 C-terminal domain (CTD) is an important site of p53 regulation but by nature is difficult to study, as it is intrinsically disordered. In this study, we performed molecular dynamics simulations on the p53 CTD and five known regulatory binding partners. We identified distinct trends in fluctuation within and around the p53 CTD binding site on each partner demonstrating a behavior that facilitates association. Further, we present evidence that the size of the hydrophobic pocket in each p53 CTD binding site governs the secondary structure of the p53 CTD when in the bound state. This information will be useful for predicting new binding partners for the p53 CTD, identifying interacting regions within other known partners, and discovering inhibitors that provide additional points of control over p53 activity.

Practical considerations for building GROMOS-compatible small-molecule topologies

Molecular dynamics simulations are being applied to increasingly complex systems, including those involving small endogenous compounds and drug molecules. In order to obtain meaningful and accurate data from these simulations, high-quality topologies for small molecules must be generated in a manner that is consistent with the derivation of the force field applied to the system. Often, force fields are designed for use with macromolecules such as proteins, making their transferability to other species challenging. Investigators are increasingly attracted to automated topology generation programs, although the quality of the resulting topologies remains unknown. Here we assess the applicability of the popular PRODRG server that generates small-molecule topologies for use with the GROMOS family of force fields. We find that PRODRG does not reproduce topologies for even the most well-characterized species in the force field due to inconsistent charges and charge groups. We assessed the effects of PRODRG-derived charges on several systems: pure liquids, amino acids at a hydrophobic-hydrophilic interface, and an enzyme-cofactor complex. We found that partial atomic charges generated by PRODRG are largely incompatible with GROMOS force fields, and the behavior of these systems deviates substantially from that of simulations using GROMOS parameters. We conclude by proposing several points as "best practices" for parametrization of small molecules under the GROMOS force fields.

Structural biology of periplasmic chaperones

Proteins often require specific helper proteins, chaperones, to assist with their correct folding and to protect them from denaturation and aggregation. The cell envelope of Gram-negative bacteria provides a particularly challenging environment for chaperones to function in as it lacks readily available energy sources such as adenosine 5' triphosphate (ATP) to power reaction cycles. Periplasmic chaperones have therefore evolved specialized mechanisms to carry out their functions without the input of external energy and in many cases to transduce energy provided by protein folding or ATP hydrolysis at the inner membrane. Structural and biochemical studies have in recent years begun to elucidate the specific functions of many important periplasmic chaperones and how these functions are carried out. This includes not only specific carrier chaperones, such as those involved in the biosynthesis of adhesive fimbriae in pathogenic bacteria, but also more general pathways including the periplasmic transport of outer membrane proteins and the extracytoplasmic stress responses. This chapter aims to provide an overview of protein chaperones so far identified in the periplasm and how structural biology has assisted with the elucidation of their functions.

An intramolecular FRET system monitors fingers subdomain opening in Klentaq1

A major goal of polymerase research is to determine the mechanism through which a nucleotide complementary to a templating DNA base is selected and delivered to the polymerase active site. Structural evidence suggests a large open-to-closed conformational change affecting the fingers subdomain as being crucial to the process. We previously designed a FRET system capable of measuring the rate of fingers subdomain closure in the presence of correct nucleotide. However, this FRET system was limited in that it could not directly measure the rate of fingers subdomain opening by FRET after polymerization or in the absence of DNA. Here we report the development of a new system capable of measuring both fingers subdomain closure and reopening by FRET, and show that the rate of fingers subdomain opening is limited only by the rate of polymerization. We anticipate that this system will scale down to the single molecule level, allowing measurement of fingers subdomain movements in the presence of incorrect nucleotide and in the absence of DNA.

Novel suture device for beating-heart mitral leaflet approximation

BACKGROUND

This investigation evaluates the potential of using a novel suturing device to achieve mitral valve repair (Alfieri type) on a beating heart without cardiopulmonary bypass.

METHODS

Eight healthy adult sheep were anesthetized and the chest was opened via a left thoracotomy. The suture device was directly inserted into the appendage of the left atrium. Suction ports on the distal tip of the device grasped and approximated the mitral leaflets while the heart was beating. Two-dimensional echocardiography and intracardiac pressure monitoring at the tip of the device were utilized to guide the procedure. The device was used to place two single sutures across the two leaflets at the center of the mitral valve. A knot pusher with integrated cutter was used to tie the sutures and cut the suture ends.

RESULTS

In all animals, the free margins of the mitral leaflets were successfully grasped and approximated by this device. Echocardiography confirmed successful deployment of the sutures in all cases, with a figure-of-eight appearance of the valve and normal valve hemodynamic function after placement of the sutures. Mid-leaflet approximation was verified at autopsy immediately after the procedure. No tissue damage was observed.

CONCLUSIONS

This study demonstrates that mitral valve repair (Alfieri type) can be performed safely and consistently on a beating heart without cardiopulmonary bypass using this new tissue approximation suture device. This technique may be applicable to the treatment of ischemic mitral regurgitation in conjunction with revascularization procedures or to mitral regurgitation in heart failure patients.

An Integrated Approach for Maximising Local and Scientific Knowledge for Land Management Decision-Making in the New Zealand High Country

This paper describes the development of a process to facilitate the identification and introduction of sustainable land management practices in the high country of New Zealand. The process was designed to gather and structure community knowledge (both local and scientific) into a single, accessible decision support system (DSS). The development and provision of appropriate, and user-friendly monitoring tools is supported. An outline is given of how this integrated system can be used to integrate monitoring with adaptive management. Special reference is made to how this process is used as a large-scale ecological 'experiment', to enhance continually the knowledge base available for land use decision-making in the South Island high country of New Zealand.

On the action of adrenaline on the blood vessels in human skeletal muscle

Images