Tributes to Professor Cyril A. Fox, Jr.

Cy Fox often forgets that he was supposed to be “only” an academic. For three decades, as he helped thousands of law students through the intricacies of the Rule in Shelley’s Case, or watched them calculate a “life in being plus twenty one years” for the Rule Against Perpetuities, he failed to appreciate that law school was supposed to be an amalgam of theory and confusion, not the place for his teaching law students about helping real people solve real life problems

Exploiting transient protein states for the design of small-molecule stabilizers of mutant p53

The destabilizing p53 cancer mutation Y220C creates an extended crevice on the surface of the protein that can be targeted by small-molecule stabilizers. Here, we identify different classes of small molecules that bind to this crevice and determine their binding modes by X-ray crystallography. These structures reveal two major conformational states of the pocket and a cryptic, transiently open hydrophobic subpocket that is modulated by Cys220. In one instance, specifically targeting this transient protein state by a pyrrole moiety resulted in a 40-fold increase in binding affinity. Molecular dynamics simulations showed that both open and closed states of this subsite were populated at comparable frequencies along the trajectories. Our data extend the framework for the design of high-affinity Y220C mutant binders for use in personalized anticancer therapy and, more generally, highlight the importance of implementing protein dynamics and hydration patterns in the drug-discovery process

Role of the Subunits Interactions in the Conformational Transitions in Adult Human Hemoglobin: an Explicit Solvent Molecular Dynamics Study

Hemoglobin exhibits allosteric structural changes upon ligand binding due to the dynamic interactions between the ligand binding sites, the amino acids residues and some other solutes present under physiological conditions. In the present study, the dynamical and quaternary structural changes occurring in two unligated (deoxy-) T structures, and two fully ligated (oxy-) R, R2 structures of adult human hemoglobin were investigated with molecular dynamics. It is shown that, in the sub-microsecond time scale, there is no marked difference in the global dynamics of the amino acids residues in both the oxy- and the deoxy- forms of the individual structures. In addition, the R, R2 are relatively stable and do not present quaternary conformational changes within the time scale of our simulations while the T structure is dynamically more flexible and exhibited the T\rightarrow R quaternary conformational transition, which is propagated by the relative rotation of the residues at the {\alpha}1{\beta}2 and {\alpha}2{\beta}1 interface.Comment: Reprinted (adapted) with permission from J. Phys. Chem. B DOI:10.1021/jp3022908. Copyright (2012) American Chemical Societ

Orbitofrontal cortex, emotional decision-making and response to cognitive behavioural therapy for psychosis

Grey matter volume (GMV) in the orbitofrontal cortex (OFC) may relate to better response to cognitive behavioural therapy for psychosis (CBTp) because of the region's role in emotional decision-making and cognitive flexibility. This study aimed to determine the relation between pre-therapy OFC GMV or asymmetry and CBTp responsiveness and emotional decision-making as measured by the Iowa Gambling Task (IGT). Thirty patients received CBTp + standard care (CBTp+SC; 25 completers) for 6-8 months. All patients (before receiving CBTp) and 25 healthy participants underwent structural magnetic resonance imaging and performed the IGT. Patients' symptoms were assessed before and after therapy. Pre-therapy OFC GMV, measured using a region-of-interest approach, and IGT performance, measured as overall learning, attention to reward, memory for past outcomes and choice consistency, were comparable between patient and healthy groups. In the CBTp+SC group, greater OFC GMV was correlated with positive symptom improvement, specifically hallucinations and persecution. Greater rightward OFC asymmetry correlated with improvement in several negative and general psychopathology symptoms. Greater left OFC GMV was associated with lower IGT attention to reward. The findings suggest that greater OFC volume and rightward asymmetry, which maintain the OFC's function in emotional decision-making and cognitive flexibility, are beneficial for CBTp responsiveness

Hidden bias in the DUD-E dataset leads to misleading performance of deep learning in structure-based virtual screening

Recently much effort has been invested in using convolutional neural network (CNN) models trained on 3D structural images of protein-ligand complexes to distinguish binding from non-binding ligands for virtual screening. However, the dearth of reliable protein-ligand x-ray structures and binding affinity data has required the use of constructed datasets for the training and evaluation of CNN molecular recognition models. Here, we outline various sources of bias in one such widely-used dataset, the Directory of Useful Decoys: Enhanced (DUDE). We have constructed and performed tests to investigate whether CNN models developed using DUD-E are properly learning the underlying physics of molecular recognition, as intended, or are instead learning biases inherent in the dataset itself. We find that superior enrichment efficiency in CNN models can be attributed to the analogue and decoy bias hidden in the DUD-E dataset rather than successful generalization of the pattern of proteinligand interactions. Comparing additional deep learning models trained on PDBbind datasets, we found that their enrichment performances using DUD-E are not superior to the performance of the docking program AutoDock Vina. Together, these results suggest that biases that could be present in constructed datasets should be thoroughly evaluated before applying them to machine learning based methodology development

Hidden bias in the DUD-E dataset leads to misleading performance of deep learning in structure-based virtual screening

Recently much effort has been invested in using convolutional neural network (CNN) models trained on 3D structural images of protein-ligand complexes to distinguish binding from non-binding ligands for virtual screening. However, the dearth of reliable protein-ligand x-ray structures and binding affinity data has required the use of constructed datasets for the training and evaluation of CNN molecular recognition models. Here, we outline various sources of bias in one such widely-used dataset, the Directory of Useful Decoys: Enhanced (DUDE). We have constructed and performed tests to investigate whether CNN models developed using DUD-E are properly learning the underlying physics of molecular recognition, as intended, or are instead learning biases inherent in the dataset itself. We find that superior enrichment efficiency in CNN models can be attributed to the analogue and decoy bias hidden in the DUD-E dataset rather than successful generalization of the pattern of proteinligand interactions. Comparing additional deep learning models trained on PDBbind datasets, we found that their enrichment performances using DUD-E are not superior to the performance of the docking program AutoDock Vina. Together, these results suggest that biases that could be present in constructed datasets should be thoroughly evaluated before applying them to machine learning based methodology development

Mechanism of 150-cavity formation in influenza neuraminidase

The recently discovered 150-cavity in the active site of group-1 influenza A neuraminidase (NA) proteins provides a target for rational structure-based drug development to counter the increasing frequency of antiviral resistance in influenza. Surprisingly, the 2009 H1N1 pandemic virus (09N1) neuramidase was crystalized without the 150-cavity characteristic of group-1 NAs. Here we demonstrate, through a total sum of 1.6 μs of biophysical simulations, that 09N1 NA exists in solution preferentially with an open 150-cavity. Comparison with simulations using avian N1, human N2 and 09N1 with a I149V mutation and an extensive bioinformatics analysis suggests that the conservation of a key salt bridge is crucial in the stabilization of the 150-cavity across both subtypes. This result provides an atomic-level structural understanding of the recent finding that antiviral compounds designed to take advantage of contacts in the 150-cavity can inactivate both 2009 H1N1 pandemic and avian H5N1 viruses

Absolute Single-Molecule Entropies from Quasi-Harmonic Analysis of Microsecond Molecular Dynamics: Correction Terms and Convergence Properties

The convergence properties of the absolute single-molecule configurational entropy and the correction terms used to estimate it are investigated using microsecond molecular dynamics simulation of a peptide test system and an improved methodology. The results are compared with previous applications for systems of diverse chemical nature. It is shown that (i) the effect of anharmonicity is small, (ii) the effect of pairwise correlation is typically large, and (iii) the latter affects to a larger extent the entropy estimate of thermodynamic states characterized by a higher motional correlation. The causes of such deviations from a quasi-harmonic behavior are explained. This improved approach provides entropies also for molecular systems undergoing conformational transitions and characterized by highly frustrated energy surfaces, thus not limited to systems sampling a single quasi-harmonic basin. Overall, this study emphasizes the need for extensive phase-space sampling in order to obtain a reliable estimation of entropic contributions