Molecular modeling as a predictive tool for the development of solid dispersions.

In this study molecular modeling is introduced as a novel approach for the development of pharmaceutical solid dispersions. A computational model based on quantum mechanical (QM) calculations was used to predict the miscibility of various drugs in various polymers by predicting the binding strength between the drug and dimeric form of the polymer. The drug/polymer miscibility was also estimated by using traditional approaches such as Van Krevelen/Hoftyzer and Bagley solubility parameters or Flory-Huggins interaction parameter in comparison to the molecular modeling approach. The molecular modeling studies predicted successfully the drug-polymer binding energies and the preferable site of interaction between the functional groups. The drug-polymer miscibility and the physical state of bulk materials, physical mixtures, and solid dispersions were determined by thermal analysis (DSC/MTDSC) and X-ray diffraction. The produced solid dispersions were analyzed by X-ray photoelectron spectroscopy (XPS), which confirmed not only the exact type of the intermolecular interactions between the drug-polymer functional groups but also the binding strength by estimating the N coefficient values. The findings demonstrate that QM-based molecular modeling is a powerful tool to predict the strength and type of intermolecular interactions in a range of drug/polymeric systems for the development of solid dispersions.

Quantum mechanics/molecular mechanics studies on the mechanism of action of cofactor pyridoxal 5'-phosphate in ornithine 4,5-aminomutase.

A computational study was performed on the experimentally elusive cyclisation step in the cofactor pyridoxal 5'-phosphate (PLP)-dependent D-ornithine 4,5-aminomutase (OAM)-catalysed reaction. Calculations using both model systems and a combined quantum mechanics/molecular mechanics approach suggest that regulation of the cyclic radical intermediate is achieved through the synergy of the intrinsic catalytic power of cofactor PLP and the active site of the enzyme. The captodative effect of PLP is balanced by an enzyme active site that controls the deprotonation of both the pyridine nitrogen atom (N1) and the Schiff-base nitrogen atom (N2). Furthermore, electrostatic interactions between the terminal carboxylate and amino groups of the substrate and Arg297 and Glu81 impose substantial "strain" energy on the orientation of the cyclic intermediate to control its trajectory. In addition the "strain" energy, which appears to be sensitive to both the number of carbon atoms in the substrate/analogue and the position of the radical intermediates, may play a key role in controlling the transition of the enzyme from the closed to the open state. Our results provide new insights into several aspects of the radical mechanism in aminomutase catalysis and broaden our understanding of cofactor PLP-dependent reactions.

Reaction of vascular adhesion protein-1 (VAP-1) with primary amines: mechanistic insights from isotope effects and quantitative structure-activity relationships.

Human vascular adhesion protein-1 (VAP-1) is an endothelial copper-dependent amine oxidase involved in the recruitment and extravasation of leukocytes at sites of inflammation. VAP-1 is an important therapeutic target for several pathological conditions. We expressed soluble VAP-1 in HEK293 EBNA1 cells at levels suitable for detailed mechanistic studies with model substrates. Using the model substrate benzylamine, we analyzed the steady-state kinetic parameters of VAP-1 as a function of solution pH. We found two macroscopic pK(a) values that defined a bell-shaped plot of turnover number k(cat,app) as a function of pH, representing ionizable groups in the enzyme-substrate complex. The dependence of (k(cat)/K(m))(app) on pH revealed a single pK(a) value (∼9) that we assigned to ionization of the amine group in free benzylamine substrate. A kinetic isotope effect (KIE) of 6 to 7.6 on (k(cat)/K(m))(app) over the pH range of 6 to 10 was observed with d(2)-benzylamine. Over the same pH range, the KIE on k(cat) was found to be close to unity. The unusual KIE values on (k(cat)/K(m))(app) were rationalized using a mechanistic scheme that includes the possibility of multiple isotopically sensitive steps. We also report the analysis of quantitative structure-activity relationships (QSAR) using para-substituted protiated and deuterated phenylethylamines. With phenylethylamines we observed a large KIE on k(cat,app) (8.01 ± 0.28 with phenylethylamine), indicating that C-H bond breakage is limiting for 2,4,5-trihydroxyphenylalanine quinone reduction. Poor correlations were observed between steady-state rate constants and QSAR parameters. We show the importance of combining KIE, QSAR, and structural studies to gain insight into the complexity of the VAP-1 steady-state mechanism.

Temperature-dependent study reveals that dynamics of hydrophobic residues plays an important functional role in the mitochondrial Tim9-Tim10 complex.

Protein-protein interaction is a fundamental process in all major biological processes. The hexameric Tim9-Tim10 (translocase of inner membrane) complex of the mitochondrial intermembrane space plays an essential chaperone-like role during import of mitochondrial membrane proteins. However, little is known about the functional mechanism of the complex because the interaction is weak and transient. This study investigates how electrostatic and hydrophobic interactions affect the conformation and function of the complex at physiological temperatures, using both experimental and computational methods. The results suggest that, first, different complex conformational states exist at equilibrium, and the major difference between these states is the degree of hydrophobic interactions. Second, the conformational change mimics the biological activity of the complex as measured by substrate binding at the same temperatures. Finally, molecular dynamics simulation and detailed energy decomposition analysis provided supporting evidence at the atomic level for the presence of an excited state of the complex, the formation of which is largely driven by the disruption of hydrophobic interactions. Taken together, this study indicates that the dynamics of the hydrophobic residues plays an important role in regulating the function of the Tim9-Tim10 complex.

Large-scale domain conformational change is coupled to the activation of the Co-C bond in the B12-dependent enzyme ornithine 4,5-aminomutase: a computational study.

We present here an energetic and atomistic description of how D-ornithine 4,5-aminomutase (OAM), an adenosylcobalamin (AdoCbl; coenzyme B(12))-dependent isomerase, employs a large-scale protein domain conformational change to orchestrate the homolytic rupture of the Co-C bond. Our results suggest that in going from the open form (catalytically inactive) to the closed form (catalytically active), the Rossmann domain of OAM effectively approaches the active site as a rigid body. It undergoes a combination of a ~52° rotation and a ~14 Å translation to bring AdoCbl-initially positioned ~25 Å away-into the active-site cavity. This process is coupled to repositioning of the Ado moiety of AdoCbl from the eastern conformation to the northern conformation. Combined quantum mechanics and molecular mechanics calculations further indicate that in the open form, the protein environment does not impact significantly on the Co-C bond homolytic rupture, rendering it unusually stable, and thus catalytically inactive. Upon formation of the closed form, the Co-C bond is activated through the synergy of steric and electrostatic effects arising from tighter interactions with the surrounding enzyme. The more pronounced effect of the protein in the closed form gives rise to an elongated Co-C bond (by 0.03 Å), puckering of the ribose and increased "strain" energy on the Ado group and to a lesser extent the corrin ring. Our computational studies reveal novel strategies employed by AdoCbl-dependent enzymes in the control of radical catalysis.

Manipulations of phenylnorbornyl palladium species for multicomponent construction of a bridged polycyclic privileged scaffold.

Hexahydromethanocarbazole is a privileged scaffold in the discovery of new drugs and photoactive organic materials due to its good balance between structural complexity and minimized entropy penalty upon receptor binding. To address the difficulty of synthesizing this highly desirable bridged polycyclic scaffold, we designed a convenient multicomponent reaction cascade as intercepted Heck addition/C-H activation/C-palladacycle formation/electrophilic attack of ANP/N-palladacycle formation/Buchwald amination. A distinguishing feature of this sophisticated strategy is the successive generation of two key phenylnorbornyl palladium species to control the reaction flow towards desired products. DFT calculations further reveal the crucial roles of Cs2CO3 and 5,6-diester substitutions on the norbornene reactant in preventing multiple side-reactions. This innovative method exhibits a broad scope with good yields, and therefore will enable the construction of natural-product-like compound libraries based on hexahydromethanocarbazole.

Delivering Antisense Oligonucleotides across the Blood-Brain Barrier by Tumor Cell-Derived Small Apoptotic Bodies.

The blood-brain barrier (BBB) is the most restrictive and complicated barrier that keeps most biomolecules and drugs from the brain. An efficient brain delivery strategy is urgently needed for the treatment of brain diseases. Based on the studies of brain-targeting extracellular vesicles (EVs), the potential of using small apoptotic bodies (sABs) from brain metastatic cancer cells for brain-targeting drug delivery is explored. It is found that anti-TNF-α antisense oligonucleotide (ASO) combined with cationic konjac glucomannan (cKGM) can be successfully loaded into sABs via a transfection/apoptosis induction process and that the sABs generated by B16F10 cells have an extraordinarily high brain delivery efficiency. Further studies suggest that ASO-loaded sABs (sCABs) are transcytosed by b. End3 (brain microvascular endothelial cells, BMECs) to penetrate the BBB, which is mediated by CD44v6, and eventually taken up by microglial cells in the brain. In a Parkinson's disease (PD) mouse model, sCABs dramatically ameliorate PD symptoms via the anti-inflammatory effect of ASO. This study suggests that sABs from brain metastatic cancer cells are excellent carriers for brain-targeted delivery, as they have not only an extraordinary delivery efficiency but also a much higher scale-up production potential than other EVs.

Assignment of the vibrational spectra of enzyme-bound tryptophan tryptophyl quinones using a combined QM/MM approach.

Fourier transform infrared (FTIR) spectroscopy can be used to provide a detailed time-resolved probe of reaction intermediates in enzyme-catalyzed systems. Accurate assignment of the respective chemical species being studied is key to the success of this approach. The plethora of signals from the protein environment, leading to complexity in the spectra, presents a particular challenge. Here we present a combined QM/MM-based approach that can be used to assign key resonances in the FTIR spectrum of tryptophan tryptophyl quinone (TTQ) in the TTQ-dependent quinoprotein aromatic amine dehydrogenase (AADH). We show that consideration of the cofactor alone is not sufficient to identify correctly the experimentally observed resonances-inclusion of the protein is required for this. However, to enable accurate peak assignment, a stepwise approach is needed that builds up increasing levels of complexity from a simple system. This study serves as a benchmark for future QM/MM-based studies to predict the spectroscopic changes during the interconversion of intermediates in the reductive half-reaction catalyzed by AADH, and more generally for using a combined QM/MM approach to calculate spectroscopic data of protein cofactors and cofactor-based adducts.

Evidence to support the hypothesis that promoting vibrations enhance the rate of an enzyme catalyzed H-tunneling reaction.

In recent years there has been a shift away from transition state theory models for H-transfer reactions. Models that incorporate tunneling as the mechanism of H-transfer are now recognized as a better description of such reactions. Central to many models of H-tunneling is the notion that specific vibrational modes of the protein and/or substrate can increase the probability of a H-tunneling reaction, modes that are termed promoting vibrations. Thus far there has been limited evidence that promoting vibrations can increase the rate of H-transfer. In the present communication we examine the single hydride transfer from both NADPH and NADH to FMN in the reductive half-reaction of pentaerythritol tetranitrate reductase (PETNR). We find that there is a significant promoting vibration with NADPH but not with NADH and that the observed rate of hydride transfer is significantly (approximately 15x) faster with NADPH. We rule out differences in rate due to variation in driving force and the donor-acceptor distance, suggesting it is the promoting vibration with NADPH that is the origin of the increased observed rate. This study therefore provides direct evidence that promoting vibrations can lead to an increase in rate.

Parallel pathways and free-energy landscapes for enzymatic hydride transfer probed by hydrostatic pressure.

Mutation of an active-site residue in morphinone reductase leads to a conformationally rich landscape that enhances the rate of hydride transfer from NADH to FMN at standard pressure (1 bar). Increasing the pressure causes interconversion between different conformational substates in the mutant enzyme. While high pressure reduces the donor-acceptor distance in the wild-type enzyme, increased conformational freedom "dampens" its effect in the mutant.We show that hydride transfer from NADH to FMN catalysed by the N189A mutant of morphinone reductase occurs along parallel "chemical" pathways in a conformationally rich free-energy landscape. We have developed experimental kinetic and spectroscopic tools by using hydrostatic pressure to explore this free-energy landscape. The crystal structure of the N189A mutant enzyme in complex with the unreactive coenzyme analogue NADH(4) indicates that the nicotinamide moiety of the analogue is conformationally less restrained than the corresponding structure of the wild-type NADH(4) complex. This increased degree of conformational freedom in the N189A enzyme gives rise to the concept of multiple reactive configurations (MRCs), and we show that the relative population of these states across the free-energy landscape can be perturbed experimentally as a function of pressure. Specifically, the amplitudes of individual kinetic phases that were observed in stopped-flow studies of the hydride transfer reaction are sensitive to pressure; this indicates that pressure drives an altered distribution across the energy landscape. We show by absorbance spectroscopy that the loss of charge-transfer character of the enzyme-coenzyme complex is attributed to the altered population of MRCs on the landscape. The existence of a conformationally rich landscape in the N189A mutant is supported by molecular dynamics simulations at low and high pressure. The work provides firm experimental and computational support for the existence of parallel pathways arising from multiple conformational states of the enzyme-coenzyme complex. Hydrostatic pressure is a powerful and general probe of multidimensional energy landscapes that can be used to analyse experimentally parallel pathways for enzyme-catalysed reactions. We suggest that this is especially the case following directed mutation of a protein, which can lead to increased population of reactant states that are essentially inaccessible in the free-energy landscape of wild-type enzyme.

Barrier compression enhances an enzymatic hydrogen-transfer reaction.

Putting the squeeze on: Hydrostatic pressure causes a shortening of the charge-transfer bond in the binary complex of morphinone reductase and NADH(4) (see diagram). Molecular dynamics simulations suggest that pressure reduces the average reaction barrier width by restricting the conformational space available to the flavin mononucleotide and NADH within the active site. The apparent rate of catalysis increases with pressure.

Development of boronic acid-functionalized mesoporous silica-coated core/shell magnetic microspheres with large pores for endotoxin removal.

Endotoxins are found almost everywhere and possess high toxicity in vivo and in vitro. Here we design a novel boronate affinity material, called boronic acid-functionalized mesoporous silica-coated core/shell magnetic microspheres (Fe3O4@nSiO2@mSiO2-BA) with large pores (pore size > 20 nm) based on the chemical structure and physical properties of endotoxins, for facile and highly efficient removal of endotoxins. Dual modes for endotoxin removal were proposed and confirmed in this work: the endotoxin aggregates with size < 20 nm were bound with boronic acid ligands chemically modified on the inner and outer surface of the large pores of Fe3O4@nSiO2@mSiO2-BA microspheres; while the larger endotoxin micelles (size >20 nm) were absorbed on the outer surface of the prepared material based on boronate affinity. Transmission electron microscopy (TEM), X-ray diffraction (XRD), nitrogen adsorption/desorption isotherms and Fourier transform infrared (FT-IR) spectroscopy confirm that Fe3O4@nSiO2@mSiO2-BA microspheres possess core/shell structure, uniform diameter (520 nm), high surface area (205.57 m2/g), large mesopores (21.8 nm) and boronic acid ligands. The purification procedures of Fe3O4@nSiO2@mSiO2-BA microspheres for endotoxin were optimized, and 50 mM NH4HCO3 (pH 8.0) and 0.05 M fructose were selected as loading/washing, elution buffers, respectively. The binding capacity of Fe3O4@nSiO2@mSiO2-BA microspheres for endotoxin was calculated to be 60.84 EU/g under the optimized conditions. Finally, the established analytical method was applied to remove endotoxins from plasmid DNA. After endotoxin removal, the endotoxin content in plasmid DNA was reduced from 0.0026 to 0.0006 EU/mL for two-fold concentration, and from 0.0088 to 0.0022 EU/mL for five-fold concentration after binding, respectively. Additional advantages of the prepared boronate affinity material include excellent stability, reusability/repeatability, and low cost. Boronate affinity materials with large pores could thus prove to be powerful adsorbents for endotoxin removal and the potential applications in the aspects of biological research, pharmaceutical industry, and life health.

Deep tunneling dominates the biologically important hydride transfer reaction from NADH to FMN in morphinone reductase.

The temperature dependence of the primary kinetic isotope effect (KIE), combined temperature-pressure studies of the primary KIE, and studies of the alpha-secondary KIE previously led us to infer that hydride transfer from nicotinamide adenine dinucleotide to flavin mononucleotide in morphinone reductase proceeds via environmentally coupled hydride tunneling. We present here a computational analysis of this hydride transfer reaction using QM/MM molecular dynamics simulations and variational transition-state theory calculations. Our calculated primary and secondary KIEs are in good agreement with the corresponding experimental values. Although the experimentally observed KIE lies below the semiclassical limit, our calculations suggest that approximately 99% of the reaction proceeds via tunneling: this is the first "deep tunneling" reaction observed for hydride transfer. We also show that the dominant tunneling mechanism is controlled by the isotope at the primary rather than the secondary position: with protium in the primary position, large-curvature tunneling dominates, whereas with deuterium in this position, small-curvature tunneling dominates. Also, our study is consistent with tunneling being preceded by reorganization: in the reactant, the rings of the nicotinamide and isoalloxazine moieties are stacked roughly parallel to each other, and as the system moves toward a "tunneling-ready" configuration, the nicotinamide ring rotates to become almost perpendicular to the isoalloxazine ring.

A half-site multimeric enzyme achieves its cooperativity without conformational changes.

Cooperativity is a feature many multimeric proteins use to control activity. Here we show that the bacterial heptose isomerase GmhA displays homotropic positive and negative cooperativity among its four protomers. Most similar proteins achieve this through conformational changes: GmhA instead employs a delicate network of hydrogen bonds, and couples pairs of active sites controlled by a unique water channel. This network apparently raises the Lewis acidity of the catalytic zinc, thus increasing the activity at one active site at the cost of preventing substrate from adopting a reactive conformation at the paired negatively cooperative site - a "half-site" behavior. Our study establishes the principle that multimeric enzymes can exploit this cooperativity without conformational changes to maximize their catalytic power and control. More broadly, this subtlety by which enzymes regulate functions could be used to explore new inhibitor design strategies.

Role of tryptophan residues of Erv1: Trp95 and Trp183 are important for its folding and oxidase function.

Erv1 is an FAD-dependent thiol oxidase of the ERV (essential for respiration and viability)/ALR (augmenter of liver regeneration) sub-family and an essential component of the mitochondrial import and assembly pathway. Erv1 contains six tryptophan residues, which are all located in the highly conserved C-terminal FAD-binding domain. Though important structural roles were predicted for the invariable Trp(95), no experimental study has been reported. In the present study, we investigated the structural and functional roles of individual tryptophan residues of Erv1. Six single tryptophan-to-phenylalanine yeast mutant strains were generated and their effects on cell viability were tested at various temperatures. Then, the mutants were purified from Escherichia coli. Their effects on folding, FAD-binding and Erv1 activity were characterized. Our results showed that Erv1(W95F) has the strongest effect on the stability and function of Erv1 and followed by Erv1(W183F). Erv1(W95F) results in a decrease in the Tm of Erv1 by 23°C, a significant loss of the oxidase activity and thus causing cell growth defects at both 30°C and 37°C. Erv1(W183F) induces changes in the oligomerization state of Erv1, along with a pronounced effect on the stability of Erv1 and its function at 37°C, whereas the other mutants had no clear effect on the function of Erv1 including the highly conserved Trp(157) mutant. Finally, computational analysis indicates that Trp(95) plays a key role in stabilizing the isoalloxazine ring to interact with Cys(133). Taken together, the present study provided important insights into the molecular mechanism of how thiol oxidases use FAD in catalysing disulfide bond formation.

Drug-polymer intermolecular interactions in hot-melt extruded solid dispersions.

The purpose of the study was to investigate and identify the interactions within solid dispersions of cationic drugs and anionic polymers processed by hot-melt extrusion (HME) technique. Propranolol HCl (PRP) and diphenhydramine HCl (DPD) were used as model cationic active substances while pH sensitive anionic methacrylic acid based methyl methacrylate copolymers Eudragit L100 (L100) and ethyl acrylate copolymer Eudragit L100-55 (Acryl EZE) (L100-55) were used as polymeric carriers. The extrudates were further characterised using various physicochemical characterisation techniques to determine the morphology, the drug state within the polymer matrices and the type of drug-polymer interactions. Molecular modelling predicted the existence of two possible H-bonding types while the X-ray photon spectroscopy (XPS) advanced surface analysis of the extrudates revealed intermolecular ionic interactions between the API amino functional groups and the polymer carboxylic groups through the formation of hydrogen bonding. The magnitude of the intermolecular interactions varied according to the drug-polymer miscibility.

New insights into the multi-step reaction pathway of the reductive half-reaction catalysed by aromatic amine dehydrogenase: a QM/MM study.

Computational insight into the multi-step reaction cycle of aromatic amine dehydrogenase is presented, identifying the energy landscape and pathway for multiple proton transfers. This atomistic picture of the reaction sequence--including short-lived reaction intermediates and a stepwise reaction mechanism--bridges the gap between a small number of crystallographic snapshots.

Atomistic insight into the origin of the temperature-dependence of kinetic isotope effects and H-tunnelling in enzyme systems is revealed through combined experimental studies and biomolecular simulation.

The physical basis of the catalytic power of enzymes remains contentious despite sustained and intensive research efforts. Knowledge of enzyme catalysis is predominantly descriptive, gained from traditional protein crystallography and solution studies. Our goal is to understand catalysis by developing a complete and quantitative picture of catalytic processes, incorporating dynamic aspects and the role of quantum tunnelling. Embracing ideas that we have spearheaded from our work on quantum mechanical tunnelling effects linked to protein dynamics for H-transfer reactions, we review our recent progress in mapping macroscopic kinetic descriptors to an atomistic understanding of dynamics linked to biological H-tunnelling reactions.

Molecular dynamics simulation of thermal unfolding of Thermatoga maritima DHFR.

Molecular dynamics simulations of the temperature-induced unfolding reaction of native dimeric dihydrofolate reductase from the hyperthermophile Thermatoga maritima (TmDHFR) and the experimentally inaccessible TmDHFR monomer were carried out at 400 K, 450 K and 500 K. The results revealed that the unfolding of TmDHFR subunits followed a similar path to that of the monomeric DHFR from the mesophile E. coli (EcDHFR). An initial collapse of the adenosine-binding domain (ABD) was followed by the loss of the N-terminal and loop domains (NDLD). Interestingly, the elements of the secondary structure of the isolated TmDHFR monomer were maintained for significantly longer periods of time for the hyperthermophilic enzyme, suggesting that subunit stability contributes to the enhanced resistance of TmDHFR to temperature-induced unfolding. The interactions between the subunits of the TmDHFR dimer led to a stabilisation of the NDLD. The hydrogen bonds between residues 140-143 in betaG of one subunit and residues 125-127 in betaF of the other subunit were retained for significant parts of the simulations at all temperatures. These intermolecular hydrogen bonds were lost after the unfolding of the individual subunits. The high stability of the dimer mediated by strong intersubunit contacts together with an intrinsically enhanced stability of the subunits compared to EcDHFR provides a molecular rational for the higher stability of the thermophilic enzyme. The computed unfolding pathways suggest that the partly folded dimer may be a genuine folding intermediate.