Principals in Programming Languages: Technical Results

This is the companion technical report for ``Principals in Programming Languages'' [20]. See that document for a more readable version of these results. In this paper, we describe two variants of the simply typed $\lambda$-calculus extended with a notion of {\em principal}. The results are languages in which intuitive statements like ``the client must call $\mathtt{open}$ to obtain a file handle'' can be phrased and proven formally. The first language is a two-agent calculus with references and recursive types, while the second language explores the possibility of multiple agents with varying amounts of type information. We use these calculi to give syntactic proofs of some type abstraction results that traditionally require semantic arguments

Altruistic Metadynamics: Multisystem Biased Simulation

We present a new simple extension of multiple walker metadynamics which makes it possible to simulate simultaneously multiple different molecular systems and to predict their free energy surfaces, named Altruistic metadynamics. Similarly to basic metadynamics, it uses a bias potential in the form of hills summed over the simulation. Each system adds a big hill to its “own” bias potential and smaller hills to bias potentials of other systems, hence, each system enhances sampling of other systems. This makes it possible to achieve either faster reaching of the stationary point or higher accuracy of the calculated free energy surfaces. This should be efficient in modeling of series of chemically similar systems, for example, in computational drug screening by metadynamics. The method was tested on model energy surfaces, alanine dipeptide modeled in different force fields and monosaccharides of d-hexopyranose series

Decoding the Role of Water Dynamics in Ligand–Protein Unbinding: CRF₁R as a Test Case

The residence time of a ligand–protein complex is a crucial aspect in determining biological effect in vivo. Despite its importance, the prediction of ligand <i>k</i>_off still remains challenging for modern computational chemistry. We have developed aMetaD, a fast and generally applicable computational protocol to predict ligand–protein unbinding events using a molecular dynamics (MD) method based on adiabatic-bias MD and metadynamics. This physics-based, fully flexible, and pose-dependent ligand scoring function evaluates the maximum energy (RTscore) required to move the ligand from the bound-state energy basin to the next. Unbinding trajectories are automatically analyzed and translated into atomic solvation factor (SF) values representing the water dynamics during the unbinding event. This novel computational protocol was initially tested on two M₃ muscarinic receptor and two adenosine A_2A receptor antagonists and then evaluated on a test set of 12 CRF₁R ligands. The resulting RTscores were used successfully to classify ligands with different residence times. Additionally, the SF analysis was used to detect key differences in the degree of accessibility to water molecules during the predicted ligand unbinding events. The protocol provides actionable working hypotheses that are applicable in a drug discovery program for the rational optimization of ligand binding kinetics

GPCR-Bench: A Benchmarking Set and Practitioners’ Guide for G Protein-Coupled Receptor Docking

Virtual screening is routinely used to discover new ligands and in particular new ligand chemotypes for G protein-coupled receptors (GPCRs). To prepare for a virtual screen, we often tailor a docking protocol that will enable us to select the best candidates for further screening. To aid this, we created GPCR-Bench, a publically available docking benchmarking set in the spirit of the DUD and DUD-E reference data sets for validation studies, containing 25 nonredundant high-resolution GPCR costructures with an accompanying set of diverse ligands and computational decoy molecules for each target. Benchmarking sets are often used to compare docking protocols; however, it is important to evaluate docking methods not by “retrospective” hit rates but by the actual likelihood that they will produce novel <i>prospective</i> hits. Therefore, docking protocols must not only rank active molecules highly but also produce good poses that a chemist will select for purchase and screening. Currently, no simple objective machine-scriptable function exists that can do this; instead, docking hit lists must be subjectively examined in a consistent way to compare between docking methods. We present here a case study highlighting considerations we feel are of importance when evaluating a method, intended to be useful as a practitioners’ guide

Making Structural Sense of Dimerization Interfaces of Delta Opioid Receptor Homodimers

Opioid receptors, like other members of the G protein-coupled receptor (GPCR) family, have been shown to associate to form dimers and/or oligomers at the plasma membrane. Whether this association is stable or transient is not known. Recent compelling evidence suggests that at least some GPCRs rapidly associate and dissociate. We have recently calculated binding affinities from free energy estimates to predict transient association between mouse delta opioid receptor (DOR) protomers at a symmetric interface involving the fourth transmembrane (TM4) helix (herein termed “4” dimer). Here we present disulfide cross-linking experiments with DOR constructs with cysteines substituted at the extracellular ends of TM4 or TM5 that confirm the formation of DOR complexes involving these helices. Our results are consistent with the involvement of TM4 and/or TM5 at the DOR homodimer interface, but possibly with differing association propensities. Coarse-grained (CG) well-tempered metadynamics simulations of two different dimeric arrangements of DOR involving TM4 alone or with TM5 (herein termed “4/5” dimer) in an explicit lipid−water environment confirmed the presence of two structurally and energetically similar configurations of the 4 dimer, as previously assessed by umbrella sampling calculations, and revealed a single energetic minimum of the 4/5 dimer. Additional CG umbrella sampling simulations of the 4/5 dimer indicated that the strength of association between DOR protomers varies depending on the protein region at the interface, with the 4 dimer being more stable than the 4/5 dimer

Controlling the Dissociation of Ligands from the Adenosine A_2A Receptor through Modulation of Salt Bridge Strength

The association and dissociation kinetics of ligands binding to proteins vary considerably, but the mechanisms behind this variability are poorly understood, limiting their utilization for drug discovery. This is particularly so for G protein-coupled receptors (GPCRs) where high resolution structural information is only beginning to emerge. Engineering the human A_2A adenosine receptor has allowed structures to be solved in complex with the reference compound ZM241385 and four related ligands at high resolution. Differences between the structures are limited, with the most pronounced being the interaction of each ligand with a salt bridge on the extracellular side of the receptor. Mutagenesis experiments confirm the role of this salt bridge in controlling the dissociation kinetics of the ligands from the receptor, while molecular dynamics simulations demonstrate the ability of ligands to modulate salt bridge stability. These results shed light on a structural determinant of ligand dissociation kinetics and identify a means by which this property may be optimized