Gaurav Chopra

KoBaMIN: a knowledge-based minimization web server for protein structure refinement

The KoBaMIN web server provides an online interface to a simple, consistent and computationally efficient protein structure refinement protocol based on minimization of a knowledge-based potential of mean force. The server can be used to refine either a single protein structure or an ensemble of proteins starting from their unrefined coordinates in PDB format. The refinement method is particularly fast and accurate due to the underlying knowledge-based potential derived from structures deposited in the PDB; as such, the energy function implicitly includes the effects of solvent and the crystal environment. Our server allows for an optional but recommended step that optimizes stereochemistry using the MESHI software. The KoBaMIN server also allows comparison of the refined structures with a provided reference structure to assess the changes brought about by the refinement protocol. The performance of KoBaMIN has been benchmarked widely on a large set of decoys, all models generated at the seventh worldwide experiments on critical assessment of techniques for protein structure prediction (CASP7) and it was also shown to produce top-ranking predictions in the refinement category at both CASP8 and CASP9, yielding consistently good results across a broad range of model quality values. The web server is fully functional and freely available at http://csb.stanford.edu/kobamin.

Solvent dramatically affects protein structure refinement

One of the most challenging problems in protein structure prediction is improvement of homology models (structures within 1–3 Å Cα rmsd of the native structure), also known as the protein structure refinement problem. It has been shown that improvement could be achieved using in vacuo energy minimization with molecular mechanics and statistically derived continuously differentiable hybrid knowledge-based (KB) potential functions. Globular proteins, however, fold and function in aqueous solution in vivo and in vitro. In this work, we study the role of solvent in protein structure refinement. Molecular dynamics in explicit solvent and energy minimization in both explicit and implicit solvent were performed on a set of 75 native proteins to test the various energy potentials. A more stringent test for refinement was performed on 729 near-native decoys for each native protein. We use a powerfully convergent energy minimization method to show that implicit solvent (GBSA) provides greater improvement for some proteins than the KB potential: 24 of 75 proteins showing an average improvement of >20% in Cα rmsd from the native structure with GBSA, compared to just 7 proteins with KB. Molecular dynamics in explicit solvent moved the structures further away from their native conformation than the initial, unrefined decoys. Implicit solvent gives rise to a deep, smooth potential energy attractor basin that pulls toward the native structure.

Divergent Phenotypes of Human Regulatory T Cells Expressing the Receptors TIGIT and CD226

Regulatory T cells (Tregs) play a central role in counteracting inflammation and autoimmunity. A more complete understanding of cellular heterogeneity and the potential for lineage plasticity in human Treg subsets may identify markers of disease pathogenesis and facilitate the development of optimized cellular therapeutics. To better elucidate human Treg subsets, we conducted direct transcriptional profiling of CD4+FOXP3+Helios+ thymic-derived Tregs and CD4+FOXP3+Helios− T cells, followed by comparison with CD4+FOXP3−Helios− T conventional cells. These analyses revealed that the coinhibitory receptor T cell Ig and ITIM domain (TIGIT) was highly expressed on thymic-derived Tregs. TIGIT and the costimulatory factor CD226 bind the common ligand CD155. Thus, we analyzed the cellular distribution and suppressive activity of isolated subsets of CD4+CD25+CD127lo/− T cells expressing CD226 and/or TIGIT. We observed TIGIT is highly expressed and upregulated on Tregs after activation and in vitro expansion, and is associated with lineage stability and suppressive capacity. Conversely, the CD226+TIGIT− population was associated with reduced Treg purity and suppressive capacity after expansion, along with a marked increase in IL-10 and effector cytokine production. These studies provide additional markers to delineate functionally distinct Treg subsets that may help direct cellular therapies and provide important phenotypic markers for assessing the role of Tregs in health and disease.

Consistent Refinement of Submitted Models at CASP using a Knowledge-based Potential

Protein structure refinement is an important but unsolved problem; it must be solved if we are to predict biological function that is very sensitive to structural details. Specifically, critical assessment of techniques for protein structure prediction (CASP) shows that the accuracy of predictions in the comparative modeling category is often worse than that of the template on which the homology model is based. Here we describe a refinement protocol that is able to consistently refine submitted predictions for all categories at CASP7. The protocol uses direct energy minimization of the knowledge‐based potential of mean force that is based on the interaction statistics of 167 atom types (Summa and Levitt, Proc Natl Acad Sci USA 2007; 104:3177–3182). Our protocol is thus computationally very efficient; it only takes a few minutes of CPU time to run typical protein models (300 residues). We observe an average structural improvement of 1% in GDT_TS, for predictions that have low and medium homology to known PDB structures (Global Distance Test score or GDT_TS between 50 and 80%). We also observe a marked improvement in the stereochemistry of the models. The level of improvement varies amongst the various participants at CASP, but we see large improvements (>10% increase in GDT_TS) even for models predicted by the best performing groups at CASP7. In addition, our protocol consistently improved the best predicted models in the refinement category at CASP7 and CASP8. These improvements in structure and stereochemistry prove the usefulness of our computationally inexpensive, powerful and automatic refinement protocol. Proteins 2010.

WeFold: A coopetition for protein structure prediction

The protein structure prediction problem continues to elude scientists. Despite the introduction of many methods, only modest gains were made over the last decade for certain classes of prediction targets. To address this challenge, a social‐media based worldwide collaborative effort, named WeFold, was undertaken by 13 labs. During the collaboration, the laboratories were simultaneously competing with each other. Here, we present the first attempt at “coopetition” in scientific research applied to the protein structure prediction and refinement problems. The coopetition was possible by allowing the participating labs to contribute different components of their protein structure prediction pipelines and create new hybrid pipelines that they tested during CASP10. This manuscript describes both successes and areas needing improvement as identified throughout the first WeFold experiment and discusses the efforts that are underway to advance this initiative. A footprint of all contributions and structures are publicly accessible at http://www.wefold.org. Proteins 2014; 82:1850–1868.

Oxido-reductive regulation of vascular remodeling by receptor tyrosine kinase ROS1

Angioplasty and stenting is the primary treatment for flow-limiting atherosclerosis; however, this strategy is limited by pathological vascular remodeling. Using a systems approach, we identified a role for the network hub gene glutathione peroxidase-1 (GPX1) in pathological remodeling following human blood vessel stenting. Constitutive deletion of Gpx1 in atherosclerotic mice recapitulated this phenotype of increased vascular smooth muscle cell (VSMC) proliferation and plaque formation. In an independent patient cohort, gene variant pair analysis identified an interaction of GPX1 with the orphan protooncogene receptor tyrosine kinase ROS1. A meta-analysis of the only genome-wide association studies of human neointima-induced in-stent stenosis confirmed the association of the ROS1 variant with pathological remodeling. Decreased GPX1 expression in atherosclerotic mice led to reductive stress via a time-dependent increase in glutathione, corresponding to phosphorylation of the ROS1 kinase activation site Y2274. Loss of GPX1 function was associated with both oxidative and reductive stress, the latter driving ROS1 activity via s-glutathiolation of critical residues of the ROS1 tyrosine phosphatase SHP-2. ROS1 inhibition with crizotinib and deglutathiolation of SHP-2 abolished GPX1-mediated increases in VSMC proliferation while leaving endothelialization intact. Our results indicate that GPX1-dependent alterations in oxido-reductive stress promote ROS1 activation and mediate vascular remodeling.

Remarkable patterns of surface water ordering around polarized buckminsterfullerene

Accurate description of water structure affects simulation of protein folding, substrate binding, macromolecular recognition, and complex formation. We study the hydration of buckminsterfullerene, the smallest hydrophobic nanosphere, by molecular dynamics simulations using a state-of-the-art quantum mechanical polarizable force field (QMPFF3), derived from quantum mechanical data at the MP2/aug-cc-pVTZ(-hp) level augmented by CCSD(T). QMPFF3 calculation of the hydrophobic effect is compared to that obtained with empirical force fields. Using a novel and highly sensitive method, we see polarization increases ordered water structure so that the imprint of the hydrophobic surface atoms on the surrounding waters is stronger and extends to long-range. We see less water order for empirical force fields. The greater order seen with QMPFF3 will affect biological processes through a stronger hydrophobic effect.

Strategic Protein Target Analysis for Developing Drugs to Stop Dental Caries

Dental caries is the most common disease to cause irreversible damage in humans. Several therapeutic agents are available to treat or prevent dental caries, but none besides fluoride has significantly influenced the disease burden globally. Etiologic mechanisms of the mutans group streptococci and specific Lactobacillus species have been characterized to various degrees of detail, from identification of physiologic processes to specific proteins. Here, we analyze the entire Streptococcus mutans proteome for potential drug targets by investigating their uniqueness with respect to non-cariogenic dental plaque bacteria, quality of protein structure models, and the likelihood of finding a drug for the active site. Our results suggest specific targets for rational drug discovery, including 15 known virulence factors, 16 proteins for which crystallographic structures are available, and 84 previously uncharacterized proteins, with various levels of similarity to homologs in dental plaque bacteria. This analysis provides a map to streamline the process of clinical development of effective multispecies pharmacologic interventions for dental caries.

Distal Effect of Amino Acid Substitutions in CYP2C9 Polymorphic Variants Causes Differences in Interatomic Interactions against (S)-Warfarin

Cytochrome P450 2C9 (CYP2C9) is crucial in excretion of commonly prescribed drugs. However, changes in metabolic activity caused by CYP2C9 polymorphisms inevitably result in adverse drug effects. CYP2C9*2 and *3 are prevalent in Caucasian populations whereas CYP2C9*13 is remarkable in Asian populations. Single amino acid substitutions caused by these mutations are located outside catalytic cavity but affect kinetic activities of mutants compared to wild-type enzyme. To relate distal effects of these mutations and defective drug metabolisms, simulations of CYP2C9 binding to anti-coagulant (S)-warfarin were performed as a system model. Representative (S)-warfarin-bound forms of wild-type and mutants were sorted and assessed through knowledge-based scoring function. Interatomic interactions towards (S)-warfarin were predicted to be less favorable in mutant structures in correlation with larger distance between hydroxylation site of (S)-warfarin and reactive oxyferryl heme than wild-type structure. Using computational approach could delineate complication of CYP polymorphism in management of drug therapy.

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.