Mutations on RBD of SARS-CoV-2 Omicron variant result in stronger binding to human ACE2 receptor.

The COVID-19 pandemic caused by the SARS-CoV-2 virus has led to more than 270 million infections and 5.3 million of deaths worldwide. Several major variants of SARS-CoV-2 have emerged and posed challenges in controlling the pandemic. The recently occurred Omicron variant raised serious concerns about reducing the efficacy of vaccines and neutralization antibodies due to its vast mutations. We have modelled the complex structure of the human ACE2 protein and the receptor binding domain (RBD) of Omicron Spike protein (S-protein), and conducted atomistic molecular dynamics simulations to study the binding interactions. The analysis shows that the Omicron RBD binds more strongly to the human ACE2 protein than the original strain. The mutations at the ACE2-RBD interface enhance the tight binding by increasing hydrogen bonding interaction and enlarging buried solvent accessible surface area.

Using molecular dynamics for the refinement of atomistic models of GPCRs by homology modeling.

Despite GPCRs sharing a common seven helix bundle, analysis of the diverse crystallographic structures available reveal specific features that might be relevant for ligand design. Despite the number of crystallographic structures of GPCRs steadily increasing, there are still challenges that hamper the availability of new structures. In the absence of a crystallographic structure, homology modeling remains one of the important techniques for constructing 3D models of proteins. In the present study we investigated the use of molecular dynamics simulations for the refinement of GPCRs models constructed by homology modeling. Specifically, we investigated the relevance of template selection, ligand inclusion as well as the length of the simulation on the quality of the GPCRs models constructed. For this purpose we chose the crystallographic structure of the rat muscarinic M3 receptor as reference and constructed diverse atomistic models by homology modeling, using different templates. Specifically, templates used in the present work include the human muscarinic M2; the more distant human histamine H1 and the even more distant bovine rhodopsin as shown in the GPCRs phylogenetic tree. We also investigated the use or not of a ligand in the refinement process. Hence, we conducted the refinement process of the M3 model using the M2 muscarinic as template with tiotropium or NMS docked in the orthosteric site and compared with the results obtained with a model refined without any ligand bound.

Flavonoid allosteric modulation of mutated visual rhodopsin associated with retinitis pigmentosa.

Dietary flavonoids exhibit many biologically-relevant functions and can potentially have beneficial effects in the treatment of pathological conditions. In spite of its well known antioxidant properties, scarce structural information is available on the interaction of flavonoids with membrane receptors. Advances in the structural biology of a specific class of membrane receptors, the G protein-coupled receptors, have significantly increased our understanding of drug action and paved the way for developing improved therapeutic approaches. We have analyzed the effect of the flavonoid quercetin on the conformation, stability and function of the G protein-coupled receptor rhodopsin, and the G90V mutant associated with the retinal degenerative disease retinitis pigmentosa. By using a combination of experimental and computational methods, we suggest that quercetin can act as an allosteric modulator of opsin regenerated with 9-cis-retinal and more importantly, that this binding has a positive effect on the stability and conformational properties of the G90V mutant associated with retinitis pigmentosa. These results open new possibilities to use quercetin and other flavonoids, in combination with specific retinoids like 9-cis-retinal, for the treatment of retinal degeneration associated with retinitis pigmentosa. Moreover, the use of flavonoids as allosteric modulators may also be applicable to other members of the G protein-coupled receptors superfamily.

New insights into the stereochemical requirements of the bradykinin B1 receptor antagonists binding.

Bradykinin (BK) is a nonapeptide involved in several pathophysiological conditions including among others, septic and haemorrhagic shock, anaphylaxis, arthritis, rhinitis, asthma, inflammatory bowel disease. Accordingly, BK antagonists have long been sought after for therapeutic intervention. Action of BK is mediated through two different G-protein coupled receptors known as B1 and B2. Although there are several B1 antagonists reported in literature, their pharmacological profile is not yet optimal so that new molecules need to be discovered. In the present work we have constructed an atomistic model of the B1 receptor and docked diverse available non-peptide antagonists in order to get a deeper insight into the structure-activity relationships involving binding to this receptor. The model was constructed by homology modeling using the chemokine CXC4 and bovine rhodopsin receptors as template. The model was further refined using molecular dynamics for 600ns with the protein embedded in a POPC bilayer. From the refinement process we obtained an average structure that was used for docking studies using the Glide software. Antagonists selected for the docking studies include Compound 11, Compound 12, Chroman28, SSR240612, NPV-SAA164 and PS020990. The results of the docking study underline the role of specific receptor residues in ligand binding. The results of this study permitted to define a pharmacophore that describes the stereochemical requirements of antagonist binding, and can be used for the discovery of new compounds.

New insights into the stereochemical requirements of the bradykinin B2 receptor antagonists binding.

Bradykinin (BK) is a member of the kinin family, released in response to inflammation, trauma, burns, shock, allergy and some cardiovascular diseases, provoking vasodilatation and increased vascular permeability among other effects. Their actions are mediated through at least two G-protein coupled receptors, B1 a receptor up-regulated during inflammation episodes or tissue trauma and B2 that is constitutively expressed in a variety of cell types. The goal of the present work is to carry out a structure-activity study of BK B2 antagonism, taking into account the stereochemical features of diverse non-peptide antagonists and the way these features translate into ligand anchoring points to complementary regions of the receptor, through the analysis of the respective ligand-receptor complex. For this purpose an atomistic model of the BK B2 receptor was built by homology modeling and subsequently refined embedded in a lipid bilayer by means of a 600 ns molecular dynamics trajectory. The average structure from the last hundred nanoseconds of the molecular dynamics trajectory was energy minimized and used as model of the receptor for docking studies. For this purpose, a set of compounds with antagonistic profile, covering maximal diversity were selected from the literature. Specifically, the set of compounds include Fasitibant, FR173657, Anatibant, WIN64338, Bradyzide, CHEMBL442294, and JSM10292. Molecules were docked into the BK B2 receptor model and the corresponding complexes analyzed to understand ligand-receptor interactions. The outcome of this study is summarized in a 3D pharmacophore that explains the observed structure-activity results and provides insight into the design of novel molecules with antagonistic profile. To prove the validity of the pharmacophore hypothesized a virtual screening process was also carried out. The pharmacophore was used as query to identify new hits using diverse databases of molecules. The results of this study revealed a set of new hits with structures not connected to the molecules used for pharmacophore development. A few of these structures were purchased and tested. The results of the binding studies show about a 33% success rate with a correlation between the number of pharmacophore points fulfilled and their antagonistic potency. Some of these structures are disclosed in the present work.

Zinc Is Involved in Depression by Modulating G Protein-Coupled Receptor Heterodimerization.

5-Hydroxytryptamine 1A receptor and galanin receptor 1 belong to the G protein-coupled receptors superfamily, and they have been described to heterodimerize triggering an anomalous physiological state that would underlie depression. Zinc supplementation has been widely reported to improve treatment against major depressive disorder. Our work has focused on the study and characterization of these receptors and its relationships with zinc both under purified conditions and in cell culture. To this aim, we have designed a strategy to purify the receptors in a conformationally active state. We have used receptors tagged with the monoclonal Rho-1D4 antibody and employed ligand-assisted purification in order to successfully purify both receptors in a properly folded and active state. The interaction between both purified receptors has been analyzed by surface plasmon resonance in order to determine the kinetics of dimerization. Zinc effect on heteromer has also been tested using the same methodology but exposing the 5-hydroxytryptamine 1A receptor to zinc before the binding experiment. These results, combined with Förster resonance energy transfer (FRET) measurements, in the absence and presence of zinc, suggest that this ion is capable of disrupting this interaction. Moreover, molecular modeling suggests that there is a coincidence between zinc-binding sites and heterodimerization interfaces for the serotonin receptor. Our results establish a rational explanation for the role of zinc in the molecular processes associated with receptor-receptor interactions and its relationship with depression, in agreement with previously reported evidence for the positive effects of zinc in depression treatment, and the involvement of our target dimer in the same disease.

Designing type II topoisomerase inhibitors: A molecular modeling approach.

Nosocomial infections are produced by pathogens with the ability to persist in hospital environments and with the propensity to develop resistance to diverse antimicrobials. In order to tackle resistance, it has been pointed as good strategy to select resilient drug targets that are evolutionally constrained to design drugs less susceptible to develop resistance. Molecular modeling can help to fulfill this goal by providing a rationalization of the observed resistance at the molecular level and, suggesting modifications on existing drugs or in the design of new ones to overcome the problem. The present report focus on type II topoisomerases, a clinical validated target for antibacterials and describe diverse modes of intervention including, inhibition of their ATPase function, stabilization of the cleavage complex or prevention of DNA strand hydrolysis. Moreover, the origin of resistance is also rationalized on the base of ligand-target interactions. Finally, efforts are described to circumvent the effect of non-susceptible strains by the design of new drugs based on existing ones, like the case of diones that act through the same mechanism as quinolones or the newly released quinole-carbonitrile derivatives that inhibit type II topoisomerases through a new mechanism.

Molecular determinants of the bacterial resistance to fluoroquinolones: a computational study.

Quinolones constitute a large class of antibacterial agents whose action is mediated through the formation of a ternary complex with DNA and either, DNA Gyrase or topoisomerase IV, resulting in the inhibition of DNA replication. In order to get a deeper insight into the features of the complex formation, we carried out docking studies of fifteen diverse quinolones to the cleaved topoisomerase IV-DNA complex. Docking studies were performed using the crystal structures of the cleaved complex with levofloxacin and moxifloxacin (pdb entries 3K9F and 2XKK, respectively) using the GOLD software. Ligands dock in positions similar to those of the crystal structures. Analysis of the results reveals that bound quinolones appear intercalated between the two nucleotides that are involved in the DNA cleavage and exhibit hydrogen bonds with Arg(117) and, the latter mediated though a water molecule. Arg(117) has not been described to be involved in resistance, since it is putatively involved in the enzymatic reaction and its mutation would be lethal for the organism. Mutants of Ser(79) exhibit resistance to quinolones which can be explained by the loss of an important anchoring point. Interestingly, quinolone resistance observed in Asp(83) mutants cannot be explained directly on the basis of the loss of a direct interaction, but could be explained on the basis of its involvement at the entrance of the ligands to their binding pocket since the residue is located at the mouth of the pocket. The results of the present study suggest that the 4-keto and 3-carboxyl groups of the fluoroquinolones bind a Mg(2+) before binding to the cleaved topoisomarase IV-DNA complex and use Asp(83) for entry into the binding pocket. Accordingly, mutations that do not conserve the binding capacity for the quinolone-Mg(2+) complex will prevent the binding of this class of ligands. The results we present here are also compared with the structure of PD0305970 a 2,4-dione active against the Ser(79) and Asp(83) mutants.