Visualizing and quantifying molecular goodness-of-fit: small-probe contact dots with explicit hydrogen atoms

Unknotting RNA: A method to resolve computational artifacts

RNA 3D structure prediction often encounters entanglements, computational artifacts that complicate structural models, resulting in their exclusion from further studies despite the potentially accurate prediction of regions outside the entanglement. This study presents a protocol aimed at resolving such issues in RNA models while preserving the overall 3D fold and structural integrity. By employing the SPQR coarse-grained model and short Molecular Dynamics simulations, the protocol imposes energy terms that enable selective modifications to disentangle structures without causing significant distortions. The method was validated on 195 entangled RNA models from CASP15 and RNA-Puzzles, successfully resolving over 70% of interlaces and approximately 40% of lassos, with minimal impact on the original geometry but notable improvement in ClashScore. The efficiency of untangling conformations that are unequivocally classified as artifacts is 81%. Certain cases, particularly those involving dense packing of atoms or complex secondary structures, posed challenges that limited the efficiency of the method. In this paper, we present quantitative results from the application of the protocol and discuss examples of both successfully disentangled and unresolved structures. We show a viable approach for refining models previously deemed unsuitable due to topological artifacts.

Quick-and-Easy Validation of Protein-Ligand Binding Models Using Fragment-Based Semiempirical Quantum Chemistry

Electronic structure calculations in enzymes converge very slowly with respect to the size of the model region that is described using quantum mechanics (QM), requiring hundreds of atoms to obtain converged results and exhibiting substantial sensitivity (at least in smaller models) to which amino acids are included in the QM region. As such, there is considerable interest in developing automated procedures to construct a QM model region based on well-defined criteria. However, testing such procedures is burdensome due to the cost of large-scale electronic structure calculations. Here, we show that semiempirical methods can be used as alternatives to density functional theory (DFT) to assess convergence in sequences of models generated by various automated protocols. The cost of these convergence tests is reduced even further by means of a many-body expansion. We use this approach to examine convergence (with respect to model size) of protein-ligand binding energies. Fragment-based semiempirical calculations afford well-converged interaction energies in a tiny fraction of the cost required for DFT calculations. Two-body interactions between the ligand and single-residue amino acid fragments afford a low-cost way to construct a "QM-informed" enzyme model of reduced size, furnishing an automatable active-site model-building procedure. This provides a streamlined, user-friendly approach for constructing ligand binding-site models that requires neither a priori information nor manual adjustments. Extension to model-building for thermochemical calculations should be straightforward.

Effect of Direct and Water-Mediated Interactions on the Identification of Hotspots in Biomolecular Complexes with Multiple Subsystems

Identification of important residues in biochemical complexes is often a crucial step for many problems in molecular biology and biochemistry. A method is proposed to identify hotspots in biomolecular complexes based on a new metric, derived from networks representing molecular subunits (residues, bridging solvent molecules, ligands etc.) connected by interactions. A singular value decomposition of the weighted adjacency matrix is used to construct a scalar rank for each subunit that reflects its importance in the residue interaction network. This metric is called the singular value centrality. In addition, a new formalism is proposed to account for water-mediated interactions in the ranking of residues. Interactions for a residue network can be provided by various computational methods. In this work interactions are obtained from full quantum-mechanical calculations of protein-protein complexes using the fragment molecular orbital method. The ranking results are shown to be in good agreement with earlier computational and experimental studies. The developed method can be used to gain a deeper insight into the role of subunits in complex biomolecular systems.

Protocol for Designing De Novo Noncanonical Peptide Binders in OSPREY

D-peptides, the mirror image of canonical L-peptides, offer numerous biological advantages that make them effective therapeutics. This article details how to use DexDesign, the newest OSPREY-based algorithm, for designing these D-peptides de novo. OSPREY physics-based models precisely mimic energy-equivariant reflection operations, enabling the generation of D-peptide scaffolds from L-peptide templates. Due to the scarcity of D-peptide:L-protein structural data, DexDesign calls a geometric hashing algorithm, Method of Accelerated Search for Tertiary Ensemble Representatives, as a subroutine to produce a synthetic structural dataset. DexDesign enables mixed-chirality designs with a new user interface and also reduces the conformation and sequence search space using three new design techniques: Minimum Flexible Set, Inverse Alanine Scanning, and K*-based Mutational Scanning.

Building molecular model series from heterogeneous CryoEM structures using Gaussian mixture models and deep neural networks

Cryogenic electron microscopy (CryoEM) produces structures of macromolecules at near-atomic resolution. However, building molecular models with good stereochemical geometry from those structures can be challenging and time-consuming, especially when many structures are obtained from datasets with conformational heterogeneity. Here we present a model refinement protocol that automatically generates series of molecular models from CryoEM datasets, which describe the dynamics of the macromolecular system and have near-perfect geometry scores.

The influence of model building schemes and molecular dynamics sampling on QM-cluster models: the chorismate mutase case study

Most QM-cluster models of enzymes are constructed based on X-ray crystal structures, which limits comparison to in vivo structure and mechanism. The active site of chorismate mutase from Bacillus subtilis and the enzymatic transformation of chorismate to prephenate is used as a case study to guide construction of QM-cluster models built first from the X-ray crystal structure, then from molecular dynamics (MD) simulation snapshots. The Residue Interaction Network ResidUe Selector (RINRUS) software toolkit, developed by our group to simplify and automate the construction of QM-cluster models, is expanded to handle MD to QM-cluster model workflows. Several options, some employing novel topological clustering from residue interaction network (RIN) information, are evaluated for generating conformational clustering from MD simulation. RINRUS then generates a statistical thermodynamic framework for QM-cluster modeling of the chorismate mutase mechanism via refining 250 MD frames with density functional theory (DFT). The 250 QM-cluster models sampled provide a mean ΔG‡ of 10.3 ± 2.6 kcal mol-1 compared to the experimental value of 15.4 kcal mol-1 at 25 °C. While the difference between theory and experiment is consequential, the level of theory used is modest and therefore "chemical" accuracy is unexpected. More important are the comparisons made between QM-cluster models designed from the X-ray crystal structure versus those from MD frames. The large variations in kinetic and thermodynamic properties arise from geometric changes in the ensemble of QM-cluster models, rather from the composition of the QM-cluster models or from the active site-solvent interface. The findings open the way for further quantitative and reproducible calibration in the field of computational enzymology using the model construction framework afforded with the RINRUS software toolkit.

Xplor-NIH: Better parameters and protocols for NMR protein structure determination

The present work describes an update to the protein covalent geometry and atomic radii parameters in the Xplor-NIH biomolecular structure determination package. In combination with an improved treatment of selected non-bonded interactions between atoms three bonds apart, such as those involving methyl hydrogens, and a previously developed term that affects the system's gyration volume, the new parameters are tested using structure calculations on 30 proteins with restraints derived from nuclear magnetic resonance data. Using modern structure validation criteria, including several formally adopted by the Protein Data Bank, and a clear measure of structural accuracy, the results show superior performance relative to previous Xplor-NIH implementations. Additionally, the Xplor-NIH structures compare favorably against originally determined NMR models.

DexDesign: A new OSPREY-based algorithm for designing de novo D-peptide inhibitors

With over 270 unique occurrences in the human genome, peptide-recognizing PDZ domains play a central role in modulating polarization, signaling, and trafficking pathways. Mutations in PDZ domains lead to diseases such as cancer and cystic fibrosis, making PDZ domains attractive targets for therapeutic intervention. D-peptide inhibitors offer unique advantages as therapeutics, including increased metabolic stability and low immunogenicity. Here, we introduce DexDesign, a novel OSPREY-based algorithm for computationally designing de novo D-peptide inhibitors. DexDesign leverages three novel techniques that are broadly applicable to computational protein design: the Minimum Flexible Set, K*-based Mutational Scan, and Inverse Alanine Scan, which enable exponential reductions in the size of the peptide sequence search space. We apply these techniques and DexDesign to generate novel D-peptide inhibitors of two biomedically important PDZ domain targets: CAL and MAST2. We introduce a new framework for analyzing de novo peptides-evaluation along a replication/restitution axis-and apply it to the DexDesign-generated D-peptides. Notably, the peptides we generated are predicted to bind their targets tighter than their targets' endogenous ligands, validating the peptides' potential as lead therapeutic candidates. We provide an implementation of DexDesign in the free and open source computational protein design software OSPREY.

DexDesign: an OSPREY-based algorithm for designing de novo D-peptide inhibitors

With over 270 unique occurrences in the human genome, peptide-recognizing PDZ domains play a central role in modulating polarization, signaling, and trafficking pathways. Mutations in PDZ domains lead to diseases such as cancer and cystic fibrosis, making PDZ domains attractive targets for therapeutic intervention. D-peptide inhibitors offer unique advantages as therapeutics, including increased metabolic stability and low immunogenicity. Here, we introduce DexDesign, a novel OSPREY-based algorithm for computationally designing de novo D-peptide inhibitors. DexDesign leverages three novel techniques that are broadly applicable to computational protein design: the Minimum Flexible Set, K*-based Mutational Scan, and Inverse Alanine Scan. We apply these techniques and DexDesign to generate novel D-peptide inhibitors of two biomedically important PDZ domain targets: CAL and MAST2. We introduce a framework for analyzing de novo peptides-evaluation along a replication/restitution axis-and apply it to the DexDesign-generated D-peptides. Notably, the peptides we generated are predicted to bind their targets tighter than their targets' endogenous ligands, validating the peptides' potential as lead inhibitors. We also provide an implementation of DexDesign in the free and open source computational protein design software OSPREY.

Assessment of three-dimensional RNA structure prediction in CASP15

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty-two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and x-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as noncanonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

The Glycine N-Methyltransferase Case Study: Another Challenge for QM-Cluster Models?

The methyl transfer reaction between SAM and glycine catalyzed by glycine N-methyltransferase (GNMT) was examined using QM-cluster models generated by Residue Interaction Network ResidUe Selector (RINRUS). RINRUS is a Python-based tool that can build QM-cluster models with rules-based processing of the active site residue interaction network. This way of enzyme model-building allows quantitative analysis of residue and fragment contributions to kinetic and thermodynamic properties of the enzyme. Many residue fragments are important for the GNMT catalytic reaction, such as Gly137, Asn138, and Arg175, which interact with the glycine substrate, and Trp30, Asp85, and Tyr242, which interact with the SAM cofactor. Our study shows that active site fragments that interact with the glycine substrate and the SAM cofactor must both be included in the QM-cluster models. Even though the proposed mechanism is a simple one-step reaction, GNMT may be a rather challenging case study for QM-cluster models because convergence in energetics requires models with >350 atoms. "Maximal" QM-cluster models built with either qualitative contact count ranking or quantitative interaction energies from functional group symmetry adapted perturbation theory provide acceptable results. Hence, important residue fragments that contribute to the energetics of the methyl-transfer reaction in GNMT are correctly identified in the RIN. Observations from this work suggest new directions to better establish an effective approach for constructing atomic-level enzyme models.

Assessment of three-dimensional RNA structure prediction in CASP15

The prediction of RNA three-dimensional structures remains an unsolved problem. Here, we report assessments of RNA structure predictions in CASP15, the first CASP exercise that involved RNA structure modeling. Forty two predictor groups submitted models for at least one of twelve RNA-containing targets. These models were evaluated by the RNA-Puzzles organizers and, separately, by a CASP-recruited team using metrics (GDT, lDDT) and approaches (Z-score rankings) initially developed for assessment of proteins and generalized here for RNA assessment. The two assessments independently ranked the same predictor groups as first (AIchemy_RNA2), second (Chen), and third (RNAPolis and GeneSilico, tied); predictions from deep learning approaches were significantly worse than these top ranked groups, which did not use deep learning. Further analyses based on direct comparison of predicted models to cryogenic electron microscopy (cryo-EM) maps and X-ray diffraction data support these rankings. With the exception of two RNA-protein complexes, models submitted by CASP15 groups correctly predicted the global fold of the RNA targets. Comparisons of CASP15 submissions to designed RNA nanostructures as well as molecular replacement trials highlight the potential utility of current RNA modeling approaches for RNA nanotechnology and structural biology, respectively. Nevertheless, challenges remain in modeling fine details such as non-canonical pairs, in ranking among submitted models, and in prediction of multiple structures resolved by cryo-EM or crystallography.

Cross-protein transfer learning substantially improves disease variant prediction

BACKGROUND

Genetic variation in the human genome is a major determinant of individual disease risk, but the vast majority of missense variants have unknown etiological effects. Here, we present a robust learning framework for leveraging saturation mutagenesis experiments to construct accurate computational predictors of proteome-wide missense variant pathogenicity.

RESULTS

We train cross-protein transfer (CPT) models using deep mutational scanning (DMS) data from only five proteins and achieve state-of-the-art performance on clinical variant interpretation for unseen proteins across the human proteome. We also improve predictive accuracy on DMS data from held-out proteins. High sensitivity is crucial for clinical applications and our model CPT-1 particularly excels in this regime. For instance, at 95% sensitivity of detecting human disease variants annotated in ClinVar, CPT-1 improves specificity to 68%, from 27% for ESM-1v and 55% for EVE. Furthermore, for genes not used to train REVEL, a supervised method widely used by clinicians, we show that CPT-1 compares favorably with REVEL. Our framework combines predictive features derived from general protein sequence models, vertebrate sequence alignments, and AlphaFold structures, and it is adaptable to the future inclusion of other sources of information. We find that vertebrate alignments, albeit rather shallow with only 100 genomes, provide a strong signal for variant pathogenicity prediction that is complementary to recent deep learning-based models trained on massive amounts of protein sequence data. We release predictions for all possible missense variants in 90% of human genes.

CONCLUSIONS

Our results demonstrate the utility of mutational scanning data for learning properties of variants that transfer to unseen proteins.

A graph-based machine learning framework identifies critical properties of FVIII that lead to hemophilia A

Introduction: Blood coagulation is an essential process to cease bleeding in humans and other species. This mechanism is characterized by a molecular cascade of more than a dozen components activated after an injury to a blood vessel. In this process, the coagulation factor VIII (FVIII) is a master regulator, enhancing the activity of other components by thousands of times. In this sense, it is unsurprising that even single amino acid substitutions result in hemophilia A (HA)-a disease marked by uncontrolled bleeding and that leaves patients at permanent risk of hemorrhagic complications. Methods: Despite recent advances in the diagnosis and treatment of HA, the precise role of each residue of the FVIII protein remains unclear. In this study, we developed a graph-based machine learning framework that explores in detail the network formed by the residues of the FVIII protein, where each residue is a node, and two nodes are connected if they are in close proximity on the FVIII 3D structure. Results: Using this system, we identified the properties that lead to severe and mild forms of the disease. Finally, in an effort to advance the development of novel recombinant therapeutic FVIII proteins, we adapted our framework to predict the activity and expression of more than 300 in vitro alanine mutations, once more observing a close agreement between the in silico and the in vitro results. Discussion: Together, the results derived from this study demonstrate how graph-based classifiers can leverage the diagnostic and treatment of a rare disease.

Room-temperature crystallography reveals altered binding of small-molecule fragments to PTP1B

Much of our current understanding of how small-molecule ligands interact with proteins stems from X-ray crystal structures determined at cryogenic (cryo) temperature. For proteins alone, room-temperature (RT) crystallography can reveal previously hidden, biologically relevant alternate conformations. However, less is understood about how RT crystallography may impact the conformational landscapes of protein-ligand complexes. Previously, we showed that small-molecule fragments cluster in putative allosteric sites using a cryo crystallographic screen of the therapeutic target PTP1B (Keedy et al., 2018). Here, we have performed two RT crystallographic screens of PTP1B using many of the same fragments, representing the largest RT crystallographic screens of a diverse library of ligands to date, and enabling a direct interrogation of the effect of data collection temperature on protein-ligand interactions. We show that at RT, fewer ligands bind, and often more weakly - but with a variety of temperature-dependent differences, including unique binding poses, changes in solvation, new binding sites, and distinct protein allosteric conformational responses. Overall, this work suggests that the vast body of existing cryo-temperature protein-ligand structures may provide an incomplete picture, and highlights the potential of RT crystallography to help complete this picture by revealing distinct conformational modes of protein-ligand systems. Our results may inspire future use of RT crystallography to interrogate the roles of protein-ligand conformational ensembles in biological function.

Benchmarking applicability of medium-resolution cryo-EM protein structures for structure-based drug design

Cryo-electron microscopy (cryo-EM) is gaining large attention for high-resolution protein structure determination in solutions. However, a very high percentage of cryo-EM structures correspond to resolutions of 3-5 Å, making the structures difficult to be used in in silico drug design. In this study, we analyze how useful cryo-EM protein structures are for in silico drug design by evaluating ligand docking accuracy. From realistic cross-docking scenarios using medium resolution (3-5 Å) cryo-EM structures and a popular docking tool Autodock-Vina, only 20% of docking succeeded, when the success rate doubles in the same kind of cross-docking but using high-resolution (<2 Å) crystal structures instead. We decipher the reason for failures by decomposing the contribution from resolution-dependent and independent factors. The heterogeneity in the protein side-chain and backbone conformations is identified as the major resolution-dependent factor causing docking difficulty from our analysis, while intrinsic receptor flexibility mainly comprises the resolution-independent factor. We demonstrate the flexibility implementation in current ligand docking tools is able to rescue only a portion of failures (10%), and the limited performance was majorly due to potential structural errors than conformational changes. Our work suggests the strong necessity of more robust method developments on ligand docking and EM modeling techniques in order to fully utilize cryo-EM structures for in silico drug design.

A pathogenic human Orai1 mutation unmasks STIM1-independent rapid inactivation of Orai1 channels

Ca²⁺ release-activated Ca²⁺ (CRAC) channels are activated by direct physical interactions between Orai1, the channel protein, and STIM1, the endoplasmic reticulum Ca²⁺ sensor. A hallmark of CRAC channels is fast Ca²⁺-dependent inactivation (CDI) which provides negative feedback to limit Ca²⁺ entry through CRAC channels. Although STIM1 is thought to be essential for CDI, its molecular mechanism remains largely unknown. Here, we examined a poorly understood gain-of-function (GOF) human Orai1 disease mutation, L138F, that causes tubular aggregate myopathy. Through pairwise mutational analysis, we determine that large amino acid substitutions at either L138 or the neighboring T92 locus located on the pore helix evoke highly Ca²⁺-selective currents in the absence of STIM1. We find that the GOF phenotype of the L138 pathogenic mutation arises due to steric clash between L138 and T92. Surprisingly, strongly activating L138 and T92 mutations showed CDI in the absence of STIM1, contradicting prevailing views that STIM1 is required for CDI. CDI of constitutively open T92W and L138F mutants showed enhanced intracellular Ca²⁺ sensitivity, which was normalized by re-adding STIM1 to the cells. Truncation of the Orai1 C-terminus reduced T92W CDI, indicating a key role for the Orai1 C-terminus for CDI. Overall, these results identify the molecular basis of a disease phenotype with broad implications for activation and inactivation of Orai1 channels.

Evaluating the active site-substrate interplay between x-ray crystal structure and molecular dynamics in chorismate mutase

Designing realistic quantum mechanical (QM) models of enzymes is dependent on reliably discerning and modeling residues, solvents, and cofactors important in crafting the active site microenvironment. Interatomic van der Waals contacts have previously demonstrated usefulness toward designing QM-models, but their measured values (and subsequent residue importance rankings) are expected to be influenceable by subtle changes in protein structure. Using chorismate mutase as a case study, this work examines the differences in ligand-residue interatomic contacts between an x-ray crystal structure and structures from a molecular dynamics simulation. Select structures are further analyzed using symmetry adapted perturbation theory to compute ab initio ligand-residue interaction energies. The findings of this study show that ligand-residue interatomic contacts measured for an x-ray crystal structure are not predictive of active site contacts from a sampling of molecular dynamics frames. In addition, the variability in interatomic contacts among structures is not correlated with variability in interaction energies. However, the results spotlight using interaction energies to characterize and rank residue importance in future computational enzymology workflows.

Targeting Protein Pockets with Halogen Bonds: The Role of the Halogen Environment

Protein pockets that form a halogen bond (X-bond) with a halogenated ligand molecule simultaneously form other (mainly hydrophobic) interactions with the halogen atom that can be considered as its "X-bond environment" (XBenv). Most studies in the field have focused on the X-bond, with the properties of the XBenv usually overlooked. In this work, we derived a protocol that evaluates the XBenv strength as a measure of the propensity of a protein pocket to host an X-bond. The charge density-based topological descriptors in combination with machine learning tools were employed to predict formation and strength of the interactions that conform the XBenv as a function of their geometrical parameters. On the basis of these results, we propose that the XBenv can be used as a footprint to judge the chance of a protein pocket to form an X-bond.

Activation of p53: How phosphorylated Ser15 triggers sequential phosphorylation of p53 at Thr18 by CK1δ

The N-terminal transactivation domain (TAD) of p53 is a disordered region with multiple phosphorylation sites. Phosphorylation at Thr18 is crucial for the release of p53 from its negative regulator, MDM2. In stressed cells, CK1δ is responsible for phosphorylating Thr18, but requires Ser15 to be phosphorylated. To understand the mechanistic underpinnings of this sequential phosphorylation, molecular modeling and molecular dynamics simulation studies of these phosphorylation events were carried out. Our models suggest that a positively charged region on CK1δ near the adenosine triphosphate (ATP) binding pocket, which is conserved across species, sequesters the negatively charged pSer15, thereby constraining the positioning of the rest of the peptide, such that the side chain of Thr18 is positioned close to the γ-phosphate of ATP. Furthermore, our studies show that the phosphorylated p53 TAD1 (p53pSer15) peptide binds more strongly to CK1δ than does p53. p53 adopts a helical structure when bound to CK1δ, which is lost upon phosphorylation at Ser15, thus gaining higher flexibility and ability to morph into the binding site. We propose that upon phosphorylation at Ser15 the p53 TAD1 peptide binds to CK1δ through an electrostatically driven induced fit mechanism resulting in a flanking fuzzy complex.

Metabolic and Structural Insights into Hydrogen Sulfide Mis-Regulation in Enterococcus faecalis

Hydrogen sulfide (H2S) is implicated as a cytoprotective agent that bacteria employ in response to host-induced stressors, such as oxidative stress and antibiotics. The physiological benefits often attributed to H2S, however, are likely a result of downstream, more oxidized forms of sulfur, collectively termed reactive sulfur species (RSS) and including the organic persulfide (RSSH). Here, we investigated the metabolic response of the commensal gut microorganism Enterococcus faecalis to exogenous Na2S as a proxy for H2S/RSS toxicity. We found that exogenous sulfide increases protein abundance for enzymes responsible for the biosynthesis of coenzyme A (CoA). Proteome S-sulfuration (persulfidation), a posttranslational modification implicated in H2S signal transduction, is also widespread in this organism and is significantly elevated by exogenous sulfide in CstR, the RSS sensor, coenzyme A persulfide (CoASSH) reductase (CoAPR) and enzymes associated with de novo fatty acid biosynthesis and acetyl-CoA synthesis. Exogenous sulfide significantly impacts the speciation of fatty acids as well as cellular concentrations of acetyl-CoA, suggesting that protein persulfidation may impact flux through these pathways. Indeed, CoASSH is an inhibitor of E. faecalis phosphotransacetylase (Pta), suggesting that an important metabolic consequence of increased levels of H2S/RSS may be over-persulfidation of this key metabolite, which, in turn, inhibits CoA and acyl-CoA-utilizing enzymes. Our 2.05 Å crystallographic structure of CoA-bound CoAPR provides new structural insights into CoASSH clearance in E. faecalis.

Mebendazole's Conformational Space and Its Predicted Binding to Human Heat-Shock Protein 90

Recent experimental evidence suggests that mebendazole, a popular antiparasitic drug, binds to heat shock protein 90 (Hsp90) and inhibits acute myeloid leukemia cell growth. In this study we use quantum mechanics (QM), molecular similarity, and molecular dynamics (MD) calculations to predict possible binding poses of mebendazole to the adenosine triphosphate (ATP) binding site of Hsp90. Extensive conformational searches and minimization of the five mebendazole tautomers using the MP2/aug-cc-pVTZ theory level resulted in 152 minima. Mebendazole-Hsp90 complex models were subsequently created using the QM optimized conformations and protein coordinates obtained from experimental crystal structures that were chosen through similarity calculations. Nine different poses were identified from a total of 600 ns of explicit solvent, all-atom MD simulations using two different force fields. All simulations support the hypothesis that mebendazole is able to bind to the ATP binding site of Hsp90.

A Machine Learning Framework Predicts the Clinical Severity of Hemophilia B Caused by Point-Mutations

Blood coagulation is a vital physiological mechanism to stop blood loss following an injury to a blood vessel. This process starts immediately upon damage to the endothelium lining a blood vessel, and results in the formation of a platelet plug that closes the site of injury. In this repair operation, an essential component is the coagulation factor IX (FIX), a serine protease encoded by the F9 gene and whose deficiency causes hemophilia B. If not treated by prophylaxis or gene therapy, patients with this condition are at risk of life-threatening bleeding episodes. In this sense, a deep understanding of the FIX protein and its activated form (FIXa) is essential to develop efficient therapeutics. In this study, we used well-studied structural analysis techniques to create a residue interaction network of the FIXa protein. Here, the nodes are the amino acids of FIXa, and two nodes are connected by an edge if the two residues are in close proximity in the FIXa 3D structure. This representation accurately captured fundamental properties of each amino acid of the FIXa structure, as we found by validating our findings against hundreds of clinical reports about the severity of HB. Finally, we established a machine learning framework named HemB-Class to predict the effect of mutations of all FIXa residues to all other amino acids and used it to disambiguate several conflicting medical reports. Together, these methods provide a comprehensive map of the FIXa protein architecture and establish a robust platform for the rational design of FIX therapeutics.

A Case Study of the Glycoside Hydrolase Enzyme Mechanism Using an Automated QM-Cluster Model Building Toolkit

Glycoside hydrolase enzymes are important for hydrolyzing the β-1,4 glycosidic bond in polysaccharides for deconstruction of carbohydrates. The two-step retaining reaction mechanism of Glycoside Hydrolase Family 7 (GH7) was explored with different sized QM-cluster models built by the Residue Interaction Network ResidUe Selector (RINRUS) software using both the wild-type protein and its E217Q mutant. The first step is the glycosylation, in which the acidic residue 217 donates a proton to the glycosidic oxygen leading to bond cleavage. In the subsequent deglycosylation step, one water molecule migrates into the active site and attacks the anomeric carbon. Residue interaction-based QM-cluster models lead to reliable structural and energetic results for proposed glycoside hydrolase mechanisms. The free energies of activation for glycosylation in the largest QM-cluster models were predicted to be 19.5 and 31.4 kcal mol⁻¹ for the wild-type protein and its E217Q mutant, which agree with experimental trends that mutation of the acidic residue Glu217 to Gln will slow down the reaction; and are higher in free energy than the deglycosylation transition states (13.8 and 25.5 kcal mol⁻¹ for the wild-type protein and its mutant, respectively). For the mutated protein, glycosylation led to a low-energy product. This thermodynamic sink may correspond to the intermediate state which was isolated in the X-ray crystal structure. Hence, the glycosylation is validated to be the rate-limiting step in both the wild-type and mutated enzyme.

Interface refinement of low- to medium-resolution Cryo-EM complexes using HADDOCK2.4

A wide range of cellular processes requires the formation of multimeric protein complexes. The rise of cryo-electron microscopy (cryo-EM) has enabled the structural characterization of these protein assemblies. The density maps produced can, however, still suffer from limited resolution, impeding the process of resolving structures at atomic resolution. In order to solve this issue, monomers can be fitted into low- to medium-resolution maps. Unfortunately, the models produced frequently contain atomic clashes at the protein-protein interfaces (PPIs), as intermolecular interactions are typically not considered during monomer fitting. Here, we present a refinement approach based on HADDOCK2.4 to remove intermolecular clashes and optimize PPIs. A dataset of 14 cryo-EM complexes was used to test eight protocols. The best-performing protocol, consisting of a semi-flexible simulated annealing refinement with centroid restraints on the monomers, was able to decrease intermolecular atomic clashes by 98% without significantly deteriorating the quality of the cryo-EM density fit.

Specific Human ATR and ATM Inhibitors Modulate Single Strand DNA Formation in Leishmania major Exposed to Oxidative Agent

Leishmania parasites are the causative agents of a group of neglected tropical diseases known as leishmaniasis. The molecular mechanisms employed by these parasites to adapt to the adverse conditions found in their hosts are not yet completely understood. DNA repair pathways can be used by Leishmania to enable survival in the interior of macrophages, where the parasite is constantly exposed to oxygen reactive species. In higher eukaryotes, DNA repair pathways are coordinated by the central protein kinases ataxia telangiectasia mutated (ATM) and ataxia telangiectasia and Rad3 related (ATR). The enzyme Exonuclease-1 (EXO1) plays important roles in DNA replication, repair, and recombination, and it can be regulated by ATM- and ATR-mediated signaling pathways. In this study, the DNA damage response pathways in promastigote forms of L. major were investigated using bioinformatics tools, exposure of lineages to oxidizing agents and radiation damage, treatment of cells with ATM and ATR inhibitors, and flow cytometry analysis. We demonstrated high structural and important residue conservation for the catalytic activity of the putative LmjEXO1. The overexpression of putative LmjEXO1 made L. major cells more susceptible to genotoxic damage, most likely due to the nuclease activity of this enzyme and the occurrence of hyper-resection of DNA strands. These cells could be rescued by the addition of caffeine or a selective ATM inhibitor. In contrast, ATR-specific inhibition made the control cells more susceptible to oxidative damage in an LmjEXO1 overexpression-like manner. We demonstrated that ATR-specific inhibition results in the formation of extended single-stranded DNA, most likely due to EXO1 nucleasic activity. Antagonistically, ATM inhibition prevented single-strand DNA formation, which could explain the survival phenotype of lineages overexpressing LmjEXO1. These results suggest that an ATM homolog in Leishmania could act to promote end resection by putative LmjEXO1, and an ATR homologue could prevent hyper-resection, ensuring adequate repair of the parasite DNA.

The importance of residue-level filtering and the Top2018 best-parts dataset of high-quality protein residues

We have curated a high-quality, "best-parts" reference dataset of about 3 million protein residues in about 15,000 PDB-format coordinate files, each containing only residues with good electron density support for a physically acceptable model conformation. The resulting prefiltered data typically contain the entire core of each chain, in quite long continuous fragments. Each reference file is a single protein chain, and the total set of files were selected for low redundancy, high resolution, good MolProbity score, and other chain-level criteria. Then each residue was critically tested for adequate local map quality to firmly support its conformation, which must also be free of serious clashes or covalent-geometry outliers. The resulting Top2018 prefiltered datasets have been released on the Zenodo online web service and are freely available for all uses under a Creative Commons license. Currently, one dataset is residue filtered on main chain plus Cβ atoms, and a second dataset is full-residue filtered; each is available at four different sequence-identity levels. Here, we illustrate both statistics and examples that show the beneficial consequences of residue-level filtering. That process is necessary because even the best of structures contain a few highly disordered local regions with poor density and low-confidence conformations that should not be included in reference data. Therefore, the open distribution of these very large, prefiltered reference datasets constitutes a notable advance for structural bioinformatics and the fields that depend upon it.

Cheminformatic quantum mechanical enzyme model design: A catechol-O-methyltransferase case study

To accurately simulate the inner workings of an enzyme active site with quantum mechanics (QM), not only must the reactive species be included in the model but also important surrounding residues, solvent, or coenzymes involved in crafting the microenvironment. Our lab has been developing the Residue Interaction Network Residue Selector (RINRUS) toolkit to utilize interatomic contact network information for automated, rational residue selection and QM-cluster model generation. Starting from an x-ray crystal structure of catechol-O-methyltransferase, RINRUS was used to construct a series of QM-cluster models. The reactant, product, and transition state of the methyl transfer reaction were computed for a total of 550 models, and the resulting free energies of activation and reaction were used to evaluate model convergence. RINRUS-designed models with only 200-300 atoms are shown to converge. RINRUS will serve as a cornerstone for improved and automated cheminformatics-based enzyme model design.

RNA 3D Structure Prediction Using Coarse-Grained Models

The three-dimensional (3D) structures of Ribonucleic acid (RNA) molecules are essential to understanding their various and important biological functions. However, experimental determination of the atomic structures is laborious and technically difficult. The large gap between the number of sequences and the experimentally determined structures enables the thriving development of computational approaches to modeling RNAs. However, computational methods based on all-atom simulations are intractable for large RNA systems, which demand long time simulations. Facing such a challenge, many coarse-grained (CG) models have been developed. Here, we provide a review of CG models for modeling RNA 3D structures, compare the performance of the different models, and offer insights into potential future developments.

The Serological Cross-Detection of Bat-Borne Hantaviruses: A Valid Strategy or Taking Chances?

Bats are hosts of a range of viruses, and their great diversity and unique characteristics that distinguish them from all other mammals have been related to the maintenance, evolution, and dissemination of these pathogens. Recently, very divergent hantaviruses have been discovered in distinct species of bats worldwide, but their association with human disease remains unclear. Considering the low success rates of detecting hantavirus RNA in bat tissues and that to date no hantaviruses have been isolated from bat samples, immunodiagnostic tools could be very helpful to understand pathogenesis, epidemiology, and geographic range of bat-borne hantaviruses. In this sense, we aimed to identify in silico immunogenic B-cell epitopes present on bat-borne hantaviruses nucleoprotein (NP) and verify if they are conserved among them and other selected members of Mammantavirinae, using a combination of (the three most used) different prediction algorithms, ELLIPRO, Discotope 2.0, and PEPITO server. To support our data, we in silico modeled 3D structures of NPs from representative members of bat-borne hantaviruses, using comparative and ab initio methods due to the absence of crystallographic structures of studied proteins or similar models in the Protein Data Bank. Our analysis demonstrated the antigenic complexity of the bat-borne hantaviruses group, showing a low sequence conservation of epitopes among members of its own group and a minor conservation degree in comparison to Orthohantavirus, with a recognized importance to public health. Our data suggest that the use of recombinant rodent-borne hantavirus NPs to cross-detect antibodies against bat- or shrew-borne viruses could underestimate the real impact of this virus in nature.

Protein residue network analysis reveals fundamental properties of the human coagulation factor VIII

Hemophilia A is an X-linked inherited blood coagulation disorder caused by the production and circulation of defective coagulation factor VIII protein. People living with this condition receive either prophylaxis or on-demand treatment, and approximately 30% of patients develop inhibitor antibodies, a serious complication that limits treatment options. Although previous studies performed targeted mutations to identify important residues of FVIII, a detailed understanding of the role of each amino acid and their neighboring residues is still lacking. Here, we addressed this issue by creating a residue interaction network (RIN) where the nodes are the FVIII residues, and two nodes are connected if their corresponding residues are in close proximity in the FVIII protein structure. We studied the characteristics of all residues in this network and found important properties related to disease severity, interaction to other proteins and structural stability. Importantly, we found that the RIN-derived properties were in close agreement with in vitro and clinical reports, corroborating the observation that the patterns derived from this detailed map of the FVIII protein architecture accurately capture the biological properties of FVIII.

VoroContacts: a tool for the analysis of interatomic contacts in macromolecular structures

SUMMARY

VoroContacts is a versatile tool for computing and analyzing contact surface areas (CSAs) and solvent accessible surface areas (SASAs) for 3 D structures of proteins, nucleic acids and their complexes at the atomic resolution. CSAs and SASAs are derived using Voronoi tessellation of 3 D structure, represented as a collection of atomic balls. VoroContacts web server features a highly configurable query interface, which enables on-the-fly analysis of contacts for selected set of atoms and allows filtering interatomic contacts by their type, surface areas, distance between contacting atoms and sequence separation between contacting residues. The VoroContacts functionality is also implemented as part of the standalone Voronota package, enabling batch processing.

AVAILABILITY AND IMPLEMENTATION

https://bioinformatics.lt/wtsam/vorocontacts.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A novel measure to analyze protein structures: Aspect ratio in protein alpha shapes

Proteins' three-dimensional (3D) structures are analyzed traditionally using geometric descriptors such as torsional angles and inter-atomic distances. In this study a measure that is borrowed from computational geometry, aspect ratio of each tetrahedron in alpha shapes of proteins, is utilized. This geometric descriptor differentiates alpha and beta structural classes of proteins when combined with principal components analysis. The method converts the structures of individual proteins, 3D coordinates of the atoms, to points on a plane. It has a high degree of accuracy in differentiating R and T structures of hemoglobin. Therefore, it is anticipated that the geometric measure can be used successfully in a method that is extended to solve classification problems in machine learning.

Novel mutations in malonyl-CoA-acyl carrier protein transacylase provoke autosomal recessive optic neuropathy

Inherited optic neuropathies are rare eye diseases of optic nerve dysfunction that present in various genetic forms. Previously, mutation in three genes encoding mitochondrial proteins has been implicated in autosomal recessive forms of optic atrophy that involve progressive degeneration of optic nerve and retinal ganglion cells (RGC). Using whole exome analysis, a novel double homozygous mutation p.L81R and pR212W in malonyl CoA-acyl carrier protein transacylase (MCAT), a mitochondrial protein involved in fatty acid biosynthesis, has now been identified as responsible for an autosomal recessive optic neuropathy from a Chinese consanguineous family. MCAT is expressed in RGC that are rich in mitochondria. The disease variants lead to structurally unstable MCAT protein with significantly reduced intracellular expression. RGC-specific knockdown of Mcat in mice, lead to an attenuated retinal neurofiber layer, that resembles the phenotype of optic neuropathy. These results indicated that MCAT plays an essential role in mitochondrial function and maintenance of RGC axons, while novel MCAT p.L81R and p.R212W mutations can lead to optic neuropathy.

The missing link: covalent linkages in structural models

Covalent linkages between constituent blocks of macromolecules and ligands have been subject to inconsistent treatment during the model-building, refinement and deposition process. This may stem from a number of sources, including difficulties with initially detecting the covalent linkage, identifying the correct chemistry, obtaining an appropriate restraint dictionary and ensuring its correct application. The analysis presented herein assesses the extent of problems involving covalent linkages in the Protein Data Bank (PDB). Not only will this facilitate the remediation of existing models, but also, more importantly, it will inform and thus improve the quality of future linkages. By considering linkages of known type in the CCP4 Monomer Library (CCP4-ML), failure to model a covalent linkage is identified to result in inaccurate (systematically longer) interatomic distances. Scanning the PDB for proximal atom pairs that do not have a corresponding type in the CCP4-ML reveals a large number of commonly occurring types of unannotated potential linkages; in general, these may or may not be covalently linked. Manual consideration of the most commonly occurring cases identifies a number of genuine classes of covalent linkages. The recent expansion of the CCP4-ML is discussed, which has involved the addition of over 16 000 and the replacement of over 11 000 component dictionaries using AceDRG. As part of this effort, the CCP4-ML has also been extended using AceDRG link dictionaries for the aforementioned linkage types identified in this analysis. This will facilitate the identification of such linkage types in future modelling efforts, whilst concurrently easing the process involved in their application. The need for a universal standard for maintaining link records corresponding to covalent linkages, and references to the associated dictionaries used during modelling and refinement, following deposition to the PDB is emphasized. The importance of correctly modelling covalent linkages is demonstrated using a case study, which involves the covalent linkage of an inhibitor to the main protease in various viral species, including SARS-CoV-2. This example demonstrates the importance of properly modelling covalent linkages using a comprehensive restraint dictionary, as opposed to just using a single interatomic distance restraint or failing to model the covalent linkage at all.

Protein-nanoparticle interactions and a new insight

The study of protein-nanoparticle interactions provides knowledge about the bio-reactivity of nanoparticles, and creates a database of nanoparticles for applications in nanomedicine, nanodiagnosis, and nanotherapy. The problem arises when nanoparticles come in contact with physiological fluids such as plasma or serum, wherein they interact with the proteins (or other biomolecules). This interaction leads to the coating of proteins on the nanoparticle surface, mostly due to the electrostatic interaction, called 'corona'. These proteins are usually partially unfolded. The protein corona can deter nanoparticles from their targeted functionalities, such as drug/DNA delivery at the site and fluorescence tagging of diseased tissues. The protein corona also has many repercussions on cellular intake, inflammation, accumulation, degradation, and clearance of the nanoparticles from the body depending on the exposed part of the proteins. Hence, the protein-nanoparticle interaction and the configuration of the bound-proteins on the nanosurface need thorough investigation and understanding. Several techniques such as DLS and zeta potential measurement, UV-vis spectroscopy, fluorescence spectroscopy, circular dichroism, FTIR, and DSC provide valuable information in the protein-nanoparticle interaction study. Besides, theoretical simulations also provide additional understanding. Despite a lot of research publications, the fundamental question remained unresolved. Can we aim for the application of functional nanoparticles in medicine? A new insight, given by us, in this article assumes a reasonable solution to this crucial question.

QM-Cluster Model Study of the Guaiacol Hydrogen Atom Transfer and Oxygen Rebound with Cytochrome P450 Enzyme GcoA

The key step of the O-demethylation of guaiacol by GcoA of the cytochrome P450-reductase pair was studied with DFT using two 10-residue and three 15-residue QM-cluster models. For each model, two reaction pathways were examined, beginning with a different guaiacol orientation. Based on this study, His354, Phe349, Glu249, and Pro250 residues were found to be important for keeping the heme in a planar geometry throughout the reaction. Val241 and Gly245 residues were needed in the QM-cluster models to provide the hydrophobic pocket for an appropriate guaiacol pose in the reaction. The aromatic triad Phe75, Phe169, and Phe395 may be necessary to facilitate guaiacol migrating into the enzyme active site, but it does not qualitatively affect kinetics and thermodynamics of the proposed mechanism. All QM-cluster models created by RINRUS agree very well with previous experimental work. This study provides details for better understanding enzymatic O-demethylation of lignins to form catechol derivatives by GcoA.

The critical role of QM/MM X-ray refinement and accurate tautomer/protomer determination in structure-based drug design

Conventional protein:ligand crystallographic refinement uses stereochemistry restraints coupled with a rudimentary energy functional to ensure the correct geometry of the model of the macromolecule-along with any bound ligand(s)-within the context of the experimental, X-ray density. These methods generally lack explicit terms for electrostatics, polarization, dispersion, hydrogen bonds, and other key interactions, and instead they use pre-determined parameters (e.g. bond lengths, angles, and torsions) to drive structural refinement. In order to address this deficiency and obtain a more complete and ultimately more accurate structure, we have developed an automated approach for macromolecular refinement based on a two layer, QM/MM (ONIOM) scheme as implemented within our DivCon Discovery Suite and "plugged in" to two mainstream crystallographic packages: PHENIX and BUSTER. This implementation is able to use one or more region layer(s), which is(are) characterized using linear-scaling, semi-empirical quantum mechanics, followed by a system layer which includes the balance of the model and which is described using a molecular mechanics functional. In this work, we applied our Phenix/DivCon refinement method-coupled with our XModeScore method for experimental tautomer/protomer state determination-to the characterization of structure sets relevant to structure-based drug design (SBDD). We then use these newly refined structures to show the impact of QM/MM X-ray refined structure on our understanding of function by exploring the influence of these improved structures on protein:ligand binding affinity prediction (and we likewise show how we use post-refinement scoring outliers to inform subsequent X-ray crystallographic efforts). Through this endeavor, we demonstrate a computational chemistry ↔ structural biology (X-ray crystallography) "feedback loop" which has utility in industrial and academic pharmaceutical research as well as other allied fields.

Improving SARS-CoV-2 structures: Peer review by early coordinate release

This work builds upon the record-breaking speed and generous immediate release of new experimental three-dimensional structures of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins and complexes, which are crucial to downstream vaccine and drug development. We have surveyed those structures to catch the occasional errors that could be significant for those important uses and for which we were able to provide demonstrably higher-accuracy corrections. This process relied on new validation and correction methods such as CaBLAM and ISOLDE, which are not yet in routine use. We found such important and correctable problems in seven early SARS-CoV-2 structures. Two of the structures were soon superseded by new higher-resolution data, confirming our proposed changes. For the other five, we emailed the depositors a documented and illustrated report and encouraged them to make the model corrections themselves and use the new option at the worldwide Protein Data Bank for depositors to re-version their coordinates without changing the Protein Data Bank code. This quickly and easily makes the better-accuracy coordinates available to anyone who examines or downloads their structure, even before formal publication. The changes have involved sequence misalignments, incorrect RNA conformations near a bound inhibitor, incorrect metal ligands, and cis-trans or peptide flips that prevent good contact at interaction sites. These improvements have propagated into nearly all related structures done afterward. This process constitutes a new form of highly rigorous peer review, which is actually faster and more strict than standard publication review because it has access to coordinates and maps; journal peer review would also be strengthened by such access.

Cryo-EM structure of amyloid fibrils formed by the entire low complexity domain of TDP-43

Amyotrophic lateral sclerosis and several other neurodegenerative diseases are associated with brain deposits of amyloid-like aggregates formed by the C-terminal fragments of TDP-43 that contain the low complexity domain of the protein. Here, we report the cryo-EM structure of amyloid formed from the entire TDP-43 low complexity domain in vitro at pH 4. This structure reveals single protofilament fibrils containing a large (139-residue), tightly packed core. While the C-terminal part of this core region is largely planar and characterized by a small proportion of hydrophobic amino acids, the N-terminal region contains numerous hydrophobic residues and has a non-planar backbone conformation, resulting in rugged surfaces of fibril ends. The structural features found in these fibrils differ from those previously found for fibrils generated from short protein fragments. The present atomic model for TDP-43 LCD fibrils provides insight into potential structural perturbations caused by phosphorylation and disease-related mutations.

How does α1Histidine102 affect the binding of modulators to α1β2γ2 GABAA receptors? molecular insights from in silico experiments

The activation of GABAA receptors by the neurotransmitter gamma-aminobutyric acid mediates the rapid inhibition response in the central nervous system of mammals. Many neurological and mental health disorders arise from alterations in the structure or function of these pentameric ion channels. GABAA receptors are targets for numerous drugs, including benzodiazepines, which bind to α1β2γ2 GABAA receptors with high affinity to a site in the extracellular domain, between subunits α1 and γ2. It has been established experimentally that the binding of these drugs depends on the presence of one particular amino acid in the α1 subunit: histidine 102. However, the specific role it plays in the intermolecular interaction has not been elucidated. In this work, we applied in silico methods to understand whether certain protonation and rotamer states of α1His102 facilitate the binding of modulators. We analysed diazepam binding, a benzodiazepine, and the antagonist flumazenil to the GABAA receptor using molecular dynamics simulations and adaptive biasing force simulations. The binding free energy follows changes in the protonation state for both ligands, and rotameric states of α1His102 were specific for the different compounds, suggesting distinct preferences for positive allosteric modulators and antagonists. Moreover, in the presence of diazepam and favoured by a neutral tautomer, we identified a water molecule that links loops A, B, and C and may be relevant to the modulation mechanism.

Computationally-obtained structural insights into the molecular interactions between Pidilizumab and binding partners DLL1 and PD-1

Pidilizumab is a monoclonal antibody tested against several types of malignancies, such as lymphoma and metastatic melanoma, showing promising results. In 2016, the FDA put Pidilizumab's clinical studies on partial hold due to emerging evidence pointing to the antibody target uncertainty. Although initial studies indicated an interaction with the PD-1 checkpoint receptor, recent updates assert that Pidilizumab binds primarily to Notch ligand DLL1. However, a detailed description of which interactions coordinate antibody-antigen complex formation is lacking. Therefore, this study uses computational tools to identify molecular interactions between Pidilizumab and its reported targets PD-1 and DLL1. A docking methodology was validated and applied to determine the binding modes between modeled Pidilizumab scFvs and the two antigens. We used Molecular Dynamics (MD) simulations to verify the complexes' stability and submitted the resulting trajectory files to MM/PBSA and Principal Component Analysis. A set of different prediction tools determined scFv interface hot-spots. Whereas docking and MD simulations revealed that the antibody fragments do not interact straightforwardly with PD-1, ten scFv hot-spots, including Met93 and Leu112, mediated the interaction with the DLL1 C2 domain. The interaction triggered a conformational selection-like effect on DLL1, allowing new hydrogen bonds on the β3-β4 interface loop. The unprecedented structural data on Pidilizumab's interactions provided novel evidence that its legitimate target is the DLL1 protein and offered structural insight on how these molecules interact, shedding light on the pathways that could be affected by the use of this essential immunobiological.Communicated by Ramaswamy H. Sarma.

Observation of an Unusually Large IR Red-Shift in an Unconventional S-H···S Hydrogen-Bond

The S-H···S non-covalent interaction is generally known as an extremely unconventional weak hydrogen-bond in the literature. The present gas-phase spectroscopic investigation shows that the S-H···S hydrogen-bond can be as strong as any conventional hydrogen-bond in terms of the IR red-shift in the stretching frequency of the hydrogen-bond donor group. Herein, the strength of the S-H···S hydrogen-bond has been determined by measuring the red-shift (∼150 cm^-1) of the S-H stretching frequency in a model complex of 2-chlorothiophenol and dimethyl sulfide using isolated gas-phase IR spectroscopy coupled with quantum chemistry calculations. The observation of an unusually large IR red-shift in the S-H···S hydrogen-bond is explained in terms of the presence of a significant amount of charge-transfer interactions in addition to the usual electrostatic interactions. The existence of ∼750 S-H···S interactions between the cysteine and methionine residues in 642 protein structures determined from an extensive Protein Data Bank analysis also indicates that this interaction is important for the structures of proteins.

Seeing the PDB

Ever since the first structures of proteins were determined in the 1960s, structural biologists have required methods to visualize biomolecular structures, both as an essential tool for their research and also to promote 3D comprehension of structural results by a wide audience of researchers, students, and the general public. In this review to celebrate the 50th anniversary of the Protein Data Bank, we present our own experiences in developing and applying methods of visualization and analysis to the ever-expanding archive of protein and nucleic acid structures in the worldwide Protein Data Bank. Across that timespan, Jane and David Richardson have concentrated on the organization inside and between the macromolecules, with ribbons to show the overall backbone "fold" and contact dots to show how the all-atom details fit together locally. David Goodsell has explored surface-based representations to present and explore biological subjects that range from molecules to cells. This review concludes with some ideas about the current challenges being addressed by the field of biomolecular visualization.

CERES: a cryo-EM re-refinement system for continuous improvement of deposited models

The field of electron cryomicroscopy (cryo-EM) has advanced quickly in recent years as the result of numerous technological and methodological developments. This has led to an increase in the number of atomic structures determined using this method. Recently, several tools for the analysis of cryo-EM data and models have been developed within the Phenix software package, such as phenix.real_space_refine for the refinement of atomic models against real-space maps. Also, new validation metrics have been developed for low-resolution cryo-EM models. To understand the quality of deposited cryo-EM structures and how they might be improved, models deposited in the Protein Data Bank that have map resolutions of better than 5 Å were automatically re-refined using current versions of Phenix tools. The results are available on a publicly accessible web page (https://cci.lbl.gov/ceres). The implementation of a Cryo-EM Re-refinement System (CERES) for the improvement of models deposited in the wwPDB, and the results of the re-refinements, are described. Based on these results, contents are proposed for a `cryo-EM Table 1', which summarizes experimental details and validation metrics in a similar way to `Table 1' in crystallography. The consistent use of robust metrics for the evaluation of cryo-EM models and data should accompany every structure deposition and be reported in scientific publications.

Exploring the candidates for a new protein folding - cross-α amyloid - in available protein databases

Amyloid fibrils are formed from the assembly of soluble proteins and are responsible for many diseases. They are known to have a cross-β structure, where the fibril runs perpendicular to the β-sheets. A new type of tertiary structure formed by the aggregation of peptides in their α-helical form, in naturally occurring as well as synthetic peptides, termed cross-α amyloid has been reported recently. We have studied the interactions responsible for the formation of these cross-α amyloids and proposed a model to determine the peptides that could form these structures. Eight such peptides obtained using the model have been shown to form a cross-α structure using molecular dynamics simulations. The formation of a cross-α structure from eight copies of a randomly chosen peptide and its stability over a microsecond simulation have been demonstrated. A software named Cross-Alpha-Det has been developed that can determine whether a protein can form a cross-α structure from its secondary structure.

An allosteric hot spot in the tandem-SH2 domain of ZAP-70 regulates T-cell signaling

T-cell receptor (TCR) signaling is initiated by recruiting ZAP-70 to the cytosolic part of TCR. ZAP-70, a non-receptor tyrosine kinase, is composed of an N-terminal tandem SH2 (tSH2) domain connected to the C-terminal kinase domain. The ZAP-70 is recruited to the membrane through binding of tSH2 domain and the doubly phosphorylated ITAM motifs of CD3 chains in the TCR complex. Our results show that the tSH2 domain undergoes a biphasic structural transition while binding to the doubly phosphorylated ITAM-ζ1 peptide. The C-terminal SH2 domain binds first to the phosphotyrosine residue of ITAM peptide to form an encounter complex leading to subsequent binding of second phosphotyrosine residue to the N-SH2 domain. We decipher a network of noncovalent interactions that allosterically couple the two SH2 domains during binding to doubly phosphorylated ITAMs. Mutation in the allosteric network residues, for example, W165C, uncouples the formation of encounter complex to the subsequent ITAM binding thus explaining the altered recruitment of ZAP-70 to the plasma membrane causing autoimmune arthritis in mice. The proposed mechanism of allosteric coupling is unique to ZAP-70, which is fundamentally different from Syk, a close homolog of ZAP-70 expressed in B-cells.

Potential therapeutic target identification in the novel 2019 coronavirus: insight from homology modeling and blind docking study

Background

The 2019-nCoV which is regarded as a novel coronavirus is a positive-sense single-stranded RNA virus. It is infectious to humans and is the cause of the ongoing coronavirus outbreak which has elicited an emergency in public health and a call for immediate international concern has been linked to it. The coronavirus main proteinase which is also known as the 3C-like protease (3CLpro) is a very important protein in all coronaviruses for the role it plays in the replication of the virus and the proteolytic processing of the viral polyproteins. The resultant cytotoxic effect which is a product of consistent viral replication and proteolytic processing of polyproteins can be greatly reduced through the inhibition of the viral main proteinase activities. This makes the 3C-like protease of the coronavirus a potential and promising target for therapeutic agents against the viral infection.

Results

This study describes the detailed computational process by which the 2019-nCoV main proteinase coding sequence was mapped out from the viral full genome, translated and the resultant amino acid sequence used in modeling the protein 3D structure. Comparative physiochemical studies were carried out on the resultant target protein and its template while selected HIV protease inhibitors were docked against the protein binding sites which contained no co-crystallized ligand.

Conclusion

In line with results from this study which has shown great consistency with other scientific findings on coronaviruses, we recommend the administration of the selected HIV protease inhibitors as first-line therapeutic agents for the treatment of the current coronavirus epidemic.

A defined structural unit enables de novo design of small-molecule-binding proteins

The de novo design of proteins that bind highly functionalized small molecules represents a great challenge. To enable computational design of binders, we developed a unit of protein structure-a van der Mer (vdM)-that maps the backbone of each amino acid to statistically preferred positions of interacting chemical groups. Using vdMs, we designed six de novo proteins to bind the drug apixaban; two bound with low and submicromolar affinity. X-ray crystallography and mutagenesis confirmed a structure with a precisely designed cavity that forms favorable interactions in the drug-protein complex. vdMs may enable design of functional proteins for applications in sensing, medicine, and catalysis.

Using an integrative machine learning approach utilising homology modelling to clinically interpret genetic variants: CACNA1F as an exemplar

Advances in DNA sequencing technologies have revolutionised rare disease diagnostics and have led to a dramatic increase in the volume of available genomic data. A key challenge that needs to be overcome to realise the full potential of these technologies is that of precisely predicting the effect of genetic variants on molecular and organismal phenotypes. Notably, despite recent progress, there is still a lack of robust in silico tools that accurately assign clinical significance to variants. Genetic alterations in the CACNA1F gene are the commonest cause of X-linked incomplete Congenital Stationary Night Blindness (iCSNB), a condition associated with non-progressive visual impairment. We combined genetic and homology modelling data to produce CACNA1F-vp, an in silico model that differentiates disease-implicated from benign missense CACNA1F changes. CACNA1F-vp predicts variant effects on the structure of the CACNA1F encoded protein (a calcium channel) using parameters based upon changes in amino acid properties; these include size, charge, hydrophobicity, and position. The model produces an overall score for each variant that can be used to predict its pathogenicity. CACNA1F-vp outperformed four other tools in identifying disease-implicated variants (area under receiver operating characteristic and precision recall curves = 0.84; Matthews correlation coefficient = 0.52) using a tenfold cross-validation technique. We consider this protein-specific model to be a robust stand-alone diagnostic classifier that could be replicated in other proteins and could enable precise and timely diagnosis.