Intrinsic α‐helical and β‐sheet conformational preferences: A computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle (-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (N–C_α–C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the “bridge” region of-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between -ψ conformations and τ.     

The power of hard-sphere models: explaining side-chain dihedral angle distributions of Thr and Val

The energy functions used to predict protein structures typically include both molecular-mechanics and knowledge-based terms. In contrast, our approach is to develop robust physics- and geometry-based methods. Here, we investigate to what extent simple hard-sphere models can be used to predict side-chain conformations. The distributions of the side-chain dihedral angle χ₁ of Val and Thr in proteins of known structure show distinctive features: Val side chains predominantly adopt χ₁ = 180°, whereas Thr side chains typically adopt χ₁ = 60° and 300° (i.e., χ₁ = ±60° or g− and g⁺ configurations). Several hypotheses have been proposed to explain these differences, including interresidue steric clashes and hydrogen-bonding interactions. In contrast, we show that the observed side-chain dihedral angle distributions for both Val and Thr can be explained using only local steric interactions in a dipeptide mimetic. Our results emphasize the power of simple physical approaches and their importance for future advances in protein engineering and design.     

Improving consumption rate estimates by incorporating wild activity into a bioenergetics model

Consumption is the basis of metabolic and trophic ecology and is used to assess an animal''s trophic impact. The contribution of activity to an animal''s energy budget is an important parameter when estimating consumption, yet activity is usually measured in captive animals. Developments in telemetry have allowed the energetic costs of activity to be measured for wild animals; however, wild activity is seldom incorporated into estimates of consumption rates. We calculated the consumption rate of a free‐ranging marine predator (yellowtail kingfish, Seriola lalandi) by integrating the energetic cost of free‐ranging activity into a bioenergetics model. Accelerometry transmitters were used in conjunction with laboratory respirometry trials to estimate kingfish active metabolic rate in the wild. These field‐derived consumption rate estimates were compared with those estimated by two traditional bioenergetics methods. The first method derived routine swimming speed from fish morphology as an index of activity (a “morphometric” method), and the second considered activity as a fixed proportion of standard metabolic rate (a “physiological” method). The mean consumption rate for free‐ranging kingfish measured by accelerometry was 152 J·g⁻¹·day⁻¹, which lay between the estimates from the morphometric method (μ = 134 J·g⁻¹·day⁻¹) and the physiological method (μ = 181 J·g⁻¹·day⁻¹). Incorporating field‐derived activity values resulted in the smallest variance in log‐normally distributed consumption rates (σ = 0.31), compared with the morphometric (σ = 0.57) and physiological (σ = 0.78) methods. Incorporating field‐derived activity into bioenergetics models probably provided more realistic estimates of consumption rate compared with the traditional methods, which may further our understanding of trophic interactions that underpin ecosystem‐based fisheries management. The general methods used to estimate active metabolic rates of free‐ranging fish could be extended to examine ecological energetics and trophic interactions across aquatic and terrestrial ecosystems.     

Crystal structures of a therapeutic single chain antibody in complex with two drugs of abuse—Methamphetamine and 3,4‐methylenedioxymethamphetamine

Methamphetamine (METH) is a major drug threat in the United States and worldwide. Monoclonal antibody (mAb) therapy for treating METH abuse is showing exciting promise and the understanding of how mAb structure relates to function will be essential for future development of these important therapies. We have determined crystal structures of a high affinity anti-(+)-METH therapeutic single chain antibody fragment (scFv6H4, K_D= 10 nM) derived from one of our candidate mAb in complex with METH and the (+) stereoisomer of another abused drug, 3,4-methylenedioxymethamphetamine (MDMA), known by the street name “ecstasy.” The crystal structures revealed that scFv6H4 binds to METH and MDMA in a deep pocket that almost completely encases the drugs mostly through aromatic interactions. In addition, the cationic nitrogen of METH and MDMA forms a salt bridge with the carboxylate group of a glutamic acid residue and a hydrogen bond with a histidine side chain. Interestingly, there are two water molecules in the binding pocket and one of them is positioned for a C—H⋯O interaction with the aromatic ring of METH. These first crystal structures of a high affinity therapeutic antibody fragment against METH and MDMA (resolution = 1.9 Å, and 2.4 Å, respectively) provide a structural basis for designing the next generation of higher affinity antibodies and also for carrying out rational humanization.     

Crystallization and preliminary analysis of crystals of the 24-meric hemocyanin of the emperor scorpion (Pandinus imperator)

Hemocyanins are giant oxygen transport proteins found in the hemolymph of several invertebrate phyla. They constitute giant multimeric molecules whose size range up to that of cell organelles such as ribosomes or even small viruses. Oxygen is reversibly bound by hemocyanins at binuclear copper centers. Subunit interactions within the multisubunit hemocyanin complex lead to diverse allosteric effects such as the highest cooperativity for oxygen binding found in nature. Crystal structures of a native hemocyanin oligomer larger than a hexameric substructure have not been published until now. We report for the first time growth and preliminary analysis of crystals of the 24-meric hemocyanin (M_W = 1.8 MDa) of emperor scorpion (Pandinus imperator), which diffract to a resolution of 6.5 Å. The crystals are monoclinc with space group C 1 2 1 and cell dimensions a = 311.61 Å, b = 246.58 Å and c = 251.10 Å (α = 90.00°, β = 90.02°, γ = 90.00°). The asymmetric unit contains one molecule of the 24-meric hemocyanin and the solvent content of the crystals is 56%. A preliminary analysis of the hemocyanin structure reveals that emperor scorpion hemocyanin crystallizes in the same oxygenated conformation, which is also present in solution as previously shown by cryo-EM reconstruction and small angle x-ray scattering experiments.     

Crystal structure of a DNA containing the planar, phenoxazine-derived bi-functional spectroscopic probe C

Previously, we developed the deoxycytosine analog Ç (C-spin) as a bi-functional spectroscopic probe for the study of nucleic acid structure and dynamics using electron paramagnetic resonance (EPR) and fluorescence spectroscopy. To understand the effect of Ç on nucleic acid structure, we undertook a detailed crystallographic analysis. A 1.7 Å resolution crystal structure of Ç within a decamer duplex A-form DNA confirmed that Ç forms a non-perturbing base pair with deoxyguanosine, as designed. In the context of double-stranded DNA Ç adopted a planar conformation. In contrast, a crystal structure of the free spin-labeled base ç displayed a ∼20° bend at the oxazine linkage. Density function theory calculations revealed that the bent and planar conformations are close in energy and exhibit the same frequency for bending. These results indicate a small degree of flexibility around the oxazine linkage, which may be a consequence of the antiaromaticity of a 16-π electron ring system. Within DNA, the amplitude of the bending motion is restricted, presumably due to base-stacking interactions. This structural analysis shows that the Ç forms a planar, structurally non-perturbing base pair with G indicating it can be used with high confidence in EPR- or fluorescence-based structural and dynamics studies.     

Revisiting the Ramachandran plot from a new angle

The pioneering work of Ramachandran and colleagues emphasized the dominance of steric constraints in specifying the structure of polypeptides. The ubiquitous Ramachandran plot of backbone dihedral angles (φ and ψ) defined the allowed regions of conformational space. These predictions were subsequently confirmed in proteins of known structure. Ramachandran and colleagues also investigated the influence of the backbone angle τ on the distribution of allowed φ/ψ combinations. The “bridge region” (φ ≤ 0° and −20° ≤ ψ ≤ 40°) was predicted to be particularly sensitive to the value of τ. Here we present an analysis of the distribution of φ/ψ angles in 850 non-homologous proteins whose structures are known to a resolution of 1.7 Å or less and sidechain B-factor less than 30 Ų. We show that the distribution of φ/ψ angles for all 87,000 residues in these proteins shows the same dependence on τ as predicted by Ramachandran and colleagues. Our results are important because they make clear that steric constraints alone are sufficient to explain the backbone dihedral angle distributions observed in proteins. Contrary to recent suggestions, no additional energetic contributions, such as hydrogen bonding, need be invoked.     

基于稳定同位素技术的塔里木河下游不同林龄胡杨的水分利用来源

水分是制约很多陆地生态系统植物生长和繁殖的重要因素, 在干旱地区尤为明显。利用稳定同位素技术探究塔里木河下游不同林龄胡杨(Populus euphratica)的水分来源情况, 了解生态输水背景下荒漠河岸林的水分利用循环与利用策略, 可为生态输水提供科学依据, 同时也可对同类地区的生态恢复提供借鉴。本研究通过测定塔里木河下游胡杨茎干水和各潜在水源(土壤水、地下水)的稳定氢氧同位素值(δD、δ¹⁸O), 应用多源线性混合模型(IsoSource)分析了各潜在水源对不同林龄胡杨的贡献比例, 并结合3种林龄胡杨不同土壤深度含水量的变化, 分析了胡杨的主要吸水层位。结果表明: (1)不同林龄胡杨样地的不同深度区间上的土壤水δ¹⁸O值存在显著差异(P < 0.05): 胡杨幼龄木、成熟木、过熟木木质部δ¹⁸O分别为-7.83 ± 0.07‰、-8.53 ± 0.11‰、-9.36 ± 0.21‰; 而δD值不存在显著差异(P > 0.05)。可据此来推断胡杨的主要吸水层位。(2)总体上, 三种林龄胡杨土壤水δ¹⁸O值随土壤深度增加而减小, 并趋于接近地下水的δ¹⁸O值。其中, 0-60 cm土壤水受蒸发影响比较大, 其同位素组成经历了强烈的蒸发分馏过程, 土壤含水量极少, 土壤水δ¹⁸O值偏正。(3)不同林龄胡杨所利用的水分来源不同: 胡杨幼龄木对于地表80 cm以下的土壤水以及地下水均有一定程度的利用, 对80-140 cm、140-220 cm和220-340 cm的土壤水平均利用比率依次为16.2%、21.4%和24.6%, 对地下水平均利用比率为24.5%; 成熟木主要利用220-340 cm的土壤水及地下水, 平均利用比率分别为36.9%和42.3%; 过熟木主要利用140-340 cm的土壤水及地下水, 平均利用比率分别为32.8%和49.3%。     

Binding of Sulfamethazine to β-cyclodextrin and Methyl-β-cyclodextrin

β-cyclodextrin (βCD) and methyl-β-cyclodextrin (MβCD) complexes with sulfamethazine (SMT) were prepared and characterized by different experimental techniques, and the effects of βCD and MβCD on drug solubility were assessed via phase-solubility analysis. The phase-solubility diagram for the drug showed an increase in water solubility, with the following affinity constants calculated: 40.4 ± 0.4 (pH 2.0) and 29.4 ± 0.4 (pH 8.0) M⁻¹ with βCD and 56 ± 1 (water), 39 ± 3 (pH 2.0) and 39 ± 5 (pH 8.0) M⁻¹ with MβCD. According to ¹H NMR and 2D NMR spectroscopy, the complexation mode involved the aromatic ring of SMT included in the MβCD cavity. The complexes obtained in solid state by freeze drying were characterized by Fourier transform infrared spectroscopy, scanning electron microscopy, and thermal analysis. The amorphous complexes obtained in this study may be useful in the preparation of pharmaceutical dosage forms of SMT.     

The three-dimensional structure of alpha1-purothionin in solution: combined use of nuclear magnetic resonance, distance geometry and restrained molecular dynamics

The determination of the three-dimensional solution structure of α1-purothionin using a combination of metric matrix distance geometry and restrained molecular dynamics calculations based on n.m.r. data is presented. The experimental data comprise complete sequence-specific proton resonance assignments, a set of 310 approximate interproton distance restraints derived from nuclear Overhauser effects, 27 Ø backbone torsion angle restraints derived from vicinal coupling constants, 4 distance restraints from hydrogen bonds and 12 distance restraints from disulphide bridges. The average atomic rms difference between the final nine converged structures and the mean structure obtained by averaging their coordinates is 1.5 ± 0.1 å for the backbone atoms and 2.0 ± 0.1 å for all atoms. The overall shape of α1-purothionin is that of the capital letter L, similar to that of crambin, with the longer arm comprising two approximately parallel α-helices and the shorter arm a strand and a mini anti-parallel β sheet.     

Structural and Functional Characterization of an Anesthetic Binding Site in the Second Cysteine-Rich Domain of Protein Kinase Cδ?

Elucidating the principles governing anesthetic-protein interactions requires structural determinations at high resolutions not yet achieved with ion channels. Protein kinase C (PKC) activity is modulated by general anesthetics. We solved the structure of the phorbol-binding domain (C1B) of PKCδ complexed with an ether (methoxymethylcycloprane) and with an alcohol (cyclopropylmethanol) at 1.36-Å resolution. The cyclopropane rings of both agents displace a single water molecule in a surface pocket adjacent to the phorbol-binding site, making van der Waals contacts with the backbone and/or side chains of residues Asn-237 to Ser-240. Surprisingly, two water molecules anchored in a hydrogen-bonded chain between Thr-242 and Lys-260 impart elasticity to one side of the binding pocket. The cyclopropane ring takes part in π-acceptor hydrogen bonds with the amide of Met-239. There is a crucial hydrogen bond between the oxygen atoms of the anesthetics and the hydroxyl of Tyr-236. A Tyr-236-Phe mutation results in loss of binding. Thus, both van der Waals interactions and hydrogen-bonding are essential for binding to occur. Ethanol failed to bind because it is too short to benefit from both interactions. Cyclopropylmethanol inhibited phorbol-ester-induced PKCδ activity, but failed to do so in PKCδ containing the Tyr-236-Phe mutation.     

The Binding Mode of ATP Revealed by the Solution Structure of the N-domain of Human ATP7A

We report the solution NMR structures of the N-domain of the Menkes protein (ATP7A) in the ATP-free and ATP-bound forms. The structures consist of a twisted antiparallel six-stranded β-sheet flanked by two pairs of α-helices. A protein loop of 50 amino acids located between β3 and β4 is disordered and mobile on the subnanosecond time scale. ATP binds with an affinity constant of (1.2 ± 0.1) × 10⁴ m⁻¹ and exchanges with a rate of the order of 1 × 10³ s⁻¹. The ATP-binding cavity is considerably affected by the presence of the ligand, resulting in a more compact conformation in the ATP-bound than in the ATP-free form. This structural variation is due to the movement of the α1-α2 and β2-β3 loops, both of which are highly conserved in copper(I)-transporting P_IB-type ATPases. The present structure reveals a characteristic binding mode of ATP within the protein scaffold of the copper(I)-transporting P_IB-type ATPases with respect to the other P-type ATPases. In particular, the binding cavity contains mainly hydrophobic aliphatic residues, which are involved in van der Waal''s interactions with the adenine ring of ATP, and a Glu side chain, which forms a crucial hydrogen bond to the amino group of ATP.     

Interactions between thiamine and large anions. Crystal structures of (H-thiamine)[Pt(SCN)4] · 3H2O and (H-thiamine)[Pt(SCN)6] · H2O

The reaction of thiamine with K₂Pt^IICl₄ and with Pt^IVCl₄ in the presence of excess NaSCN in aqueous solution gave thiamine salts, (H-thiamine)[Pt(SCN)₄] · 3H₂O (1) and (H-thiamine)[Pt(SCN)₆] · H₂O (2), respectively, structures of which have been determined by X-ray diffraction. The thiamine molecule adopts the usual F conformation in each salt. In 1, [Pt(SCN)₄]²⁻ ions act as large planar spacers in the crystal lattice and interact scarcely with thiamine, except for a hydrogen bonding with the terminal hydroxy O(5γ). Instead, water molecules form two types of host–guest-like interactions with the pyrimidine and the thiazolium moieties of a thiamine molecule, one being a C(2)–Hwaterpyrimidine bridge and the other being an N(4′)–Hwaterthiazolium bridge. In 2, despite the much larger ion size, octahedral [Pt(SCN)₆]²⁻ ions form a C(2)–Hanionpyrimidine bridge and an N(4′)–Hanionthiazolium bridge. An additional hydrogen bonding between the anion and the terminal O(5γ) of thiamine creates a hydrogen-bonded macrocyclic ring {thiaminium–[Pt(SCN)₆]²⁻}₂, a supramolecule.     

Consistent picture of the reversible thermal unfolding of hen egg-white lysozyme from experiment and molecular dynamics

Synchrotron radiation circular dichroism, Fourier transform infrared, and nuclear magnetic resonance spectroscopies, and small-angle x-ray scattering were used to monitor the reversible thermal unfolding of hen egg white lysozyme. The results were compared with crystal structures and high- and low-temperature structures derived from molecular-dynamics calculations. The results of both experimental and computational methods indicate that the unfolding process starts with the loss of β-structures followed by the reversible loss of helix content from ∼40% at 20°C to 27% at 70°C and ∼20% at 77°C, beyond which unfolding becomes irreversible. Concomitantly there is a reversible increase in the radius of gyration of the protein from 15 Å to 18 Å. The reversible decrease in forward x-ray scattering demonstrates a lack of aggregation upon unfolding, suggesting the change is due to a larger dilation of hydration water than of bulk water. Molecular-dynamics simulations suggest a similar sequence of events and are in good agreement with the ¹H^N chemical shift differences in nuclear magnetic resonance. This study demonstrates the power of complementary methods for elucidating unfolding/refolding processes and the nature of both the unfolded structure, for which there is no crystallographic data, and the partially unfolded forms of the protein that can lead to fibril formation and disease.     

Changes in temperature have opposing effects on current amplitude in α7 and α4β2 nicotinic acetylcholine receptors

We have examined the effect of temperature on the electrophysiological properties of three neuronal nicotinic acetylcholine receptor (nAChR) subtypes: the rapidly desensitizing homomeric α7 nAChR, the more slowly desensitizing heteromeric α4β2 nAChR and on α7 nAChRs containing a transmembrane mutation (L247T) that results in dramatically reduced desensitization. In all cases, the functional properties of receptors expressed in Xenopus oocytes at room temperature (RT; 21°C) were compared to those recorded at either physiological temperature (37°C) or at lower temperature (4°C). Alterations in temperature had dramatically differing effects on the amplitude of whole-cell responses detected with these three nAChR subtypes. Compared to responses at RT, the amplitude of agonist-evoked responses with α4β2 nAChRs was increased at high temperature (125±9%, n = 6, P<0.01) and reduced at low temperature (47±5%, n = 6, P<0.01), whereas the amplitude of α7 responses was reduced at high temperature (27±7%, n = 11, P<0.001) and increased at low temperatures (224±16%, n = 10, P<0.001). In contrast to the effects of temperature on α4β2 and wild type α7 nAChRs, the amplitude of α7 nAChRs containing the L247T mutation was unaffected by changes in temperature. In addition, changes in temperature had little or no effect on current amplitude when α7 nAChRs were activated by the largely non-desensitizing allosteric agonist 4BP-TQS. Despite these differing effects of temperature on the amplitude of agonist-evoked responses in different nAChRs, changes in temperature had a consistent effect on the rate of receptor desensitization on all subtypes examined. In all cases, higher temperature resulted in increased rates of desensitization. Thus, it appears that the differing effects of temperature on the amplitudes of whole-cell responses cannot be explained by temperature-induced changes in receptor desensitization rates.     

Structural origin of weakly ordered nitroxide motion in spin‐labeled proteins

A disulfide-linked nitroxide side chain (R1) used in site-directed spin labeling of proteins often exhibits an EPR spectrum characteristic of a weakly ordered z-axis anisotropic motion at topographically diverse surface sites, including those on helices, loops and edge strands of β-sheets. To elucidate the origin of this motion, the first crystal structures of R1 that display simple z-axis anisotropic motion at solvent-exposed helical sites (131 and 151) and a loop site (82) in T4 lysozyme have been determined. Structures of 131R1 and 151R1 determined at cryogenic or ambient temperature reveal an intraresidue C_α—H···S_δ interaction that immobilizes the disulfide group, consistent with a model in which the internal motions of R1 are dominated by rotations about the two terminal bonds (Columbus, Kálai, Jeko, Hideg, and Hubbell, Biochemistry 2001;40:3828–3846). Remarkably, the 131R1 side chain populates two rotamers equally, but the EPR spectrum reflects a single dominant dynamic population, showing that the two rotamers have similar internal motion determined by the common disulfide-backbone interaction. The anisotropic motion for loop residue 82R1 is also accounted for by a common disulfide-backbone interaction, showing that the interaction does not require a specific secondary structure. If the above observations prove to be general, then significant variations in order and rate for R1 at noninteracting solvent-exposed helical and loop sites can be assigned to backbone motion because the internal motion is essentially constant.     

Aromatic-aromatic interactions: analysis of π-π interactions in interleukins and TNF proteins

Aromatic-aromatic hydrogen bonds are important in many areas of chemistry, biology and materials science. In this study we have analyzed the roles played by the π-π interactions in interleukins (ILs) and tumor necrosis factor (TNF) proteins. Majority of π-π interacting residues are conserved in ILs and TNF proteins. The accessible surface area calculations in these proteins reveal that these interactions might be important in stabilizing the inner core regions of these proteins. In addition to π-π interactions, the aromatic residues also form π-networks in ILs and TNF proteins. The results obtained in the present study indicate that π-π interactions and π-π networks play important roles in the structural stability of ILs and TNF proteins.     

P2X7 receptor activation mediates organic cation uptake into human myeloid leukaemic KG-1 cells

The P2X7 purinergic receptor is an ATP-gated cation channel with an emerging role in neoplasia. In this study we demonstrate that the human KG-1 cell line, a model of acute myelogenous leukaemia, expresses functional P2X7. RT-PCR and immunochemical techniques demonstrated the presence of P2X7 mRNA and protein respectively in KG-l cells, as well as in positive control multiple myeloma RPMI 8226 cells. Flow cytometric measurements demonstrated that ATP induced ethidium⁺ uptake into KG-l cells suspended in sucrose medium (EC₅₀ of ∼3 μM), but not into cells in NaCl medium. In contrast, ATP induced ethidium⁺ uptake into RPMI 8226 cells suspended in either sucrose or NaCl medium (EC₅₀ of ∼3 or ∼99 μM, respectively), as well as into RPMI 8226 cells in KCl medium (EC₅₀ of ∼18 μM). BzATP and to a lesser extent ATPγS and αβ-methylene ATP, but not ADP or UTP, also induced ethidium⁺ uptake into KG-1 cells. ATP-induced ethidium⁺ uptake was completely impaired by the P2X7 antagonists, AZ10606120 and A-438079. ATP-induced ethidium⁺ uptake was also impaired by probenecid but not by carbenoxolone, both pannexin-1 antagonists. ATP induced YO-PRO-1²⁺ and propidium²⁺ uptake into KG-1 cells. Finally, sequencing of full-length P2X7 cDNA identified several single nucleotide polymorphisms (SNPs) in KG-1 cells including H155Y, A348T, T357S and Q460R. RPMI 8226 cells contained A348T, A433V and H521Q SNPs. In conclusion, the KG-1 cell line expresses functional P2X7. This cell line may help elucidate the signalling pathways involved in P2X7-induced survival and invasiveness of myeloid leukaemic cells.     

Structural basis of phosphodiesterase 6 inhibition by the C‐terminal region of the γ‐subunit

The inhibitory interaction of phosphodiesterase-6 (PDE6) with its γ-subunit (Pγ) is pivotal in vertebrate phototransduction. Here, crystal structures of a chimaeric PDE5/PDE6 catalytic domain (PDE5/6cd) complexed with sildenafil or 3-isobutyl-1-methylxanthine and the Pγ-inhibitory peptide Pγ₇₀₋₈₇ have been determined at 2.9 and 3.0 Å, respectively. These structures show the determinants and the mechanism of the PDE6 inhibition by Pγ and suggest the conformational change of Pγ on transducin activation. Two variable H- and M-loops of PDE5/6cd form a distinct interface that contributes to the Pγ-binding site. This allows the Pγ C-terminus to fit into the opening of the catalytic pocket, blocking cGMP access to the active site. Our analysis suggests that disruption of the H–M loop interface and Pγ-binding site is a molecular cause of retinal degeneration in atrd3 mice. Comparison of the two PDE5/6cd structures shows an overlap between the sildenafil and Pγ₇₀₋₈₇-binding sites, thereby providing critical insights into the side effects of PDE5 inhibitors on vision.     

A new era of the prototropic tautomerism of the quercetin molecule: A QM/QTAIM computational advances

Abstract

In this study for the first time we have revealed and investigated in details 123 different prototropic tautomers of the most stable conformer of the quercetin molecule using quantum-mechanical calculations at the MP2/6-311++G(2df,pd)//B3LYP/6-311++G(d,p) level of QM theory. We have found that in the most energetically favorable prototropic tautomer mobile hydrogen atoms are localized at the О3, О3′, О4′, О5, and О7 exocyclic oxygen atoms. Molecular tautomers are in the range of the Gibbs free energies from 0.0 to 69.8?kcal·mol^?1, while zwitterionic ones – from 30.1 до 172.8?kcal·mol^?1 at normal conditions. It was also reliably established that the weakest point causing the decyclization of the molecule is its C ring – this reaction is launched by the transition of the proton from the C8H group to the endocyclic O1 oxygen atom. All prototropic tautomers, except two cases, are joined by the intramolecular cooperative specific interactions (from 1 to 5) – H-bonds and attractive van der Waals contacts, which have been revealed and characterized by QTAIM analysis.

Communicated by Ramaswamy H. Sarma