Binding free energy calculations that make use of alchemical pathways are becoming increasingly feasible thanks to advances in hardware and algorithms. Although relative binding free energy (RBFE) ...calculations are starting to find widespread use, absolute binding free energy (ABFE) calculations are still being explored mainly in academic settings due to the high computational requirements and still uncertain predictive value. However, in some drug design scenarios, RBFE calculations are not applicable and ABFE calculations could provide an alternative. Computationally cheaper end-point calculations in implicit solvent, such as molecular mechanics Poisson–Boltzmann surface area (MMPBSA) calculations, could too be used if one is primarily interested in a relative ranking of affinities. Here, we compare MMPBSA calculations to previously performed absolute alchemical free energy calculations in their ability to correlate with experimental binding free energies for three sets of bromodomain–inhibitor pairs. Different MMPBSA approaches have been considered, including a standard single-trajectory protocol, a protocol that includes a binding entropy estimate, and protocols that take into account the ligand hydration shell. Despite the improvements observed with the latter two MMPBSA approaches, ABFE calculations were found to be overall superior in obtaining correlation with experimental affinities for the test cases considered. A difference in weighted average Pearson ( r p ̅ ) and Spearman ( r s ̅ ) correlations of 0.25 and 0.31 was observed when using a standard single-trajectory MMPBSA setup ( r p ̅ = 0.64 and r s ̅ = 0.66 for ABFE; r p ̅ = 0.39 and r s ̅ = 0.35 for MMPBSA). The best performing MMPBSA protocols returned weighted average Pearson and Spearman correlations that were about 0.1 inferior to ABFE calculations: r p ̅ = 0.55 and r s ̅ = 0.56 when including an entropy estimate, and r p ̅ = 0.53 and r s ̅ = 0.55 when including explicit water molecules. Overall, the study suggests that ABFE calculations are indeed the more accurate approach, yet there is also value in MMPBSA calculations considering the lower compute requirements, and if agreement to experimental affinities in absolute terms is not of interest. Moreover, for the specific protein–ligand systems considered in this study, we find that including an explicit ligand hydration shell or a binding entropy estimate in the MMPBSA calculations resulted in significant performance improvements at a negligible computational cost.
Cysteine plays an essential role in cellular redox homoeostasis as a key constituent of the tripeptide glutathione (GSH). A rate limiting step in cellular GSH synthesis is the availability of ...cysteine. However, circulating cysteine exists in the blood as the oxidised di-peptide cystine, requiring specialised transport systems for its import into the cell. System xc
is a dedicated cystine transporter, importing cystine in exchange for intracellular glutamate. To counteract elevated levels of reactive oxygen species in cancerous cells system xc
is frequently upregulated, making it an attractive target for anticancer therapies. However, the molecular basis for ligand recognition remains elusive, hampering efforts to specifically target this transport system. Here we present the cryo-EM structure of system xc
in both the apo and glutamate bound states. Structural comparisons reveal an allosteric mechanism for ligand discrimination, supported by molecular dynamics and cell-based assays, establishing a mechanism for cystine transport in human cells.
Water plays a critical role in ligand-protein interactions. However, it is still challenging to predict accurately not only where water molecules prefer to bind, but also which of those water ...molecules might be displaceable. The latter is often seen as a route to optimizing affinity of potential drug candidates. Using a protocol we call WaterDock, we show that the freely available AutoDock Vina tool can be used to predict accurately the binding sites of water molecules. WaterDock was validated using data from X-ray crystallography, neutron diffraction and molecular dynamics simulations and correctly predicted 97% of the water molecules in the test set. In addition, we combined data-mining, heuristic and machine learning techniques to develop probabilistic water molecule classifiers. When applied to WaterDock predictions in the Astex Diverse Set of protein ligand complexes, we could identify whether a water molecule was conserved or displaced to an accuracy of 75%. A second model predicted whether water molecules were displaced by polar groups or by non-polar groups to an accuracy of 80%. These results should prove useful for anyone wishing to undertake rational design of new compounds where the displacement of water molecules is being considered as a route to improved affinity.
Accurate prediction of binding affinities has been a central goal of computational chemistry for decades, yet remains elusive. Despite good progress, the required accuracy for use in a drug-discovery ...context has not been consistently achieved for drug-like molecules. Here, we perform absolute free energy calculations based on a thermodynamic cycle for a set of diverse inhibitors binding to bromodomain-containing protein 4 (BRD4) and demonstrate that a mean absolute error of 0.6 kcal mol
can be achieved. We also show a similar level of accuracy (1.0 kcal mol
) can be achieved in pseudo prospective approach. Bromodomains are epigenetic mark readers that recognize acetylation motifs and regulate gene transcription, and are currently being investigated as therapeutic targets for cancer and inflammation. The unprecedented accuracy offers the exciting prospect that the binding free energy of drug-like compounds can be predicted for pharmacologically relevant targets.
There are increasing numbers of ion channel structures featuring heteromeric subunit assembly, exemplified by synaptic α1βB Glycine and α4β2 Nicotinic receptors. These structures exhibit inherent ...pore asymmetry, but the relevance of this to function is unknown. Furthermore, molecular dynamics simulations performed on symmetrical homomeric channels often leads to thermal distortion whereby conformations of the resulting ensemble are also asymmetrical. When functionally annotating ion channels, researchers often rely on minimal constrictions determined via radius-profile calculations performed with computer programs, such as HOLE or CHAP, coupled with an assessment of pore hydrophobicity. However, such tools typically employ spherical probe particles, limiting their ability to accurately capture pore asymmetry. Here, we introduce an algorithm that employs ellipsoidal probe particles, enabling a more comprehensive representation of the pore geometry. Our analysis reveals that the use of non-spherical ellipsoids for pore characterization, provides a more accurate and easily interpretable depiction of conductance. To quantify the implications of pore asymmetry on conductance, we systematically investigated carbon nanotubes (CNTs) with varying degrees of pore asymmetry as model systems. The conductance through these channels shows surprising effects that would otherwise not be predicted with spherical probes. The results have broad implications not only for the functional annotation of biological ion channels, but also for the design of synthetic channel systems for use in areas such as water filtration. Furthermore, we make use of the more accurate characterization of channel pores to refine a physical conductance model to obtain a heuristic estimate for single channel conductance. The code is freely available, obtainable as pip-installable python package and provided as a webservice.There are increasing numbers of ion channel structures featuring heteromeric subunit assembly, exemplified by synaptic α1βB Glycine and α4β2 Nicotinic receptors. These structures exhibit inherent pore asymmetry, but the relevance of this to function is unknown. Furthermore, molecular dynamics simulations performed on symmetrical homomeric channels often leads to thermal distortion whereby conformations of the resulting ensemble are also asymmetrical. When functionally annotating ion channels, researchers often rely on minimal constrictions determined via radius-profile calculations performed with computer programs, such as HOLE or CHAP, coupled with an assessment of pore hydrophobicity. However, such tools typically employ spherical probe particles, limiting their ability to accurately capture pore asymmetry. Here, we introduce an algorithm that employs ellipsoidal probe particles, enabling a more comprehensive representation of the pore geometry. Our analysis reveals that the use of non-spherical ellipsoids for pore characterization, provides a more accurate and easily interpretable depiction of conductance. To quantify the implications of pore asymmetry on conductance, we systematically investigated carbon nanotubes (CNTs) with varying degrees of pore asymmetry as model systems. The conductance through these channels shows surprising effects that would otherwise not be predicted with spherical probes. The results have broad implications not only for the functional annotation of biological ion channels, but also for the design of synthetic channel systems for use in areas such as water filtration. Furthermore, we make use of the more accurate characterization of channel pores to refine a physical conductance model to obtain a heuristic estimate for single channel conductance. The code is freely available, obtainable as pip-installable python package and provided as a webservice.
Binding selectivity is a requirement for the development of a safe drug, and it is a critical property for chemical probes used in preclinical target validation. Engineering selectivity adds ...considerable complexity to the rational design of new drugs, as it involves the optimization of multiple binding affinities. Computationally, the prediction of binding selectivity is a challenge, and generally applicable methodologies are still not available to the computational and medicinal chemistry communities. Absolute binding free energy calculations based on alchemical pathways provide a rigorous framework for affinity predictions and could thus offer a general approach to the problem. We evaluated the performance of free energy calculations based on molecular dynamics for the prediction of selectivity by estimating the affinity profile of three bromodomain inhibitors across multiple bromodomain families, and by comparing the results to isothermal titration calorimetry data. Two case studies were considered. In the first one, the affinities of two similar ligands for seven bromodomains were calculated and returned excellent agreement with experiment (mean unsigned error of 0.81 kcal/mol and Pearson correlation of 0.75). In this test case, we also show how the preferred binding orientation of a ligand for different proteins can be estimated via free energy calculations. In the second case, the affinities of a broad-spectrum inhibitor for 22 bromodomains were calculated and returned a more modest accuracy (mean unsigned error of 1.76 kcal/mol and Pearson correlation of 0.48); however, the reparametrization of a sulfonamide moiety improved the agreement with experiment.
Water is often found to mediate interactions between a ligand and a protein. It can play a significant role in orientating the ligand within a binding pocket and contribute to the free energy of ...binding. It would thus be extremely useful to be able to accurately predict the position and orientation of water molecules within a binding pocket. Recently, we developed the WaterDock protocol that was able to predict 97% of the water molecules in a test set. However, this approach generated false positives at a rate of over 20% in most cases and whilst this might be acceptable for some applications, in high throughput scenarios this is not desirable. Here we tackle this problem via the inclusion of knowledge regarding the solvation structure of ligand functional groups. We call this new protocol WaterDock2 and demonstrate that this protocol maintains a similar true positive rate to the original implementation but is capable of reducing the false-positive rate by over 50%. To improve the usability of the method, we have also developed a plugin for the popular graphics program PyMOL. The plugin also contains an implementation of the original WaterDock.
Pentameric ligand-gated ion channels are key players in mediating fast neurotransmission. Glycine receptors are chloride-selective members of this receptor family that mediate inhibitory synaptic ...transmission and are implicated in neurological disorders including autism and hyperekplexia. They have been structurally characterized by both X-ray crystallography and cryoelectron microscopy (cryo-EM) studies, with the latter giving rise to what was proposed as a possible open state. However, recent work has questioned the physiological relevance of this open state structure, since it rapidly collapses in molecular dynamics simulations. Here, we show that the collapse can be avoided by a careful equilibration protocol that reconciles the more problematic regions of the original density map and gives a stable open state that shows frequent selective chloride permeation. The protocol developed in this work provides a means to refine open-like structures of the whole pentameric ligand-gated ion channel superfamily and reconciles the previous issues with the cryo-EM structure.
Display omitted
•MD is used to refine problematic regions of the open state of the glycine receptor•We functionally annotate it as open and selective for chloride ions•The open state is stabilized by the 9′ residues entering conserved hydrophobic pockets•The protocol can be more broadly applied to all members of the Cys-loop family
Dmgen and Biggin use a molecular dynamics protocol to refine problematic regions of the open state cryo-EM structure of the glycine receptor. They demonstrate that a stable open state is achieved when the leucine gate residues at the 9′ position move into well-conserved hydrophobic pockets at the interface between subunits.