Brian McClain - Academia.edu (original) (raw)

Papers by Brian McClain

J. Chem. Theory Comput., 2014

The solubility of organic molecules is of critical importance to the pharmaceutical industry; how... more The solubility of organic molecules is of critical importance to the pharmaceutical industry; however, robust computational methods to predict this quantity from firstprinciples are lacking. Solubility can be computed from a thermodynamic cycle that decomposes standard state solubility into the sum of solid−vapor sublimation and vapor−liquid solvation free energies ΔG solubility°= ΔG sub°+ ΔG solv°. Over the past few decades, alchemical simulation methods to compute solvation free energy using classical force fields have become widely used. However, analogous methods for determining the free energy of the sublimation/deposition phase transition are currently limited by the necessity of a priori knowledge of the atomic coordinates of the crystal. Here, we describe progress toward an alternative scheme based on growth of the asymmetric unit into a crystal via alchemy (GAUCHE). GAUCHE computes deposition free energy ΔG dep°= −ΔG sub°= −k B T ln(V c /V g ) + ΔG AU + ΔG AU→UC as the sum of an entropic term to account for compressing a vapor at 1 M standard state (V g ) into the molar volume of the crystal (V c ), where k B is Boltzmann's constant and T is temperature in degrees Kelvin, plus two simulation steps. In the first simulation step, the deposition free energy ΔG AU for a system composed of only N AU asymmetric unit (AU) molecule(s) is computed beginning from an arbitrary conformation in vacuum. In the second simulation step, the change in free energy ΔG AU→UC to expand the asymmetric unit degrees of freedom into a unit cell (UC) composed of N UC independent molecules is computed. This latter step accounts for the favorable free energy of removing the constraint that every symmetry mate of the asymmetric unit has an identical conformation and intermolecular interactions. The current work is based on NVT simulations, which requires knowledge of the crystal space group and unit cell parameters from experiment, but not a priori knowledge of crystalline atomic coordinates. GAUCHE was applied to 5 organic molecules whose sublimation free energy has been measured experimentally, based on the polarizable AMOEBA force field and more than a microsecond of sampling per compound in the program Force Field X. The mean unsigned and RMS errors were only 1.6 and 1.7 kcal/mol, respectively, which indicates that GAUCHE is capable of accurate prediction of absolute sublimation thermodynamics.

Journal of Molecular Biology, 2013

Double-stranded RNA (dsRNA) viruses transcribe and replicate RNA within an assembled, inner capsi... more Double-stranded RNA (dsRNA) viruses transcribe and replicate RNA within an assembled, inner capsid particle; only plus-sense mRNA emerges into the intracellular milieu. During infectious entry of a rotavirus particle, the outer layer of its three-layer structure dissociates, delivering the inner double-layered particle (DLP) into the cytosol. DLP structures determined by X-ray crystallography and electron cryomicroscopy (cryoEM) show that the RNA coils uniformly into the particle interior, avoiding a "fivefold hub" of more structured density projecting inward from the VP2 shell of the DLP along each of the twelve 5-fold axes. Analysis of the X-ray crystallographic electron density map suggested that principal contributors to the hub are the N-terminal arms of VP2, but reexamination of the cryoEM map has shown that many features come from a molecule of VP1, randomly occupying five equivalent and partly overlapping positions. We confirm here that the electron density in the X-ray map leads to the same conclusion, and we describe the functional implications of the orientation and position of the polymerase. The exit channel for the nascent transcript directs the nascent transcript toward an opening along the 5-fold axis. The template strand enters from within the particle, and the dsRNA product of the initial replication step exits in a direction tangential to the inner surface of the VP2 shell, allowing it to coil optimally within the DLP. The polymerases of reoviruses appear to have similar positions and functional orientations.

Journal of Medicinal Chemistry, 2008

Small molecule protein kinase inhibitors are widely employed as biological reagents and as leads ... more Small molecule protein kinase inhibitors are widely employed as biological reagents and as leads in the design of drugs for a variety of diseases. We investigated the phenomenon of kinase-likeness, i.e., the propensity of ligands to inhibit protein kinases, in the context of kinase-specific substructural fragments. The frequency of occurrence of multiple structural fragments in kinase inhibitor libraries relative to nonkinase compounds has been analyzed. A combination of structural fragment counts, termed the "2-0" kinaselikeness rule, provides approximately 5-fold enrichment in kinase active compounds. This rule has been validated using in-house kinase counterscreening data and applied prospectively to uncover kinase activities in marketed drugs. In addition, the role of discriminating fragments in kinase recognition was interrogated using available structural data, providing an insight into their effect on inhibitor potency and selectivity. One of these fragments, bisarylaniline, has been characterized as a kinase-privileged fragment with specific binding preferences and a link to increased activity within kinases.

Biophysical Chemistry, 2015

As compounds are optimized for greater potency during pharmaceutical discovery, their aqueous sol... more As compounds are optimized for greater potency during pharmaceutical discovery, their aqueous solubility often decreases, making them less viable as orally-administered drugs. To investigate whether potency and insolubility share a common origin, we examined the structural and thermodynamic properties of telaprevir, a sparingly soluble inhibitor of hepatitis C virus protease. Comparison of the hydrogen bond motifs in crystalline telaprevir with those present in the protease-telaprevir complex revealed striking similarities. Additionally, the thermodynamics of telaprevir dissolution closely resembles those of protein-ligand dissociation. Together, these findings point to a common origin of potency and insolubility rooted in particular amide-amide hydrogen bond patterns. The insolubility of telaprevir is shown by computational analysis to be caused by interactions in the crystal, not unfavorable hydrophobic hydration. Accordingly, competing out the particular amide-amide hydrogen bond motifs in crystalline telaprevir with 4-hydroxybenzoic acid yielded a co-crystalline solid with excellent aqueous dissolution and oral absorption. The analysis suggests a generalizable approach for identifying drug candidate compounds that either can or cannot be rendered orally bioavailable by alteration of their crystalline solid phases, in an approach that provides a pragmatic way to attain substantial enhancements in the success rate of drug discovery and development.

J. Chem. Theory Comput., 2014

The solubility of organic molecules is of critical importance to the pharmaceutical industry; how... more The solubility of organic molecules is of critical importance to the pharmaceutical industry; however, robust computational methods to predict this quantity from firstprinciples are lacking. Solubility can be computed from a thermodynamic cycle that decomposes standard state solubility into the sum of solid−vapor sublimation and vapor−liquid solvation free energies ΔG solubility°= ΔG sub°+ ΔG solv°. Over the past few decades, alchemical simulation methods to compute solvation free energy using classical force fields have become widely used. However, analogous methods for determining the free energy of the sublimation/deposition phase transition are currently limited by the necessity of a priori knowledge of the atomic coordinates of the crystal. Here, we describe progress toward an alternative scheme based on growth of the asymmetric unit into a crystal via alchemy (GAUCHE). GAUCHE computes deposition free energy ΔG dep°= −ΔG sub°= −k B T ln(V c /V g ) + ΔG AU + ΔG AU→UC as the sum of an entropic term to account for compressing a vapor at 1 M standard state (V g ) into the molar volume of the crystal (V c ), where k B is Boltzmann's constant and T is temperature in degrees Kelvin, plus two simulation steps. In the first simulation step, the deposition free energy ΔG AU for a system composed of only N AU asymmetric unit (AU) molecule(s) is computed beginning from an arbitrary conformation in vacuum. In the second simulation step, the change in free energy ΔG AU→UC to expand the asymmetric unit degrees of freedom into a unit cell (UC) composed of N UC independent molecules is computed. This latter step accounts for the favorable free energy of removing the constraint that every symmetry mate of the asymmetric unit has an identical conformation and intermolecular interactions. The current work is based on NVT simulations, which requires knowledge of the crystal space group and unit cell parameters from experiment, but not a priori knowledge of crystalline atomic coordinates. GAUCHE was applied to 5 organic molecules whose sublimation free energy has been measured experimentally, based on the polarizable AMOEBA force field and more than a microsecond of sampling per compound in the program Force Field X. The mean unsigned and RMS errors were only 1.6 and 1.7 kcal/mol, respectively, which indicates that GAUCHE is capable of accurate prediction of absolute sublimation thermodynamics.

Journal of Molecular Biology, 2013

Double-stranded RNA (dsRNA) viruses transcribe and replicate RNA within an assembled, inner capsi... more Double-stranded RNA (dsRNA) viruses transcribe and replicate RNA within an assembled, inner capsid particle; only plus-sense mRNA emerges into the intracellular milieu. During infectious entry of a rotavirus particle, the outer layer of its three-layer structure dissociates, delivering the inner double-layered particle (DLP) into the cytosol. DLP structures determined by X-ray crystallography and electron cryomicroscopy (cryoEM) show that the RNA coils uniformly into the particle interior, avoiding a "fivefold hub" of more structured density projecting inward from the VP2 shell of the DLP along each of the twelve 5-fold axes. Analysis of the X-ray crystallographic electron density map suggested that principal contributors to the hub are the N-terminal arms of VP2, but reexamination of the cryoEM map has shown that many features come from a molecule of VP1, randomly occupying five equivalent and partly overlapping positions. We confirm here that the electron density in the X-ray map leads to the same conclusion, and we describe the functional implications of the orientation and position of the polymerase. The exit channel for the nascent transcript directs the nascent transcript toward an opening along the 5-fold axis. The template strand enters from within the particle, and the dsRNA product of the initial replication step exits in a direction tangential to the inner surface of the VP2 shell, allowing it to coil optimally within the DLP. The polymerases of reoviruses appear to have similar positions and functional orientations.

Journal of Medicinal Chemistry, 2008

Small molecule protein kinase inhibitors are widely employed as biological reagents and as leads ... more Small molecule protein kinase inhibitors are widely employed as biological reagents and as leads in the design of drugs for a variety of diseases. We investigated the phenomenon of kinase-likeness, i.e., the propensity of ligands to inhibit protein kinases, in the context of kinase-specific substructural fragments. The frequency of occurrence of multiple structural fragments in kinase inhibitor libraries relative to nonkinase compounds has been analyzed. A combination of structural fragment counts, termed the "2-0" kinaselikeness rule, provides approximately 5-fold enrichment in kinase active compounds. This rule has been validated using in-house kinase counterscreening data and applied prospectively to uncover kinase activities in marketed drugs. In addition, the role of discriminating fragments in kinase recognition was interrogated using available structural data, providing an insight into their effect on inhibitor potency and selectivity. One of these fragments, bisarylaniline, has been characterized as a kinase-privileged fragment with specific binding preferences and a link to increased activity within kinases.

Biophysical Chemistry, 2015

As compounds are optimized for greater potency during pharmaceutical discovery, their aqueous sol... more As compounds are optimized for greater potency during pharmaceutical discovery, their aqueous solubility often decreases, making them less viable as orally-administered drugs. To investigate whether potency and insolubility share a common origin, we examined the structural and thermodynamic properties of telaprevir, a sparingly soluble inhibitor of hepatitis C virus protease. Comparison of the hydrogen bond motifs in crystalline telaprevir with those present in the protease-telaprevir complex revealed striking similarities. Additionally, the thermodynamics of telaprevir dissolution closely resembles those of protein-ligand dissociation. Together, these findings point to a common origin of potency and insolubility rooted in particular amide-amide hydrogen bond patterns. The insolubility of telaprevir is shown by computational analysis to be caused by interactions in the crystal, not unfavorable hydrophobic hydration. Accordingly, competing out the particular amide-amide hydrogen bond motifs in crystalline telaprevir with 4-hydroxybenzoic acid yielded a co-crystalline solid with excellent aqueous dissolution and oral absorption. The analysis suggests a generalizable approach for identifying drug candidate compounds that either can or cannot be rendered orally bioavailable by alteration of their crystalline solid phases, in an approach that provides a pragmatic way to attain substantial enhancements in the success rate of drug discovery and development.