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 Å 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.

A series of acylguanidine beta secretase 1 (BACE1) inhibitors with modified scaffold and P3 pocket substituent was synthesized and studied with free energy perturbation (FEP) calculations. The resulting molecules showed potencies in enzymatic BACE1 inhibition assays up to 1 nM. The correlation between the predicted activity from the FEP calculations and the experimental activity was good for the P3 pocket substituents. The average mean unsigned error (MUE) between prediction and experiment was 0.68 ± 0.17 kcal/mol for the default 5 ns lambda window simulation time improving to 0.35 ± 0.13 kcal/mol for 40 ns. FEP calculations for the P2' pocket substituents on the same acylguanidine scaffold also showed good agreement with experiment and the results remained stable with repeated simulations and increased simulation time. It proved more difficult to use FEP calculations to study the scaffold modification from increasing 5 to 6 and 7 membered-rings. Although prediction and experiment were in agreement for short 2 ns simulations, as the simulation time increased the results diverged. This was improved by the use of a newly developed "Core Hopping FEP+" approach, which also showed improved stability in repeat calculations. The origins of these differences along with the value of repeat and longer simulation times are discussed. This work provides a further example of the use of FEP as a computational tool for molecular design.

A congeneric series of 21 phosphodiesterase 2 (PDE2) inhibitors are reported. Crystal structures show how the molecules can occupy a ‘top-pocket’ of the active site. Molecules with small substituents do not enter the pocket, a critical leucine (Leu770) is closed and water molecules are present. Large substituents enter the pocket, opening the Leu770 conformation and displacing the waters. We also report an X-ray structure revealing a new conformation of the PDE2 active site domain. The relative binding affinities of these compounds were studied with free energy perturbation (FEP) methods and it represents an attractive real-world test case. In general, the calculations could predict the energy of small-to-small, or large-to-large molecule perturbations. However, accurately capturing the transition from small-to-large proved challenging. Only when using alternative protein conformations did results improve. The new X-ray structure, along with a modelled dimer, conferred stability to the catalytic domain during the FEP molecular dynamics (MD) simulations, increasing the convergence and thereby improving the prediction of ΔΔG of binding for some small-to-large transitions. In summary, we found the most significant improvement in results when using different protein structures, and this data set is useful for future free energy validation studies.

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 A2A adenosine receptor, and is found to accurately explain how point mutations of different nature affect the binding affinity of a potent agonist.