RESUMO
CNDOL is an a priori, approximate Fockian for molecular wave functions. In this study, we employ several modes of singly excited configuration interaction (CIS) to model molecular excitation properties by using four combinations of the one electron operator terms. Those options are compared to the experimental and theoretical data for a carefully selected set of molecules. The resulting excitons are represented by CIS wave functions that encompass all valence electrons in the system for each excited state energy. The Coulomb-exchange term associated to the calculated excitation energies is rationalized to evaluate theoretical exciton binding energies. This property is shown to be useful for discriminating the charge donation ability of molecular and supermolecular systems. Multielectronic 3D maps of exciton formal charges are showcased, demonstrating the applicability of these approximate wave functions for modeling properties of large molecules and clusters at nanoscales. This modeling proves useful in designing molecular photovoltaic devices. Our methodology holds potential applications in systematic evaluations of such systems and the development of fundamental artificial intelligence databases for predicting related properties.
RESUMO
The 1,3 dipolar cycloaddition reactions of münchnones and alkenes provide an expedite synthetic way to substituted pyrroles, an exceedingly important structural motif in the pharmaceutical and material science fields of research. The factors governing their regioselectivity rationalization are not well understood. Using several approaches, we investigate a set of 14 reactions (featuring two münchnones, 12 different alkenes, and two alkynes). The Natural Bond Theory and the Non-Covalent Interaction Index analyses of the noncovalent interaction energies fail to predict the experimental major regioisomer. Employing global cDFT descriptors or local ones such as the Fukui function and dual descriptor yields similarly inaccurate predictions. Only the local softness pairing, within Pearson's Hard and Soft Acids and Bases principle, constitutes a reliable predictor for the major reaction product. By taking into account an estimator for the steric effects, the correct regioisomer is predicted. Steric effects play a major role in driving the regioselectivity, as was corroborated by energy decomposition analysis of the transition states.
RESUMO
MOTIVATION: The term clustering designates a comprehensive family of unsupervised learning methods allowing to group similar elements into sets called clusters. Geometrical clustering of molecular dynamics (MD) trajectories is a well-established analysis to gain insights into the conformational behavior of simulated systems. However, popular variants collapse when processing relatively long trajectories because of their quadratic memory or time complexity. From the arsenal of clustering algorithms, HDBSCAN stands out as a hierarchical density-based alternative that provides robust differentiation of intimately related elements from noise data. Although a very efficient implementation of this algorithm is available for programming-skilled users (HDBSCAN*), it cannot treat long trajectories under the de facto molecular similarity metric RMSD. RESULTS: Here, we propose MDSCAN, an HDBSCAN-inspired software specifically conceived for non-programmers users to perform memory-efficient RMSD-based clustering of long MD trajectories. Methodological improvements over the original version include the encoding of trajectories as a particular class of vantage-point tree (decreasing time complexity), and a dual-heap approach to construct a quasi-minimum spanning tree (reducing memory complexity). MDSCAN was able to process a trajectory of 1 million frames using the RMSD metric in about 21 h with <8 GB of RAM, a task that would have taken a similar time but more than 32 TB of RAM with the accelerated HDBSCAN* implementation generally used. AVAILABILITY AND IMPLEMENTATION: The source code and documentation of MDSCAN are free and publicly available on GitHub (https://github.com/LQCT/MDScan.git) and as a PyPI package (https://pypi.org/project/mdscan/). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Assuntos
Simulação de Dinâmica Molecular , Software , Análise por Conglomerados , AlgoritmosRESUMO
MOTIVATION: Density Peaks is a widely spread clustering algorithm that has been previously applied to Molecular Dynamics (MD) simulations. Its conception of cluster centers as elements displaying both a high density of neighbors and a large distance to other elements of high density, particularly fits the nature of a geometrical converged MD simulation. Despite its theoretical convenience, implementations of Density Peaks carry a quadratic memory complexity that only permits the analysis of relatively short trajectories. RESULTS: Here, we describe DP+, an exact novel implementation of Density Peaks that drastically reduces the RAM consumption in comparison to the scarcely available alternatives designed for MD. Based on DP+, we developed RCDPeaks, a refined variant of the original Density Peaks algorithm. Through the use of DP+, RCDPeaks was able to cluster a one-million frames trajectory using less than 4.5 GB of RAM, a task that would have taken more than 2 TB and about 3× more time with the fastest and less memory-hunger alternative currently available. Other key features of RCDPeaks include the automatic selection of parameters, the screening of center candidates and the geometrical refining of returned clusters. AVAILABILITY AND IMPLEMENTATION: The source code and documentation of RCDPeaks are free and publicly available on GitHub (https://github.com/LQCT/RCDPeaks.git). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Assuntos
Algoritmos , Simulação de Dinâmica Molecular , Software , Análise por Conglomerados , DocumentaçãoRESUMO
MOTIVATION: Classical Molecular Dynamics (MD) is a standard computational approach to model time-dependent processes at the atomic level. The inherent sparsity of increasingly huge generated trajectories demands clustering algorithms to reduce other post-simulation analysis complexity. The Quality Threshold (QT) variant is an appealing one from the vast number of available clustering methods. It guarantees that all members of a particular cluster will maintain a collective similarity established by a user-defined threshold. Unfortunately, its high computational cost for processing big data limits its application in the molecular simulation field. RESULTS: In this work, we propose a methodological parallel between QT clustering and another well-known algorithm in the field of Graph Theory, the Maximum Clique Problem. Molecular trajectories are represented as graphs whose nodes designate conformations, while unweighted edges indicate mutual similarity between nodes. The use of a binary-encoded RMSD matrix coupled to the exploitation of bitwise operations to extract clusters significantly contributes to reaching a very affordable algorithm compared to the few implementations of QT for MD available in the literature. Our alternative provides results in good agreement with the exact one while strictly preserving the collective similarity of clusters. AVAILABILITY AND IMPLEMENTATION: The source code and documentation of BitQT are free and publicly available on GitHub (https://github.com/LQCT/BitQT.git) and ReadTheDocs (https://bitqt.readthedocs.io/en/latest/), respectively. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Assuntos
Simulação de Dinâmica Molecular , Software , Algoritmos , Análise por Conglomerados , DocumentaçãoRESUMO
Ghrelin is a peptide hormone involved in multiple functions, including growth hormone release stimulation, food intake regulation, and metabolic and cytoprotective effect. A novel family of peptides with internal cycles was designed as ghrelin analogs and the biological activity of two of them (A228 and A233) was experimentally studied in-depth. In this work, an in silico strategy was developed for describing and assessing the binding modes of A228 and A233 to GHS-R1a (ghrelin receptor) comparing it with ghrelin and GHRP-6 peptides. Several reported structures of different G protein coupled receptors were used as templates, to obtain a good quality model of GHS-R1a. The best model was selected by preliminary molecular docking with ghrelin and GHRP-6. Docking was used to estimate peptide orientations in the binding site of the best model, observing a superposition of its N-terminal and its first aromatic residue. To test the complex stability in time, the C-terminal fragments of each peptide were added and the complexes were inserted a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane, performing a molecular dynamic simulation for 100 ns using the CHARMM36 force field. Despite of the structural differences, the studied peptides share a common binding mode; the N-terminal interacts with E124 and the aromatic residue close to it, with the aromatic cluster (F279, F309, and F312). A preliminary pharmacophore model, consisting in a positive charged amine and an aromatic ring at an approximate distance of 0.79 nm, can be proposed. The results here described could represent a step forward in the efficient search of new ghrelin analogs.
Assuntos
Simulação por Computador , Peptídeos/metabolismo , Peptídeos/farmacologia , Receptores de Grelina/agonistas , Sequência de Aminoácidos , Animais , Humanos , Simulação de Acoplamento Molecular , Simulação de Dinâmica Molecular , Peptídeos/química , Ligação Proteica/efeitos dos fármacos , Receptores de Grelina/química , Receptores de Grelina/metabolismoRESUMO
We followed a comprehensive computational strategy to understand and eventually predict the structure-activity relationship of thirty-three 1,3-disubstituted imidazole [1,5-α] pyrazine derivatives described as ATP competitive inhibitors of the IGF-1 receptor related to Ewing sarcoma. The quantitative structure-activity relationship model showed that the inhibitory potency is correlated with the molar volume, a steric descriptor and the net charge calculated value on atom C1 (q1) and N4 (q4) of the pharmacophore, all of them appearing to give a positive contribution to the inhibitory activity. According to experimental and calculated values, the most potent compound would be 3-[4-(azetidin-2-ylmethyl) cyclohexyl]-1-[3-(benzyloxy) phenyl] imidazo [1,5-α]pyrazin-8-amine (compound 23). Docking was used to guess important residues involved in the ATP-competitive inhibitory activity. It was validated by 200 ns of molecular dynamics (MD) simulation using improved linear interaction energy (LIE) method. MD of previously preferred structures by docking shows that the most potent ligand could establish hydrogen bonds with the ATP-binding site of the receptor, and the Ser979 and Ser1059 residues contribute favourably to the binding stability of compound 23. MD simulation also gave arguments about the chemical structure of the compound 23 being able to fit in the ATP-binding pocket, expecting to remain stable into it during the entire simulation and allowing us to hint the significant contribution expected to be given by electrostatic and hydrophobic interactions to the ligand-receptor complex stability. This computational combined strategy here described could represent a useful and effective prime approach to guide the identification of tyrosine kinase inhibitors as new lead compounds.
Assuntos
Trifosfato de Adenosina/química , Antineoplásicos/química , Imidazóis/química , Modelos Moleculares , Pirazinas/química , Relação Quantitativa Estrutura-Atividade , Receptor IGF Tipo 1/química , Animais , Antineoplásicos/farmacologia , Sítios de Ligação , Ligação Competitiva , Linhagem Celular , Humanos , Interações Hidrofóbicas e Hidrofílicas , Ligantes , Camundongos , Simulação de Acoplamento Molecular , Simulação de Dinâmica Molecular , Estrutura Molecular , Ligação Proteica , Pirazinas/farmacologia , Receptor IGF Tipo 1/antagonistas & inibidores , Reprodutibilidade dos TestesRESUMO
Alzheimer's disease is a progressive neurodegenerative disorder characterized by the abnormal processing of the Tau and the amyloid precursor proteins. The unusual aggregation of Tau is based on the formation of intermolecular ß-sheets through two motifs: 275 VQIINK280 and 306 VQIVYK311 . Phenylthiazolyl-hydrazides (PTHs) are capable of inhibiting/disassembling Tau aggregates. However, the disaggregation mechanism of Tau oligomers by PTHs is still unknown. In this work, we studied the disruption of the oligomeric form of the Tau motif 306 VQIVYK311 by PTHs through molecular docking, molecular dynamics, and free energy calculations. We predicted hydrophobic interactions as the major driving forces for the stabilization of Tau oligomer, with V306 and I308 being the major contributors. Nonpolar component of the binding free energy is essential to stabilize Tau-PTH complexes. PTHs disrupted mainly the van der Waals interactions between the monomers, leading to oligomer destabilization. Destabilization of full Tau filament by PTHs and emodin was not observed in the sampled 20 ns; however, in all cases, the nonpolar component of the binding free energy is essential for the formation of Tau filament-PTH and Tau filament-emodin. These results provide useful clues for the design of more effective Tau-aggregation inhibitors.
Assuntos
Hidrazinas/farmacologia , Agregados Proteicos , Tiazóis/farmacologia , Proteínas tau/antagonistas & inibidores , Proteínas tau/química , Motivos de Aminoácidos , Emodina/farmacologia , Hidrazinas/química , Ligação de Hidrogênio , Modelos Moleculares , Conformação Molecular , Agregados Proteicos/efeitos dos fármacos , Termodinâmica , Tiazóis/químicaRESUMO
The aim of this work is to describe the molecular inclusion of chlordecone with α-, ß-, and γ-cyclodextrin in aqueous solution using quantum mechanics. The guest-host complexes of chlordecone and cyclodextrins are modeled in aqueous solution using the multiple minima hypersurface methodology with a PM6-D3H4X semiempirical Hamiltonian, and the lowest energy minima obtained are reoptimized using the M06-2X density functional and the intermolecular interactions described using quantum theory of atoms in molecules (QTAIM). The studied complexes are classified according to the degree of inclusion, namely, total occlusion, partial occlusion, and external interaction. More stable complexes are obtained when γ-CD is used as the host molecule. The interactions characterized through QTAIM analysis are all of electrostatic nature, predominantly of dispersive type. In this work, a method based on the counterpoise correction is also discussed to mitigate the basis set superposition error in density functional theory calculations when using an implicit solvation model.
Assuntos
Clordecona , Ciclodextrinas , Teoria Quântica , Eletricidade Estática , ÁguaRESUMO
Clustering Molecular Dynamics trajectories is a common analysis that allows grouping together similar conformations. Several algorithms have been designed and optimized to perform this routine task, and among them, Quality Threshold stands as a very attractive option. This algorithm guarantees that in retrieved clusters no pair of frames will have a similarity value greater than a specified threshold, and hence, a set of strongly correlated frames are obtained for each cluster. In this work, it is shown that various commonly used software implementations are flawed by confusing Quality Threshold with another simplistic well-known clustering algorithm published by Daura et al. (Daura, X.; van Gunsteren, W. F.; Jaun, B.; Mark, A. E.; Gademann, K.; Seebach, D. Peptide Folding: When Simulation Meets Experiment. Angew. Chemie Int. Ed. 1999, 38 (1/2), 236-240). Daura's algorithm does not impose any quality threshold for the frames contained in retrieved clusters, bringing unrelated structural configurations altogether. The advantages of using Quality Threshold whenever possible to explore Molecular Dynamic trajectories is exemplified. An in-house implementation of the original Quality Threshold algorithm has been developed in order to illustrate our comments, and its code is freely available for further use by the scientific community.
Assuntos
Simulação de Dinâmica Molecular , Algoritmos , Análise por Conglomerados , Teoria QuânticaRESUMO
The growing computational capacity allows the investigation of large biomolecular systems by increasingly extensive molecular dynamics simulations. The resulting huge trajectories demand efficient partition methods to discern relevant structural dissimilarity. Clustering algorithms are available to address this task, but their implementations still need to be improved to gain in computational speed and to reduce the consumption of random access memory. We propose the BitClust code which, based on a combination of Python and C programming languages, performs fast structural clustering of long molecular trajectories. BitClust takes advantage of bitwise operations applied to a bit-encoded pairwise similarity matrix. Our approach allowed us to process a half-million frame trajectory in 6 h using less than 35 GB, a task that is not affordable with any of the similar alternatives.
Assuntos
Algoritmos , Simulação de Dinâmica Molecular , Análise por Conglomerados , Linguagens de Programação , Fatores de TempoRESUMO
The regioselectivity of the 1,3-dipolar cycloaddition of a model nitrone with a set of dipolarophiles, presenting diverse electronic effects, is analyzed using conceptual density functional theory (DFT) methods. We deviate from standard approaches based on frontier molecular orbitals and formulations of the local hard/soft acid/base principle and use instead the dual descriptor. A detailed analysis is carried out to determine the influence of the way to calculate the dual descriptor, the computational procedure, basis set and choice of method to condensate the values of this descriptor. We show that the qualitative regioselectivity predictions depend on the choice of "computational conditions", something that indicates the danger of using black-box computational set-ups in conceptual DFT studies.
RESUMO
A novel heuristic using an iterative select-and-purge strategy is proposed. It combines statistical techniques for sampling and classification by rigid molecular docking through an inverse virtual screening scheme. This approach aims to the de novo discovery of short peptides that may act as docking receptors for small target molecules when there are no data available about known association complexes between them. The algorithm performs an unbiased stochastic exploration of the sample space, acting as a binary classifier when analyzing the entire peptides population. It uses a novel and effective criterion for weighting the likelihood of a given peptide to form an association complex with a particular ligand molecule based on amino acid sequences. The exploratory analysis relies on chemical information of peptides composition, sequence patterns, and association free energies (docking scores) in order to converge to those peptides forming the association complexes with higher affinities. Statistical estimations support these results providing an association probability by improving predictions accuracy even in cases where only a fraction of all possible combinations are sampled. False positives/false negatives ratio was also improved with this method. A simple rigid-body docking approach together with the proper information about amino acid sequences was used. The methodology was applied in a retrospective docking study to all 8000 possible tripeptide combinations using the 20 natural amino acids, screened against a training set of 77 different ligands with diverse functional groups. Afterward, all tripeptides were screened against a test set of 82 ligands, also containing different functional groups. Results show that our integrated methodology is capable of finding a representative group of the top-scoring tripeptides. The associated probability of identifying the best receptor or a group of the top-ranked receptors is more than double and about 10 times higher, respectively, when compared to classical random sampling methods.
Assuntos
Biologia Computacional/métodos , Avaliação Pré-Clínica de Medicamentos/métodos , Peptídeos/metabolismo , Receptores Artificiais/metabolismo , Bibliotecas de Moléculas Pequenas/metabolismo , Bibliotecas de Moléculas Pequenas/farmacologia , Algoritmos , Ligantes , Simulação de Acoplamento Molecular , Curva ROC , Processos Estocásticos , Termodinâmica , Interface Usuário-ComputadorRESUMO
Pyrrolyl-silicon compounds were investigated by different theoretical approaches. Model monomers consisted of a pyrrole ring N-substituted with silylmethoxy and silylhydroxy end groups through a propyl chain spacer, designated as PySi and PySiOH. Geometrical, vibrational, and electronic properties, as well as chemical reactivity, are discussed and compared with pyrrole (Py) and N-propylpyrrole (N-PrPy) that were studied in parallel for reference purposes and methods validation. The electronic distribution between PySi and PySiOH differs importantly, the former being an electron donor, as Py and N-PrPy. Conversely, PySiOH presents donor-acceptor character with the LUMO energy level localized on the silanol end group. Global and local reactivity descriptors predict PySiOH more reactive than PySi with two preferential reactive sites: electron-rich Py ring and electron-deficient silanol group. On the basis of experimental studies, oligomers of PySiOH linked α-α' via Py rings (α-α'PynSiOH, n = 2, 3) were considered as model molecules of hydrolyzed PySi. The most stable structures were derived from randomly generated α-α'PynSiOH that were optimized at semiempirical AM1 and refined with M05-2X/6-31G(d,p). Conformational analysis of dimer and trimer structures points to stability enhanced by molecular packing. Nonetheless, NBO and RDG results indicate that oligomer stability is dictated by the cooperative contribution of hydrogen bonding between silanol end groups and dispersive vdW interactions between silanol and the π system of the Py ring. The latter interaction resulting from electron delocalization induced by an electron-deficient silanol group seems to determine the smaller gap energy of T-shaped OH-π arrangements. The theoretical findings support the peculiar chemical behavior revealed by experiment.
RESUMO
Spectral shifts of rhodopsin, which are related to variations of the electron distribution in 11-cis-retinal, are investigated here using the method of deformed atoms in molecules. We found that systems carrying the M207R and S186W mutations display large perturbations of the π-conjugated system with respect to wild-type rhodopsins. These changes agree with the predicted behavior of the bond length alternation (BLA) and the blue shifts of vertical excitation energies of these systems. The effect of the planarity of the central and Schiff-base regions of retinal chain on the electronic structure of the chromophore is also investigated. By establishing nonlinear polynomial relations between BLA, chain distortions, and vertical excitation energies, we are also able to provide a semiquantitative approach for the understanding of the mechanisms regulating spectral shifts in rhodopsin and its mutants.
Assuntos
Elétrons , Retinaldeído/química , Rodopsina/química , Animais , Bovinos , Humanos , Modelos Moleculares , Simulação de Dinâmica Molecular , Mutação , Rodopsina/genética , Eletricidade EstáticaRESUMO
Quantitative structure-activity relationship (QSAR) study of 19-nor-testosterone steroids family was performed using quantum and physicochemical molecular descriptors. The quantum-chemical descriptors were calculated using semiempirical calculations. The descriptor values were statistically correlated using multi-linear regression analysis. The QSAR study indicated that the electronic properties of these derivatives have significant relationship with observed biological activities. The found QSAR equations explain that the energy difference between the LUMO and HOMO, the total dipole moment, the chemical potential and the value of the net charge of different carbon atoms in the steroid nucleus showed key interaction of these steroids with their anabolic-androgenic receptor binding site. The calculated values predict that the 17α-cyclopropyl-17ß, 3ß-hydroxy-4-estrene compound presents the highest anabolic-androgenic ratio (AAR) and the 7α-methyl-17ß-acetoxy-estr-4-en-3-one compound the lowest AAR. This study might be helpful in the future successful identification of "real" or "virtual" anabolic-androgenic steroids.
Assuntos
Anabolizantes/química , Anabolizantes/farmacologia , Androgênios/química , Androgênios/farmacologia , Nandrolona/análogos & derivados , Nandrolona/química , Animais , Desenho de Fármacos , Masculino , Modelos Químicos , Estrutura Molecular , Músculo Esquelético/efeitos dos fármacos , Nandrolona/farmacologia , Próstata/efeitos dos fármacos , Relação Quantitativa Estrutura-Atividade , Ratos , Ratos Wistar , Glândulas Seminais/efeitos dos fármacosRESUMO
The interaction of the most active natural brassinosteroid, brassinolide, with the twenty natural amino acids is studied applying the multiple minima hypersurface method to model the molecular interactions explicitly. The resulting thermodynamic data gives useful information about the amino acids with the greatest association for brassinolide and the stabilities of such complexes. Density functional theory (DFT) optimizations were further carried out to test the performance of semiempirical calculations. Additional calculations with a more accurate DFT method were performed to explore the formation of this type of molecular complexes. The semiempirical geometries and stability order of these complexes are in good agreement with the DFT calculations. Each group of amino acids possesses a preferential zone of interaction with brassinolide, forming the polar-charged amino acids the most stable complexes. This study could contribute to future investigations of the interaction of brassinosteroids with the receptor protein in plants.
Assuntos
Aminoácidos Essenciais/química , Colestanóis/química , Modelos Moleculares , Esteroides Heterocíclicos/química , Brassinosteroides , TermodinâmicaRESUMO
The interaction of three different brassinosteroids with water was studied by the Multiple Minima Hypersurface (MMH) procedure to model molecular interactions explicitly. The resulting thermodynamic data give useful information on properties of molecular association with water. This application can serve as a tool for future investigations and modelling concerning interactions of brassinosteroids with receptor proteins in plants. DFT/B3LYP calculations were also made in order to correlate and test the performance of the current AM1 Hamiltonian calculations of these complexes, which are inherent to MMH routine. Diol functionalities located in ring A and lateral chain appears as the sites that show the highest affinity to water. The oxalactone group does not appear to be a key structural requirement in the association with water. Parallel calculations with a "polarizable continuum method" (PCM) agreed with the reported experimental order of biological activities, where Brassinolide exhibited the best solubility features.
Assuntos
Colestanóis/química , Modelos Teóricos , Reguladores de Crescimento de Plantas/química , Esteroides Heterocíclicos/química , Água/química , Brassinosteroides , Estrutura Molecular , Solubilidade , TermodinâmicaRESUMO
Predictive quantitative structure-activity relationship (QSAR) models of anabolic and androgenic activities for the testosterone and dihydrotestosterone steroid analogues were obtained by means of multiple linear regression using quantum and physicochemical molecular descriptors (MD) as well as a genetic algorithm for the selection of the best subset of variables. Quantitative models found for describing the anabolic (androgenic) activity are significant from a statistical point of view: R(2) of 0.84 (0.72 and 0.70). A leave-one-out cross-validation procedure revealed that the regression models had a fairly good predictability [q(2) of 0.80 (0.60 and 0.59)]. In addition, other QSAR models were developed to predict anabolic/androgenic (A/A) ratios and the best regression equation explains 68% of the variance for the experimental values of AA ratio and has a rather adequate q(2) of 0.51. External validation, by using test sets, was also used in each experiment in order to evaluate the predictive power of the obtained models. The result shows that these QSARs have quite good predictive abilities (R(2) of 0.90, 0.72 (0.55), and 0.53) for anabolic activity, androgenic activity, and A/A ratios, respectively. Last, a Williams plot was used in order to define the domain of applicability of the models as a squared area within +/-2 band for residuals and a leverage threshold of h=0.16. No apparent outliers were detected and the models can be used with high accuracy in this applicability domain. MDs included in our QSAR models allow the structural interpretation of the biological process, evidencing the main role of the shape of molecules, hydrophobicity, and electronic properties. Attempts were made to include lipophilicity (octanol-water partition coefficient (logP)) and electronic (hardness (eta)) values of the whole molecules in the multivariate relations. It was found from the study that the logP of molecules has positive contribution to the anabolic and androgenic activities and high values of eta produce unfavorable effects. The found MDs can also be efficiently used in similarity studies based on cluster analysis. Our model for the anabolic/androgenic ratio (expressed by weight of levator ani muscle, LA, and seminal vesicle, SV, in mice) predicts that the 2-aminomethylene-17alpha-methyl-17beta-hydroxy-5alpha-androstan-3-one (43) compound is the most potent anabolic steroid, and the 17alpha-methyl-2beta,17beta-dihydroxy-5alpha-androstane (31) compound is the least potent one of this series. The approach described in this report is an alternative for the discovery and optimization of leading anabolic compounds among steroids and analogues. It also gives an important role to electron exchange terms of molecular interactions to this kind of steroid activity.
Assuntos
Anabolizantes/química , Anabolizantes/farmacologia , Androgênios/química , Androgênios/farmacologia , Di-Hidrotestosterona/análogos & derivados , Modelos Químicos , Testosterona/análogos & derivados , Algoritmos , Androgênios/genética , Análise por Conglomerados , Simulação por Computador , Humanos , Masculino , Relação Quantitativa Estrutura-Atividade , Testosterona/genéticaRESUMO
We propose two models of the human S-arrestin/rhodopsin complex in the inactive dark adapted rhodopsin and meta rhodopsin II form, obtained by homology modeling and knowledge based docking. First, a homology model for the human S-arrestin was built and validated by molecular dynamics, showing an average root mean square deviation difference from the pattern behavior of 0.76 A. Then, combining the human S-arrestin model and the modeled structure of the two human rhodopsin forms, we propose two models of interaction for the human S-arrestin/rhodopsin complex. The models involve two S-arrestin regions related to the N domain (residues 68-78; 170-182) and a third constituent of the C domain (248-253), with the rhodopsin C terminus (330-348). Of the 22 single point mutations related to retinitis pigmentosa and congenital night blindness located in the cytoplasmatic portion of rhodopsin or in S-arrestin, our models locate 16 in the interaction region and relate two others to possible dimer formation. Our calculations also predict that the light activated complex is more stable than the dark adapted rhodopsin and, therefore, of higher affinity to S-arrestin.