Cisplatin uptake and release in pH sensitive zeolitic imidazole frameworks.

Cancer remains hard to treat, partially due to the non-specificity of chemotherapeutics. Metal-organic frameworks (MOFs) are promising carriers for targeted chemotherapy, yet, to date, there have been few detailed studies to systematically enhance drug loading while maintaining controlled release. In this work, we investigate which molecular simulation methods best capture the experimental uptake and release of cisplatin from UiO-66 and UiO-66(NH2). We then screen a series of biocompatible, pH-sensitive zeolitic imidazolate frameworks (ZIFs) for their ability to retain cisplatin in healthy parts of the patient and release it in the vicinity of a tumor. Pure-component GCMC simulations show that the maximum cisplatin loading depends on the pore volume. To achieve this maximum loading in the presence of water, either the pore size needs to be large enough to occupy both cisplatin and its solvation shell or the MOF-cisplatin interaction must be more favorable than the cisplatin-shell interaction. Both solvated and non-solvated simulations show that cisplatin release rates can be controlled by either decreasing the pore limiting diameters or by manipulating framework-cisplatin interaction energies to create strong, dispersed adsorption sites. The latter method is preferable if cisplatin loading is performed from solution into a pre-synthesized framework as weak interaction energies and small pore window diameters will hinder cisplatin uptake. Here, ZIF-82 is most promising. If it is possible to load cisplatin during crystallization, ZIF-11 would outcompete the other MOFs screened as cisplatin cannot pass through its pore windows; therefore, release rates would be purely driven by the pH triggered framework degradation.

Monoclonal antibody stability can be usefully monitored using the excitation-energy-dependent fluorescence edge-shift.

Among the major challenges in the development of biopharmaceuticals are structural heterogeneity and aggregation. The development of a successful therapeutic monoclonal antibody (mAb) requires both a highly active and also stable molecule. Whilst a range of experimental (biophysical) approaches exist to track changes in stability of proteins, routine prediction of stability remains challenging. The fluorescence red edge excitation shift (REES) phenomenon is sensitive to a range of changes in protein structure. Based on recent work, we have found that quantifying the REES effect is extremely sensitive to changes in protein conformational state and dynamics. Given the extreme sensitivity, potentially this tool could provide a 'fingerprint' of the structure and stability of a protein. Such a tool would be useful in the discovery and development of biopharamceuticals and so we have explored our hypothesis with a panel of therapeutic mAbs. We demonstrate that the quantified REES data show remarkable sensitivity, being able to discern between structurally identical antibodies and showing sensitivity to unfolding and aggregation. The approach works across a broad concentration range (µg-mg/ml) and is highly consistent. We show that the approach can be applied alongside traditional characterisation testing within the context of a forced degradation study (FDS). Most importantly, we demonstrate the approach is able to predict the stability of mAbs both in the short (hours), medium (days) and long-term (months). The quantified REES data will find immediate use in the biopharmaceutical industry in quality assurance, formulation and development. The approach benefits from low technical complexity, is rapid and uses instrumentation which exists in most biochemistry laboratories without modification.

TCR-induced alteration of primary MHC peptide anchor residue.

The HLA-A*02:01-restricted decapeptide EAAGIGILTV, derived from melanoma antigen recognized by T-cells-1 (MART-1) protein, represents one of the best-studied tumor associated T-cell epitopes, but clinical results targeting this peptide have been disappointing. This limitation may reflect the dominance of the nonapeptide, AAGIGILTV, at the melanoma cell surface. The decapeptide and nonapeptide are presented in distinct conformations by HLA-A*02:01 and TCRs from clinically relevant T-cell clones recognize the nonapeptide poorly. Here, we studied the MEL5 TCR that potently recognizes the nonapeptide. The structure of the MEL5-HLA-A*02:01-AAGIGILTV complex revealed an induced fit mechanism of antigen recognition involving altered peptide-MHC anchoring. This "flexing" at the TCR-peptide-MHC interface to accommodate the peptide antigen explains previously observed incongruences in this well-studied system and has important implications for future therapeutic approaches. Finally, this study expands upon the mechanisms by which molecular plasticity can influence antigen recognition by T cells.

Exposing the Interplay Between Enzyme Turnover, Protein Dynamics, and the Membrane Environment in Monoamine Oxidase B.

There is an increasing realization that structure-based drug design may show improved success by understanding the ensemble of conformations accessible to an enzyme and how the environment affects this ensemble. Human monoamine oxidase B (MAO-B) catalyzes the oxidation of amines and is inhibited for the treatment of both Parkinson's disease and depression. Despite its clinical importance, its catalytic mechanism remains unclear, and routes to drugging this target would be valuable. Evidence of a radical in either the transition state or the resting state of MAO-B is present throughout the literature and is suggested to be a flavin semiquinone, a tyrosyl radical, or both. Here we see evidence of a resting-state flavin semiquinone, via absorption redox studies and electron paramagnetic resonance, suggesting that the anionic semiquinone is biologically relevant. On the basis of enzyme kinetic studies, enzyme variants, and molecular dynamics simulations, we find evidence for the importance of the membrane environment in mediating the activity of MAO-B and that this mediation is related to the protein dynamics of MAO-B. Further, our MD simulations identify a hitherto undescribed entrance for substrate binding, membrane modulated substrate access, and indications for half-site reactivity: only one active site is accessible to binding at a time. Our study combines both experimental and computational evidence to illustrate the subtle interplay between enzyme activity and protein dynamics and the immediate membrane environment. Understanding key biomedical enzymes to this level of detail will be crucial to inform strategies (and binding sites) for rational drug design for these targets.

Intrinsic Flexibility of the EMT Zeolite Framework under Pressure.

The roles of organic additives in the assembly and crystallisation of zeolites are still not fully understood. This is important when attempting to prepare novel frameworks to produce new zeolites. We consider 18-crown-6 ether (18C6) as an additive, which has previously been shown to differentiate between the zeolite EMC-2 (EMT) and faujasite (FAU) frameworks. However, it is unclear whether this distinction is dictated by influences on the metastable free-energy landscape or geometric templating. Using high-pressure synchrotron X-ray diffraction, we have observed that the presence of 18C6 does not impact the EMT framework flexibility-agreeing with our previous geometric simulations and suggesting that 18C6 does not behave as a geometric template. This was further studied by computational modelling using solid-state density-functional theory and lattice dynamics calculations. It is shown that the lattice energy of FAU is lower than EMT, but is strongly impacted by the presence of solvent/guest molecules in the framework. Furthermore, the EMT topology possesses a greater vibrational entropy and is stabilised by free energy at a finite temperature. Overall, these findings demonstrate that the role of the 18C6 additive is to influence the free energy of crystallisation to assemble the EMT framework as opposed to FAU.

'Something in the way she moves': The functional significance of flexibility in the multiple roles of protein disulfide isomerase (PDI).

Protein disulfide isomerase (PDI) has diverse functions in the endoplasmic reticulum as catalyst of redox transfer, disulfide isomerization and oxidative protein folding, as molecular chaperone and in multi-subunit complexes. It interacts with an extraordinarily wide range of substrate and partner proteins, but there is only limited structural information on these interactions. Extensive evidence on the flexibility of PDI in solution is not matched by any detailed picture of the scope of its motion. A new rapid method for simulating the motion of large proteins provides detailed molecular trajectories for PDI demonstrating extensive changes in the relative orientation of its four domains, great variation in the distances between key sites and internal motion within the core ligand-binding domain. The review shows that these simulations are consistent with experimental evidence and provide insight into the functional capabilities conferred by the extensive flexible motion of PDI.

The flexibility and dynamics of protein disulfide isomerase.

We have studied the mobility of the multidomain folding catalyst, protein disulfide isomerase (PDI), by a coarse-graining approach based on flexibility. We analyze our simulations of yeast PDI (yPDI) using measures of backbone movement, relative positions and orientations of domains, and distances between functional sites. We find that there is interdomain flexibility at every interdomain junction but these show very different characteristics. The extent of interdomain flexibility is such that yPDI's two active sites can approach much more closely than is found in crystal structures-and indeed hinge motion to bring these sites into proximity is the lowest energy normal mode of motion of the protein. The flexibility predicted for yPDI (based on one structure) includes the other known conformation of yPDI and is consistent with (i) the mobility observed experimentally for mammalian PDI and (ii) molecular dynamics. We also observe intradomain flexibility and clear differences between the domains in their propensity for internal motion. Our results suggest that PDI flexibility enables it to interact with many different partner molecules of widely different sizes and shapes, and highlights considerable similarities of yPDI and mammalian PDI. Proteins 2016; 84:1776-1785. © 2016 Wiley Periodicals, Inc.

Characterization of folding cores in the cyclophilin A-cyclosporin A complex.

Determining the folding core of a protein yields information about its folding process and dynamics. The experimental procedures for identifying the amino acids that make up the folding core include hydrogen-deuterium exchange and Φ-value analysis and can be expensive and time consuming. Because of this, there is a desire to improve upon existing methods for determining protein folding cores theoretically. We have obtained HDX data for the complex of cyclophilin A with the immunosuppressant cyclosporin A. We compare these data, as well as literature values for uncomplexed cyclophilin A, to theoretical predictions using a combination of rigidity analysis and coarse-grained simulations of protein motion. We find that in this case, the most specific prediction of folding cores comes from a combined approach that models the rigidity of the protein using the first software suite and the dynamics of the protein using the froda tool.

Structures of mesophilic and extremophilic citrate synthases reveal rigidity and flexibility for function.

Citrate synthase (CS) catalyses the entry of carbon into the citric acid cycle and is highly-conserved structurally across the tree of life. Crystal structures of dimeric CSs are known in both "open" and "closed" forms, which differ by a substantial domain motion that closes the substrate-binding clefts. We explore both the static rigidity and the dynamic flexibility of CS structures from mesophilic and extremophilic organisms from all three evolutionary domains. The computational expense of this wide-ranging exploration is kept to a minimum by the use of rigidity analysis and rapid all-atom simulations of flexible motion, combining geometric simulation and elastic network modeling. CS structures from thermophiles display increased structural rigidity compared with the mesophilic enzyme. A CS structure from a psychrophile, stabilized by strong ionic interactions, appears to display likewise increased rigidity in conventional rigidity analysis; however, a novel modified analysis, taking into account the weakening of the hydrophobic effect at low temperatures, shows a more appropriate decreased rigidity. These rigidity variations do not, however, affect the character of the flexible dynamics, which are well conserved across all the structures studied. Simulation trajectories not only duplicate the crystallographically observed symmetric open-to-closed transitions, but also identify motions describing a previously unidentified antisymmetric functional motion. This antisymmetric motion would not be directly observed in crystallography but is revealed as an intrinsic property of the CS structure by modeling of flexible motion. This suggests that the functional motion closing the binding clefts in CS may be independent rather than symmetric and cooperative.

High-resolution NMR studies of structure and dynamics of human ERp27 indicate extensive interdomain flexibility.

ERp27 (endoplasmic reticulum protein 27.7 kDa) is a homologue of PDI (protein disulfide-isomerase) localized to the endoplasmic reticulum. ERp27 is predicted to consist of two thioredoxin-fold domains homologous with the non-catalytic b and b' domains of PDI. The structure in solution of the N-terminal b-like domain of ERp27 was solved using high-resolution NMR data. The structure confirms that it has the thioredoxin fold and that ERp27 is a member of the PDI family. (15)N-NMR relaxation data were obtained and ModelFree analysis highlighted limited exchange contributions and slow internal motions, and indicated that the domain has an average order parameter S(2) of 0.79. Comparison of the single-domain structure determined in the present study with the equivalent domain within full-length ERp27, determined independently by X-ray diffraction, indicated very close agreement. The domain interface inferred from NMR data in solution was much more extensive than that observed in the X-ray structure, suggesting that the domains flex independently and that crystallization selects one specific interdomain orientation. This led us to apply a new rapid method to simulate the flexibility of the full-length protein, establishing that the domains show considerable freedom to flex (tilt and twist) about the interdomain linker, consistent with the NMR data.

Exploring the energy landscapes of protein folding simulations with Bayesian computation.

Nested sampling is a Bayesian sampling technique developed to explore probability distributions localized in an exponentially small area of the parameter space. The algorithm provides both posterior samples and an estimate of the evidence (marginal likelihood) of the model. The nested sampling algorithm also provides an efficient way to calculate free energies and the expectation value of thermodynamic observables at any temperature, through a simple post processing of the output. Previous applications of the algorithm have yielded large efficiency gains over other sampling techniques, including parallel tempering. In this article, we describe a parallel implementation of the nested sampling algorithm and its application to the problem of protein folding in a Go-like force field of empirical potentials that were designed to stabilize secondary structure elements in room-temperature simulations. We demonstrate the method by conducting folding simulations on a number of small proteins that are commonly used for testing protein-folding procedures. A topological analysis of the posterior samples is performed to produce energy landscape charts, which give a high-level description of the potential energy surface for the protein folding simulations. These charts provide qualitative insights into both the folding process and the nature of the model and force field used.

Structure and in silico simulations of a cold-active esterase reveals its prime cold-adaptation mechanism.

Here we determined the structure of a cold active family IV esterase (EstN7) cloned from Bacillus cohnii strain N1. EstN7 is a dimer with a classical α/ß hydrolase fold. It has an acidic surface that is thought to play a role in cold-adaption by retaining solvation under changed water solvent entropy at lower temperatures. The conformation of the functionally important cap region is significantly different to EstN7's closest relatives, forming a bridge-like structure with reduced helical content providing greater access to the active site through more than one substrate access tunnel. However, dynamics do not appear to play a major role in cold adaption. Molecular dynamics at different temperatures, rigidity analysis, normal mode analysis and geometric simulations of motion confirm the flexibility of the cap region but suggest that the rest of the protein is largely rigid. Rigidity analysis indicates the distribution of hydrophobic tethers is appropriate to colder conditions, where the hydrophobic effect is weaker than in mesophilic conditions due to reduced water entropy. Thus, it is likely that increased substrate accessibility and tolerance to changes in water entropy are important for of EstN7's cold adaptation rather than changes in dynamics.

Early stages of phase selection in MOF formation observed in molecular Monte Carlo simulations.

Metal-organic frameworks (MOF) comprising metal nodes bridged by organic linkers show great promise because of their guest-specific gas sorption, separation, drug-delivery, and catalytic properties. The selection of metal node, organic linker, and synthesis conditions in principle offers engineered control over both structure and function. For MOFs to realise their potential and to become more than just promising materials, a degree of predictability in the synthesis and a better understanding of the self-assembly or initial growth processes is of paramount importance. Using cobalt succinate, a MOF that exhibits a variety of phases depending on synthesis temperature and ligand to metal ratio, as proof of concept, we present a molecular Monte Carlo approach that allows us to simulate the early stage of MOF assembly. We introduce a new Contact Cluster Monte Carlo (CCMC) algorithm which uses a system of overlapping "virtual sites" to represent the coordination environment of the cobalt and both metal-metal and metal-ligand associations. Our simulations capture the experimentally observed synthesis phase distinction in cobalt succinate at 348 K. To the best of our knowledge this is the first case in which the formation of different MOF phases as a function of composition is captured by unbiased molecular simulations. The CCMC algorithm is equally applicable to any system in which short-range attractive interactions are a dominant feature, including hydrogen-bonding networks, metal-ligand coordination networks, or the assembly of particles with "sticky" patches, such as colloidal systems or the formation of protein complexes.

Pressure-induced symmetry changes in body-centred cubic zeolites.

Previous work has shown a strong correlation between zeolite framework flexibility and the nature of structural symmetry and phase transitions. However, there is little experimental data regarding this relationship, in addition to how flexibility can be connected to the synthesis of these open-framework materials. This is of interest for the synthesis of novel zeolites, which require organic additives to permutate the resulting geometry and symmetry of the framework. Here, we have used high-pressure powder X-ray diffraction to study the three zeolites: Na-X, RHO and ZK-5, which can all be prepared using 18-crown-6 ether as an organic additive. We observe significant differences in how the occluded 18-crown-6 ether influences the framework flexibility-this being dependent on the geometry of the framework. We use these differences as an indicator to define the role of 18-crown-6 ether during zeolite crystallization. Furthermore, in conjunction with previous work, we predict that pressure-induced symmetry transitions are intrinsic to body-centred cubic zeolites. The high symmetry yields fewer degrees of freedom, meaning it is energetically favourable to lower the symmetry to facilitate further compression.

Fitting low-resolution cryo-EM maps of proteins using constrained geometric simulations.

Recent experimental advances in producing density maps from cryo-electron microscopy (cryo-EM) have challenged theorists to develop improved techniques to provide structural models that are consistent with the data and that preserve all the local stereochemistry associated with the biomolecule. We develop a new technique that maintains the local geometry and chemistry at each stage of the fitting procedure. A geometric simulation is used to drive the structure from some appropriate starting point (a nearby experimental structure or a modeled structure) toward the experimental density, via a set of small incremental motions. Structural motifs such as alpha-helices can be held rigid during the fitting procedure as the starting structure is brought into alignment with the experimental density. After validating this procedure on simulated data for adenylate kinase and lactoferrin, we show how cryo-EM data for two different GroEL structures can be fit using a starting x-ray crystal structure. We show that by incorporating the correct local stereochemistry in the modeling, structures can be obtained with effective resolution that is significantly higher than might be expected from the nominal cryo-EM resolution.

Structure and function of L-threonine-3-dehydrogenase from the parasitic protozoan Trypanosoma brucei revealed by X-ray crystallography and geometric simulations.

Two of the world's most neglected tropical diseases, human African trypanosomiasis (HAT) and Chagas disease, are caused by protozoan parasites of the genus Trypanosoma. These organisms possess specialized metabolic pathways, frequently distinct from those in humans, which have potential to be exploited as novel drug targets. This study elucidates the structure and function of L-threonine-3-dehydrogenase (TDH) from T. brucei, the causative pathogen of HAT. TDH is a key enzyme in the metabolism of L-threonine, and an inhibitor of TDH has been shown to have trypanocidal activity in the procyclic form of T. brucei. TDH is a nonfunctional pseudogene in humans, suggesting that it may be possible to rationally design safe and specific therapies for trypanosomiasis by targeting this parasite enzyme. As an initial step, the TDH gene from T. brucei was expressed and the three-dimensional structure of the enzyme was solved by X-ray crystallography. In multiple crystallographic structures, T. brucei TDH is revealed to be a dimeric short-chain dehydrogenase that displays a considerable degree of conformational variation in its ligand-binding regions. Geometric simulations of the structure have provided insight into the dynamic behaviour of this enzyme. Furthermore, structures of TDH bound to its natural substrates and known inhibitors have been determined, giving an indication of the mechanism of catalysis of the enzyme. Collectively, these results provide vital details for future drug design to target TDH or related enzymes.

Defining the flexibility window in ordered aluminosilicate zeolites.

The flexibility window in zeolites was originally identified using geometric simulation as a hypothetical property of SiO2 systems. The existence of the flexibility window in hypothetical structures may help us to identify those we might be able to synthesize in the future. We have previously found that the flexibility window in silicates is connected to phase transitions under pressure, structure amorphization and other physical behaviours and phenomena. We here extend the concept to ordered aluminosilicate systems using softer 'bar' constraints that permit additional flexibility around aluminium centres. Our experimental investigation of pressure-induced amorphization in sodalites is consistent with the results of our modelling. The softer constraints allow us to identify a flexibility window in the anomalous case of goosecreekite.

Evidence of Correlated Static Disorder in the Fenna-Matthews-Olson Complex.

Observation of excitonic quantum beats in photosynthetic antennae has prompted wide debate regarding the function of excitonic coherence in pigment-protein complexes. Much of this work focuses on the interactions of excitons with the femto-to-picosecond dynamical fluctuations of their environment. However, in experiments these effects can be masked by static disorder of the excited-state energies across ensembles, whose microscopic origins are challenging to predict. Here the excited-state properties of ∼2000 atom clusters of the Fenna-Matthews-Olson complex are simulated using a unique combination of linear-scaling density functional theory and constrained geometric dynamics. While slow, large amplitude protein motion leads to large variations in the Qy transitions of two pigments, we identify pigment-protein correlations that greatly reduce variations in the energy gap across the ensemble, which is consistent with experimental observations of suppressed inhomogeneous dephasing of quantum beats.

A complete thermodynamic analysis of enzyme turnover links the free energy landscape to enzyme catalysis.

Our understanding of how enzymes work is coloured by static structure depictions where the enzyme scaffold is presented as either immobile, or in equilibrium between well-defined static conformations. Proteins, however, exhibit a large degree of motion over a broad range of timescales and magnitudes and this is defined thermodynamically by the enzyme free energy landscape (FEL). The role and importance of enzyme motion is extremely contentious. Much of the challenge is in the experimental detection of so called 'conformational sampling' involved in enzyme turnover. Herein we apply combined pressure and temperature kinetics studies to elucidate the full suite of thermodynamic parameters defining an enzyme FEL as it relates to enzyme turnover. We find that the key thermodynamic parameters governing vibrational modes related to enzyme turnover are the isobaric expansivity term and the change in heat capacity for enzyme catalysis. Variation in the enzyme FEL affects these terms. Our analysis is supported by a range of biophysical and computational approaches that specifically capture information on protein vibrational modes and the FEL (all atom flexibility calculations, red edge excitation shift spectroscopy and viscosity studies) that provide independent evidence for our findings. Our data suggest that restricting the enzyme FEL may be a powerful strategy when attempting to rationally engineer enzymes, particularly to alter thermal activity. Moreover, we demonstrate how rational predictions can be made with a rapid computational approach.

Thermal stability, storage and release of proteins with tailored fit in silica.

Biological substances based on proteins, including vaccines, antibodies, and enzymes, typically degrade at room temperature over time due to denaturation, as proteins unfold with loss of secondary and tertiary structure. Their storage and distribution therefore relies on a "cold chain" of continuous refrigeration; this is costly and not always effective, as any break in the chain leads to rapid loss of effectiveness and potency. Efforts have been made to make vaccines thermally stable using treatments including freeze-drying (lyophilisation), biomineralisation, and encapsulation in sugar glass and organic polymers. Here for the first time we show that proteins can be enclosed in a deposited silica "cage", rendering them stable against denaturing thermal treatment and long-term ambient-temperature storage, and subsequently released into solution with their structure and function intact. This "ensilication" method produces a storable solid protein-loaded material without the need for desiccation or freeze-drying. Ensilication offers the prospect of a solution to the "cold chain" problem for biological materials, in particular for vaccines.