Sílvia G. Estácio

The Kohn-Sham density of states and band gap of water: from small clusters to liquid water.

Electronic properties of water clusters (H2O)(n), with n=2, 4, 8, 10, 15, 20, and 30 molecules were investigated by sequential Monte Carlo/density-functional theory (DFT) calculations. DFT calculations were carried out over uncorrelated configurations generated by Monte Carlo simulations of liquid water with a reparametrized exchange-correlation functional that reproduces the experimental information on the electronic properties (first ionization energy and highest occupied molecular orbital-lowest unoccupied molecular orbital gap) of the water dimer. The dependence of electronic properties on the cluster size (n) shows that the density of states (DOS) of small water clusters (n>10) exhibits the same basic features that are typical of larger aggregates, such as the mixing of the 3a1 and 1b1 valence bands. When long-ranged polarization effects are taken into account by the introduction of embedding charges, the DOS associated with 3a1 orbitals is significantly enhanced. In agreement with valence-band photoelectron spectra of liquid water, the 1b1, 3a1, and 1b2 electron binding energies in water aggregates are redshifted by approximately 1 eV relative to the isolated molecule. By extrapolating the results for larger clusters the threshold energy for photoelectron emission is 9.6+/-0.15 eV (free clusters) and 10.58+/-0.10 eV (embedded clusters). Our results for the electron affinity (V0=-0.17+/-0.05 eV) and adiabatic band gap (E(G,Ad)=6.83+/-0.05 eV) of liquid water are in excellent agreement with recent information from theoretical and experimental works.

Identification of a Conserved Aggregation-Prone Intermediate State in the Folding Pathways of Spc-SH3 Amyloidogenic Variants

We compared the folding pathways of selected mutational variants of the α-spectrin SH3 domain (Spc-SH3) by using a continuum model that combines a full atomistic protein representation with the Gō potential. Experimental data show that the N47G mutant shows very little tendency to aggregate while the N47A and triple mutant D48G(2Y) are both amyloidogenic, with the latter being clearly more aggregation prone. We identified a strikingly similar native-like folding intermediate across the three mutants, in which strand β(1) is totally unstructured and more than half of the major hydrophobic core residues are highly solvent exposed. Results from extensive docking simulations show that the ability of the intermediates to dimerize is largely driven by strand β(1) and is consistent with the in vitro aggregation behavior reported for the corresponding mutants. They further suggest that residues 44 and 53, which are key players in the nucleation-condensation mechanism of folding, are also important triggers of the aggregation process.

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.

Robustness of atomistic Gō models in predicting native-like folding intermediates.

Gō models are exceedingly popular tools in computer simulations of protein folding. These models are native-centric, i.e., they are directly constructed from the proteins native structure. Therefore, it is important to understand up to which extent the atomistic details of the native structure dictate the folding behavior exhibited by Gō models. Here we address this challenge by performing exhaustive discrete molecular dynamics simulations of a Gō potential combined with a full atomistic protein representation. In particular, we investigate the robustness of this particular type of Gō models in predicting the existence of intermediate states in protein folding. We focus on the N47G mutational form of the Spc-SH3 folding domain (x-ray structure) and compare its folding pathway with that of alternative native structures produced in silico. Our methodological strategy comprises equilibrium folding simulations, structural clustering, and principal component analysis.

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 Gō 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.

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.

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.

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.

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.

Energetics of tert-Butoxyl Addition Reaction to Norbornadiene: A Method for Estimating the π-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.