RESUMO

The correct evaluation of ligand binding free energies by computational methods is still a very challenging active area of research. The most employed methods for these calculations can be roughly classified into four groups: (i) the fastest and less accurate methods, such as molecular docking, designed to sample a large number of molecules and rapidly rank them according to the potential binding energy; (ii) the second class of methods use a thermodynamic ensemble, typically generated by molecular dynamics, to analyze the endpoints of the thermodynamic cycle for binding and extract differences, in the so-called 'end-point' methods; (iii) the third class of methods is based on the Zwanzig relationship and computes the free energy difference after a chemical change of the system (alchemical methods); and (iv) methods based on biased simulations, such as metadynamics, for example. These methods require increased computational power and as expected, result in increased accuracy for the determination of the strength of binding. Here, we describe an intermediate approach, based on the Monte Carlo Recursion (MCR) method first developed by Harold Scheraga. In this method, the system is sampled at increasing effective temperatures, and the free energy of the system is assessed from a series of terms W(b,T), computed from Monte Carlo (MC) averages at each iteration. We show the application of the MCR for ligand binding with datasets of guest-hosts systems (N = 75) and we observed that a good correlation is obtained between experimental data and the binding energies computed with MCR. We also compared the experimental data with an end-point calculation from equilibrium Monte Carlo calculations that allowed us to conclude that the lower-energy (lower-temperature) terms in the calculation are the most relevant to the estimation of the binding energies, resulting in similar correlations between MCR and MC data and the experimental values. On the other hand, the MCR method provides a reasonable view of the binding energy funnel, with possible connections with the ligand binding kinetics, as well. The codes developed for this analysis are publicly available on GitHub as a part of the LiBELa/MCLiBELa project (https://github.com/alessandronascimento/LiBELa).

RESUMO

Septins possess a conserved guanine nucleotide-binding (G) domain that participates in the stabilization of organized hetero-oligomeric complexes which assemble into filaments, rings and network-like structures. The fruit fly, Drosophila melanogaster, has five such septin genes encoding Sep1, Sep2, Sep4, Sep5 and Pnut. Here, we report the crystal structure of the heterodimer formed between the G-domains of Sep1 and Sep2, the first from an insect to be described to date. A G-interface stabilizes the dimer (in agreement with the expected arrangement for the Drosophila hexameric particle) and this bears significant resemblance to its human counterparts, even down to the level of individual amino acid interactions. On the other hand, a model for the G-interface formed between the two copies of Pnut which occupy the centre of the hexamer, shows important structural differences, including the loss of a highly favourable bifurcated salt-bridge network. Whereas wild-type Pnut purifies as a monomer, the reintroduction of the salt-bridge network results in stabilizing the dimeric interface in solution as shown by size exclusion chromatography and thermal stability measurements. Adaptive steered molecular dynamics reveals an unzipping mechanism for dimer dissociation which initiates at a point of electrostatic repulsion within the switch II region. Overall, the data contribute to a better understanding of the molecular interactions involved in septin assembly/disassembly.

RESUMO

Structural biology projects are highly dependent on the large-scale expression of soluble protein and, for this purpose, heterologous expression using bacteria or yeast as host systems is usually employed. In this scenario, some of the parameters to be optimized include (i) those related to the protein construct, such as the use of a fusion protein, the choice of an N-terminus fusion/tag or a C-terminus fusion/tag; (ii) those related to the expression stage, such as the concentration and selection of inducer agent and temperature expression and (iii) the choice of the host system, which includes the selection of a prokaryotic or eukaryotic cell and the adoption of a strain. The optimization of some of the parameters related to protein expression, stage (ii), is straightforward. On the other hand, the determination of the most suitable parameters related to protein construction requires a new cycle of gene cloning, while the optimization of the host cell is less straightforward. Here, we evaluated a scalable approach for the screening of host cells for protein expression in a structural biology pipeline. We evaluated four Escherichia coli strains looking for the best yield of soluble heterologous protein expression using the same strategy for protein construction and gene cloning and comparing it to our standard strain, Rosetta 2 (DE3). Using a liquid handling device (robot), E. coli pT-GroE, Lemo21(DE3), Arctic Express (DE3), and Rosetta Gami 2 (DE3) strains were screened for the maximal yield of soluble heterologous protein recovery. For the genes used in this experiment, the Arctic Express (DE3) strain resulted in better yields of soluble heterologous proteins. We propose that screening of host cell/strain is feasible, even for smaller laboratories and the experiment as proposed can easily be scalable to a high-throughput approach.

Assuntos

RESUMO

Studies have shown that the level of ascorbic acid (AA) is reduced in the brain of Alzheimer's disease (AD) patients. However, its effect on amyloid-ß 1-42 (Aß42) aggregation has not yet been elucidated. Here we investigated for the first time the effect of AA on Aß42 aggregation using fluorescence assay, circular dichroism, atomic force microscopy, isothermal titration calorimetry, ligand docking, and molecular dynamics. Our results showed that the fibril content decreases in the growth phase when the peptides are co-incubated with AA. AA molecules bind to Aß42 peptides with high binding affinity and a binding site for AA between the ß-strands of Aß42 oligomers prevents the stack of adjacent strands. We demonstrate the inhibitory effect of AA on the aggregation of Aß42 and its molecular interactions, which can contribute to the development of an accessible therapy for AD and also to the design of novel drugs for other amyloidogenic diseases.

Assuntos

RESUMO

Infectious diseases account for 25% of the causes of death worldwide and this rate is expected to increase due to antibiotic resistance. Among the bacteria associated with healthcare infections, Staphylococcus aureus is a prevalent pathogen and about 50% of the isolates are found to be methicillin-resistant. Here we describe the identification of ticarcillin as a weak binder of the S. aureus UDP-N-acetylglucosamine 2-epimerase. After a docking screening, ticarcillin was identified as a ligand in using the recently proposed isothermal analysis of differential scanning fluorimetry data. Finally, an equilibrium MD simulation confirmed the docking binding mode as a stable pose, with large contributions to the binding energy coming from interactions between Arg206 and Arg207 and the carboxylate groups in ticarcillin.

Assuntos

RESUMO

Malaria is still today one of the most concerning diseases, with 219 million infections in 2019, most of them in Sub-Saharan Africa and Latin America, causing approx. 409,000 deaths per year. Despite the tremendous advances in malaria treatment and prevention, there is still no vaccine for this disease yet available and the increasing parasite resistance to already existing drugs is becoming an alarming issue globally. In this context, several potential targets for the development of new drug candidates have been proposed and, among those, the de novo biosynthesis pathway for the B6 vitamin was identified to be a promising candidate. The reason behind its significance is the absence of the pathway in humans and its essential presence in the metabolism of major pathogenic organisms. The pathway consists of two enzymes i.e. Pdx1 (PLP synthase domain) and Pdx2 (glutaminase domain), the last constituting a transient and dynamic complex with Pdx1 as the prime player and harboring the catalytic center. In this review, we discuss the structural biology of Pdx1 and Pdx2, together with and the understanding of the PLP biosynthesis provided by the crystallographic data. We also highlight the existing evidence of the effect of PLP synthesis inhibition on parasite proliferation. The existing data provide a flourishing environment for the structure-based design and optimization of new substrate analogs that could serve as inhibitors or even suicide inhibitors.

Assuntos

RESUMO

Diabetes is an important chronic disease affecting about 10% of the adult population in the US and over 420 million people worldwide, resulting in 1.6 million deaths every year, according to the World Health Organization. The most common type of the disease, type 2 diabetes, can be pharmacologically managed using oral hypoglycemic agents or thiazolidinediones (TZDs), such as pioglitazone, which act by activating the Peroxisome Proliferated-Activated Receptor γ. Despite their beneficial effects in diabetes treatment, TZDs like rosiglitazone and troglitazone were withdrawn due to safety reasons, creating a void in the pharmacological options for the treatment of this important disease. Here, we explored a structure-based approach in the screening for new chemical probes for a deeper investigation of the effects of PPARγ activation. A class of tetrazole compounds was identified and the compounds named T1, T2 and T3 were purchased and evaluated for their ability to interact with the PPARγ ligand binding domain (LBD). The compounds were binders with micromolar range affinity, as determined by their IC50 values. A Monte Carlo simulation of the compound T2 revealed that the tetrazole ring makes favorable interaction with the polar arm of the receptor binding pocket. Finally, the crystal structure of the PPARγ-LBD-T2 complex was solved at 2.3 Å, confirming the binding mode for this compound. The structure also revealed that, when the helix H12 is mispositioned, an alternative binding conformation is observed for the ligand suggesting an H12-dependent binding conformation for the tetrazole compound.

Assuntos

RESUMO

Glycoside hydrolases (GHs) are essential for plant biomass deconstruction. GH11 family consist of endo-ß-1,4-xylanases which hydrolyze xylan, the second most abundant cell wall biopolymer after cellulose, into small bioavailable oligomers. Structural requirements for enzymatic mechanism of xylan hydrolysis is well described for GH11 members. However, over the last years, it has been discovered that some enzymes from GH11 family have a secondary binding sites (SBS), which modulate the enzymes activities, but mechanistic details of the molecular communication between the active site and SBS of the enzymes remain a conundrum. In the present work we structurally characterized GH11 xylanase from Paenibacillus xylanivorans A57 (PxXyn11B), a microorganism of agricultural importance, using protein crystallography and molecular dynamics simulations. The PxXyn11B structure was solved to 2.5 Å resolution and different substrates (xylo-oligosaccharides from X3 to X6), were modelled in its active and SBS sites. Molecular Dynamics (MD) simulations revealed an important role of SBS in the activity and conformational mobility of PxXyn11B, demonstrating that binding of the reaction products to the SBS of the enzyme stabilizes the N-terminal region and, consequently, the active site. Furthermore, MD simulations showed that the longer the ligand, the better is the stabilization within active site, and the positive subsites contribute less to the stabilization of the substrates than the negative ones. These findings provide rationale for the observed enzyme kinetics, shedding light on the conformational modulation of the GH11 enzymes via their SBS mediated by the positive molecular feedback loop which involve the products of the enzymatic reaction.

RESUMO

Antibiotic resistance is a worldwide concern that requires a concerted action from physicians, patients, governmental agencies, and academia to prevent infections and the spread of resistance, track resistant bacteria, improve the use of current antibiotics, and develop new antibiotics. Despite the efforts spent so far, the current antibiotics in the market are restricted to only five general targets/pathways highlighting the need for basic research focusing on the discovery and evaluation of new potential targets. Here we interrogate two biosynthetic pathways as potentially druggable pathways in bacteria. The biosynthesis pathway for thiamine (vitamin B1), absent in humans, but found in many bacteria, including organisms in the group of the ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumanii, Pseudomonas aeruginosa, and Enterobacter sp.) and the biosynthesis pathway for pyridoxal 5'-phosphate and its vitamers (vitamin B6), found in S. aureus. Using current genomic data, we discuss the possibilities of inhibition of enzymes in the pathway and review the current state of the art in the scientific literature.

Assuntos

RESUMO

The precise modeling of molecular interactions remains an important goal among molecular modeling techniques. Some of the challenges in the field include the precise definition of a Hamiltonian for biomolecular systems, together with precise parameters derived from Molecular Mechanics Force Fields, for example. The problem is even more challenging when interaction energies from different species are computed, such as the interaction energy involving a ligand and a protein, given that small differences must be computed from large energies. Here we evaluated the effects of the electrostatic model for ligand binding energy evaluation in the context of ligand docking. For this purpose, a classical Coulomb potential with distance-dependent dielectrics was compared with a Poisson-Boltzmann (PB) model for electrostatic potential computation, based on DelPhi calculations. We found that, although the electrostatic energies were highly correlated for the Coulomb and PB models, the ligand pose and the enrichment of actual ligands against decoy compounds, were improved when binding energies were computed using PB as compared to the Coulomb model. We observed that the electrostatic energies computed with the Coulomb model were, on average, ten times larger than the energies computed with the PB model, suggesting a strong overestimation of the polar interactions in the Coulomb model. We also found that a slightly smoothed Lennard-Jones potential combined with the PB model resulted in a good compromise between ligand sampling and energetic scoring.