Axel Bidon-Chanal

MDpocket: open-source cavity detection and characterization on molecular dynamics trajectories

MOTIVATION A variety of pocket detection algorithms are now freely or commercially available to the scientific community for the analysis of static protein structures. However, since proteins are dynamic entities, enhancing the capabilities of these programs for the straightforward detection and characterization of cavities taking into account protein conformational ensembles should be valuable for capturing the plasticity of pockets, and therefore allow gaining insight into structure-function relationships. RESULTS This article describes a new method, called MDpocket, providing a fast, free and open-source tool for tracking small molecule binding sites and gas migration pathways on molecular dynamics (MDs) trajectories or other conformational ensembles. MDpocket is based on the fpocket cavity detection algorithm and a valuable contribution to existing analysis tools. The capabilities of MDpocket are illustrated for three relevant cases: (i) the detection of transient subpockets using an ensemble of crystal structures of HSP90; (ii) the detection of known xenon binding sites and migration pathways in myoglobin; and (iii) the identification of suitable pockets for molecular docking in P38 Map kinase. AVAILABILITY MDpocket is free and open-source software and can be downloaded at http://fpocket.sourceforge.net. CONTACT pschmidtke@ub.edu SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Ligand‐induced dynamical regulation of NO conversion in Mycobacterium tuberculosis truncated hemoglobin‐N

Mycobacterium tuberculosis, the causative agent of human tuberculosis, is forced into latency by nitric oxide produced by macrophages during infection. In response to nitrosative stress M. tuberculosis has evolved a defense mechanism that relies on the oxygenated form of “truncated hemoglobin” N (trHbN), formally acting as NO‐dioxygenase, yielding the harmless nitrate ion. X‐ray crystal structures have shown that trHbN hosts a two‐branched protein matrix tunnel system, proposed to control diatomic ligand migration to the heme, as the rate‐limiting step in NO conversion to nitrate. Extended molecular dynamics simulations (0.1 μs), employed here to characterize the factors controlling diatomic ligand diffusion through the apolar tunnel system, suggest that O2 migration in deoxy‐trHbN is restricted to a short branch of the tunnel, and that O2 binding to the heme drives conformational and dynamical fluctuations promoting NO migration through the long tunnel branch. The simulation results suggest that trHbN has evolved a dual‐path mechanism for migration of O2 and NO to the heme, to achieve the most efficient NO detoxification. Proteins 2006.

Essential Dynamics: A Tool for Efficient Trajectory Compression and Management

We present a simple method for compression and management of very large molecular dynamics trajectories. The approach is based on the projection of the Cartesian snapshots collected along the trajectory into an orthogonal space defined by the eigenvectors obtained by diagonalization of the covariance matrix. The transformation is mathematically exact when the number of eigenvectors equals 3N-6 (N being the number of atoms), and in practice very accurate even when the number of eigenvectors is much smaller, permitting a dramatic reduction in the size of trajectory files. In addition, we have examined the ability of the method, when combined with interpolation, to recover dense samplings (snapshots collected at a high frequency) from more sparse (lower frequency) data as a method for further data compression. Finally, we have investigated the possibility of using the approach when extrapolating the behavior of the system to times longer than the original simulation period. Overall our results suggest that the method is an attractive alternative to current approaches for including dynamic information in static structure files such as those deposited in the Protein Data Bank.

Protein flexibility and ligand recognition: challenges for molecular modeling.

The intrinsic dynamics of macromolecules is an essential property to relate the structure of biomolecular systems with their function in the cell. In the field of ligand-receptor recognition, numerous evidences have revealed the limitations of the lock-and-key theory, and the need to elaborate models that take into account the inherent plasticity of biomolecules, such as the induced-fit model or the existence of an ensemble of pre-equilibrated conformations. Depending on the nature of the target system, ligand binding can be associated with small local adjustments in side chains or even the backbone to large-scale motions of structural fragments, domains or even subunits. Reproducing the inherent flexibility of biomolecules has thus become one of the most challenging issues in molecular modeling and simulation studies, as it has direct implications in our understanding of the structure-function relationships, but even in areas such as virtual screening and structure-based drug discovery. Given the intrinsic limitation of conventional simulation tools, only events occurring in short time scales can be reproduced at a high accuracy level through all-atom techniques such as Molecular Dynamics simulations. However, larger structural rearrangements demand the use of enhanced sampling methods relying on modified descriptions of the biomolecular system or the potential surface. This review illustrates the crucial role that structural plasticity plays in mediating ligand recognition through representative examples. In addition, it discusses some of the most powerful computational tools developed to characterize the conformational flexibility in ligand-receptor complexes.

Role of Pre-A Motif in Nitric Oxide Scavenging by Truncated Hemoglobin, HbN, of Mycobacterium tuberculosis

Mycobacterium tuberculosis truncated hemoglobin, HbN, is endowed with a potent nitric-oxide dioxygenase activity and has been found to relieve nitrosative stress and enhance in vivo survival of a heterologous host, Salmonella enterica Typhimurium, within the macrophages. These findings implicate involvement of HbN in the defense of M. tuberculosis against nitrosative stress. The protein carries a tunnel system composed of a short and a long tunnel branch that has been proposed to facilitate diatomic ligand migration to the heme and an unusual Pre-A motif at the N terminus, which does not contribute significantly to the structural integrity of the protein, as it protrudes out of the compact globin fold. Strikingly, deletion of Pre-A region from the M. tuberculosis HbN drastically reduces its ability to scavenge nitric oxide (NO), whereas its insertion at the N terminus of Pre-A lacking HbN of Mycobacterium smegmatis improved its nitric-oxide dioxygenase activity. Titration of the oxygenated adduct of HbN and its mutants with NO indicated that the stoichiometric oxidation of protein is severalfold slower when the Pre-A region is deleted in HbN. Molecular dynamics simulations show that the excision of Pre-A motif results in distinct changes in the protein dynamics, which cause the gate of the tunnel long branch to be trapped into a closed conformation, thus impeding migration of diatomic ligands toward the heme active site. The present study, thus, unequivocally demonstrates vital function of Pre-A region in NO scavenging and unravels its unique role by which HbN might attain its efficient NO-detoxification ability.

Huprine–Tacrine Heterodimers as Anti-Amyloidogenic Compounds of Potential Interest against Alzheimer’s and Prion Diseases

A family of huprine-tacrine heterodimers has been developed to simultaneously block the active and peripheral sites of acetylcholinesterase (AChE). Their dual site binding for AChE, supported by kinetic and molecular modeling studies, results in a highly potent inhibition of the catalytic activity of human AChE and, more importantly, in the in vitro neutralization of the pathological chaperoning effect of AChE toward the aggregation of both the β-amyloid peptide (Aβ) and a prion peptide with a key role in the aggregation of the prion protein. Huprine-tacrine heterodimers take on added value in that they display a potent in vitro inhibitory activity toward human butyrylcholinesterase, self-induced Aβ aggregation, and β-secretase. Finally, they are able to cross the blood-brain barrier, as predicted in an artificial membrane model assay and demonstrated in ex vivo experiments with OF1 mice, reaching their multiple biological targets in the central nervous system. Overall, these compounds are promising lead compounds for the treatment of Alzheimers and prion diseases.

Continuum solvation models: Dissecting the free energy of solvation

The most usual self-consistent reaction field (SCRF) continuum models for the description of solvation within the quantum mechanical (QM) framework are reviewed, trying to emphasize their common roots as well as the inherent approximations assumed in the calculation of the free energy of solvation. Particular attention is also paid to the specific features involved in the development of current state-of-the-art QM SCRF continuum models. This is used to discuss the need to maintain a close correspondence between each SCRF formalism and the specific details entailing its parametrization, as well as the need to be cautious in analyzing the balance between electrostatic and non-electrostatic contributions to the solvation free energy between different SCRF models. Finally, special emphasis is given to the post-processing of the free energy of solvation to derive parameters providing a compact picture of the ability of a molecule to interact with different solvents, which can be of particular interest in biopharmaceutical studies.

Tacrine-based dual binding site acetylcholinesterase inhibitors as potential disease-modifying anti-Alzheimer drug candidates.

Two novel families of dual binding site acetylcholinesterase (AChE) inhibitors have been developed, consisting of a tacrine or 6-chlorotacrine unit as the active site interacting moiety, either the 5,6-dimethoxy-2-[(4-piperidinyl)methyl]-1-indanone fragment of donepezil (or the indane derivative thereof) or a 5-phenylpyrano[3,2-c]quinoline system, reminiscent to the tryciclic core of propidium, as the peripheral site interacting unit, and a linker of suitable length as to allow the simultaneous binding at both sites. These hybrid compounds are all potent and selective inhibitors of human AChE, and more interestingly, are able to interfere in vitro both formation and aggregation of the beta-amyloid peptide, the latter effects endowing these compounds with the potential to modify Alzheimers disease progression.

Novel huprine derivatives with inhibitory activity toward β-amyloid aggregation and formation as disease-modifying anti-Alzheimer drug candidates.

A new family of dual binding site acetylcholinesterase (AChE) inhibitors has been designed, synthesized, and tested for their ability to inhibit AChE, butyrylcholinesterase (BChE), AChE‐induced and self‐induced β‐amyloid (Aβ) aggregation and β‐secretase (BACE‐1), and to cross the blood–brain barrier. The new heterodimers consist of a unit of racemic or enantiopure huprine Y or X and a donepezil‐related 5,6‐dimethoxy‐2‐[(4‐piperidinyl)methyl]indane moiety as the active site and peripheral site to mid‐gorge‐interacting moieties, respectively, connected through a short oligomethylene linker. Molecular dynamics simulations and kinetics studies support the dual site binding to AChE. The new heterodimers are potent inhibitors of human AChE and moderately potent inhibitors of human BChE, AChE‐induced and self‐induced Aβ aggregation, and BACE‐1, and are predicted to be able to enter the central nervous system (CNS), thus constituting promising multitarget anti‐Alzheimer drug candidates with the potential to modify the natural course of this disease.

High pressure reveals structural determinants for globin hexacoordination: Neuroglobin and myoglobin cases

The influence of pressure on the equilibrium between five‐(5c) and six‐coordination (6c) forms in neuroglobin (Ngb) and myoglobin (Mb) has been examined by means of molecular dynamics (MD) simulations at normal and high pressure. The results show that the main effect of high pressure is to reduce the protein mobility without altering the structure in a significant manner. Moreover, our data suggest that the equilibrium between 5c and 6c states in globins is largely controlled by the structure and dynamics of the C‐D region. Finally, in agreement with the available experimental data, the free energy profiles obtained from steered MD for both proteins indicate that high pressure enhances hexacoordination. In Ngb, the shift in equilibrium is mainly related to an increase in the 6c→5c transition barrier, whereas in Mb such a shift is primarily due to a destabilization of the 5c state. Proteins 2009.