Structure-based virtual screening toward the discovery of novel inhibitors of the DNA repair activity of the human apurinic/apyrimidinic endonuclease 1.

The DNA repair activity of human apurinic/apyrimidinic endonuclease 1 (APE1) has been recognized as a promising target for the development of small-molecule inhibitors to be used in combination with anticancer agents. In an attempt to identify novel inhibitors of APE1, we present a structure-based virtual screening (SBVS) study based on molecular docking analysis of the compounds of NCI database using the GOLD 5.1.0 (Genetic Optimization for Ligand Docking) suite of programs. Compounds selected in this screening were tested with a fluorescence-based APE1 endonuclease activity assay. Two compounds (37 and 41) were able to inhibit the multifunctional enzyme APE1 in the micromolar range, while compound 22 showed inhibitory effects at nanomolar concentrations. These results were confirmed by a plasmid DNA nicking assay. In addition, the potential APE1 inhibitors did not affect the cell viability of non-tumor MCF10A cells. Overall, compounds 22, 37, and 41 appear to be important scaffolds for the design of novel APE1 inhibitors and this study highlights the relevance of in silico-based approaches as valuable tools in drug discovery.

Thermal unfolding simulations of NBD1 domain variants reveal structural motifs associated with the impaired folding of F508del-CFTR.

We employed high-temperature classical molecular dynamics (MD) simulations to investigate the unfolding process of the wild-type (WT) and F508del-NBD1 domains of CFTR protein, with and without second-site mutations. To rationalize the in vitro behavior of F508del-NBD1, namely its lower folding yield and higher aggregation propensity, we focused our analysis of the MD data on the existence of intermediate states with aggregation potential and/or stabilized by a significant number of non-native interactions (i.e. misfolded states). We find that the deletion of phenylalanine 508 is able to forcefully reshape the conformational space of the NBD1 domain to the extent that it uniquely populates intermediate states whose structural traits provide important insights into the molecular events that underlie the impaired folding of F508del-NBD1. In particular, our simulations predict the formation of a misfolded intermediate whose population is highly enhanced by deletion of residue 508. The stabilization of this intermediate is a direct consequence of the enhanced non-native coupling between various key regions of the α-helical subdomain and ATP-binding subdomain; it is singularly characterized by a set of non-native interactions within the ATP-binding subdomain and between that domain and the α-helical subdomain region. The formation of this intermediate is not blocked by second-site suppressor mutations which indicates a limited role of the latter in correcting the rather complex folding process of the CFTR protein missing residue 508.

Calcium binding to gatekeeper residues flanking aggregation-prone segments underlies non-fibrillar amyloid traits in superoxide dismutase 1 (SOD1).

Calcium deregulation is a central feature among neurodegenerative diseases, including amyotrophic lateral sclerosis (ALS). Calcium accumulates in the spinal and brain stem motor neurons of ALS patients triggering multiple pathophysiological processes which have been recently shown to include direct effects on the aggregation cascade of superoxide dismutase 1 (SOD1). SOD1 is a Cu/Zn enzyme whose demetallated form is implicated in ALS protein deposits, contributing to toxic gain of function phenotypes. Here we undertake a combined experimental and computational study aimed at establishing the molecular details underlying the regulatory effects of Ca(2+) over SOD1 aggregation potential. Isothermal titration calorimetry indicates entropy driven low affinity association of Ca(2+) ions to apo SOD1, at pH7.5 and 37°C. Molecular dynamics simulations denote a noticeable loss of native structure upon Ca(2+) association that is especially prominent at the zinc-binding and electrostatic loops, whose decoupling is known to expose the central SOD1 ß-barrel triggering aggregation. Structural mapping of the preferential apo SOD1 Ca(2+) binding locations reveals that among the most frequent ligands for Ca(2+) are negatively-charged gatekeeper residues located in boundary positions with respect to segments highly prone to edge-to-edge aggregation. Calcium interactions thus diminish gatekeeping roles of these residues, by shielding repulsive interactions via stacking between aggregating ß-sheets, partly blocking fibril formation and promoting amyloidogenic oligomers such as those found in ALS inclusions. Interestingly, many fALS mutations occur at these positions, disclosing how Ca(2+) interactions recreate effects similar to those of genetic defects, a finding with relevance to understand sporadic ALS pathomechanisms.

A simulated intermediate state for folding and aggregation provides insights into ΔN6 ß2-microglobulin amyloidogenic behavior.

A major component of ex vivo amyloid plaques of patients with dialysis-related amyloidosis (DRA) is a cleaved variant of ß2-microglobulin (ΔN6) lacking the first six N-terminal residues. Here we perform a computational study on ΔN6, which provides clues to understand the amyloidogenicity of the full-length ß2-microglobulin. Contrary to the wild-type form, ΔN6 is able to efficiently nucleate fibrillogenesis in vitro at physiological pH. This behavior is enhanced by a mild acidification of the medium such as that occurring in the synovial fluid of DRA patients. Results reported in this work, based on molecular simulations, indicate that deletion of the N-terminal hexapeptide triggers the formation of an intermediate state for folding and aggregation with an unstructured strand A and a native-like core. Strand A plays a pivotal role in aggregation by acting as a sticky hook in dimer assembly. This study further predicts that the detachment of strand A from the core is maximized at pH 6.2 resulting into higher aggregation efficiency. The structural mapping of the dimerization interface suggests that Tyr10, His13, Phe30 and His84 are hot-spot residues in ΔN6 amyloidogenesis.

Assessing the effect of loop mutations in the folding space of ß2-microglobulin with molecular dynamics simulations.

We use molecular dynamics simulations of a full atomistic Go model to explore the impact of selected DE-loop mutations (D59P and W60C) on the folding space of protein human ß2-microglobulin (Hß2m), the causing agent of dialysis-related amyloidosis, a conformational disorder characterized by the deposition of insoluble amyloid fibrils in the osteoarticular system. Our simulations replicate the effect of mutations on the thermal stability that is observed in experiments in vitro. Furthermore, they predict the population of a partially folded state, with 60% of native internal free energy, which is akin to a molten globule. In the intermediate state, the solvent accessible surface area increases up to 40 times relative to the native state in 38% of the hydrophobic core residues, indicating that the identified species has aggregation potential. The intermediate state preserves the disulfide bond established between residue Cys25 and residue Cys80, which helps maintain the integrity of the core region, and is characterized by having two unstructured termini. The movements of the termini dominate the essential modes of the intermediate state, and exhibit the largest displacements in the D59P mutant, which is the most aggregation prone variant. PROPKA predictions of pKa suggest that the population of the intermediate state may be enhanced at acidic pH explaining the larger amyloidogenic potential observed in vitro at low pH for the WT protein and mutant forms.

In silico strategies toward enzyme function and dynamics.

Enzymes are outstanding biocatalysts involved in a plethora of chemical reactions occurring in the cell. Despite their incommensurable importance, a comprehensive understanding of enzyme catalysis is still missing. This task becomes more laborious given the unavoidability of including the inherent dynamic nature of enzymes into that description. As such, it is essential to ascertain the nature and contribution of enzyme conformational changes to catalysis and to evaluate the adequacy of the proposal associating protein internal motions to the rate enhancement achieved. Dynamic events in enzymes span a wide range of time- and length-scales which have led to a surge in multiscale methodologies targeting enzyme function and dynamics. Computational strategies assume a preponderant role in such studies by allowing the atomic detail investigation of the fundamental mechanisms of enzyme catalysis thus surpassing what is achievable through experiments. While high-accuracy quantum mechanical methods are indicated to uncover the details of the chemical reaction occurring at the active site, molecular mechanical force fields and molecular dynamics approaches provide powerful means to access the conformational energy landscape accessible to enzymes. This review outlines some of the most important in silico methodologies in this area, highlighting examples of problems tackled and the insights obtained.

Characterizing the dynamics and ligand-specific interactions in the human leukocyte elastase through molecular dynamics simulations.

The human leukocyte elastase (HLE), a neutrophil serine protease of the chymotrypsin superfamily, is a major therapeutic target for a number of inflammatory diseases, such as chronic obstructive pulmonary disease (COPD). In this work, we present a comparative explicit water molecular dynamics (MD) study on the free and inhibitor-bound HLE. Knowledge of the flexibility and conformational changes induced by this irreversible inhibitor, whether in a prebound stage or covalently bound at the enzyme binding site, encases fundamental biological interest and is particularly relevant to ongoing structure-based drug design studies. Our results suggest that HLE operates by an induced-fit mechanism with direct intervention of a surface loop which is open toward the solvent in the free enzyme and closed while in the presence of the ligand. MM-PBSA free energy calculations furthermore elucidate the energetic contributions to the distinct conformations adopted by this loop. Additionally, a survey of the major contributions to the inhibitor binding free energies was attained. Our findings enforce the need to account for HLE flexibility, whether through the use of MD-generated ensembles of HLE conformations as targets for molecular docking or via sophisticated flexible-docking algorithms. We anticipate that inclusion of the observed HLE dynamic behavior into future drug design methodologies will have a relevant impact in the development of novel, more efficient, inhibitors.

Energetics of tert-butoxyl addition reaction to norbornadiene: a method for estimating the pi-bond strength of a carbon-carbon double bond.

The energetics of tert-butoxyl radical addition reaction to norbornadiene was investigated by time-resolved photoacoustic calorimetry (TR-PAC). The result, together with the C-O bond dissociation enthalpy (BDE) in the addition product, allowed us to calculate the pi-bond dissociation enthalpy in norbornadiene. Quantum chemistry (QC) methods were also used to obtain several enthalpies of reaction of the addition of oxygen-centered radicals to alkenes. The pi-bond dissociation enthalpies in these molecules were calculated by a procedure similar to that used in the case of norbornadiene and were compared with the pi-BDE values obtained by the method proposed by Benson. These two different approaches yield similar values for the pi-BDEs in alkenes, indicating that the addition method proposed in the present study is a valid way to derive that quantity. The influence of strain in the pi-BDEs of cyclic alkenes was investigated and allowed us to justify the difference between the pi-BDE in norbornene and norbornadiene. Finally, the thermochemistry of the addition and abstraction reactions involving these two molecules and tert-butoxyl radical was analyzed.

C-H bond dissociation enthalpies in norbornane. An experimental and computational study.

Gas-phase C-H bond dissociation enthalpies (BDEs) in norbornane were determined by quantum chemistry calculations and the C2-H BDE was experimentally obtained for the first time by time-resolved photoacoustic calorimetry. CBS-Q and CBS-QB3 methods were used to derive the values DH degrees (C1-H) = 449 kJ mol-1, DH degrees (C7-H) = 439 kJ mol-1, and DH degrees (C2-H) = 413 kJ mol-1. The experimental result DH degrees (C2-H) = 414.6 +/- 5.4 kJ mol-1 is in excellent agreement with the theoretical value. The trend DH degrees (C1-H) > DH degrees (C7-H) > DH degrees (C2-H) is discussed.