Lucas A. Defelipe

Preclinical development of novel Rac1-GEF signaling inhibitors using a rational design approach in highly aggressive breast cancer cell lines.

Rho GTPases play a key role in the regulation of multiple essential cellular processes, including actin dynamics, gene transcription and cell cycle progression. Aberrant activation of Rac1, a member of Rho family of small GTPases, is associated with tumorigenesis, cancer progression, invasion and metastasis. Particularly, Rac1 is overexpressed and hyperactivated in highly aggressive breast cancer. Thus, Rac1 appears to be a promising and relevant target for the development of novel anticancer drugs. We identified the novel Rac1 inhibitor ZINC69391 through a docking-based virtual library screening targeting Rac1 activation by GEFs. This compound was able to block Rac1 interaction with its GEF Tiam1, prevented EGF-induced Rac1 activation and inhibited cell proliferation, cell migration and cell cycle progression in highly aggressive breast cancer cell lines. Moreover, ZINC69391 showed an in vivo antimetastatic effect in a syngeneic animal model. We further developed the novel analog 1A-116 by rational design and showed to be specific and more potent than the parental compound in vitro and interfered Rac1-P-Rex1 interaction. We also showed an enhanced in vivo potency of 1A-116 analog. These results show that we have developed novel Rac1 inhibitors that may be used as a novel anticancer therapy.

Juvenile Hormone Synthesis: “esterify then epoxidize” or “epoxidize then esterify”? Insights from the Structural Characterization of Juvenile Hormone Acid Methyltransferase

Juvenile hormones (JHs) play key roles in regulating metamorphosis and reproduction in insects. The last two steps of JH synthesis diverge depending on the insect order. In Lepidoptera, epoxidation by a P450 monooxygenase precedes esterification by a juvenile hormone acid methyltransferase (JHAMT). In Orthoptera, Dictyoptera, Coleoptera and Diptera epoxidation follows methylation. The aim of our study was to gain insight into the structural basis of JHAMTs substrate recognition as a means to understand the divergence of these pathways. Homology modeling was used to build the structure of Aedes aegypti JHAMT. The substrate binding site was identified, as well as the residues that interact with the methyl donor (S-adenosylmethionine) and the carboxylic acid of the substrate methyl acceptors, farnesoic acid (FA) and juvenile hormone acid (JHA). To gain further insight we generated the structures of Anopheles gambiae, Bombyx mori, Drosophila melanogaster and Tribolium castaneum JHAMTs. The modeling results were compared with previous experimental studies using recombinant proteins, whole insects, corpora allata or tissue extracts. The computational study helps explain the selectivity toward the (10R)-JHA isomer and the reduced activity for palmitic and lauric acids. The analysis of our results supports the hypothesis that all insect JHAMTs are able to recognize both FA and JHA as substrates. Therefore, the order of the methylation/epoxidation reactions may be primarily imposed by the epoxidases substrate specificity. In Lepidoptera, epoxidase might have higher affinity than JHAMT for FA, so epoxidation precedes methylation, while in most other insects there is no epoxidation of FA, but esterification of FA to form MF, followed by epoxidation to JH III.

Protein Topology Determines Cysteine Oxidation Fate: The Case of Sulfenyl Amide Formation among Protein Families

Cysteine residues have a rich chemistry and play a critical role in the catalytic activity of a plethora of enzymes. However, cysteines are susceptible to oxidation by Reactive Oxygen and Nitrogen Species, leading to a loss of their catalytic function. Therefore, cysteine oxidation is emerging as a relevant physiological regulatory mechanism. Formation of a cyclic sulfenyl amide residue at the active site of redox-regulated proteins has been proposed as a protection mechanism against irreversible oxidation as the sulfenyl amide intermediate has been identified in several proteins. However, how and why only some specific cysteine residues in particular proteins react to form this intermediate is still unknown. In the present work using in-silico based tools, we have identified a constrained conformation that accelerates sulfenyl amide formation. By means of combined MD and QM/MM calculation we show that this conformation positions the NH backbone towards the sulfenic acid and promotes the reaction to yield the sulfenyl amide intermediate, in one step with the concomitant release of a water molecule. Moreover, in a large subset of the proteins we found a conserved beta sheet-loop-helix motif, which is present across different protein folds, that is key for sulfenyl amide production as it promotes the previous formation of sulfenic acid. For catalytic activity, in several cases, proteins need the Cysteine to be in the cysteinate form, i.e. a low pKa Cys. We found that the conserved motif stabilizes the cysteinate by hydrogen bonding to several NH backbone moieties. As cysteinate is also more reactive toward ROS we propose that the sheet-loop-helix motif and the constraint conformation have been selected by evolution for proteins that need a reactive Cys protected from irreversible oxidation. Our results also highlight how fold conservation can be correlated to redox chemistry regulation of protein function.

TuberQ: a Mycobacterium tuberculosis protein druggability database

In 2012 an estimated 8.6 million people developed tuberculosis (TB) and 1.3 million died from the disease [including 320 000 deaths among human immunodeficiency virus (HIV)-positive people]. There is an urgent need for new anti-TB drugs owing to the following: the fact that current treatments have severe side effects, the increasing emergence of multidrug-resistant strains of Mycobacterium tuberculosis (Mtb), the negative drug–drug interactions with certain HIV (or other disease) treatments and the ineffectiveness against dormant Mtb. In this context we present here the TuberQ database, a novel resource for all researchers working in the field of drug development in TB. The main feature of TuberQ is to provide a druggability analysis of Mtb proteins in a consistent and effective manner, contributing to a better selection of potential drug targets for screening campaigns and the analysis of targets for structure-based drug design projects. The structural druggability analysis is combined with features related to the characteristics of putative inhibitor binding pockets and with functional and biological data of proteins. The structural analysis is performed on all available unique Mtb structures and high-quality structural homology-based models. This information is shown in an interactive manner, depicting the protein structure, the pockets and the associated characteristics for each protein. TuberQ also provides information about gene essentiality information, as determined from whole cell–based knockout experiments, and expression information obtained from microarray experiments done in different stress-related conditions. We hope that TuberQ will be a powerful tool for researchers working in TB and eventually will lead to the identification of novel putative targets and progresses in therapeutic activities. Database URL: http://tuberq.proteinq.com.ar/

Thiol redox biochemistry: insights from computer simulations

Thiol redox chemical reactions play a key role in a variety of physiological processes, mainly due to the presence of low-molecular-weight thiols and cysteine residues in proteins involved in catalysis and regulation. Specifically, the subtle sensitivity of thiol reactivity to the environment makes the use of simulation techniques extremely valuable for obtaining microscopic insights. In this work we review the application of classical and quantum–mechanical atomistic simulation tools to the investigation of selected relevant issues in thiol redox biochemistry, such as investigations on (1) the protonation state of cysteine in protein, (2) two-electron oxidation of thiols by hydroperoxides, chloramines, and hypochlorous acid, (3) mechanistic and kinetics aspects of the de novo formation of disulfide bonds and thiol−disulfide exchange, (4) formation of sulfenamides, (5) formation of nitrosothiols and transnitrosation reactions, and (6) one-electron oxidation pathways.

QM/MM study of the C—C coupling reaction mechanism of CYP121, an essential Cytochrome p450 of Mycobacterium tuberculosis

Among 20 p450s of Mycobacterium tuberculosis (Mt), CYP121 has received an outstanding interest, not only due to its essentiality for bacterial viability but also because it catalyzes an unusual carbon–carbon coupling reaction. Based on the structure of the substrate bound enzyme, several reaction mechanisms were proposed involving first Tyr radical formation, second Tyr radical formation, and C—C coupling. Key and unknown features, being the nature of the species that generate the first and second radicals, and the role played by the protein scaffold each step. In the present work we have used classical and quantum based computer simulation methods to study in detail its reaction mechanism. Our results show that substrate binding promotes formation of the initial oxy complex, Compound I is the responsible for first Tyr radical formation, and that the second Tyr radical is formed subsequently, through a PCET reaction, promoted by the presence of key residue Arg386. The final C—C coupling reaction possibly occurs in bulk solution, thus yielding the product in one oxygen reduction cycle. Our results thus contribute to a better comprehension of MtCYP121 reaction mechanism, with direct implications for inhibitor design, and also contribute to our general understanding of these type of enzymes. Proteins 2014; 82:1004–1021.

Molecular Dynamics Simulations Provide Atomistic Insight into Hydrogen Exchange Mass Spectrometry Experiments

It is now clear that proteins are flexible entities that in solution switch between conformations to achieve their function. Hydrogen/Deuterium Exchange Mass Spectrometry (HX/MS) is an invaluable tool to understand dynamic changes in proteins modulated by cofactor binding, post-transductional modifications, or protein-protein interactions. ERK2MAPK, a protein involved in highly conserved signal transduction pathways of paramount importance for normal cellular function, has been extensively studied by HX/MS. Experiments of the ERK2MAPK in the inactive and active states (in the presence or absence of bound ATP) have provided valuable information on the plasticity of the MAPK domain. However, interpretation of the HX/MS data is difficult, and changes are mostly explained in relation to available X-ray structures, precluding a complete atomic picture of protein dynamics. In the present work, we have used all atom Molecular Dynamics simulations (MD) to provide a theoretical framework for the interpretation of HX/MS data. Our results show that detailed analysis of protein-solvent interaction along the MD simulations allows (i) prediction of the number of protons exchanged for each peptide in the HX/MS experiments, (ii) rationalization of the experimentally observed changes in exchange rates in different protein conditions at the residue level, and (iii) that at least for ERK2MAPK, most of the functionally observed differences in protein dynamics are related to what can be considered the native state conformational ensemble. In summary, the combination of HX/MS experiments with all atom MD simulations emerges as a powerful approach to study protein native state dynamics with atomic resolution.

FUNCTIONAL ANALYSIS OF A MOSQUITO SHORT-CHAIN DEHYDROGENASE CLUSTER

The short-chain dehydrogenases (SDR) constitute one of the oldest and largest families of enzymes with over 46,000 members in sequence databases. About 25% of all known dehydrogenases belong to the SDR family. SDR enzymes have critical roles in lipid, amino acid, carbohydrate, hormone, and xenobiotic metabolism as well as in redox sensor mechanisms. This family is present in archaea, bacteria, and eukaryota, emphasizing their versatility and fundamental importance for metabolic processes. We identified a cluster of eight SDRs in the mosquito Aedes aegypti (AaSDRs). Members of the cluster differ in tissue specificity and developmental expression. Heterologous expression produced recombinant proteins that had diverse substrate specificities, but distinct from the conventional insect alcohol (ethanol) dehydrogenases. They are all NADP⁺-dependent and they have S-enantioselectivity and preference for secondary alcohols with 8-15 carbons. Homology modeling was used to build the structure of AaSDR1 and two additional cluster members. The computational study helped explain the selectivity toward the (10S)-isomers as well as the reduced activity of AaSDR4 and AaSDR9 for longer isoprenoid substrates. Similar clusters of SDRs are present in other species of insects, suggesting similar selection mechanisms causing duplication and diversification of this family of enzymes.

Molecular Dynamics in Mixed Solvents Reveals Protein–Ligand Interactions, Improves Docking, and Allows Accurate Binding Free Energy Predictions

One of the most important biological processes at the molecular level is the formation of protein-ligand complexes. Therefore, determining their structure and underlying key interactions is of paramount relevance and has direct applications in drug development. Because of its low cost relative to its experimental sibling, molecular dynamics (MD) simulations in the presence of different solvent probes mimicking specific types of interactions have been increasingly used to analyze protein binding sites and reveal protein-ligand interaction hot spots. However, a systematic comparison of different probes and their real predictive power from a quantitative and thermodynamic point of view is still missing. In the present work, we have performed MD simulations of 18 different proteins in pure water as well as water mixtures of ethanol, acetamide, acetonitrile and methylammonium acetate, leading to a total of 5.4 μs simulation time. For each system, we determined the corresponding solvent sites, defined as space regions adjacent to the protein surface where the probability of finding a probe atom is higher than that in the bulk solvent. Finally, we compared the identified solvent sites with 121 different protein-ligand complexes and used them to perform molecular docking and ligand binding free energy estimates. Our results show that combining solely water and ethanol sites allows sampling over 70% of all possible protein-ligand interactions, especially those that coincide with ligand-based pharmacophoric points. Most important, we also show how the solvent sites can be used to significantly improve ligand docking in terms of both accuracy and precision, and that accurate predictions of ligand binding free energies, along with relative ranking of ligand affinity, can be performed.

Two novel DNA variants associated with glucose-6-phosphate dehydrogenase deficiency found in Argentine pediatric patients

OBJECTIVE The enzyme glucose-6-phosphate dehydrogenase (G6PD) catalyses the first step in the pentose phosphate pathway, producing nicotinamide adenine dinucleotide phosphate (NADPH). NADPH plays a crucial role in preventing oxidative damage to proteins and other molecules in cells, mostly red blood cells. G6PD deficiency has an x-linked pattern of inheritance in which hemizygous males are deficient, while females may or may not be deficient depending on the number of affected alleles. We report two novel DNA variants in the G6PD gene detected in two male probands with chronic nonspherocytic hemolytic anemia (CNSHA), who were referred for hematological evaluation. METHOD Probands and their relatives underwent clinical, biochemical, and molecular assessment. RESULTS Two novel DNA variants, c.995C>T and c.1226C>A, were found in this study. At the protein level, they produce the substitution of Ser332Phe and Pro409Gln, respectively. These DNA variants were analyzed in the female relatives of probands for genetic counseling. CONCLUSIONS The novel DNA variants were classified as class I based on the clinical, biochemical, and molecular evaluations performed.