Christopher C. Roberts

DNA bending by a phantom protein

BACKGROUND Despite its stiffness, duplex DNA is extensively bent and folded during packaging and gene expression in biological systems. Modulation of the electrostatic repulsion between phosphates in the DNA backbone may be important in the bending of DNA by proteins. Here, we analyze the shape of DNA molecules that have been modified chemically to mimic the electrostatic consequences of a bound protein. RESULTS We have simulated salt bridges between DNA phosphates and cationic amino acid sidechains of a phantom protein by tethering ammonium cations to one face of the DNA helix. Tethered ammonium cations, but not neutral acetylated controls, induce DNA to bend toward its neutralized surface. CONCLUSIONS The shape of DNA molecules bearing a laterally-asymmetric distribution of tethered cations agrees qualitatively with theoretical predictions and with results previously obtained using neutral phosphate analogs. These data suggest principles that might be applied to the design of artificial DNA-bending proteins.

SYNTHESIS OF OLIGONUCLEOTIDES BEARING THE NON-STANDARD BASES ISO-C AND ISO-G. COMPARISON OF ISO-C-ISO-G, C-G AND U-A BASE-PAIR STABILITIES IN RNA/DNA DUPLEXES

Abstract The synthesis of iso -cytidine and 2′-deoxy- iso -guanosine phosphoramidites in protected form is described. The synthesis of complementary dodecanucleotides containing a central iso -C- iso -G base-pair was achieved. It is found that the iso -C- iso -G base-pair has comparable stability to a C-G base-pair.

Modeling of Enhanced Catalysis in Multienzyme Nanostructures: Effect of Molecular Scaffolds, Spatial Organization, and Concentration

Colocalized multistep enzymatic reaction pathways within biological catabolic and metabolic processes occur with high yield and specificity. Spatial organization on membranes or surfaces may be associated with increased efficiency of intermediate substrate transfer. Using a new Brownian dynamics package, GeomBD, we explored the geometric features of a surface-anchored enzyme system by parallel coarse-grained Brownian dynamics simulations of substrate diffusion over microsecond (μs) to millisecond (ms) time scales. We focused on a recently developed glucose oxidase (GOx), horseradish peroxidase (HRP), and DNA origami-scaffold enzyme system, where the H2O2 substrate of HRP is produced by GOx. The results revealed and explained a significant advantage in catalytic enhancement by optimizing interenzyme distance and orientation in the presence of the scaffold model. The planar scaffold colocalized the enzymes and provided a diffusive barrier that enhanced substrate transfer probability, becoming more relevant with increasing interenzyme distance. The results highlight the importance of protein geometry in the proper assessment of distance and orientation dependence on the probability of substrate transfer. They shed light on strategies for engineering multienzyme complexes and further investigation of enhanced catalytic efficiency for substrate diffusion between membrane-anchoring proteins.

Etched Glass Microarrays with Differential Resonance for Enhanced Contrast and Sensitivity of Surface Plasmon Resonance Imaging Analysis

We report the fabrication and characterization of gold-coated etched glass array substrates for surface plasmon resonance imaging (SPRi) analysis with significantly enhanced performance, in particular image contrast and sensitivity. The etching of the glass substrate induces a variation in the resonance condition and thus in the resonance angle between the etched wells and the surrounding area, leading to the isolation of the array spot resonance with a significant reduction of the background signal. FDTD simulations show arrays with large spots and minimal spot-to-spot spacing yield ideal differential resonance conditions, which are verified by experimental results. Simulations also indicate the etched well structure exhibits enhanced SPR electric field intensity by 3-fold as compared to standard planar gold chips. Changes in the bulk sensitivity of the etched arrays have been obtained at the 10(-4) RIU level based on image intensity difference. The strong image contrast allows for improved microarray imaging analysis with easily distinguished signals from background resonance. The etched array chips are demonstrated for SPRi detection of bacterial toxins through the coating of an ultrathin SiO(2) film for direct vesicle fusion that establishes a supported membrane-based biosensing interface. Protein detection with cholera toxin (CT) at 5 nM is obtained, making this chip one of the most sensitive SPR imaging substrates ever reported without a postbinding amplification scheme. Furthermore, the surface can be regenerated by Triton X-100 for repeated cycles of membrane formation, protein binding, and biomolecular removal. The reusability and enhanced performance of the etched glass array chips should find a broad range of applications, opening up new avenues for high-throughput SPR imaging detection with convenience and marked surface sensitivity.

Mechanisms of enhanced catalysis in enzyme‐DNA nanostructures revealed through molecular simulations and experimental analysis

Understanding and controlling the molecular interactions between enzyme substrates and DNA nanostructures has important implications in the advancement of enzyme–DNA technologies as solutions in biocatalysis. Such hybrid nanostructures can be used to create enzyme systems with enhanced catalysis by controlling the local chemical and physical environments and the spatial organization of enzymes. Here we have used molecular simulations with corresponding experiments to describe a mechanism of enhanced catalysis due to locally increased substrate concentrations. With a series of DNA nanostructures conjugated to horseradish peroxidase, we show that binding interactions between substrates and the DNA structures can increase local substrate concentrations. Increased local substrate concentrations in HRP(DNA) nanostructures resulted in 2.9‐ and 2.4‐fold decreases in the apparent Michaelis constants of tetramethylbenzidine and 4‐aminophenol, substrates of HRP with tunable binding interactions to DNA nanostructures with dissociation constants in the micromolar range. Molecular simulations and kinetic analysis also revealed that increased local substrate concentrations enhanced the rates of substrate association. Identification of the mechanism of increased local concentration of substrates in close proximity to enzymes and their active sites adds to our understanding of nanostructured biocatalysis from which we can develop guidelines for enhancing catalysis in rationally designed systems.

Agonist and antagonist binding to the nuclear vitamin D receptor: dynamics, mutation effects and functional implications

PurposeThe thermodynamically favored complex between the nuclear vitamin D receptor (VDR) and 1α,25(OH)2-vitamin D3 (1,25D3) triggers a shift in equilibrium to favor VDR binding to DNA, heterodimerization with the nuclear retinoid x receptor (RXR) and subsequent regulation of gene transcription. The key amino acids and structural requirements governing VDR binding to nuclear coactivators (NCoA) are well defined. Yet very little is understood about the internal changes in amino acid flexibility underpinning the control of ligand affinity, helix 12 conformation and function. Herein, we use molecular dynamics (MD) to study how the backbone and side-chain flexibility of the VDR differs when a) complexed to 1α,25(OH)2-vitamin D3 (1,25D3, agonist) and (23S),25-dehydro-1α(OH)-vitamin D3-26,23-lactone (MK, antagonist); b) residues that form hydrogen bonds with the C25-OH (H305 and H397) of 1,25D3 are mutated to phenylalanine; c) helix 12 conformation is changed and ligand is removed; and d) x-ray water near the C1- and C3-OH groups of 1,25D3 are present or replaced with explicit solvent.MethodsWe performed molecular dynamic simulations on the apo- and holo-VDRs and used T-Analyst to monitor the changes in the backbone and side-chain flexibility of residues that form regions of the VDR ligand binding pocket (LBP), NCoA surface and control helix 12 conformation.ResultsThe VDR-1,25D3 and VDR-MK MD simulations demonstrate that 1,25D3 and MK induce highly similar changes in backbone and side-chain flexibility in residues that form the LBP. MK however did increase the backbone and side-chain flexibility of L404 and R274 respectively. MK also induced expansion of the VDR charge clamp (i.e. NCoA surface) and weakened the intramolecular interaction between H305---V418 (helix 12) and TYR401 (helix 11). In VDR_FF, MK induced a generally more rigid LBP and stronger interaction between F397 and F422 than 1,25D3, and reduced the flexibility of the R274 side-chain. Lastly the VDR MD simulations indicate that R274 can sample multiple conformations in the presence of ligand. When the R274 is extended, the β-OH group of 1,25D3 lies proximal to the backbone carbonyl oxygen of R274 and the side-chain forms H-bonds with hinge domain residues. This differs from the x-ray, kinked geometry, where the side-chain forms an H-bond with the 1α-OH group. Furthermore, 1,25D3, but not MK was observed to stabilize the x-ray geometry of R274 during the > 30 ns MD runs.ConclusionsThe MD methodology applied herein provides an in silico foundation to be expanded upon to better understand the intrinsic flexibility of the VDR and better understand key side-chain and backbone movements involved in the bimolecular interaction between the VDR and its’ ligands.

Immunoproteasome inhibition and bioactivity of thiasyrbactins

A family of macrodilactam natural products, the syrbactins, are known proteasome inhibitors. A small group of syrbactin analogs was prepared with a sulfur-for-carbon substitution to enhance synthetic accessibility and facilitate modulation of their solubility. Two of these compounds surprisingly proved to be inhibitors of the trypsin-like catalytic site, including of the immunoproteasome. Their bound and free conformations suggest special properties of the thiasyrbactin ring are responsible for this unusual preference, which may be exploited to develop drug-like immunoproteasome inhibitors. These compounds show greater selectivity than earlier compounds used to infer phenotypes of immunoproteasome inhibition, like ONX-0914.

Analysis of Ligand–Receptor Association and Intermediate Transfer Rates in Multienzyme Nanostructures with All-Atom Brownian Dynamics Simulations

We present the second-generation GeomBD Brownian dynamics software for determining interenzyme intermediate transfer rates and substrate association rates in biomolecular complexes. Substrate and intermediate association rates for a series of enzymes or biomolecules can be compared between the freely diffusing disorganized configuration and various colocalized or complexed arrangements for kinetic investigation of enhanced intermediate transfer. In addition, enzyme engineering techniques, such as synthetic protein conjugation, can be computationally modeled and analyzed to better understand changes in substrate association relative to native enzymes. Tools are provided to determine nonspecific ligand-receptor association residence times, and to visualize common sites of nonspecific association of substrates on receptor surfaces. To demonstrate features of the software, interenzyme intermediate substrate transfer rate constants are calculated and compared for all-atom models of DNA origami scaffold-bound bienzyme systems of glucose oxidase and horseradish peroxidase. Also, a DNA conjugated horseradish peroxidase enzyme was analyzed for its propensity to increase substrate association rates and substrate local residence times relative to the unmodified enzyme. We also demonstrate the rapid determination and visualization of common sites of nonspecific ligand-receptor association by using HIV-1 protease and an inhibitor, XK263. GeomBD2 accelerates simulations by precomputing van der Waals potential energy grids and electrostatic potential grid maps, and has a flexible and extensible support for all-atom and coarse-grained force fields. Simulation software is written in C++ and utilizes modern parallelization techniques for potential grid preparation and Brownian dynamics simulation processes. Analysis scripts, written in the Python scripting language, are provided for quantitative simulation analysis. GeomBD2 is applicable to the fields of biophysics, bioengineering, and enzymology in both predictive and explanatory roles.

Syrbactin Structural Analog TIR-199 Blocks Proteasome Activity and Induces Tumor Cell Death

Multiple myeloma is an aggressive hematopoietic cancer of plasma cells. The recent emergence of three effective FDA-approved proteasome-inhibiting drugs, bortezomib (Velcade®), carfilzomib (Kyprolis®), and ixazomib (Ninlaro®), confirms that proteasome inhibitors are therapeutically useful against neoplastic disease, in particular refractory multiple myeloma and mantle cell lymphoma. This study describes the synthesis, computational affinity assessment, and preclinical evaluation of TIR-199, a natural product-derived syrbactin structural analog. Molecular modeling and simulation suggested that TIR-199 covalently binds each of the three catalytic subunits (β1, β2, and β5) and revealed key interaction sites. In vitro and cell culture-based proteasome activity measurements confirmed that TIR-199 inhibits the proteasome in a dose-dependent manner and induces tumor cell death in multiple myeloma and neuroblastoma cells as well as other cancer types in the NCI-60 cell panel. It is particularly effective against kidney tumor cell lines, with >250-fold higher anti-tumor activities than observed with the natural product syringolin A. In vivo studies in mice revealed a maximum tolerated dose of TIR-199 at 25 mg/kg. The anti-tumor activity of TIR-199 was confirmed in hollow fiber assays in mice. Adverse drug reaction screens in a kidney panel revealed no off-targets of concern. This is the first study to examine the efficacy of a syrbactin in animals. Taken together, the results suggest that TIR-199 is a potent new proteasome inhibitor with promise for further development into a clinical drug for the treatment of multiple myeloma and other forms of cancer.

Ligand Binding Pathway Elucidation for Cryptophane Host-Guest Complexes.

Modeling binding pathways can provide insight into molecular recognition, including kinetic mechanisms, barriers to binding, and gating effects. This work represents a novel computational approach, Hopping Minima, for the determination of conformational transitions of single molecules as well as binding pathways for molecular complexes. The method begins by thoroughly sampling a set of conformational minima for a molecular system. The natural motions of the system are modeled using the normal modes of the sampled minima. The natural motions are utilized to connect conformational minima and are finally combined to form association/binding pathways in the case of molecular complexes. We provide an implementation and example application of the method using alanine dipeptide and a set of chemical host-guest systems: two cryptophane hosts with two guest cations, trimethylammonium and tetramethylammonium. Our results demonstrate that conformational transitions can be modeled and extended to find binding pathways as well as energetic information relevant to the minimum conformations involved. This approach has advantages over simulation-based methods for studying systems with slow binding processes and can help design molecules with preferred binding kinetics.