Cation-π versus OH-π Interactions in Proteins: A Density Functional Study

Switching a HO···π Interaction to a Nonconventional OH···π Hydrogen Bond: A Completed Crystallographic Puzzle

In this article, we present crystallographic and spectroscopic evidence of a tunable system wherein a HO···π interaction switches incrementally to a nonconventional OH···π hydrogen bonding (HB) interaction. More specifically, we report the synthesis of substituted forms of model system 1 to study the effects of aryl ring electronic density on the qualitative characteristics of OH···π hydrogen bonds therein. The OH stretch in experimental infrared data, in agreement with density-functional theory (DFT) calculations, shows continuous red-shifts as the adjacent ring becomes more electron rich. For example, the OH stretch of an amino-substituted analogue is red-shifted by roughly 50 cm^-1 compared to the same stretch in the CF₃ analogue, indicating a significantly stronger HB interaction in the former. Moreover, DFT calculations (ωB97XD/6-311+G**) predict that increasing electronic density on the adjacent top ring reduces the aryl π-OH σ* energy gap with a concomitant enhancement of the OH n-π* energy gap. Consequently, a dominant π-σ* interaction in the amino substituted analogue locks the system in the in-form while a favorable n-π* interaction reverses the orientation of the oxygen-bound hydrogen in its protonated form. Additionally, the ¹H NMR data of various analogues reveal that stronger OH···π interactions in systems with electron-rich aromatic rings slow exchange of the alcoholic proton, thereby revealing coupling with the geminal proton. Finally, X-ray crystallographic analyses of a spectrum of analogues clearly visualize the three distinct stages of "switch"-starting with exclusive HO···π, to partitioned HO···π/OH···π, and finally to achieving exclusive OH···π forms. Furthermore, the crystal structure of the amino analogue reveals an interesting feature in which an extended HB network, involving two conventional (NH···O) and two nonconventional (OH···π) HBs, dimerizes and anchors the molecule in the unit cell.

Serine and Cysteine π-Interactions in Nature: A Comparison of the Frequency, Structure, and Stability of Contacts Involving Oxygen and Sulfur

Despite many DNA–protein π-interactions in high-resolution crystal structures, only four X–H···π or X···π interactions were found between serine (Ser) or cysteine (Cys) and DNA nucleobase π-systems in over 100 DNA–protein complexes (where X = O for Ser and X = S for Cys). Nevertheless, 126 non-covalent contacts occur between Ser or Cys and the aromatic amino acids in many binding arrangements within proteins. Furthermore, Ser and Cys protein–protein π-interactions occur with similar frequencies and strengths. Most importantly, due to the great stability that can be provided to biological macromolecules (up to –20 kJ mol–1 for neutral π-systems or –40 kJ mol–1 for cationic π-systems), Ser and Cys π-interactions should be considered when analyzing protein stability and function.

Benzophenone-based derivatives: a novel series of potent and selective dual inhibitors of acetylcholinesterase and acetylcholinesterase-induced beta-amyloid aggregation

The leading mechanistic theory of Alzheimer's disease (AD) is the "amyloid hypothesis" which states that the accumulation of the amyloid β protein (Aβ), and its subsequent aggregation into plaques, is responsible for the initiation of a cascade of events resulting in neurodegeneration and dementia. The anti-amyloid disease-modifying approach, based on the decrease in the production of Aβ, gained thus a paramount importance. The aim of this study was the design and synthesis of a new series of acetylcholinesterase inhibitors (AChEIs) endowed with anti-Aβ aggregating capability. These dual binding inhibitors, being able to interact both with the peripheral anionic site (PAS) of AChE and the catalytic subsite, proved to be able to inhibit the AChE-induced Aβ aggregation. Thus, starting from the lead compound 1, an AChEI composed by a benzophenone scaffold and a N,N'-methylbenzylamino group, a substantial modification aimed at targeting the PAS was performed. To this aim, different amino-terminal side chains were incorporated into this main framework, in order to mimic the diethylmethylammonium alkyl moiety of the pure PAS ligand propidium. The synthesized compounds proved to effectively and selectively inhibit AChE. Moreover, compounds 16a-c and 18a,b, with a propoxy and a hexyloxy tether respectively, showed a good activity against the AChE-induced Aβ aggregation. In particular, molecular modeling studies confirmed that compounds carrying the diethylaminopropoxy and the diethylaminohexyloxy side chains (compounds 16a and 19a, respectively) could suitably contact the PAS pocket of the enzyme.

Differential bonding interactions of inverse agonists of angiotensin II type 1 receptor in stabilizing the inactive state

Although the sartan family of angiotensin II type 1 (AT(1)) receptor blockers (ARBs), which includes valsartan, olmesartan, and losartan, have a common pharmacophore structure, their effectiveness in therapy differs. Although their efficacy may be related to their binding strength, this notion has changed with a better understanding of the molecular mechanism. Therefore, we hypothesized that each ARB differs with regard to its molecular interactions with AT(1) receptor in inducing inverse agonism. Interactions between valsartan and residues Ser(105), Ser(109), and Lys(199) were important for binding. Valsartan is a strong inverse agonist of constitutive inositol phosphate production by the wild-type and N111G mutant receptors. Substituted cysteine accessibility mapping studies indicated that valsartan, but not losartan, which has only weak inverse agonism, may stabilize the N111G receptor in an inactive state upon binding. In addition, the inverse agonism by valsatan was mostly abolished with S105A/S109A/K199Q substitutions in the N111G background. Molecular modeling suggested that Ser(109) and Lys(199) bind to phenyl and tetrazole groups of valsartan, respectively. Ser(105) is a candidate for binding to the carboxyl group of valsartan. Thus, the most critical interaction for inducing inverse agonism involves transmembrane (TM) V (Lys(199)) of AT(1) receptor although its inverse agonist potency is comparable to olmesartan, which bonds with TM III (Tyr(113)) and TM VI (His(256)). These results provide new insights into improving ARBs and development of new G protein-coupled receptor antagonists.

MCH-R1 antagonists based on an arginine scaffold: SAR studies on the amino-terminus

We have identified a novel series of potent MCH-R1 antagonists based on l-arginine. As predicted by computational methods, there was an activity dependence on the pi-electronic character of the aromatic systems corresponding to the amino-terminus of these molecules. These results have enhanced our understanding of the MCH-R1 receptor and the potential for a predictive homology model.

Structural requirements for factor Xa inhibition by 3-oxybenzamides with neutral P1 substituents: combining X-ray crystallography, 3D-QSAR, and tailored scoring functions

The design, synthesis, and structure-activity relationship of 3-oxybenzamides as potent inhibitors of the coagulation protease factor Xa are described on the basis of X-ray structures, privileged structure motifs, and SAR information. A total of six X-ray structures of fXa/inhibitor complexes led us to identify the major protein-ligand interactions. The binding mode is characterized by a lipophilic dichlorophenyl substituent interacting with Tyr228 in the protease S1 pocket, while polar parts are accommodated in S4. This alignment in combination with docking allowed derivation of 3D-QSAR models and tailored scoring functions to rationalize biological affinity and provide guidelines for optimization. The resulting models showed good correlation coefficients and predictions of external test sets. Furthermore, they correspond to binding site topologies in terms of steric, electrostatic, and hydrophobic complementarity. Two approaches to derive tailored scoring functions combining binding site and ligand information led to predictive models with acceptable predictions of the external set. Good correlations to experimental affinities were obtained for both AFMoC (adaptation of fields for molecular comparison) and the novel TScore function. The SAR information from 3D-QSAR and tailored scoring functions agrees with all experimental data and provides guidelines and reasonable activity estimations for novel fXa inhibitors.

3-(4-[[Benzyl(methyl)amino]methyl]phenyl)-6,7-dimethoxy-2H-2-chromenone (AP2238) inhibits both acetylcholinesterase and acetylcholinesterase-induced beta-amyloid aggregation: a dual function lead for Alzheimer's disease therapy

In recent years, the investigation of acetylcholinesterase (AChE) inhibitors has gained further interest, because the involvement of the peripheral site of the enzyme in the beta-amyloid (Abeta) aggregation process has been disclosed. We present here, for the first time, a direct evidence of the Abeta antiaggregating action of an AChE inhibitor (AP2238) purposely designed to bind at both the catalytic and the peripheral sites of the human enzyme.

Structurally conserved amino Acid w501 is required for RNA helicase activity but is not essential for DNA helicase activity of hepatitis C virus NS3 protein

Hepatitis C virus (HCV) is a positive-strand RNA virus that encodes a helicase required for viral replication. Although HCV does not replicate through a DNA intermediate, HCV helicase unwinds both RNA and DNA duplexes. An X-ray crystal structure of the HCV helicase complexed with (dU)(8) has been solved, and the substrate-amino acids interactions within the catalytic pocket were shown. Among these, residues W501 and V432 were reported to have base stacking interactions and to be important for the unwinding function of HCV helicase. It has been hypothesized that specific interactions between the enzyme and substrate in the catalytic pocket are responsible for the substrate specificity phenotype. We therefore mutagenized W501 and V432 to investigate their role in substrate specificity in HCV helicase. Replacement of W501, but not V432, with nonaromatic residues resulted in complete loss of RNA unwinding activity, whereas DNA unwinding activity was largely unaffected. The loss of unwinding activity was fully restored in the W501F mutant, indicating that the aromatic ring is crucial for RNA helicase function. Analysis of ATPase and nucleic acid binding activities in the W501 mutant enzymes revealed that these activities are not directly responsible for the substrate specificity phenotype. Molecular modeling of the enzyme-substrate interaction at W501 revealed a putative pi-facial hydrogen bond between the 2'-OH of ribose and the aromatic tryptophan ring. This evidence correlates with biochemical results suggesting that the pi-facial bond may play an important role in the RNA unwinding activity of the HCV NS3 protein.

Design and quantitative structure-activity relationship of 3-amidinobenzyl-1H-indole-2-carboxamides as potent, nonchiral, and selective inhibitors of blood coagulation factor Xa

A series of 138 nonchiral 3-amidinobenzyl-1H-indole-2-carboxamides and analogues as inhibitors of the blood coagulation enzyme factor Xa (fXa) were designed, synthesized, and investigated by X-ray structure analysis and 3D quantitative structure-activity relationship (QSAR) studies (CoMFA, CoMSIA) in order to identify important protein-ligand interactions responsible for biological affinity and selectivity. Several compounds from this series are highly potent and selective inhibitors of this important enzyme linking extrinsic and intrinsic coagulation pathways. To rationalize biological affinity and to provide guidelines for further design, all compounds were docked into the factor Xa binding site. Those docking studies were based on X-ray structures of factor Xa in complex with literature-known inhibitors. It was possible to validate those binding modes by four X-ray crystal structures of representative ligands in factor Xa, while one ligand was additionally crystallized in trypsin to rationalize requirements for selective factor Xa inhibition. The 3D-QSAR models based on a superposition rule derived from these docking studies were validated using conventional and cross-validated r(2) values using the leave-one-out method and repeated analyses using two randomly chosen cross-validation groups plus randomization of biological activities. This led to consistent and highly predictive 3D-QSAR models with good correlation coefficients for both CoMFA and CoMSIA, which were found to correspond to experimentally determined factor Xa binding site topology in terms of steric, electrostatic, and hydrophobic complementarity. Subsets selected as smaller training sets using 2D fingerprints and maximum dissimilarity methods resulted in 3D-QSAR models with remarkable correlation coefficients and a high predictive power. The final quantitative SAR information agrees with all experimental data for the binding topology and thus provides reasonable activity predictions for novel factor Xa inhibitors.

Enzymatic GTP hydrolysis: insights from an ab initio molecular dynamics study

Ab initio methods were used to shed light on fundamental aspects of the enzymatic mechanism of guanosine triphosphate hydrolysis in the Cdc42/Cdc42GAP complex. The calculations focused on the nucleophilic addition of the catalytic water molecule to the gamma-phosphate phosphorus atom. A large model system was required to correctly reproduce the electrostatic properties on the active site. The model turned out to reproduce most of the electrostatic field of the biological complex at the reactants. Our calculations established the H-bond pattern of the catalytic water (WAT), which turned out to interact with Q61 and T35, in the most stable conformation. This ruled out the possibility that the catalytic water transferred its proton directly to the gamma-phosphate. Furthermore, the calculations suggested that the electronic structure of WAT was very different from that in the bulk. Finally, this study showed that during the reaction, WAT transferred a proton to Gln61, consistent with the available X-ray data on a transition-state analogue/enzyme complex(19) and with the decrease of activity in the Q61E mutant.

Investigation of synthetic hosts that model cation-pi sites found at protein binding domains

Small cyclophanes containing aromatic groups and dialkyl ammonium ions were created as model systems of the cation-pi complexes found at some protein binding domains. The hosts had different shapes in order to investigate the effect the arrangement of ammonium ions to aromatic surfaces has on their reactivity. pK(a) values of the hosts were substantially different in DMSO or (95/5) DMSO/D(2)O solutions, which showed that the ions existed in different environments of the hosts. Electrostatic charges, as determined by density functional calculations, revealed that the magnitude of a cationic charge depends on its position relative to an aromatic ring. Association constants of the hosts bound to the sodium salt of N-acetyl phenylalanine in d(6)-DMSO and in (95/5) d(6)-DMSO/D(2)O solutions were inversely proportional to the magnitude of the hosts' acidity constants. These results suggest that the magnitude of the positive charge for cationic groups of cation-pi complexes is reduced by being associated with electron-rich faces of aromatic rings. The aromatic rings, however, lessen the desolvation penalty that must be overcome for ligand binding, giving an overall more favorable association.

Tetraethylammonium binding to the outer mouth of the KcsA potassium channel: implications for ion permeation

Extracellular tetraethylammonium (TEA+) inhibits the current carried out by K+ ions in potassium channels. Structural models of wild-type (WT) and Y82C KcsA K+ channel/TEA+ complexes are here built using docking procedures, electrostatics calculations and molecular dynamics simulations. The calculations are based on the structure determined by Doyle et al. (11) Our calculations suggest that in WT, the TEA+ cation turns binds at the outer mouth of the selectivity filter, stabilized by electrostatic and hydrophobic interactions with the four Tyr82 side chains. Replacement of Tyr82 with Cys causes a decrease of the affinity of the cation for the channel, consistently with the available site-directed mutagenesis data (16). An MD simulation in which K+ replaces TEA+ provides evidence that TEA+ binding site can accommodate a potassium ion, in agreement with the high-resolution structure recently reported by Zhou et al. (20)

A classical and ab initio study of the interaction of the myosin triphosphate binding domain with ATP

We used classical molecular mechanics (MM) simulations and quantum mechanical (QM) structural relaxations to examine the active site of myosin when bound to ATP. Two conformations of myosin have been determined by x-ray crystallography. In one, there is no direct interaction between switch 2 and the nucleotide (open state). In the other (closed state), the universally conserved switch 2 glycine forms a hydrogen bond with a gamma-phosphate oxygen. MM simulations indicate that the two states are thermodynamically stable and allow us to investigate the extent to which the P-loop, switch 1, and switch 2 are involved in hydrolysis. We find that the open structure has a higher affinity for ATP than the closed structure, and that ATP is distorted toward a transition state by interactions with the protein. We also examine how the structure of the binding site changes with either MgATP or CaATP as the nucleotide in myosin in the open conformer. Our analyses suggest that higher CaATPase rates occur because the leaving phosphate (P(i)) group is more weakly bound and dissociation occurs faster. Finally, we validate the use of a particular formulation of a QM methodology (Car-Parrinello) to further refine the structures of the active site.