Search

The cytoplasm is compartmentalized into different translation environments. mRNAs use their 3′UTRs to localize to distinct cytoplasmic compartments, including TIS granules (TGs). Many transcription factors, including MYC, are translated in TGs. It was shown that translation of proteins in TGs enables the formation of protein complexes that cannot be established when these proteins are translated in the cytosol, but the mechanism is poorly understood. Here we show that MYC protein complexes that involve binding to the intrinsically disordered region (IDR) of MYC are only formed when MYC is translated in TGs. TG-dependent protein complexes require TG-enriched mRNAs for assembly. These mRNAs bind to a new and widespread RNA-binding domain in neutral or negatively charged IDRs in several transcription factors, including MYC. RNA-IDR interaction changes the conformational ensemble of the IDR, enabling the formation of MYC protein complexes that act in the nucleus and control functions that cannot be accomplished by cytosolically-translated MYC. We propose that certain mRNAs have IDR chaperone activity as they control IDR conformations. In addition to post-translational modifications, we found a novel mode of protein activity regulation. Since RNA-IDR interactions are prevalent, we suggest that mRNA-dependent control of protein functional states is widespread.

The cowpea chlorotic mottle virus (CCMV) is a plant virus explored as a nanotechnological platform. The robust self-assembly mechanism of its capsid protein allows for drug encapsulation and targeted delivery. Additionally, the capsid nanoparticle can be used as a programmable platform to display different molecular moieties. In view of future applications, efficient production and purification of plant viruses are key steps. In established protocols, the need for ultracentrifugation is a significant limitation due to cost, difficult scalability, and safety issues. In addition, the purity of the final virus isolate often remains unclear. Here, an advanced protocol for the purification of the CCMV from infected plant tissue was developed, focusing on efficiency, economy, and final purity. The protocol involves precipitation with PEG 8000, followed by affinity extraction using a novel peptide aptamer. The efficiency of the protocol was validated using size exclusion chromatography, MALDI-TOF mass spectrometry, reversed-phase HPLC, and sandwich immunoassay. Furthermore, it was demonstrated that the final eluate of the affinity column is of exceptional purity (98.4%) determined by HPLC and detection at 220 nm. The scale-up of our proposed method seems to be straightforward, which opens the way to the large-scale production of such nanomaterials. This highly improved protocol may facilitate the use and implementation of plant viruses as nanotechnological platforms for in vitro and in vivo applications.

The mechanisms by which a protein’s 3D structure can be determined based on its amino acid sequence have long been one of the key mysteries of biophysics. Often simplistic models, such as those derived from geometric constraints, capture bulk real-world 3D protein-protein properties well. One approach is using protein contact maps (PCMs) to better understand proteins’ properties. In this study, we explore the emergent behaviour of contact maps for different geometrically constrained models and compare them to real-world protein systems. Specifically, we derive an analytical approximation for the distribution of amino acid distances, denoted as P(s), using a mean-field approach based on a geometric constraint model. This approximation is then validated for amino acid distance distributions generated from a 2D and 3D version of the geometrically constrained random interaction model. For real protein data, we show how the analytical approximation can be used to fit amino acid distance distributions of protein chain lengths of L ≈ 100, L ≈ 200, and L ≈ 300 generated from two different methods of evaluating a PCM, a simple cutoff based method and a shadow map based method. We present evidence that geometric constraints are sufficient to model the amino acid distance distributions of protein chains in bulk and amino acid sequences only play a secondary role, regardless of the definition of the PCM.

The ratio of the nucleon F2 structure functions, F2n/F2p, is determined by the MARATHON experiment from measurements of deep inelastic scattering of electrons from H3 and He3 nuclei. The experiment was performed in the Hall A Facility of Jefferson Lab using two high-resolution spectrometers for electron detection, and a cryogenic target system which included a low-activity tritium cell. The data analysis used a novel technique exploiting the mirror symmetry of the two nuclei, which essentially eliminates many theoretical uncertainties in the extraction of the ratio. The results, which cover the Bjorken scaling variable range 0.19

Nepenthes diabolica, a new species of tropical pitcher plant from Central Sulawesi, is described and illustrated. It is a high-montane, terrestrial to epiphytic taxon, distinguished by its elaborately flanged peristome and rufous, woolly pitcher indumentum. N. diabolica is a member of section Tentaculatae, showing morphological affinities with N. hamata. Nevertheless, the two species can be distinguished by morphological characters as well as by ecology and geography. A diagnostic key to the species of section Tentaculatae is provided. Owing to its highly restricted range, small number of recorded individuals, and desirability among plant collectors, the conservation status of N. diabolica is assessed as Critically Endangered under the IUCN criteria. Finally, N. hamata, until now based on lost types and mixed paratypes, is lectotypified using recently uncovered original material.

The ratio of the nucleon F_{2} structure functions, F_{2}^{n}/F_{2}^{p}, is determined by the MARATHON experiment from measurements of deep inelastic scattering of electrons from ^{3}H and ^{3}He nuclei. The experiment was performed in the Hall A Facility of Jefferson Lab using two high-resolution spectrometers for electron detection, and a cryogenic target system which included a low-activity tritium cell. The data analysis used a novel technique exploiting the mirror symmetry of the two nuclei, which essentially eliminates many theoretical uncertainties in the extraction of the ratio. The results, which cover the Bjorken scaling variable range 0.19x0.83, represent a significant improvement compared to previous SLAC and Jefferson Lab measurements for the ratio. They are compared to recent theoretical calculations and empirical determinations of the F_{2}^{n}/F_{2}^{p} ratio.

When protons and neutrons (nucleons) are bound into atomic nuclei, they are close enough to feel significant attraction, or repulsion, from the strong, short-distance part of the nucleon-nucleon interaction. These strong interactions lead to hard collisions between nucleons, generating pairs of highly energetic nucleons referred to as short-range correlations (SRCs). SRCs are an important but relatively poorly understood part of nuclear structure

Data availability statement: The trajectories of Mpro simulations and models of the metastable states can be downloaded from 10.5281/zenodo.4782284. β-coronavirus (CoVs) alone has been responsible for three major global outbreaks in the 21st century. The current crisis has led to an urgent requirement to develop therapeutics. Even though a number of vaccines are available, alternative strategies targeting essential viral components are required as a backup against the emergence of lethal viral variants. One such target is the main protease (Mpro) that plays an indispensable role in viral replication. The availability of over 270 Mpro X-ray structures in complex with inhibitors provides unique insights into ligand–protein interactions. Herein, we provide a comprehensive comparison of all nonredundant ligand-binding sites available for SARS-CoV2, SARS-CoV, and MERS-CoV Mpro. Extensive adaptive sampling has been used to investigate structural conservation of ligand-binding sites using Markov state models (MSMs) and compare conformational dynamics employing convolutional variational auto-encoder-based deep learning. Our results indicate that not all ligand-binding sites are dynamically conserved despite high sequence and structural conservation across β-CoV homologs. This highlights the complexity in targeting all three Mpro enzymes with a single pan inhibitor. There was no funding for this work

This work investigates the connection between stochastic protein dynamics and function for the enzyme cyclophilin A (CypA) in wild-type form, and three variants that feature several mutations distal from the active site. Previous biophysical studies have suggested that conformational exchange between a ‘major’ active and a ‘minor’ inactive state on millisecond timescales plays a key role in catalysis for CypA. Here this hypothesis is addressed by a variety of molecular dynamics simulation techniques. Strikingly we show that exchange between major and minor active site conformations occurs at a rate that is 5 to 6 orders of magnitude faster than previously proposed. The minor active site conformation is found to be catalytically impaired, and decreased catalytic activity of the mutants is caused by changes in Phe113 motions on a ns-μs timescale. Therefore millisecond timescale motions may not be necessary to explain allosteric effects in cyclophilins. The relationship between molecular motion and catalysis in enzymes is debated. Here, simulations of cyclophilin A and three catalytically-impaired mutants reveal a nanosecond-scale interconversion between active and inactive conformations, orders of magnitude faster than previously suggested.

The quantum mechanical bespoke (QUBE) force-field approach has been developed to facilitate the automated derivation of potential energy function parameters for modeling protein-ligand binding. To date, the approach has been validated in the context of Monte Carlo simulations of protein-ligand complexes. We describe here the implementation of the QUBE force field in the alchemical free-energy calculation molecular dynamics simulation package SOMD. The implementation is validated by demonstrating the reproducibility of absolute hydration free energies computed with the QUBE force field across the SOMD and GROMACS software packages. We further demonstrate, by way of a case study involving two series of non-nucleoside inhibitors of HIV-1 reverse transcriptase, that the availability of QUBE in a modern simulation package that makes efficient use of graphics processing unit acceleration will facilitate high-throughput alchemical free-energy calculations.

This article is a preprint and has not been certified by peer review. It was eventually published online on 14 June 2021 by American Chemical Society (ACS Publications) as: Cho, E. et al. (2021) 'Dynamic Profiling of β-Coronavirus 3CL Mpro Protease Ligand-Binding Sites', Journal of Chemical Information and Modeling, 61 (6), pp. 3058 - 3073. doi: 10.1021/acs.jcim.1c00449. Data Availability Statement: The trajectories of Mpro simulations and models of the metastable states can be obtained from the corresponding author. Jupyter-notebooks to generate Markov State Models can be downloaded from 10.6084/m9.figshare.14343725 bioRxiv posts many COVID19-related papers. A reminder: they have not been formally peer-reviewed and should not guide health-related behavior or be reported in the press as conclusive. β-coronaviruses alone have been responsible for three major global outbreaks in the 21st century. The current crisis has led to an urgent requirement to develop therapeutics. Even though a number of vaccines are available, alternative strategies targeting essential viral components are required as a back-up against the emergence of lethal viral variants. One such target is the main protease (Mpro) that plays an indispensible role in viral replication. The availability of over 270 Mpro X-ray structures in complex with inhibitors provides unique insights into ligand-protein interactions. Herein, we provide a comprehensive comparison of all non-redundant ligand-binding sites available for SARS-CoV2, SARS-CoV and MERS-CoV Mpro. Extensive adaptive sampling has been used to explore conformational dynamics employing convolutional variational auto encoder-based deep learning, and investigates structural conservation of the ligand binding sites using Markov state models across β-coronavirus homologs. Our results indicate that not all ligand-binding sites are dynamically conserved despite high sequence and structural conservation across β-coronavirus homologs. This highlights the complexity in targeting all three Mpro enzymes with a single pan inhibitor. U.S. Department of Energy, Office of Science, through the Advanced Scientific Computing Research (ASCR), under contract number DEAC05-00OR22725 and the Exascale Computing Project (ECP) (17-SC-20-SC). BI would like to acknowledge the COVID-19 pump-priming grant from the University of Huddersfield for funding computing resources for analysis.

A methodology that combines alchemical free energy calculations (FEP) with machine learning (ML) has been developed to compute accurate absolute hydration free energies. The hybrid FEP/ML methodology was trained on a subset of the FreeSolv database and retrospectively shown to outperform most submissions from the SAMPL4 competition. Compared to pure machine-learning approaches, FEP/ML yields more precise estimates of free energies of hydration and requires a fraction of the training set size to outperform standalone FEP calculations. The ML-derived correction terms are further shown to be transferable to a range of related FEP simulation protocols. The approach may be used to inexpensively improve the accuracy of FEP calculations and to flag molecules which will benefit the most from bespoke force field parametrization efforts.

A methodology that combines alchemical free energy calculations (FEP) with machine learning (ML) has been developed to compute accurate absolute hydration free energies. The hybrid FEP/ML methodology was trained on a subset of the FreeSolv database, and retrospectively shown to outperform most submissions from the SAMPL4 competition. Compared to pure machine-learning approaches, FEP/ML yields more precise estimates of free energies of hydration, and requires a fraction of the training set size to outperform standalone FEP calculations. The ML-derived correction terms are further shown to be transferable to a range of related FEP simulation protocols. The approach may be used to inexpensively improve the accuracy of FEP calculations, and to flag molecules which will benefit the most from bespoke forcefield parameterisation efforts.

We report the first measurement of the (e,e^{'}p) three-body breakup reaction cross sections in helium-3 (^{3}He) and tritium (^{3}H) at large momentum transfer [⟨Q^{2}⟩≈1.9 (GeV/c)^{2}] and x_{B}>1 kinematics, where the cross section should be sensitive to quasielastic (QE) scattering from single nucleons. The data cover missing momenta 40≤p_{miss}≤500 MeV/c that, in the QE limit with no rescattering, equals the initial momentum of the probed nucleon. The measured cross sections are compared with state-of-the-art ab initio calculations. Overall good agreement, within ±20%, is observed between data and calculations for the full p_{miss} range for ^{3}H and for 100≤p_{miss}≤350 MeV/c for ^{3}He. Including the effects of rescattering of the outgoing nucleon improves agreement with the data at p_{miss}>250 MeV/c and suggests contributions from charge-exchange (SCX) rescattering. The isoscalar sum of ^{3}He plus ^{3}H, which is largely insensitive to SCX, is described by calculations to within the accuracy of the data over the entire p_{miss} range. This validates current models of the ground state of the three-nucleon system up to very high initial nucleon momenta of 500 MeV/c.

Free-energy calculations have seen increased usage in structure-based drug design. Despite the rising interest, automation of the complex calculations and subsequent analysis of their results are still hampered by the restricted choice of available tools. In this work, an application for automated setup and processing of free-energy calculations is presented. Several sanity checks for assessing the reliability of the calculations were implemented, constituting a distinct advantage over existing open-source tools. The underlying workflow is built on top of the software Sire, SOMD, BioSimSpace, and OpenMM and uses the AMBER 14SB and GAFF2.1 force fields. It was validated on two datasets originally composed by Schrodinger, consisting of 14 protein structures and 220 ligands. Predicted binding affinities were in good agreement with experimental values. For the larger dataset, the average correlation coefficient Rp was 0.70 ± 0.05 and average Kendall's τ was 0.53 ± 0.05, which are broadly comparable to or better than previously reported results using other methods.

Free energy calculations have seen increased usage in structure-based drug design. Despite the rising interest, automation of the complex calculations and subsequent analysis of their results are still hampered by the restricted choice of available tools. In this work, an application for automated setup and processing of free energy calculations is presented. Several sanity checks for assessing the reliability of the calculations were implemented, constituting a distinct advantage over existing open-source tools. The underlying workflow is built on top of the software Sire, SOMD, BioSimSpace and OpenMM and uses the AMBER14SB and GAFF2.1 force fields. It was validated on two datasets originally composed by Schrödinger, consisting of 14 protein structures and 220 ligands. Predicted binding affinities were in good agreement with experimental values. For the larger dataset the average correlation coefficient Rp was 0.70 ± 0.05 and average Kendall’s τ was 0.53 ± 0.05 which is broadly comparable to or better than previously reported results using other methods.

Proteins need to interconvert between many conformations in order to function, many of which are formed transiently, and sparsely populated. Particularly when the lifetimes of these states approach the millisecond timescale, identifying the relevant structures and the mechanism by which they interconvert remains a tremendous challenge. Here we introduce a novel combination of accelerated MD (aMD) simulations and Markov state modelling (MSM) to explore these ‘excited’ conformational states. Applying this to the highly dynamic protein CypA, a protein involved in immune response and associated with HIV infection, we identify five principally populated conformational states and the atomistic mechanism by which they interconvert. A rational design strategy predicted that the mutant D66A should stabilise the minor conformations and substantially alter the dynamics, whereas the similar mutant H70A should leave the landscape broadly unchanged. These predictions are confirmed using CPMG and R1ρ solution state NMR measurements. By efficiently exploring functionally relevant, but sparsely populated conformations with millisecond lifetimes in silico, our aMD/MSM method has tremendous promise for the design of dynamic protein free energy landscapes for both protein engineering and drug discovery., Molecular simulations were used to design large scale loop motions in the enzyme cyclophilin A and NMR and biophysical methods were employed to validate the models.

We report the first measurement of the (e,e^{'}p) three-body breakup reaction cross sections in helium-3 (^{3}He) and tritium (^{3}H) at large momentum transfer [⟨Q^{2}⟩≈1.9 (GeV/c)^{2}] and x_{B}1 kinematics, where the cross section should be sensitive to quasielastic (QE) scattering from single nucleons. The data cover missing momenta 40≤p_{miss}≤500 MeV/c that, in the QE limit with no rescattering, equals the initial momentum of the probed nucleon. The measured cross sections are compared with state-of-the-art ab initio calculations. Overall good agreement, within ±20%, is observed between data and calculations for the full p_{miss} range for ^{3}H and for 100≤p_{miss}≤350 MeV/c for ^{3}He. Including the effects of rescattering of the outgoing nucleon improves agreement with the data at p_{miss}250 MeV/c and suggests contributions from charge-exchange (SCX) rescattering. The isoscalar sum of ^{3}He plus ^{3}H, which is largely insensitive to SCX, is described by calculations to within the accuracy of the data over the entire p_{miss} range. This validates current models of the ground state of the three-nucleon system up to very high initial nucleon momenta of 500 MeV/c.

Alchemical free energy calculations are a useful tool for predicting free energy differences associated with the transfer of molecules from one environment to another. The hallmark of these methods is the use of “bridging” potential energy functions representing alchemical intermediate states that cannot exist as real chemical species. The data collected from these bridging alchemical thermodynamic states allows the efficient computation of transfer free energies (or differences in transfer free energies) with orders of magnitude less simulation time than simulating the transfer process directly. While these methods are highly flexible, care must be taken in avoiding common pitfalls to ensure that computed free energy differences can be robust and reproducible for the chosen force field, and that appropriate corrections are included to permit direct comparison with experimental data. In this paper, we review current best practices for several popular application domains of alchemical free energy calculations performed with equilibrium simulations, in particular relative and absolute small molecule binding free energy calculations to biomolecular targets.

High-momentum configurations of nucleon pairs at short-distance are probed using measurements of the $^{12}$C$(e,e'p)$ and $^{12}$C$(e,e'pN)$ reactions (where $N$ is either $n$ or $p$), at high-$Q^2$ and $x_B>1.1$. The data span a missing-momentum range of 300--1000 MeV/c and are predominantly sensitive to the transition region of the strong nuclear interaction from a Tensor to Scalar interaction. The data are well reproduced by theoretical calculations using the Generalized Contact Formalism with both chiral and phenomenological nucleon-nucleon ($NN$) interaction models. This agreement suggests that the measured high missing-momentum protons up to $1000$ MeV/c predominantly belong to short-ranged correlated (SRC) pairs. The measured $^{12}$C$(e,e'pN)$ / $^{12}$C$(e,e'p)$ and $^{12}$C$(e,e'pp)$ / $^{12}$C$(e,e'pn)$ cross-section ratios are consistent with a decrease in the fraction of proton-neutron SRC pairs and increase in the fraction of proton-proton SRC pairs with increasing missing momentum. This confirms the transition from an isospin-dependent tensor $NN$ interaction at $\sim 400$ MeV/c to an isospin-independent scalar interaction at high-momentum around $\sim 800$ MeV/c as predicted by theoretical calculation., Accepted for publication in Physics Letters B. 7 pages, 3 figures, and online supplementary materials

The Drug Design Data Resource (D3R) consortium organises blinded challenges to address the latest advances in computational methods for ligand pose prediction, affinity ranking, and free energy calculations. Within the context of the second D3R Grand Challenge several blinded binding free energies predictions were made for two congeneric series of Farsenoid X Receptor (FXR) inhibitors with a semi-automated alchemical free energy calculation workflow featuring FESetup and SOMD software tools. Reasonable performance was observed in retrospective analyses of literature datasets. Nevertheless, blinded predictions on the full D3R datasets were poor due to difficulties encountered with the ranking of compounds that vary in their net-charge. Performance increased for predictions that were restricted to subsets of compounds carrying the same net-charge. Disclosure of X-ray crystallography derived binding modes maintained or improved the correlation with experiment in a subsequent rounds of predictions. The best performing protocols on D3R set1 and set2 were comparable or superior to predictions made on the basis of analysis of literature structure activity relationships (SAR)s only, and comparable or slightly inferior, to the best submissions from other groups. Electronic supplementary material The online version of this article (10.1007/s10822-017-0083-9) contains supplementary material, which is available to authorized users.

We report the first measurement of the $(e,e'p)$ reaction cross-section ratios for Helium-3 ($^3$He), Tritium ($^3$H), and Deuterium ($d$). The measurement covered a missing momentum range of $40 \le p_{miss} \le 550$ MeV$/c$, at large momentum transfer ($\langle Q^2 \rangle \approx 1.9$ (GeV$/c$)$^2$) and $x_B>1$, which minimized contributions from non quasi-elastic (QE) reaction mechanisms. The data is compared with plane-wave impulse approximation (PWIA) calculations using realistic spectral functions and momentum distributions. The measured and PWIA-calculated cross-section ratios for $^3$He$/d$ and $^3$H$/d$ extend to just above the typical nucleon Fermi-momentum ($k_F \approx 250$ MeV$/c$) and differ from each other by $\sim 20\%$, while for $^3$He/$^3$H they agree within the measurement accuracy of about 3\%. At momenta above $k_F$, the measured $^3$He/$^3$H ratios differ from the calculation by $20\% - 50\%$. Final state interaction (FSI) calculations using the generalized Eikonal Approximation indicate that FSI should change the $^3$He/$^3$H cross-section ratio for this measurement by less than 5\%. If these calculations are correct, then the differences at large missing momenta between the $^3$He/$^3$H experimental and calculated ratios could be due to the underlying $NN$ interaction, and thus could provide new constraints on the previously loosely-constrained short-distance parts of the $NN$ interaction., 8 pages, 3 figures (4 panels)

We measured the triple coincidence A(e,e'np) and A(e,e'pp) reactions on carbon, aluminum, iron, and lead targets at Q2 > 1.5 (GeV/c)2, xB > 1.1 and missing momentum > 400 MeV/c. This was the first direct measurement of both proton-proton (pp) and neutron-proton (np) short-range correlated (SRC) pair knockout from heavy asymmetric nuclei. For all measured nuclei, the average proton-proton (pp) to neutron-proton (np) reduced cross-section ratio is about 6%, in agreement with previous indirect measurements. Correcting for Single-Charge Exchange effects decreased the SRC pairs ratio to ~ 3%, which is lower than previous results. Comparisons to theoretical Generalized Contact Formalism (GCF) cross-section calculations show good agreement using both phenomenological and chiral nucleon-nucleon potentials, favoring a lower pp to np pair ratio. The ability of the GCF calculation to describe the experimental data using either phenomenological or chiral potentials suggests possible reduction of scale- and scheme-dependence in cross section ratios. Our results also support the high-resolution description of high-momentum states being predominantly due to nucleons in SRC pairs., Accepted for publication in PRL. 9 pages, 3 figures, 3 tables and supplementary materials

The intricate three-dimensional geometries of protein tertiary structures underlie protein function and emerge through a folding process from one-dimensional chains of amino acids. The exact spatial sequence and configuration of amino acids, the biochemical environment and the temporal sequence of distinct interactions yield a complex folding process that cannot yet be easily tracked for all proteins. To gain qualitative insights into the fundamental mechanisms behind the folding dynamics and generic features of the folded structure, we propose a simple model of structure formation that takes into account only fundamental geometric constraints and otherwise assumes randomly paired connections. We find that despite its simplicity, the model results in a network ensemble consistent with key overall features of the ensemble of Protein Residue Networks we obtained from more than 1000 biological protein geometries as available through the Protein Data Base. Specifically, the distribution of the number of interaction neighbors a unit (amino acid) has, the scaling of the structure's spatial extent with chain length, the eigenvalue spectrum and the scaling of the smallest relaxation time with chain length are all consistent between model and real proteins. These results indicate that geometric constraints alone may already account for a number of generic features of protein tertiary structures.

In the framework of the 2015 D3R inaugural grand challenge, blind binding pose and affinity predictions were performed for a set of 180 ligands of the Heat Shock Protein HSP90-α protein, a relevant cancer target. Spectral clustering was used to rapidly identify alternative binding site conformations in publicly available crystallographic HSP90-α structures. Subsequently, multiple docking and scoring protocols employing the software Autodock Vina and rDock were applied to predict binding modes and rank order ligands. Alchemical free energy calculations were performed with the software FESetup and Sire/OpenMM to predict binding affinities for three congeneric series subsets. Some of the protocols used here were ranked among the top submissions according to most of the evaluation metrics. Docking performance was excellent, but the scoring results were disappointing. A critical assessment of the results is reported, as well as suggestions for future similar competitions.

Accurate predictions of free energies of binding of ligands to proteins are challenging partly because of the nonadditivity of protein-ligand interactions; i.e., the free energy of binding is the sum of numerous enthalpic and entropic contributions that cannot be separated into functional group contributions. In principle, molecular simulations methodologies that compute free energies of binding do capture nonadditivity of protein-ligand interactions, but efficient protocols are necessary to compute well-converged free energies of binding that clearly resolve nonadditive effects. To this end, an efficient GPU-accelerated implementation of alchemical free energy calculations has been developed and applied to two congeneric series of ligands of the enzyme thrombin. The results show that accurate binding affinities are computed across the two congeneric series and positive coupling between nonpolar R1 substituents and a X = NH3+ substituent is reproduced, albeit with a weaker trend than experimentally observed. By contrast, a docking methodology completely fails to capture nonadditive effects. Further analysis shows that the nonadditive effects are partly due to variations in the strength of a hydrogen-bond between the X = NH3+ ligands family and thrombin residue Gly216. However, other partially compensating interactions occur across the entire binding site, and no single interaction dictates the magnitude of the nonadditive effects for all the analyzed protein-ligand complexes.

The intricate three-dimensional geometries of protein tertiary structures underlie protein function and emerge through a folding process from one-dimensional chains of amino acids. The exact spatial sequence and configuration of amino acids, the biochemical environment and the temporal sequence of distinct interactions yield a complex folding process that cannot yet be easily tracked for all proteins. To gain qualitative insights into the fundamental mechanisms behind the folding dynamics and generic features of the folded structure, we propose a simple model of structure formation that takes into account only fundamental geometric constraints and otherwise assumes randomly paired connections. We find that despite its simplicity, the model results in a network ensemble consistent with key overall features of the ensemble of Protein Residue Networks we obtained from more than 1000 biological protein geometries as available through the Protein Data Base. Specifically, the distribution of the number of interaction neighbors a unit (amino acid) has, the scaling of the structure’s spatial extent with chain length, the eigenvalue spectrum and the scaling of the smallest relaxation time with chain length are all consistent between model and real proteins. These results indicate that geometric constraints alone may already account for a number of generic features of protein tertiary structures.

Hit-to-lead virtual screening frequently relies on a cascade of computational methods that starts with rapid calculations applied to a large number of compounds and ends with more expensive computations restricted to a subset of compounds that passed initial filters. This work focuses on set up protocols for alchemical free energy (AFE) scoring in the context of a Docking–MM/PBSA–AFE cascade. A dataset of 15 congeneric inhibitors of the ACK1 protein was used to evaluate the performance of AFE set up protocols that varied in the steps taken to prepare input files (using previously docked and best scored poses, manual selection of poses, manual placement of binding site water molecules). The main finding is that use of knowledge derived from X-ray structures to model binding modes, together with the manual placement of a bridging water molecule, improves the R2 from 0.45 ± 0.06 to 0.76 ± 0.02 and decreases the mean unsigned error from 2.11 ± 0.08 to 1.24 ± 0.04 kcal mol-1. By contrast a brute force automated protocol that increased the sampling time ten-fold lead to little improvements in accuracy. Besides, it is shown that for the present dataset hysteresis can be used to flag poses that need further attention even without prior knowledge of experimental binding affinities.