Conformational preference of dipeptide zwitterions in aqueous solvents.

Proper description of solvent effects is challenging for theoretical methods, particularly if the solute is a zwitterion. Here, a series of theoretical procedures are used to determine the preferred solvated conformations of twelve hydrophobic dipeptides (Leu-Leu, Leu-Phe, Phe-Leu, Ile-Leu, Phe-Phe, Ala-Val, Val-Ala, Ala-Ile, Ile-Ala, Ile-Val, Val-Ile and Val-Val) in the zwitterionic state. First, the accuracy of density functional theory (DFT), combined with different implicit solvent models, for describing zwitterions in aqueous solvent is assessed by comparing the predicted against the experimental glycine tautomerization energy, i.e., the energetic difference between canonical and zwitterionic glycine in aqueous solvents. It is found that among the tested solvation schemes, the charge-asymmetric nonlocally determined local-electric solvation model (CANDLE) predicts an energetic difference in excellent agreement with the experimental value. Next, DFT-CANDLE is used to determine the most favorable solvated conformation for each of the investigated dipeptide zwitterions. The CANDLE-solvated structures are obtained by exploring the conformational space of each dipeptide zwitterion concatenating DFT calculations, in vacuum, with classical molecular dynamics simulations, in explicit solvents, and DFT calculations including explicit water molecules. It is found that the energetically most favorable conformations are similar to those of the dipeptide zwitterions in their respective crystal structures. Such structural agreement is indicative of the DFT-CANDLE accomplishment of the description of solvated zwitterions, and suggests that these biomolecules self-assemble as quasi-rigid objects.

Effect of strand register in the stability and reactivity of crystals from peptides forming amyloid fibrils.

Amyloids are cytotoxic protein aggregates that deposit in human tissues, leading to several health disorders. Their aggregates can also exhibit catalytic properties, and they have been used as candidates for the development of functional biomaterials. Despite being polymorphic, amyloids often assemble as cross-ß fibrils formed by in-register ß sheet layers. Recent studies of some amyloidogenic protein segments revealed that they crystallize as antiparallel out-of-register ß sheets. Such arrangement has been proposed to be responsible for the cytotoxicity in amyloid diseases, however, there is still no consensus on the molecular mechanism. Interestingly, two amyloidogenic peptide segments, NFGAILS and FGAILSS, arrange into out-of-register and in-register ß sheets, respectively, even though they solely differ by one aminoacid residue at both termini. In this work, we used density functional theory (DFT) to address how the strand register contributes into the packing and molecular properties of the NFGAILS and FGAILSS crystals. Our results show that the out-of-register structure is substantially more stable, at 0 K, than the in-register one due to stronger inter-strand contacts. Based on an analysis of the electrostatic potential of the crystal slabs, it is suggested that the out-of-register may potentially interact with negatively charged groups, like those found in cell membranes. Moreover, calculated reactivity descriptors indicate a similar outcome, where only the out-of-register peptide exhibits intrinsic reactive surface sites at the exposed amine and carboxylic groups. It is therefore suggested that the out-of-register arrangement may indeed be crucial for amyloid cytotoxicity. The findings presented here could help to further our understanding of amyloid aggregation, function, and toxicity.

PSIQUE: Protein Secondary Structure Identification on the Basis of Quaternions and Electronic Structure Calculations.

The secondary structure is important in protein structure analysis, classification, and modeling. We have developed a novel method for secondary structure assignment, termed PSIQUE, based on the potential energy surface (PES) of polyalanine obtained using an infinitely long chain model and density functional theory calculations. First, uniform protein segments are determined in terms of a difference of quaternions between neighboring amino acids along the protein backbone. Then, the identification of the secondary structure motifs is carried out based on the minima found in the PES. PSIQUE shows good agreement with other secondary structure assignment methods. However, it provides better discrimination of subtle secondary structures (e.g., helix types) and termini and produces more uniform segments while also accounting for local distortions. Overall, PSIQUE provides a precise and reliable assignment of secondary structures, so it should be helpful for the detailed characterization of the protein structure.

Tetrahydroquinoline-Isoxazole/Isoxazoline Hybrid Compounds as Potential Cholinesterases Inhibitors: Synthesis, Enzyme Inhibition Assays, and Molecular Modeling Studies.

A series of 44 hybrid compounds that included in their structure tetrahydroquinoline (THQ) and isoxazole/isoxazoline moieties were synthesized through the 1,3-dipolar cycloaddition reaction (1,3-DC) from the corresponding N-allyl/propargyl THQs, previously obtained via cationic Povarov reaction. In vitro cholinergic enzymes inhibition potential of all compounds was tested. Enzyme inhibition assays showed that some hybrids exhibited significant potency to inhibit acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). Especially, the hybrid compound 5n presented the more effective inhibition against AChE (4.24 µM) with an acceptable selectivity index versus BChE (SI: 5.19), while compound 6aa exhibited the greatest inhibition activity on BChE (3.97 µM) and a significant selectivity index against AChE (SI: 0.04). Kinetic studies were carried out for compounds with greater inhibitory activity of cholinesterases. Structure-activity relationships of the molecular hybrids were analyzed, through computational models using a molecular cross-docking algorithm and Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) binding free energy approach, which indicated a good correlation between the experimental inhibition values and the predicted free binding energy.

Docking and quantitative structure-activity relationship of bi-cyclic heteroaromatic pyridazinone and pyrazolone derivatives as phosphodiesterase 3A (PDE3A) inhibitors.

PDE3s belong to the phosphodiesterases family, where the PDE3A isoform is the major subtype in platelets involved in the cAMP regulation pathway of platelet aggregation. PDE3A inhibitors have been designed as potential antiplatelet agents. In this work, a homology model of PDE3A was developed and used to obtain the binding modes of bicyclic heteroaromatic pyridazinones and pyrazolones. Most of the studied compounds adopted similar orientations within the PDE3A active site, establishing hydrogen bonds with catalytic amino acids. Besides, the structure-activity relationship of the studied inhibitors was described by using a field-based 3D-QSAR method. Different structure alignment strategies were employed, including template-based and docking-based alignments. Adequate correlation models were obtained according to internal and external validations. In general, QSAR models revealed that steric and hydrophobic fields describe the different inhibitory activities of the compounds, where the hydrogen bond donor and acceptor fields have minor contributions. It should be stressed that structural elements of PDE3A inhibitors are reported here, through descriptions of their binding interactions and their differential affinities. In this sense, the present results could be useful in the future design of more specific and potent PDE3A inhibitors that may be used for the treatment of cardiovascular diseases.

Modeling cooperative effects in halogen-bonded infinite linear chains.

Non-additivity in noncovalent interactions is an important aspect of complex systems that can lead to stronger (cooperative) interactions when three or more molecular units influence each other. The halogen bond (XB) is a highly-directional noncovalent interaction that has been found to be cooperative. Here the strength and nature of cooperativity arising in X-bonded infinite linear chains of cyanogen halides and 4-halopyridines are investigated by means of density functional theory calculations. It is found that cyanogen halide chains are highly cooperative (up to 77%), whereas pyridines are only slightly cooperative (below 21%). It is demonstrated that XB and its non-additivity can be modeled as the sum of a local term, which depends on first nearest-neighbors only, and long-range effective dipole-dipole attractions. It is shown that the local term in cyanogen halides primarily accounts for repulsive short-range screened Coulomb interactions, whereas in 4-halopyridines such a term also includes attractive contributions, which are particularly sizeable in some elongated XB conformations. This outcome reveals differences in the nature of the XBs formed in these molecular systems. Nevertheless, it is shown that both systems behave as effective point dipoles regarding cooperative effects, at any point of the XB dissociation path. As such, these results are useful contributions for the understanding and modeling of non-additive effects of noncovalent interactions.

Binding-affinity predictions of HSP90 in the D3R Grand Challenge 2015 with docking, MM/GBSA, QM/MM, and free-energy simulations.

We have estimated the binding affinity of three sets of ligands of the heat-shock protein 90 in the D3R grand challenge blind test competition. We have employed four different methods, based on five different crystal structures: first, we docked the ligands to the proteins with induced-fit docking with the Glide software and calculated binding affinities with three energy functions. Second, the docked structures were minimised in a continuum solvent and binding affinities were calculated with the MM/GBSA method (molecular mechanics combined with generalised Born and solvent-accessible surface area solvation). Third, the docked structures were re-optimised by combined quantum mechanics and molecular mechanics (QM/MM) calculations. Then, interaction energies were calculated with quantum mechanical calculations employing 970-1160 atoms in a continuum solvent, combined with energy corrections for dispersion, zero-point energy and entropy, ligand distortion, ligand solvation, and an increase of the basis set to quadruple-zeta quality. Fourth, relative binding affinities were estimated by free-energy simulations, using the multi-state Bennett acceptance-ratio approach. Unfortunately, the results were varying and rather poor, with only one calculation giving a correlation to the experimental affinities larger than 0.7, and with no consistent difference in the quality of the predictions from the various methods. For one set of ligands, the results could be strongly improved (after experimental data were revealed) if it was recognised that one of the ligands displaced one or two water molecules. For the other two sets, the problem is probably that the ligands bind in different modes than in the crystal structures employed or that the conformation of the ligand-binding site or the whole protein changes.

Performance of the MM/GBSA scoring using a binding site hydrogen bond network-based frame selection: the protein kinase case.

A conformational selection method, based on hydrogen bond (Hbond) network analysis, has been designed in order to rationalize the configurations sampled using molecular dynamics (MD), which are commonly used in the estimation of the relative binding free energy of ligands to macromolecules through the MM/GBSA or MM/PBSA method. This approach makes use of protein-ligand complexes obtained from X-ray crystallographic data, as well as from molecular docking calculations. The combination of several computational approaches, like long MD simulations on protein-ligand complexes, Hbond network-based selection by scripting techniques and finally MM/GBSA, provides better statistical correlations against experimental binding data than previous similar reported studies. This approach has been successfully applied in the ranking of several protein kinase inhibitors (CDK2, Aurora A and p38), which present both diverse and related chemical structures.

Neonicotinic analogues: selective antagonists for α4ß2 nicotinic acetylcholine receptors.

Nicotine is an agonist of nicotinic acetylcholine receptors (nAChRs) that has been extensively used as a template for the synthesis of α4ß2-preferring nAChRs. Here, we used the N-methyl-pyrrolidine moiety of nicotine to design and synthesise novel α4ß2-preferring neonicotinic ligands. We increased the distance between the basic nitrogen and aromatic group of nicotine by introducing an ester functionality that also mimics acetylcholine (Fig. 2). Additionally, we introduced a benzyloxy group linked to the benzoyl moiety. Although the neonicotinic compounds fully inhibited binding of both [α-(125)I]bungarotoxin to human α7 nAChRs and [(3)H]cytisine to human α4ß2 nAChRs, they were markedly more potent at displacing radioligand binding to human α4ß2 nAChRs than to α7 nAChRs. Functional assays showed that the neonicotinic compounds behave as antagonists at α4ß2 and α4ß2α5 nAChRs. Substitutions on the aromatic ring of the compounds produced compounds that displayed marked selectivity for α4ß2 or α4ß2α5 nAChRs. Docking of the compounds on homology models of the agonist binding site at the α4/ß2 subunit interfaces of α4ß2 nAChRs suggested the compounds inhibit function of this nAChR type by binding the agonist binding site.