Exploring new scaffolds for angiotensin II receptor antagonism

Nowadays, AT₁ receptor (AT₁R) antagonists (ARBs) constitute the one of the most prevalent classes of antihypertensive drugs that modulate the renin-angiotensin system (RAS). Their main uses include also treatment of diabetic nephropathy (kidney damage due to diabetes) and congestive heart failure. Towards this direction, our study has been focused on the discovery of novel agents bearing different scaffolds which may evolve as a new class of AT₁ receptor antagonists. To fulfill this aim, a combination of computational approaches and biological assays were implemented. Particularly, a pharmacophore model was established and served as a 3D search query to screen the ChEMBL15 database. The reliability and accuracy of virtual screening results were improved by using molecular docking studies. In total, 4 compounds with completely diverse chemical scaffolds from potential ARBs, were picked and tested for their binding affinity to AT₁ receptor. Results revealed high nanomolar to micromolar affinity (IC₅₀) for all the compounds. Especially, compound 4 exhibited a binding affinity of 199 nM. Molecular dynamics simulations were utilized in an effort to provide a molecular basis of their binding to AT₁R in accordance to their biological activities.     

Computational Investigations of the Chalcogen-Substituted Carboxylic Acids RC(=O)XH and RC(=X)OH and their Dimers [RC(=O)XH]2 and [RC(=X)OH]2 (X=S, Se, Te)

Semiempirical and ab initio theoretical methods have been used to investigate molecular structures of the chalcogen-substituted carboxylic acid isomers RC(=O)XH (chalcogenol acid) and RC(=X)OH (chalcogenon acid). A recent experimental report suggests that the chalcogenon isomers, although less stable at room temperature, predominate at low temperature in polar solvents and that there is only a small barrier to isomerization between the isomers. Theoretical calculations have been used to locate minimum energy structures of chalcogen-substituted carboxylic acid isomers and to calculate energy differences between pairs of isomers. Carboxylic acids are well known to dimerize, especially in the gas phase and in non-polar solvents. We have, therefore, also calculated energies of dimerization of the chalcogen-substituted acids by optimizing the geometries of the symmetric dimers. We note that the PM3 level of theory is only qualitatively correct for sulfur- and selenium-containing species but fails even qualitatively for the tellurium-containing compounds. Ab initio results confirm the experimental observations and provide good estimates of both isomerization and dimerization energies. We conclude that for many functional groups with tautomers RC(=X)YH and RC(=Y)XH, the more acidic tautomer is the one with the acid proton on the smaller, more electronegative atom, although in many cases this may not be the more stable tautomer.Electronic Supplementary Material available.     

Estimation of uncertainties in X-ray refinement results by use of perturbed structures

The uncertainties in the refined parameters for a 1.5-A X-ray structure of carbon-monoxy (FeII) myoglobin are estimated by combining energy minimization with least-squares refinement against the X-ray data. The energy minimizations, done without reference to the X-ray data, provide perturbed structures which are used to restart conventional X-ray refinement. The resulting refined structures have the same, or better, R-factor and stereochemical parameters as the original X-ray structure, but deviate from it by 0.13 A rms for the backbone atoms and 0.31 A rms for the sidechain atoms. Atoms interacting with a disordered sidechain, Arg 45 CD3, are observed to have larger positional uncertainties. The uncertainty in the B-factors, within the isotropic harmonic motion approximation, is estimated to be 15%. The resulting X-ray structures are more consistent with the energy parameters used in simulations.     

Kinetic and molecular orbital studies on the rate of oxidation of monosubstituted phenols and anilines by lactoperoxidase compound II in comparison with the case of horseradish peroxidase

The second-order rate constant (k4) for the oxidation of monosubstituted phenols and anilines by lactoperoxidase compound II was examined by Chance's method [B. Chance, Arch. Biochem. Biophys. 71 (1957), 130–136]. When the electronic states of these substrates were calculated by an ab initio molecular orbital method, it was found that the log k4 value correlates well with the highest occupied molecular orbital (HOMO) energy level but not with the net charge or frontier electron density. These results are essentially similar to those reported previously in the case of horseradish peroxidase [J. Sakurada, R. Sekiguchi, K. Sato, and T. Hosoya, Biochemistry 29 (1990), 4093–4098], showing some dissimilar features which are considered to reflect the structural difference between the two enzymes.Abbreviations HOMO highest occupied molecular orbital - HRP horseradish peroxidase - LPO lactoperoxidase (EC 1.11.1.7) - LUMO lowest unoccupied molecular orbital     

Aerodynamics of Ephedra trifurca

Computer simulations are used to predict the behavior of pollen grains with different physical properties within the acceleration field created around the ovules of the gymnosperm Ephedra trifurca. A modelling procedure is given that (1) calculates the number of pollen grains captured by an ovule's pollination-droplet and (2) gives a correlation between pollination efficiency and the physical properties (= mass, size) of different types of pollen. Based on this procedure, the number of Ephedra pollen grains captured by micropyles can be less than the number captured from other species. However, the mass and size of Ephedra pollen grains appear to coincide with those predicted to yield a local maximum of pollination efficiency, i.e. slightly larger or smaller values of either mass or size would decrease the probability of capture. In addition, the properties of Ephedra pollen grains operate synergistically in the aerodynamic environment around ovules and are focused to collide with pollination-droplets. By analogy, the properties of Ephedra pollen coincide with those predicted for a localized adaptive peak. The physical properties of pollen grain types other than E. trifurca that can maximize pollen capture are not generally represented in the aerobiology of Ephedra during the pollination season. Therefore, the phenology of pollen release, community taxonomic-composition, and the physics of particle capture play collectively important roles in the reproductive success of Ephedra trifurca.     

Generation of rhythmic and arrhythmic behaviour of Crassulacean acid metabolism in Kalanchoë daigremontiana under continuous light by varying the irradiance or temperature: Measurements in vivo and model simulations

Leaves of Kalanchoë daigremontiana Hamet et Perr. at a photon flux density (PFD) above 220 mol·m^–2s^–1 (400–700 nm) or at leaf temperatures above 27.0 °C showed a rapid loss of rhythmicity, and a more or less pronounced damping-out of the endogenous circadian rhythm of CO₂ exchange under continuous illumination. This rhythm was reinitiated after reduction of the PFD by 90–120 mol·m^–2·s^–1 or reduction of leaf temperature by 3.5–11.0 °C under otherwise unchanged external conditions. The reduction in the magnitude of the external control parameter of the Crassulacean acid metabolism (CAM) rhythm (i.e. PFD or leaf temperature) set the phase of the new rhythm. The maxima of CO₂ uptake occurred about 5, 28, 51, 75 h after the reduction. Simulations with a CAM model under comparable conditions showed a similar behaviour. The influence of temperature on the endogenous CAM rhythm observed in K. daigremontiana in vivo could be simulated by incorporating into the model temperature-dependent switch modes for passive efflux of malate from the vacuole to the cytoplasm. Thus, the model indicates that tonoplast function plays an important role in regulation of the endogenous CAM rhythm in K. daigremontiana.Abbreviations CAM Crassulacean acid metabolism - PAR photosynthetically active radiation - PFD photon flux density This work was supported by a grant to F.B. and U.L. from Teilprojekt B5 in the Sonderforschungsbereich 199 of the Deutsche Forschungsgemeinschaft (Bonn, Germany) and by a grant to T. E. E. G. from the Sudienstiftung des deutschen Volkes (Bonn, Germany). Erika Ball is thanked for processing of time-course data for the analysis of Fourier spectra.     

Prerequisite for His4 in deltorphin A for highδ opioid receptor selectivity

Summary Analysis of deltorphin A position 4 analogues included: backbone constrained N MeHis, spinacine (Spi), N MePhe and the tetrahydroisoquinoline-3-carboxylic acid (Tic); spatially confined side-chain (Phg); and imidazole alkylation ofl- andd-His⁴ enantiomers. High selectivity was lost with the following replacements: N MeHis⁴, N MePhe⁴ and Phg⁴ reduced binding and the constrained residues also increasedµ binding; ring closure between the side-chain and amino group to yield Spi⁴ or Tic⁴ increasedµ affinity. Imidazole methylation of His⁴ marginally affected opioid binding and doubled selectivity; alkylatedd-His⁴-derivatives generally maintained selectivity in spite of decreased affinities. Thus, His⁴ imidazole preserves selectivity by facilitating high binding and by repulsion at theµ receptor. Several low energy conformers of deltorphin A indicated that the His⁴ imidazole preferred a spatial orientation parallel to the phenolic side-chain of Tyr¹ suggestive that this conformation might contribute to high affinity and selectivity.     

Optimization of the electrostatic interactions in proteins of different functional and folding type

The 3-dimensional optimization of the electrostatic interactions between the charged amino acid residues was studied by Monte Carlo simulations on an extended representative set of 141 protein structures with known atomic coordinates. The proteins were classified by different functional and structural criteria, and the optimization of the electrostatic interactions was analyzed. The optimization parameters were obtained by comparison of the contribution of charge-charge interactions to the free energy of the native protein structures and for a large number of randomly distributed charge constellations obtained by the Monte Carlo technique. On the basis of the results obtained, one can conclude that the charge-charge interactions are better optimized in the enzymes than in the proteins without enzymatic functions. Proteins that belong to the mixed αβ folding type are electrostatically better optimized than pure α-helical or β-strand structures. Proteins that are stabilized by disulfide bonds show a lower degree of electrostatic optimization. The electrostatic interactions in a native protein are effectively optimized by rejection of the conformers that lead to repulsive charge-charge interactions. Particularly, the rejection of the repulsive contacts seems to be a major goal in the protein folding process. The dependence of the optimization parameters on the choice of the potential function was tested. The majority of the potential functions gave practically identical results.     

Integrating experimental data with molecular simulations to investigate RNA structural dynamics

Conformational dynamics is crucial for ribonucleic acid (RNA) function. Techniques such as nuclear magnetic resonance, cryo-electron microscopy, small- and wide-angle X-ray scattering, chemical probing, single-molecule Förster resonance energy transfer, or even thermal or mechanical denaturation experiments probe RNA dynamics at different time and space resolutions. Their combination with accurate atomistic molecular dynamics (MD) simulations paves the way for quantitative and detailed studies of RNA dynamics. First, experiments provide a quantitative validation tool for MD simulations. Second, available data can be used to refine simulated structural ensembles to match experiments. Finally, comparison with experiments allows for improving MD force fields that are transferable to new systems for which data is not available. Here we review the recent literature and provide our perspective on this field.