Chris M. W. Ho

Cavity search: An algorithm for the isolation and display of cavity-like binding regions

SummaryA set of algorithms designed to enhance the display of protein binding cavities is presented. These algorithms, collectively entitled CAVITY SEARCH, allow the user to isolate and fully define the extent of a particular cavity. Solid modeling techniques are employed to produce a detailed cast of the active site region, which can then be color-coded to show electrostatic and steric interactions between the protein cavity and a bound ligand.

FOUNDATION: A program to retrieve all possible structures containing a user-defined minimum number of matching query elements from three-dimensional databases

SummaryA program is described that searches three-dimensional, structural databases, given a user-defined query, in order to retrieve all structures that contain any combination of a user-specified minimum number of matching elements. Queries consist of three-dimensional coordinates of atoms and/or bonds. Numerous query constraints are described which allow the investigator to define the chemical nature of the desired structures as well as the environment within which these structures must reside. They include:(1)Bonded vs. isolated atom distinction;(2)Atom type designation;(3)Definition of subsets with occupancy specification (>, =, < X atoms);(4)RMS-fit;(5)Active site volume accessibility of atoms linking query elements;(6)Number, atom type, and cyclic structure constraints for atoms linking pharmacophoric elements;(7)Automatic error boundary adjustment-ad infinitum constraint. To illustrate the capabilities of this program, queries based on the crystal structure of a thermolysin-inhibitor complex were tested against a subset of the Cambridge Crystallographic Database. Several compounds were returned which satisfied various aspects of the query, including fitting within the active site. Combination of segments of compounds which satisfy partial queries should provide a method for generating unique compounds with affinity for sites of known three-dimensional structure.

SPLICE: A program to assemble partial query solutions from three-dimensional database searches into novel ligands

SummarySPLICE is a program that processes partial query solutions retrieved from 3D, structural databases to generate novel, aggregate ligands. It is designed to interface with the database searching program FOUNDATION, which retrieves fragments containing any combination of a user-specified minimum number of matching query elements. SPLICE eliminates aspects of structures that are physically incapable of binding within the active site. Then, a systematic rule-based procedure is performed upon the remaining fragments to ensure receptor complementarity. All modifications are automated and remain transparent to the user. Ligands are then assembled by linking components into composite structures through overlapping bonds. As a control experiment, FOUNDATION and SPLICE were used to reconstruct a know HIV-1 protease inhibitor after it had been fragmented, reoriented, and added to a sham database of fifty different small molecules. To illustrate the capabilities of this program, a 3D search query containing the pharmacophoric elements of an aspartic proteinase-inhibitor crystal complex was searched using FOUNDATION against a subset of the Cambridge Structural Database. One hundred thirty-one compounds were retrieved, each containing any combination of at least four query elements. Compounds were automatically screened and edited for receptor complementarity. Numerous combinations of fragments were discovered that could be linked to form novel structures, containing a greater number of pharmacophoric elements than any single retrieved fragment.

Modelling proteins’ hidden conformations to predict antibiotic resistance

TEM β-lactamase confers bacteria with resistance to many antibiotics and rapidly evolves activity against new drugs. However, functional changes are not easily explained by differences in crystal structures. We employ Markov state models to identify hidden conformations and explore their role in determining TEM’s specificity. We integrate these models with existing drug-design tools to create a new technique, called Boltzmann docking, which better predicts TEM specificity by accounting for conformational heterogeneity. Using our MSMs, we identify hidden states whose populations correlate with activity against cefotaxime. To experimentally detect our predicted hidden states, we use rapid mass spectrometric footprinting and confirm our models’ prediction that increased cefotaxime activity correlates with reduced Ω-loop flexibility. Finally, we design novel variants to stabilize the hidden cefotaximase states, and find their populations predict activity against cefotaxime in vitro and in vivo. Therefore, we expect this framework to have numerous applications in drug and protein design.

Conformationally restricted TRH analogues: constraining the pyroglutamate region.

A modified synthetic route has been developed so that the steric size of constraints added to the pyroglutamate region of TRH (pGluHisProNH(2)) can be varied. Both an analogue with a smaller ethylene bridge and a larger, more flexible propane bridge in this region have been synthesized. These analogues were synthesized in order to probe why the initial incorporation of an ethane bridge into this region of the molecule had led to an analogue with a binding constant and potency three times lower than that of an directly analogous unconstrained analogue. The data for both analogues indicated that the fall off in activity caused by the ethane bridge in the initial analogue was not caused by the size of the bridge.

Designing small molecules to target cryptic pockets yields both positive and negative allosteric modulators

Allosteric drugs, which bind to proteins in regions other than their main ligand-binding or active sites, make it possible to target proteins considered “undruggable” and to develop new therapies that circumvent existing resistance. Despite growing interest in allosteric drug discovery, rational design is limited by a lack of sufficient structural information about alternative binding sites in proteins. Previously, we used Markov State Models (MSMs) to identify such “cryptic pockets,” and here we describe a method for identifying compounds that bind in these cryptic pockets and modulate enzyme activity. Experimental tests validate our approach by revealing both an inhibitor and two activators of TEM β-lactamase (TEM). To identify hits, a library of compounds is first virtually screened against either the crystal structure of a known cryptic pocket or an ensemble of structures containing the same cryptic pocket that is extracted from an MSM. Hit compounds are then screened experimentally and characterized kinetically in individual assays. We identify three hits, one inhibitor and two activators, demonstrating that screening for binding to allosteric sites can result in both positive and negative modulation. The hit compounds have modest effects on TEM activity, but all have higher affinities than previously identified inhibitors, which bind the same cryptic pocket but were found, by chance, via a computational screen targeting the active site. Site-directed mutagenesis of key contact residues predicted by the docking models is used to confirm that the compounds bind in the cryptic pocket as intended. Because hit compounds are identified from docking against both the crystal structure and structures from the MSM, this platform should prove suitable for many proteins, particularly targets whose crystal structures lack obvious druggable pockets, and for identifying both inhibitory and activating small-molecule modulators.

Defining NADH-Driven Allostery Regulating Apoptosis-Inducing Factor.

Apoptosis-inducing factor (AIF) is critical for mitochondrial respiratory complex biogenesis and for mediating necroptotic parthanatos; these functions are seemingly regulated by enigmatic allosteric switching driven by NADH charge-transfer complex (CTC) formation. Here, we define molecular pathways linking AIFs active site to allosteric switching regions by characterizing dimer-permissive mutants using small-angle X-ray scattering (SAXS) and crystallography and by probing AIF-CTC communication networks using molecular dynamics simulations. Collective results identify two pathways propagating allostery from the CTC active site: (1) active-site H454 links to S480 of AIFs central β-strand to modulate a hydrophobic border at the dimerization interface, and (2) an interaction network links AIFs FAD cofactor, central β-strand, and Cβ-clasp whereby R529 reorientation initiates C-loop release during CTC formation. This knowledge of AIF allostery and its flavoswitch mechanism provides a foundation for biologically understanding and biomedically controlling its participation in mitochondrial homeostasis and cell death.

Automated Solid Models From Serial Section Images

A new method for creating unambiguous and complete boundary representation solid models with a hybrid polygonal/nonuniform rational B spline representation was developed and tested using computed tomography scans of the wrist. Polygon surface approximation was applied to a sequence of parallel planar outlines of individual bone elements in the wrist. An automated technique for the transformation of edge contours into solid models was implemented. This was performed using a custom batch file command sequence generator coupled to a commercially available mechanical computer-aided design and engineering software system known as I-DEAS (Structural Dynamics Research Corporation, Milford, OH). This transformation software allows the use of biomedical scan slice data with a solid modeler.

Investigating the effect of forestry on leaf-litter arthropods (Algonquin Park, Ontario, Canada)

Arthropods are the most diverse taxonomic group of terrestrial eukaryotes and are sensitive to physical alterations in their environment such as those caused by forestry. With their enormous diversity and physical omnipresence, arthropods could be powerful indicators of the effects of disturbance following forestry. When arthropods have been used to measure the effects of disturbance, the total diversity of some groups is often found to increase following forestry. However, these findings are frequently derived using a coarse taxonomic grain (family or order) to accommodate for various taxonomic impediments (including cryptic diversity and poorly resourced taxonomists). Our intent with this work was to determine the diversity of arthropods in and around Algonquin Park, and how this diversity was influenced by disturbance (in this case, forestry within the past 25 years). We used DNA barcode-derived diversity estimates (Barcode Index Number (BIN) richness) to avoid taxonomic impediments and as a source of genetic information with which we could conduct phylogenetic estimates of diversity (PD). Diversity patterns elucidated with PD are often, but not always congruent with taxonomic estimates–and departures from these expectations can help clarify disturbance effects that are hidden from richness studies alone. We found that BIN richness and PD were greater in disturbed (forested) areas, however when we controlled for the expected relationship between PD and BIN richness, we found that cut sites contained less PD than expected and that this diversity was more phylogenetically clustered than would be predicted by taxonomic richness. While disturbance may cause an evident increase in diversity, this diversity may not reflect the full evolutionary history of the assemblage within that area and thus a subtle effect of disturbance can be found decades following forestry.

Design, synthesis, and biological evaluation of stable β 6.3 -Helices: Discovery of non-hemolytic antibacterial peptides

Gramicidin A, a topical antibiotic made from alternating L and D amino acids, is characterized by its wide central pore; upon insertion into membranes, it forms channels that disrupts ion gradients. We present helical peptidomimetics with this characteristic wide central pore that have been designed to mimic gramicidin A channels. Mimetics were designed using molecular modeling focused on oligomers of heterochiral dipeptides of proline analogs, in particular azaproline (AzPro). Molecular Dynamics simulations in water confirmed the stability of the designed helices. A sixteen-residue Formyl-(AzPro-Pro)8-NHCH2CH2OH helix was synthesized as well as a full thirty-two residue Cbz-(AzPro-Pro)16-OtBu channels. No liposomal lysis activity was observed suggesting lack of channel formation, possibly due to inappropriate hydrogen-bonding interactions in the membrane. These peptidomimetics also did not hemolyze red blood cells, unlike gramicidin A.