Inclusion of Enclosed Hydration Effects in the Binding Free Energy Estimation of Dopamine D3 Receptor Complexes

Confined hydration and conformational flexibility are some of the challenges encountered for the rational design of selective antagonists of G-protein coupled receptors. We present a set of C3-substituted (-)-stepholidine derivatives as potent binders of the dopamine D3 receptor. The compounds are characterized biochemically, as well as by computer modeling using a novel molecular dynamics-based alchemical binding free energy approach which incorporates the effect of the displacement of enclosed water molecules from the binding site. The free energy of displacement of specific hydration sites is obtained using the Hydration Site Analysis method with explicit solvation. This work underscores the critical role of confined hydration and conformational reorganization in the molecular recognition mechanism of dopamine receptors and illustrates the potential of binding free energy models to represent these key phenomena.Comment: This is the first report of using enclosed hydration in estimating binding free energies of protein-ligand complexes using implicit solvatio

Hydrogen-activation mechanism of [Fe] hydrogenase revealed by multi-scale modeling

When investigating the mode of hydrogen activation by [Fe] hydrogenases, not only the chemical reactivity at the active site is of importance but also the large-scale conformational change between the so-called open and closed conformations, which leads to a special spatial arrangement of substrate and iron cofactor. To study H2 activation, a complete model of the solvated and cofactor-bound enzyme in complex with the substrate methenyl-H4MPT+ was constructed. Both the closed and open conformations were simulated with classical molecular dynamics on the 100 ns time scale. Quantum-mechanics/molecular-mechanics calculations on snapshots then revealed the features of the active site that enable the facile H2 cleavage. The hydroxyl group of the pyridinol ligand can easily be deprotonated. With the deprotonated hydroxyl group and the structural arrangement in the closed conformation, H2 coordinated to the Fe center is subject to an ionic and orbital push-pull effect and can be rapidly cleaved with a concerted hydride transfer to methenyl-H4MPT+. An intermediary hydride species is not formed

Molecular determinants of drug-specific sensitivity for epidermal growth factor receptor (EGFR) exon 19 and 20 mutants in non-small cell lung cancer.

We hypothesized that aberrations activating epidermal growth factor receptor (EGFR) via dimerization would be more sensitive to anti-dimerization agents (e.g., cetuximab). EGFR exon 19 abnormalities (L747_A750del; deletes amino acids LREA) respond to reversible EGFR kinase inhibitors (TKIs). Exon 20 in-frame insertions and/or duplications (codons 767 to 774) and T790M mutations are clinically resistant to reversible/some irreversible TKIs. Their impact on protein function/therapeutic actionability are not fully elucidated.In our study, the index patient with non-small cell lung cancer (NSCLC) harbored EGFR D770_P772del_insKG (exon 20). A twenty patient trial (NSCLC cohort) (cetuximab-based regimen) included two participants with EGFR TKI-resistant mutations ((i) exon 20 D770>GY; and (ii) exon 19 LREA plus exon 20 T790M mutations). Structural modeling predicted that EGFR exon 20 anomalies (D770_P772del_insKG and D770>GY), but not T790M mutations, stabilize the active dimer configuration by increasing the interaction between the kinase domains, hence sensitizing to an agent preventing dimerization. Consistent with predictions, the two patients harboring D770_P772del_insKG and D770>GY, respectively, responded to an EGFR antibody (cetuximab)-based regimen; the T790M-bearing patient showed no response to cetuximab combined with erlotinib. In silico modeling merits investigation of its ability to optimize therapeutic selection based on structural/functional implications of different aberrations within the same gene

Magnetism, FeS colloids, and Origins of Life

A number of features of living systems: reversible interactions and weak bonds underlying motor-dynamics; gel-sol transitions; cellular connected fractal organization; asymmetry in interactions and organization; quantum coherent phenomena; to name some, can have a natural accounting via $physical$ interactions, which we therefore seek to incorporate by expanding the horizons of `chemistry-only' approaches to the origins of life. It is suggested that the magnetic 'face' of the minerals from the inorganic world, recognized to have played a pivotal role in initiating Life, may throw light on some of these issues. A magnetic environment in the form of rocks in the Hadean Ocean could have enabled the accretion and therefore an ordered confinement of super-paramagnetic colloids within a structured phase. A moderate H-field can help magnetic nano-particles to not only overcome thermal fluctuations but also harness them. Such controlled dynamics brings in the possibility of accessing quantum effects, which together with frustrations in magnetic ordering and hysteresis (a natural mechanism for a primitive memory) could throw light on the birth of biological information which, as Abel argues, requires a combination of order and complexity. This scenario gains strength from observations of scale-free framboidal forms of the greigite mineral, with a magnetic basis of assembly. And greigite's metabolic potential plays a key role in the mound scenario of Russell and coworkers-an expansion of which is suggested for including magnetism.Comment: 42 pages, 5 figures, to be published in A.R. Memorial volume, Ed Krishnaswami Alladi, Springer 201

Optimal assignment methods for ligand-based virtual screening

<p>Abstract</p> <p>Background</p> <p>Ligand-based virtual screening experiments are an important task in the early drug discovery stage. An ambitious aim in each experiment is to disclose active structures based on new scaffolds. To perform these "scaffold-hoppings" for individual problems and targets, a plethora of different similarity methods based on diverse techniques were published in the last years. The optimal assignment approach on molecular graphs, a successful method in the field of quantitative structure-activity relationships, has not been tested as a ligand-based virtual screening method so far.</p> <p>Results</p> <p>We evaluated two already published and two new optimal assignment methods on various data sets. To emphasize the "scaffold-hopping" ability, we used the information of chemotype clustering analyses in our evaluation metrics. Comparisons with literature results show an improved early recognition performance and comparable results over the complete data set. A new method based on two different assignment steps shows an increased "scaffold-hopping" behavior together with a good early recognition performance.</p> <p>Conclusion</p> <p>The presented methods show a good combination of chemotype discovery and enrichment of active structures. Additionally, the optimal assignment on molecular graphs has the advantage to investigate and interpret the mappings, allowing precise modifications of internal parameters of the similarity measure for specific targets. All methods have low computation times which make them applicable to screen large data sets.</p

Measuring cellular traction forces on non-planar substrates

Animal cells use traction forces to sense the mechanics and geometry of their environment. Measuring these traction forces requires a workflow combining cell experiments, image processing and force reconstruction based on elasticity theory. Such procedures have been established before mainly for planar substrates, in which case one can use the Green's function formalism. Here we introduce a worksflow to measure traction forces of cardiac myofibroblasts on non-planar elastic substrates. Soft elastic substrates with a wave-like topology were micromolded from polydimethylsiloxane (PDMS) and fluorescent marker beads were distributed homogeneously in the substrate. Using feature vector based tracking of these marker beads, we first constructed a hexahedral mesh for the substrate. We then solved the direct elastic boundary volume problem on this mesh using the finite element method (FEM). Using data simulations, we show that the traction forces can be reconstructed from the substrate deformations by solving the corresponding inverse problem with a L1-norm for the residue and a L2-norm for 0th order Tikhonov regularization. Applying this procedure to the experimental data, we find that cardiac myofibroblast cells tend to align both their shapes and their forces with the long axis of the deformable wavy substrate.Comment: 34 pages, 9 figure