Surrogate Modeling of the Relative Entropy for Inverse Design Using Smolyak Sparse Grids.

Relative entropy minimization, a statistical-mechanics approach for finding potential energy functions that produce target structural ensembles, has proven to be a powerful strategy for the inverse design of nanoparticle self-assembly. For a given target structure, the gradient of the relative entropy with respect to the adjustable parameters of the potential energy function is computed by performing a simulation, and then these parameters are updated using iterative gradient-based optimization. Small parameter updates per iteration and many iterations can be required for numerical stability, but this incurs considerable computational expense because a new simulation must be performed to reevaluate the gradient at each iteration. Here, we investigate the use of surrogate modeling to decouple the process of minimizing the relative entropy from the computationally demanding process of determining its gradient. We approximate the relative-entropy gradient using Chebyshev polynomial interpolation on Smolyak sparse grids. Our approach potentially increases the robustness and computational efficiency of using the relative entropy for inverse design, primarily for physically informed potential energy functions that have a small number of adjustable parameters.

Role of Electrostatic Hotspots in the Selectivity of Complement Control Proteins Toward Human and Bovine Complement Inhibition.

Poxviruses are dangerous pathogens, which can cause fatal infection in unvaccinated individuals. The causative agent of smallpox in humans, variola virus, is closely related to the bovine vaccinia virus, yet the molecular basis of their selectivity is currently incompletely understood. Here, we examine the role of the electrostatics in the selectivity of the smallpox protein SPICE and vaccinia protein VCP toward the human and bovine complement protein C3b, a key component of the complement immune response. Electrostatic calculations, in-silico alanine-scan and electrostatic hotspot analysis, as introduced by Kieslich and Morikis (PLoS Comput. Biol. 2012), are used to assess the electrostatic complementarity and to identify sites resistant to local perturbation where the electrostatic potential is likely to be evolutionary conserved. The calculations suggest that the bovine C3b is electrostatically prone to selectively bind its VCP ligand. On the other hand, the human isoform of C3b exhibits a lower electrostatic complementarity toward its SPICE ligand. Yet, the human C3b displays a highly preserved electrostatic core, which suggests that this isoform could be less selective in binding different ligands like SPICE and the human Factor H. This is supported by experimental cofactor activity assays revealing that the human C3b is prone to bind both SPICE and Factor H, which exhibit diverse electrostatic properties. Additional investigations considering mutants of SPICE and VCP that revert their selectivity reveal an "electrostatic switch" into the central modules of the ligands, supporting the critical role of the electrostatics in the selectivity. Taken together, these evidences provide insights into the selectivity mechanism of the complement regulator proteins encoded by the variola and vaccinia viruses to circumvent the complement immunity and exert their pathogenic action. These fundamental aspects are valuable for the development of novel vaccines and therapeutic strategies.

Reassessing enzyme kinetics: Considering protease-as-substrate interactions in proteolytic networks.

Enzymes are catalysts in biochemical reactions that, by definition, increase rates of reactions without being altered or destroyed. However, when that enzyme is a protease, a subclass of enzymes that hydrolyze other proteins, and that protease is in a multiprotease system, protease-as-substrate dynamics must be included, challenging assumptions of enzyme inertness, shifting kinetic predictions of that system. Protease-on-protease inactivating hydrolysis can alter predicted protease concentrations used to determine pharmaceutical dosing strategies. Cysteine cathepsins are proteases capable of cathepsin cannibalism, where one cathepsin hydrolyzes another with substrate present, and misunderstanding of these dynamics may cause miscalculations of multiple proteases working in one proteolytic network of interactions occurring in a defined compartment. Once rates for individual protease-on-protease binding and catalysis are determined, proteolytic network dynamics can be explored using computational models of cooperative/competitive degradation by multiple proteases in one system, while simultaneously incorporating substrate cleavage. During parameter optimization, it was revealed that additional distraction reactions, where inactivated proteases become competitive inhibitors to remaining, active proteases, occurred, introducing another network reaction node. Taken together, improved predictions of substrate degradation in a multiple protease network were achieved after including reaction terms of autodigestion, inactivation, cannibalism, and distraction, altering kinetic considerations from other enzymatic systems, since enzyme can be lost to proteolytic degradation. We compiled and encoded these dynamics into an online platform (https://plattlab.shinyapps.io/catKLS/) for individual users to test hypotheses of specific perturbations to multiple cathepsins, substrates, and inhibitors, and predict shifts in proteolytic network reactions and system dynamics.

A Nonlinear Support Vector Machine-Based Feature Selection Approach for Fault Detection and Diagnosis: Application to the Tennessee Eastman Process.

In this article, we present (1) a feature selection algorithm based on nonlinear support vector machine (SVM) for fault detection and diagnosis in continuous processes and (2) results for the Tennessee Eastman benchmark process. The presented feature selection algorithm is derived from the sensitivity analysis of the dual C-SVM objective function. This enables simultaneous modeling and feature selection paving the way for simultaneous fault detection and diagnosis, where feature ranking guides fault diagnosis. We train fault-specific two-class SVM models to detect faulty operations, while using the feature selection algorithm to improve the accuracy and perform the fault diagnosis. Our results show that the developed SVM models outperform the available ones in the literature both in terms of detection accuracy and latency. Moreover, it is shown that the loss of information is minimized with the use of feature selection techniques compared to feature extraction techniques such as principal component analysis (PCA). This further facilitates a more accurate interpretation of the results.

Simultaneous Fault Detection and Identification in Continuous Processes via nonlinear Support Vector Machine based Feature Selection.

Rapid detection and identification of process faults in industrial applications is crucial to sustain a safe and profitable operation. Today, the advances in sensor technologies have facilitated large amounts of chemical process data collection in real time which subsequently broadened the use of data-driven process monitoring techniques via machine learning and multivariate statistical analysis. One of the well-known machine learning techniques is Support Vector Machines (SVM) which allows the use of high dimensional feature sets for learning problems such as classification and regression. In this paper, we present the application of a novel nonlinear (kernel-dependent) SVM-based feature selection algorithm to process monitoring and fault detection of continuous processes. The developed methodology is derived from sensitivity analysis of the dual SVM objective and utilizes existing and novel greedy algorithms to rank features that also guides fault diagnosis. Specifically, we train fault-specific two-class SVM models to detect faulty operations, while using the feature selection algorithm to improve the accuracy of the fault detection models and perform fault diagnosis. We present results for the Tennessee Eastman process as a case study and compare our approach to existing approaches for fault detection, diagnosis and identification.

Big Data Approach to Batch Process Monitoring: Simultaneous Fault Detection and Diagnosis Using Nonlinear Support Vector Machine-based Feature Selection.

This paper presents a novel data-driven framework for process monitoring in batch processes, a critical task in industry to attain a safe operability and minimize loss of productivity and profit. We exploit high dimensional process data with nonlinear Support Vector Machine-based feature selection algorithm, where we aim to retrieve the most informative process measurements for accurate and simultaneous fault detection and diagnosis. The proposed framework is applied to an extensive benchmark dataset which includes process data describing 22,200 batches with 15 faults. We train fault and time-specific models on the prealigned batch data trajectories via three distinct time horizon approaches: one-step rolling, two-step rolling, and evolving which varies the amount of data incorporation during modeling. The results show that two-step rolling and evolving time horizon approaches perform superior to the other. Regardless of the approach, proposed framework provides a promising decision support tool for online simultaneous fault detection and diagnosis for batch processes.

Virtual Screening of Chemical Compounds for Discovery of Complement C3 Ligands.

The complement system is our first line of defense against foreign pathogens, but when it is not properly regulated, complement is implicated in the pathology of several autoimmune and inflammatory disorders. Compstatin is a peptidic complement inhibitor that acts by blocking the cleavage of complement protein C3 to the proinflammatory fragment C3a and opsonin fragment C3b. In this study, we aim to identify druglike small-molecule complement inhibitors with physicochemical, geometric, and binding properties similar to those of compstatin. We employed two approaches using various high-throughput virtual screening methods, which incorporate molecular dynamics (MD) simulations, pharmacophore model design, energy calculations, and molecular docking and scoring. We have generated a library of 274 chemical compounds with computationally predicted binding affinities for the compstatin binding site of C3. We have tested subsets of these chemical compounds experimentally for complement inhibitory activity, using hemolytic assays, and for binding affinity, using microscale thermophoresis. As a result, although none of the compounds showed inhibitory activity, compound 29 was identified to exhibit weak competitive binding against a potent compstatin analogue, therefore validating our computational approaches. Additional docking and MD simulation studies suggest that compound 29 interacts with C3 residues, which have been shown to be important in binding of compstatin to the C3c fragment of C3. Compound 29 is amenable to physicochemical optimization to acquire inhibitory properties. Additionally, it is possible that some of the untested compounds will demonstrate binding and inhibition in future experimental studies.

An analysis and evaluation of the WeFold collaborative for protein structure prediction and its pipelines in CASP11 and CASP12.

Every two years groups worldwide participate in the Critical Assessment of Protein Structure Prediction (CASP) experiment to blindly test the strengths and weaknesses of their computational methods. CASP has significantly advanced the field but many hurdles still remain, which may require new ideas and collaborations. In 2012 a web-based effort called WeFold, was initiated to promote collaboration within the CASP community and attract researchers from other fields to contribute new ideas to CASP. Members of the WeFold coopetition (cooperation and competition) participated in CASP as individual teams, but also shared components of their methods to create hybrid pipelines and actively contributed to this effort. We assert that the scale and diversity of integrative prediction pipelines could not have been achieved by any individual lab or even by any collaboration among a few partners. The models contributed by the participating groups and generated by the pipelines are publicly available at the WeFold website providing a wealth of data that remains to be tapped. Here, we analyze the results of the 2014 and 2016 pipelines showing improvements according to the CASP assessment as well as areas that require further adjustments and research.

AESOP: A Python Library for Investigating Electrostatics in Protein Interactions.

Electric fields often play a role in guiding the association of protein complexes. Such interactions can be further engineered to accelerate complex association, resulting in protein systems with increased productivity. This is especially true for enzymes where reaction rates are typically diffusion limited. To facilitate quantitative comparisons of electrostatics in protein families and to describe electrostatic contributions of individual amino acids, we previously developed a computational framework called AESOP. We now implement this computational tool in Python with increased usability and the capability of performing calculations in parallel. AESOP utilizes PDB2PQR and Adaptive Poisson-Boltzmann Solver to generate grid-based electrostatic potential files for protein structures provided by the end user. There are methods within AESOP for quantitatively comparing sets of grid-based electrostatic potentials in terms of similarity or generating ensembles of electrostatic potential files for a library of mutants to quantify the effects of perturbations in protein structure and protein-protein association.

Princeton_TIGRESS 2.0: High refinement consistency and net gains through support vector machines and molecular dynamics in double-blind predictions during the CASP11 experiment.

Protein structure refinement is the challenging problem of operating on any protein structure prediction to improve its accuracy with respect to the native structure in a blind fashion. Although many approaches have been developed and tested during the last four CASP experiments, a majority of the methods continue to degrade models rather than improve them. Princeton_TIGRESS (Khoury et al., Proteins 2014;82:794-814) was developed previously and utilizes separate sampling and selection stages involving Monte Carlo and molecular dynamics simulations and classification using an SVM predictor. The initial implementation was shown to consistently refine protein structures 76% of the time in our own internal benchmarking on CASP 7-10 targets. In this work, we improved the sampling and selection stages and tested the method in blind predictions during CASP11. We added a decomposition of physics-based and hybrid energy functions, as well as a coordinate-free representation of the protein structure through distance-binning Cα-Cα distances to capture fine-grained movements. We performed parameter estimation to optimize the adjustable SVM parameters to maximize precision while balancing sensitivity and specificity across all cross-validated data sets, finding enrichment in our ability to select models from the populations of similar decoys generated for targets in CASPs 7-10. The MD stage was enhanced such that larger structures could be further refined. Among refinement methods that are currently implemented as web-servers, Princeton_TIGRESS 2.0 demonstrated the most consistent and most substantial net refinement in blind predictions during CASP11. The enhanced refinement protocol Princeton_TIGRESS 2.0 is freely available as a web server at http://atlas.engr.tamu.edu/refinement/. Proteins 2017; 85:1078-1098. © 2017 Wiley Periodicals, Inc.

conSSert: Consensus SVM Model for Accurate Prediction of Ordered Secondary Structure.

Accurate prediction of protein secondary structure remains a crucial step in most approaches to the protein-folding problem, yet the prediction of ordered secondary structure, specifically beta-strands, remains a challenge. We developed a consensus secondary structure prediction method, conSSert, which is based on support vector machines (SVM) and provides exceptional accuracy for the prediction of beta-strands with QE accuracy of over 0.82 and a Q2-EH of 0.86. conSSert uses as input probabilities for the three types of secondary structure (helix, strand, and coil) that are predicted by four top performing methods: PSSpred, PSIPRED, SPINE-X, and RAPTOR. conSSert was trained/tested using 4261 protein chains from PDBSelect25, and 8632 chains from PISCES. Further validation was performed using targets from CASP9, CASP10, and CASP11. Our data suggest that poor performance in strand prediction is likely a result of training bias and not solely due to the nonlocal nature of beta-sheet contacts. conSSert is freely available for noncommercial use as a webservice: http://ares.tamu.edu/conSSert/.

Highly Accurate Structure-Based Prediction of HIV-1 Coreceptor Usage Suggests Intermolecular Interactions Driving Tropism.

HIV-1 entry into host cells is mediated by interactions between the V3-loop of viral glycoprotein gp120 and chemokine receptor CCR5 or CXCR4, collectively known as HIV-1 coreceptors. Accurate genotypic prediction of coreceptor usage is of significant clinical interest and determination of the factors driving tropism has been the focus of extensive study. We have developed a method based on nonlinear support vector machines to elucidate the interacting residue pairs driving coreceptor usage and provide highly accurate coreceptor usage predictions. Our models utilize centroid-centroid interaction energies from computationally derived structures of the V3-loop:coreceptor complexes as primary features, while additional features based on established rules regarding V3-loop sequences are also investigated. We tested our method on 2455 V3-loop sequences of various lengths and subtypes, and produce a median area under the receiver operator curve of 0.977 based on 500 runs of 10-fold cross validation. Our study is the first to elucidate a small set of specific interacting residue pairs between the V3-loop and coreceptors capable of predicting coreceptor usage with high accuracy across major HIV-1 subtypes. The developed method has been implemented as a web tool named CRUSH, CoReceptor USage prediction for HIV-1, which is available at http://ares.tamu.edu/CRUSH/.

Derivation of ligands for the complement C3a receptor from the C-terminus of C5a.

The complement cascade is a highly sophisticated network of proteins that are well regulated and directed in response to invading pathogens or tissue injury. Complement C3a and C5a are key mediators produced by this cascade, and their dysregulation has been linked to a plethora of inflammatory and autoimmune diseases. Consequently, this has stimulated interest in the development of ligands for the receptors for these complement peptides, C3a receptor, and C5a1 (C5aR/CD88). In this study we used computational methods to design novel C5a1 receptor ligands. However, functional screening in human monocyte-derived macrophages using the xCELLigence label-free platform demonstrated altered specificity of our ligands. No agonist/antagonist activity was observed at C5a1, but we instead saw that the ligands were able to partially agonize the closely related complement receptor C3a receptor. This was verified in the presence of C3a receptor antagonist SB 290157 and in a stable cell line expressing either C5a1 or C3a receptor alone. C3a agonism has been suggested to be a potential treatment of acute neutrophil-driven traumatic pathologies, and may have great potential as a therapeutic avenue in this arena.

Insights into the mechanism of C5aR inhibition by PMX53 via implicit solvent molecular dynamics simulations and docking.

BACKGROUND: The complement protein C5a acts by primarily binding and activating the G-protein coupled C5a receptor C5aR (CD88), and is implicated in many inflammatory diseases. The cyclic hexapeptide PMX53 (sequence Ace-Phe-[Orn-Pro-dCha-Trp-Arg]) is a full C5aR antagonist of nanomolar potency, and is widely used to study C5aR function in disease. RESULTS: We construct for the first time molecular models for the C5aR:PMX53 complex without the a priori use of experimental constraints, via a computational framework of molecular dynamics (MD) simulations, docking, conformational clustering and free energy filtering. The models agree with experimental data, and are used to propose important intermolecular interactions contributing to binding, and to develop a hypothesis for the mechanism of PMX53 antagonism. CONCLUSION: This work forms the basis for the design of improved C5aR antagonists, as well as for atomic-detail mechanistic studies of complement activation and function. Our computational framework can be widely used to develop GPCR-ligand structural models in membrane environments, peptidomimetics and other chemical compounds with potential clinical use.

Forcefield_NCAA: ab initio charge parameters to aid in the discovery and design of therapeutic proteins and peptides with unnatural amino acids and their application to complement inhibitors of the compstatin family.

We describe the development and testing of ab initio derived, AMBER ff03 compatible charge parameters for a large library of 147 noncanonical amino acids including ß- and N-methylated amino acids for use in applications such as protein structure prediction and de novo protein design. The charge parameter derivation was performed using the RESP fitting approach. Studies were performed assessing the suitability of the derived charge parameters in discriminating the activity/inactivity between 63 analogs of the complement inhibitor Compstatin on the basis of previously published experimental IC50 data and a screening procedure involving short simulations and binding free energy calculations. We found that both the approximate binding affinity (K*) and the binding free energy calculated through MM-GBSA are capable of discriminating between active and inactive Compstatin analogs, with MM-GBSA performing significantly better. Key interactions between the most potent Compstatin analog that contains a noncanonical amino acid are presented and compared to the most potent analog containing only natural amino acids and native Compstatin. We make the derived parameters and an associated web interface that is capable of performing modifications on proteins using Forcefield_NCAA and outputting AMBER-ready topology and parameter files freely available for academic use at http://selene.princeton.edu/FFNCAA . The forcefield allows one to incorporate these customized amino acids into design applications with control over size, van der Waals, and electrostatic interactions.

Engineering pre-SUMO4 as efficient substrate of SENP2.

SUMOylation, one of the most important protein post-translational modifications, plays critical roles in a variety of physiological and pathological processes. SENP (Sentrin/SUMO-specific protease), a family of SUMO-specific proteases, is responsible for the processing of pre-SUMO and removal of SUMO from conjugated substrates. SUMO4, the latest discovered member in the SUMO family, has been found as a type 1 diabetes susceptibility gene and its maturation is not understood so far. Despite the 14 amino acid differences between pre-SUMO4 and SUMO2, pre-SUMO4 is not processed by SENP2 but pre-SUMO2 does. A novel interdisciplinary approach involving computational modeling and a FRET-based protease assay was taken to engineer pre-SUMO4 as a substrate of SENP2. Given the difference in net charge between pre-SUMO4 and pre-SUMO2, the computational framework analysis of electrostatic similarities of proteins was applied to determine the contribution of each ionizable amino acid in a model of SENP2-(pre-SUMO4) binding, and to propose pre-SUMO4 mutations. The specificities of the SENP2 toward different pre-SUMO4 mutants were determined using a quantitative FRET assay by characterizing the catalytic efficiencies (kcat/KM). A single amino acid mutation made pre-SUMO4 amenable to SENP2 processing and a combination of two amino acid mutations made it highly accessible as SENP2 substrate. The combination of the two approaches provides a powerful protein engineering tool for future SUMOylation studies.

Princeton_TIGRESS: protein geometry refinement using simulations and support vector machines.

Protein structure refinement aims to perform a set of operations given a predicted structure to improve model quality and accuracy with respect to the native in a blind fashion. Despite the numerous computational approaches to the protein refinement problem reported in the previous three CASPs, an overwhelming majority of methods degrade models rather than improve them. We initially developed a method tested using blind predictions during CASP10 which was officially ranked in 5th place among all methods in the refinement category. Here, we present Princeton_TIGRESS, which when benchmarked on all CASP 7,8,9, and 10 refinement targets, simultaneously increased GDT_TS 76% of the time with an average improvement of 0.83 GDT_TS points per structure. The method was additionally benchmarked on models produced by top performing three-dimensional structure prediction servers during CASP10. The robustness of the Princeton_TIGRESS protocol was also tested for different random seeds. We make the Princeton_TIGRESS refinement protocol freely available as a web server at http://atlas.princeton.edu/refinement. Using this protocol, one can consistently refine a prediction to help bridge the gap between a predicted structure and the actual native structure.

Protein folding and de novo protein design for biotechnological applications.

In the postgenomic era, the medical/biological fields are advancing faster than ever. However, before the power of full-genome sequencing can be fully realized, the connection between amino acid sequence and protein structure, known as the protein folding problem, needs to be elucidated. The protein folding problem remains elusive, with significant difficulties still arising when modeling amino acid sequences lacking an identifiable template. Understanding protein folding will allow for unforeseen advances in protein design; often referred to as the inverse protein folding problem. Despite challenges in protein folding, de novo protein design has recently demonstrated significant success via computational techniques. We review advances and challenges in protein structure prediction and de novo protein design, and highlight their interplay in successful biotechnological applications.

Novel compstatin family peptides inhibit complement activation by drusen-like deposits in human retinal pigmented epithelial cell cultures.

We have used a novel human retinal pigmented epithelial (RPE) cell-based model that mimics drusen biogenesis and the pathobiology of age-related macular degeneration to evaluate the efficacy of newly designed peptide inhibitors of the complement system. The peptides belong to the compstatin family and, compared to existing compstatin analogs, have been optimized to promote binding to their target, complement protein C3, and to enhance solubility by improving their polarity/hydrophobicity ratios. Based on analysis of molecular dynamics simulation data of peptide-C3 complexes, novel binding features were designed by introducing intermolecular salt bridge-forming arginines at the N-terminus and at position -1 of N-terminal dipeptide extensions. Our study demonstrates that the RPE cell assay has discriminatory capability for measuring the efficacy and potency of inhibitory peptides in a macular disease environment.

A predictive model for HIV type 1 coreceptor selectivity.

Despite its sequence variability and structural flexibility, the V3 loop of the HIV-1 envelope glycoprotein gp120 is capable of recognizing cell-bound coreceptors CCR5 and CXCR4 and infecting cells. Viral selection of CCR5 is associated with the early stages of infection, and transition to selection of CXCR4 indicates disease progression. We have developed a predictive statistical model for coreceptor selectivity that uses the discrete property of net charge and the binary coreceptor preference markers of the N(6)X(7)[T/S](8)X(9) glycosylation motif and 11/24/25 positive amino acid rule. The model is based on analysis of 2,054 V3 loop sequences from patient data and allows us to infer the most likely state of the disease from physicochemical characteristics of the sequences. The performance of the model is comparable to established sequence-based predictive methods, and may be used in combination with other methods as a supportive diagnostic for coreceptor selection. This model may be used for personalized medical decisions in administering coreceptor-specific therapies.