Search

We present a robust protocol based on iterations of free energy perturbation (FEP) calculations, chemical synthesis, biophysical mapping and X-ray crystallography to reveal the binding mode of an antagonist series to the A(2A) adenosine receptor (AR). Eight A(2A)AR binding site mutations from biophysical mapping experiments were initially analyzed with sidechain FEP simulations, performed on alternate binding modes. The results distinctively supported one binding mode, which was subsequently used to design new chromone derivatives. Their affinities for the A(2A)AR were experimentally determined and investigated through a cycle of ligand-FEP calculations, validating the binding orientation of the different chemical substituents proposed. Subsequent X-ray crystallography of the A(2A)AR with a low and a high affinity chromone derivative confirmed the predicted binding orientation. The new molecules and structures here reported were driven by free energy calculations, and provide new insights on antagonist binding to the A(2A)AR, an emerging target in immunooncology.

We present a robust protocol based on iterations of free energy perturbation (FEP) calculations, chemical synthesis, biophysical mapping and X-ray crystallography to reveal the binding mode of an antagonist series to the adenosine A2A receptor (AR). Eight A2AAR binding site mutations from biophysical mapping experiments were initially analysed with sidechain FEP simulations, performed on alternate binding modes. The results distinctively supported one binding mode, which was subsequently used to design new chromone derivatives. Their affinities for the A2AAR were experimentally determined and investigated through a cycle of ligandFEP calculations, validating the binding orientation of the different chemical substituents proposed. Subsequent X-ray crystallography of the A2AAR with a low and a high affinity chromone derivative confirmed the predicted binding orientation. The new molecules and structures here reported were driven by free energy calculations, and provide new insights on antagonist binding to the A2AAR, an emerging target in immuno-oncology.

We present a robust protocol based on iterations of free energy perturbation (FEP) calculations, chemical synthesis, biophysical mapping and X‐ray crystallography to reveal the binding mode of an antagonist series to the A2A adenosine receptor (AR). Eight A2AAR binding site mutations from biophysical mapping experiments were initially analyzed with sidechain FEP simulations, performed on alternate binding modes. The results distinctively supported one binding mode, which was subsequently used to design new chromone derivatives. Their affinities for the A2AAR were experimentally determined and investigated through a cycle of ligand‐FEP calculations, validating the binding orientation of the different chemical substituents proposed. Subsequent X‐ray crystallography of the A2AAR with a low and a high affinity chromone derivative confirmed the predicted binding orientation. The new molecules and structures here reported were driven by free energy calculations, and provide new insights on antagonist binding to the A2AAR, an emerging target in immuno‐oncology

A general computational scheme to evaluate the effects of single point mutations on ligand binding is reported. This scheme is applied to characterize agonist binding to the A(2A) adenosine receptor, and is found to accurately explain how point mutations of different nature affect the binding affinity of a potent agonist.

To predict structural and energetic effects of point mutations on ligand binding is of considerable interest in biochemistry and pharmacology. This is not only useful in connection with site-directed mutagenesis experiments, but could also allow interpretation and prediction of individual responses to drug treatment. For G-protein coupled receptors systematic mutagenesis has provided the major part of functional data as structural information until recently has been very limited. For the pharmacologically important A(2A) adenosine receptor, extensive site-directed mutagenesis data on agonist and antagonist binding is available and crystal structures of both types of complexes have been determined. Here, we employ a computational strategy, based on molecular dynamics free energy simulations, to rationalize and interpret available alanine-scanning experiments for both agonist and antagonist binding to this receptor. These computer simulations show excellent agreement with the experimental data and, most importantly, reveal the molecular details behind the observed effects which are often not immediately evident from the crystal structures. The work further provides a distinct validation of the computational strategy used to assess effects of point-mutations on ligand binding. It also highlights the importance of considering not only protein-ligand interactions but also those mediated by solvent water molecules, in ligand design projects.

Despite tremendous successes of GPCR crystallography, the receptors with available structures represent only a small fraction of human GPCRs. An important role of the modeling community is to maximize structural insights for the remaining receptors and complexes. The community-wide GPCR Dock assessment was established to stimulate and monitor the progress in molecular modeling and ligand docking for GPCRs. The four targets in the present third assessment round presented new and diverse challenges for modelers, including prediction of allosteric ligand interaction and activation states in 5-hydroxytryptamine receptors 1B and 2B, and modeling by extremely distant homology for smoothened receptor. Forty-four modeling groups participated in the assessment. State-of-the-art modeling approaches achieved close-to-experimental accuracy for small rigid orthosteric ligands and models built by close homology, and they correctly predicted protein fold for distant homology targets. Predictions of long loops and GPCR activation states remain unsolved problems.

BRICHOS domains are encoded in >30 human genes, which are associated with cancer, neurodegeneration, and interstitial lung disease (ILD). The BRICHOS domain from lung surfactant protein C proprotein (proSP-C) is required for membrane insertion of SP-C and has anti-amyloid activity in vitro. Here, we report the 2.1 angstrom crystal structure of the human proSP-C BRICHOS domain, which, together with molecular dynamics simulations and hydrogen-deuterium exchange mass spectrometry, reveals how BRICHOS domains may mediate chaperone activity. Observation of amyloid deposits composed of mature SP-C in lung tissue samples from ILD patients with mutations in the BRICHOS domain or in its peptide-binding linker region supports the in vivo relevance of the proposed mechanism. The results indicate that ILD mutations interfering with proSP-C BRICHOS activity cause amyloid disease secondary to intramolecular chaperone malfunction.

The protozoan parasite Trypanosoma cruzi, the etiological agent of Chagas disease, affects millions of individuals and continues to be an important global health concern. The poor efficacy and unfavorable side effects of current treatments necessitate novel therapeutics. Cruzain, the major cysteine protease of T.similar to cruzi, is one potential novel target. Recent advances in a class of vinyl sulfone inhibitors are encouraging; however, as most potential therapeutics fail in clinical trials and both disease progression and resistance call for combination therapy with several drugs, the identification of additional classes of inhibitory molecules is essential. Using an exhaustive virtual-screening and experimental validation approach, we identify several additional small-molecule cruzain inhibitors. Further optimization of these chemical scaffolds could lead to the development of novel drugs useful in the treatment of Chagas disease.

The hERG potassium channel is of major pharmaceutical importance, and its blockade by various compounds, potentially causing serious cardiac side effects, is a major problem in drug development. Despite the large amounts of existing biochemical data on blockade of hERG by drugs and druglike compounds, relatively little is known regarding the structural basis of binding of blockers to the channel. Here, we have used a recently developed homology model of hERG to conduct molecular docking experiments with a series of channel blockers, followed by molecular dynamics simulations of the complexes and evaluation of binding free energies with the linear interaction energy method. The calculations yield a remarkably good agreement with experimental binding affinities and allow for a rationalization of three-dimensional structure-activity relationships in terms of a number of key interactions. Two main interaction regions of the channel are thus identified with implications for further mutagenesis experiments and design of new compounds.

Chagas disease, caused by the unicellular parasite Trypanosoma cruzi, claims 50,000 lives annually and is the leading cause of infectious myocarditis in the world. As current antichagastic therapies like nifurtimox and benznidazole are highly toxic, ineffective at parasite eradication, and subject to increasing resistance, novel therapeutics are urgently needed. Cruzain, the major cysteine protease of Trypanosoma cruzi, is one attractive drug target. In the current work, molecular dynamics simulations and a sequence alignment of a non-redundant, unbiased set of peptidase C1 family members are used to identify uncharacterized cruzain binding sites. The two sites identified may serve as targets for future pharmacological intervention.