Art and Science of the Cellular Mesoscale.

Experimental information from microscopy, structural biology, and bioinformatics may be integrated to build structural models of entire cells with molecular detail. This integrative modeling is challenging in several ways: the intrinsic complexity of biology results in models with many closely packed and heterogeneous components; the wealth of available experimental data is scattered among multiple resources and must be gathered, reconciled, and curated; and computational infrastructure is only now gaining the capability of modeling and visualizing systems of this complexity. We present recent efforts to address these challenges, both with artistic approaches to depicting the cellular mesoscale, and development and application of methods to build quantitative models.

Icosahedral virus structures and the protein data bank.

The structural study of icosahedral viruses has a long and impactful history in both crystallographic methodology and molecular biology. The evolution of the Protein Data Bank has paralleled and supported these studies providing readily accessible formats dealing with novel features associated with viral particle symmetries and subunit interactions. This overview describes the growth in size and complexity of icosahedral viruses from the first early studies of small RNA plant viruses and human picornaviruses up to the larger and more complex bacterial phage, insect, and human disease viruses such as Zika, hepatitis B, Adeno and Polyoma virus. The analysis of icosahedral viral capsid protein domain folds has shown striking similarities, with the beta jelly roll motif observed across multiple evolutionarily divergent species. The icosahedral symmetry of viruses drove the development of noncrystallographic symmetry averaging as a powerful phasing method, and the constraints of maintaining this symmetry resulted in the concept of quasi-equivalence in viral structures. Symmetry also played an important early role in demonstrating the power of cryo-electron microscopy as an alternative to crystallography in generating atomic resolution structures of these viruses. The Protein Data Bank has been a critical resource for assembling and disseminating these structures to a wide community, and the virus particle explorer (VIPER) was developed to enable users to easily generate and view complete viral capsid structures from their asymmetric building blocks. Finally, we share a personal perspective on the early use of computer graphics to communicate the intricacies, interactions, and beauty of these virus structures.

Structure-based virtual screening workflow to identify antivirals targeting HIV-1 capsid.

We have identified novel HIV-1 capsid inhibitors targeting the PF74 binding site. Acting as the building block of the HIV-1 capsid core, the HIV-1 capsid protein plays an important role in the viral life cycle and is an attractive target for antiviral development. A structure-based virtual screening workflow for hit identification was employed, which includes docking 1.6 million commercially-available drug-like compounds from the ZINC database to the capsid dimer, followed by applying two absolute binding free energy (ABFE) filters on the 500 top-ranked molecules from docking. The first employs the Binding Energy Distribution Analysis Method (BEDAM) in implicit solvent. The top-ranked compounds are then refined using the Double Decoupling method in explicit solvent. Both docking and BEDAM refinement were carried out on the IBM World Community Grid as part of the FightAIDS@Home project. Using this virtual screening workflow, we identified 24 molecules with calculated binding free energies between - 6 and - 12 kcal/mol. We performed thermal shift assays on these molecules to examine their potential effects on the stability of HIV-1 capsid hexamer and found that two compounds, ZINC520357473 and ZINC4119064 increased the melting point of the latter by 14.8 °C and 33 °C, respectively. These results support the conclusion that the two ZINC compounds are primary hits targeting the capsid dimer interface. Our simulations also suggest that the two hit molecules may bind at the capsid dimer interface by occupying a new sub-pocket that has not been exploited by existing CA inhibitors. The possible causes for why other top-scored compounds suggested by ABFE filters failed to show measurable activity are discussed.

Proteome-wide covalent ligand discovery in native biological systems.

Small molecules are powerful tools for investigating protein function and can serve as leads for new therapeutics. Most human proteins, however, lack small-molecule ligands, and entire protein classes are considered 'undruggable'. Fragment-based ligand discovery can identify small-molecule probes for proteins that have proven difficult to target using high-throughput screening of complex compound libraries. Although reversibly binding ligands are commonly pursued, covalent fragments provide an alternative route to small-molecule probes, including those that can access regions of proteins that are difficult to target through binding affinity alone. Here we report a quantitative analysis of cysteine-reactive small-molecule fragments screened against thousands of proteins in human proteomes and cells. Covalent ligands were identified for >700 cysteines found in both druggable proteins and proteins deficient in chemical probes, including transcription factors, adaptor/scaffolding proteins, and uncharacterized proteins. Among the atypical ligand-protein interactions discovered were compounds that react preferentially with pro- (inactive) caspases. We used these ligands to distinguish extrinsic apoptosis pathways in human cell lines versus primary human T cells, showing that the former is largely mediated by caspase-8 while the latter depends on both caspase-8 and -10. Fragment-based covalent ligand discovery provides a greatly expanded portrait of the ligandable proteome and furnishes compounds that can illuminate protein functions in native biological systems.

Integrative modeling of the HIV-1 ribonucleoprotein complex.

A coarse-grain computational method integrates biophysical and structural data to generate models of HIV-1 genomic RNA, nucleocapsid and integrase condensed into a mature ribonucleoprotein complex. Several hypotheses for the initial structure of the genomic RNA and oligomeric state of integrase are tested. In these models, integrase interaction captures features of the relative distribution of gRNA in the immature virion and increases the size of the RNP globule, and exclusion of nucleocapsid from regions with RNA secondary structure drives an asymmetric placement of the dimerized 5'UTR at the surface of the RNP globule.

Massive-Scale Binding Free Energy Simulations of HIV Integrase Complexes Using Asynchronous Replica Exchange Framework Implemented on the IBM WCG Distributed Network.

To perform massive-scale replica exchange molecular dynamics (REMD) simulations for calculating binding free energies of protein-ligand complexes, we implemented the asynchronous replica exchange (AsyncRE) framework of the binding energy distribution analysis method (BEDAM) in implicit solvent on the IBM World Community Grid (WCG) and optimized the simulation parameters to reduce the overhead and improve the prediction power of the WCG AsyncRE simulations. We also performed the first massive-scale binding free energy calculations using the WCG distributed computing grid and 301 ligands from the SAMPL4 challenge for large-scale binding free energy predictions of HIV-1 integrase complexes. In total there are ∼10000 simulated complexes, ∼1 million replicas, and ∼2000 µs of aggregated MD simulations. Running AsyncRE MD simulations on the WCG requires accepting a trade-off between the number of replicas that can be run (breadth) and the number of full RE cycles that can be completed per replica (depth). As compared with synchronous Replica Exchange (SyncRE) running on tightly coupled clusters like XSEDE, on the WCG many more replicas can be launched simultaneously on heterogeneous distributed hardware, but each full RE cycle requires more overhead. We compared the WCG results with that from AutoDock and more advanced RE simulations including the use of flattening potentials to accelerate sampling of selected degrees of freedom of ligands and/or receptors related to slow dynamics due to high energy barriers. We propose a suitable strategy of RE simulations to refine high throughput docking results which can be matched to corresponding computing resources: from HPC clusters, to small or medium-size distributed campus grids, and finally to massive-scale computing networks including millions of CPUs like the resources available on the WCG.

cellPACK: a virtual mesoscope to model and visualize structural systems biology.

cellPACK assembles computational models of the biological mesoscale, an intermediate scale (10-100 nm) between molecular and cellular biology scales. cellPACK's modular architecture unites existing and novel packing algorithms to generate, visualize and analyze comprehensive three-dimensional models of complex biological environments that integrate data from multiple experimental systems biology and structural biology sources. cellPACK is available as open-source code, with tools for validation of models and with 'recipes' and models for five biological systems: blood plasma, cytoplasm, synaptic vesicles, HIV and a mycoplasma cell. We have applied cellPACK to model distributions of HIV envelope protein to test several hypotheses for consistency with experimental observations. Biologists, educators and outreach specialists can interact with cellPACK models, develop new recipes and perform packing experiments through scripting and graphical user interfaces at http://cellPACK.org/.

Dense Array of Spikes on HIV-1 Virion Particles.

HIV-1 is rare among viruses for having a low number of envelope glycoprotein (Env) spikes per virion, i.e., ∼7 to 14. This exceptional feature has been associated with avoidance of humoral immunity, i.e., B cell activation and antibody neutralization. Virus-like particles (VLPs) with increased density of Env are being pursued for vaccine development; however, these typically require protein engineering that alters Env structure. Here, we used instead a strategy that targets the producer cell. We employed fluorescence-activated cell sorting (FACS) to sort for cells that are recognized by trimer cross-reactive broadly neutralizing antibody (bnAb) and not by nonneutralizing antibodies. Following multiple iterations of FACS, cells and progeny virions were shown to display higher levels of antigenically correct Env in a manner that correlated between cells and cognate virions (P = 0.027). High-Env VLPs, or hVLPs, were shown to be monodisperse and to display more than a 10-fold increase in spikes per particle by electron microscopy (average, 127 spikes; range, 90 to 214 spikes). Sequencing revealed a partial truncation in the C-terminal tail of Env that had emerged in the sort; however, iterative rounds of "cell factory" selection were required for the high-Env phenotype. hVLPs showed greater infectivity than standard pseudovirions but largely similar neutralization sensitivity. Importantly, hVLPs also showed superior activation of Env-specific B cells. Hence, high-Env HIV-1 virions, obtained through selection of producer cells, represent an adaptable platform for vaccine design and should aid in the study of native Env.IMPORTANCE The paucity of spikes on HIV is a unique feature that has been associated with evasion of the immune system, while increasing spike density has been a goal of vaccine design. Increasing the density of Env by modifying it in various ways has met with limited success. Here, we focused instead on the producer cell. Cells that stably express HIV spikes were screened on the basis of high binding by bnAbs and low binding by nonneutralizing antibodies. Levels of spikes on cells correlated well with those on progeny virions. Importantly, high-Env virus-like particles (hVLPs) were produced with a manifest array of well-defined spikes, and these were shown to be superior in activating desirable B cells. Our study describes HIV particles that are densely coated with functional spikes, which should facilitate the study of HIV spikes and their development as immunogens.

A New Class of Allosteric HIV-1 Integrase Inhibitors Identified by Crystallographic Fragment Screening of the Catalytic Core Domain.

HIV-1 integrase (IN) is essential for virus replication and represents an important multifunctional therapeutic target. Recently discovered quinoline-based allosteric IN inhibitors (ALLINIs) potently impair HIV-1 replication and are currently in clinical trials. ALLINIs exhibit a multimodal mechanism of action by inducing aberrant IN multimerization during virion morphogenesis and by competing with IN for binding to its cognate cellular cofactor LEDGF/p75 during early steps of HIV-1 infection. However, quinoline-based ALLINIs impose a low genetic barrier for the evolution of resistant phenotypes, which highlights a need for discovery of second-generation inhibitors. Using crystallographic screening of a library of 971 fragments against the HIV-1 IN catalytic core domain (CCD) followed by a fragment expansion approach, we have identified thiophenecarboxylic acid derivatives that bind at the CCD-CCD dimer interface at the principal lens epithelium-derived growth factor (LEDGF)/p75 binding pocket. The most active derivative (5) inhibited LEDGF/p75-dependent HIV-1 IN activity in vitro with an IC50 of 72 µm and impaired HIV-1 infection of T cells at an EC50 of 36 µm The identified lead compound, with a relatively small molecular weight (221 Da), provides an optimal building block for developing a new class of inhibitors. Furthermore, although structurally distinct thiophenecarboxylic acid derivatives target a similar pocket at the IN dimer interface as the quinoline-based ALLINIs, the lead compound, 5, inhibited IN mutants that confer resistance to quinoline-based compounds. Collectively, our findings provide a plausible path for structure-based development of second-generation ALLINIs.

Computational challenges of structure-based approaches applied to HIV.

Here, we review some of the opportunities and challenges that we face in computational modeling of HIV therapeutic targets and structural biology, both in terms of methodology development and structure-based drug design (SBDD). Computational methods have provided fundamental support to HIV research since the initial structural studies, helping to unravel details of HIV biology. Computational models have proved to be a powerful tool to analyze and understand the impact of mutations and to overcome their structural and functional influence in drug resistance. With the availability of structural data, in silico experiments have been instrumental in exploiting and improving interactions between drugs and viral targets, such as HIV protease, reverse transcriptase, and integrase. Issues such as viral target dynamics and mutational variability, as well as the role of water and estimates of binding free energy in characterizing ligand interactions, are areas of active computational research. Ever-increasing computational resources and theoretical and algorithmic advances have played a significant role in progress to date, and we envision a continually expanding role for computational methods in our understanding of HIV biology and SBDD in the future.

AutoDockFR: Advances in Protein-Ligand Docking with Explicitly Specified Binding Site Flexibility.

Automated docking of drug-like molecules into receptors is an essential tool in structure-based drug design. While modeling receptor flexibility is important for correctly predicting ligand binding, it still remains challenging. This work focuses on an approach in which receptor flexibility is modeled by explicitly specifying a set of receptor side-chains a-priori. The challenges of this approach include the: 1) exponential growth of the search space, demanding more efficient search methods; and 2) increased number of false positives, calling for scoring functions tailored for flexible receptor docking. We present AutoDockFR-AutoDock for Flexible Receptors (ADFR), a new docking engine based on the AutoDock4 scoring function, which addresses the aforementioned challenges with a new Genetic Algorithm (GA) and customized scoring function. We validate ADFR using the Astex Diverse Set, demonstrating an increase in efficiency and reliability of its GA over the one implemented in AutoDock4. We demonstrate greatly increased success rates when cross-docking ligands into apo receptors that require side-chain conformational changes for ligand binding. These cross-docking experiments are based on two datasets: 1) SEQ17 -a receptor diversity set containing 17 pairs of apo-holo structures; and 2) CDK2 -a ligand diversity set composed of one CDK2 apo structure and 52 known bound inhibitors. We show that, when cross-docking ligands into the apo conformation of the receptors with up to 14 flexible side-chains, ADFR reports more correctly cross-docked ligands than AutoDock Vina on both datasets with solutions found for 70.6% vs. 35.3% systems on SEQ17, and 76.9% vs. 61.5% on CDK2. ADFR also outperforms AutoDock Vina in number of top ranking solutions on both datasets. Furthermore, we show that correctly docked CDK2 complexes re-create on average 79.8% of all pairwise atomic interactions between the ligand and moving receptor atoms in the holo complexes. Finally, we show that down-weighting the receptor internal energy improves the ranking of correctly docked poses and that runtime for AutoDockFR scales linearly when side-chain flexibility is added.

Fragment-Based Analysis of Ligand Dockings Improves Classification of Actives.

We describe ADChemCast, a method for using results from virtual screening to create a richer representation of a target binding site, which may be used to improve ranking of compounds and characterize the determinants of ligand-receptor specificity. ADChemCast clusters docked conformations of ligands based on shared pairwise receptor-ligand interactions within chemically similar structural fragments, building a set of attributes characteristic of binders and nonbinders. Machine learning is then used to build rules from the most informational attributes for use in reranking of compounds. In this report, we use ADChemCast to improve the ranking of compounds in 11 diverse proteins from the Database of Useful Decoys-Enhanced (DUD-E) and demonstrate the utility of the method for characterizing relevant binding attributes in HIV reverse transcriptase.

A virtual screen discovers novel, fragment-sized inhibitors of Mycobacterium tuberculosis InhA.

Isoniazid (INH) is usually administered to treat latent Mycobacterium tuberculosis (Mtb) infections and is used in combination therapy to treat active tuberculosis (TB). Unfortunately, resistance to this drug is hampering its clinical effectiveness. INH is a prodrug that must be activated by Mtb catalase-peroxidase (KatG) before it can inhibit InhA (Mtb enoyl-acyl-carrier-protein reductase). Isoniazid-resistant cases of TB found in clinical settings usually involve mutations in or deletion of katG, which abrogate INH activation. Compounds that inhibit InhA without requiring prior activation by KatG would not be affected by this resistance mechanism and hence would display continued potency against these drug-resistant isolates of Mtb. Virtual screening experiments versus InhA in the GO Fight Against Malaria (GO FAM) project were designed to discover new scaffolds that display base-stacking interactions with the NAD cofactor. GO FAM experiments included targets from other pathogens, including Mtb, when they had structural similarity to a malaria target. Eight of the 16 soluble compounds identified by docking against InhA plus visual inspection were modest inhibitors and did not require prior activation by KatG. The best two inhibitors discovered are both fragment-sized compounds and displayed Ki values of 54 and 59 µM, respectively. Importantly, the novel inhibitors discovered have low structural similarity to known InhA inhibitors and thus help expand the number of chemotypes on which future medicinal chemistry efforts can be focused. These new fragment hits could eventually help advance the fight against INH-resistant Mtb strains, which pose a significant global health threat.

AutoDock4(Zn): an improved AutoDock force field for small-molecule docking to zinc metalloproteins.

Zinc is present in a wide variety of proteins and is important in the metabolism of most organisms. Zinc metalloenzymes are therapeutically relevant targets in diseases such as cancer, heart disease, bacterial infection, and Alzheimer's disease. In most cases a drug molecule targeting such enzymes establishes an interaction that coordinates with the zinc ion. Thus, accurate prediction of the interaction of ligands with zinc is an important aspect of computational docking and virtual screening against zinc containing proteins. We have extended the AutoDock force field to include a specialized potential describing the interactions of zinc-coordinating ligands. This potential describes both the energetic and geometric components of the interaction. The new force field, named AutoDock4Zn, was calibrated on a data set of 292 crystal complexes containing zinc. Redocking experiments show that the force field provides significant improvement in performance in both free energy of binding estimation as well as in root-mean-square deviation from the crystal structure pose. The new force field has been implemented in AutoDock without modification to the source code.

Virtual screening with AutoDock Vina and the common pharmacophore engine of a low diversity library of fragments and hits against the three allosteric sites of HIV integrase: participation in the SAMPL4 protein-ligand binding challenge.

To rigorously assess the tools and protocols that can be used to understand and predict macromolecular recognition, and to gain more structural insight into three newly discovered allosteric binding sites on a critical drug target involved in the treatment of HIV infections, the Olson and Levy labs collaborated on the SAMPL4 challenge. This computational blind challenge involved predicting protein-ligand binding against the three allosteric sites of HIV integrase (IN), a viral enzyme for which two drugs (that target the active site) have been approved by the FDA. Positive control cross-docking experiments were utilized to select 13 receptor models out of an initial ensemble of 41 different crystal structures of HIV IN. These 13 models of the targets were selected using our new "Rank Difference Ratio" metric. The first stage of SAMPL4 involved using virtual screens to identify 62 active, allosteric IN inhibitors out of a set of 321 compounds. The second stage involved predicting the binding site(s) and crystallographic binding mode(s) for 57 of these inhibitors. Our team submitted four entries for the first stage that utilized: (1) AutoDock Vina (AD Vina) plus visual inspection; (2) a new common pharmacophore engine; (3) BEDAM replica exchange free energy simulations, and a Consensus approach that combined the predictions of all three strategies. Even with the SAMPL4's very challenging compound library that displayed a significantly lower amount of structural diversity than most libraries that are conventionally employed in prospective virtual screens, these approaches produced hit rates of 24, 25, 34, and 27 %, respectively, on a set with 19 % declared binders. Our only entry for the second stage challenge was based on the results of AD Vina plus visual inspection, and it ranked third place overall according to several different metrics provided by the SAMPL4 organizers. The successful results displayed by these approaches highlight the utility of the computational structure-based drug discovery tools and strategies that are being developed to advance the goals of the newly created, multi-institution, NIH-funded center called the "HIV Interaction and Viral Evolution Center".

Blind prediction of HIV integrase binding from the SAMPL4 challenge.

Here, we give an overview of the protein-ligand binding portion of the Statistical Assessment of Modeling of Proteins and Ligands 4 (SAMPL4) challenge, which focused on predicting binding of HIV integrase inhibitors in the catalytic core domain. The challenge encompassed three components--a small "virtual screening" challenge, a binding mode prediction component, and a small affinity prediction component. Here, we give summary results and statistics concerning the performance of all submissions at each of these challenges. Virtual screening was particularly challenging here in part because, in contrast to more typical virtual screening test sets, the inactive compounds were tested because they were thought to be likely binders, so only the very top predictions performed significantly better than random. Pose prediction was also quite challenging, in part because inhibitors in the set bind to three different sites, so even identifying the correct binding site was challenging. Still, the best methods managed low root mean squared deviation predictions in many cases. Here, we give an overview of results, highlight some features of methods which worked particularly well, and refer the interested reader to papers in this issue which describe specific submissions for additional details.

Virtual screening of integrase inhibitors by large scale binding free energy calculations: the SAMPL4 challenge.

As part of the SAMPL4 blind challenge, filtered AutoDock Vina ligand docking predictions and large scale binding energy distribution analysis method binding free energy calculations have been applied to the virtual screening of a focused library of candidate binders to the LEDGF site of the HIV integrase protein. The computational protocol leveraged docking and high level atomistic models to improve enrichment. The enrichment factor of our blind predictions ranked best among all of the computational submissions, and second best overall. This work represents to our knowledge the first example of the application of an all-atom physics-based binding free energy model to large scale virtual screening. A total of 285 parallel Hamiltonian replica exchange molecular dynamics absolute protein-ligand binding free energy simulations were conducted starting from docked poses. The setup of the simulations was fully automated, calculations were distributed on multiple computing resources and were completed in a 6-weeks period. The accuracy of the docked poses and the inclusion of intramolecular strain and entropic losses in the binding free energy estimates were the major factors behind the success of the method. Lack of sufficient time and computing resources to investigate additional protonation states of the ligands was a major cause of mispredictions. The experiment demonstrated the applicability of binding free energy modeling to improve hit rates in challenging virtual screening of focused ligand libraries during lead optimization.

Identification of clickable HIV-1 capsid-targeting probes for viral replication inhibition.

Of the targets for HIV-1 therapeutics, the capsid core is a relatively unexploited but alluring drug target due to its indispensable roles throughout virus replication. Because of this, we aimed to identify "clickable" covalent modifiers of the HIV-1 capsid protein (CA) for future functionalization. We screened a library of fluorosulfate compounds that can undergo sulfur(VI) fluoride exchange (SuFEx) reactions, and five compounds were identified as hits. These molecules were further characterized for antiviral effects. Several compounds impacted in vitro capsid assembly. One compound, BBS-103, covalently bound CA via a SuFEx reaction to Tyr145 and had antiviral activity in cell-based assays by perturbing virus production, but not uncoating. The covalent binding of compounds that target the HIV-1 capsid could aid in the future design of antiretroviral drugs or chemical probes that will help study aspects of HIV-1 replication.

Visualizing biological data-now and in the future.

Methods and tools for visualizing biological data have improved considerably over the last decades, but they are still inadequate for some high-throughput data sets. For most users, a key challenge is to benefit from the deluge of data without being overwhelmed by it. This challenge is still largely unfulfilled and will require the development of truly integrated and highly useable tools.

Visualization of macromolecular structures.

Structural biology is rapidly accumulating a wealth of detailed information about protein function, binding sites, RNA, large assemblies and molecular motions. These data are increasingly of interest to a broader community of life scientists, not just structural experts. Visualization is a primary means for accessing and using these data, yet visualization is also a stumbling block that prevents many life scientists from benefiting from three-dimensional structural data. In this review, we focus on key biological questions where visualizing three-dimensional structures can provide insight and describe available methods and tools.