Unraveling the conformational dynamics of glycerol 3-phosphate dehydrogenase, a nicotinamide adenine dinucleotide-dependent enzyme of Leishmania mexicana.

Allosteric changes modulate the enzymatic activity, leading to activation or inhibition of the molecular target. Understanding the induced fit accommodation mechanism of a ligand in its lowest-free energy state and the subsequent conformational changes induced in the protein are important questions for drug design. In the present study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were applied to analyze the glycerol-3-phosphate dehydrogenase of Leishmania mexicana (LmGPDH) conformational changes induced by its cofactor and substrate binding. GPDH is a nicotinamide adenine dinucleotide (NAD)-dependent enzyme, which has been reported as an interesting target for drug discovery and development against leishmaniasis. Despite its relevance for glycolysis and pentose phosphate pathways, the structural flexibility and conformational motions of LmGPDH in complex with NADH and dihydroxyacetone phosphate (DHAP) remain unexplored. Here, we analyzed the conformational dynamics of the enzyme-NADH complex (cofactor), and the enzyme-NADH-DHAP complex (adduct), mapped the hydrogen-bond interactions for the complexes and pointed some structural determinants of the enzyme that emerge from these contacts to NADH and DHAP. Finally, we proposed a consistent mechanism for the conformational changes on the first step of the reversible redox conversion of dihydroxyacetone phosphate to glycerol 3-phosphate, indicating key residues and interactions that could be further explored in drug discovery.

Computational Investigation of Bisphosphate Inhibitors of 3-Deoxy-d-manno-octulosonate 8-phosphate Synthase.

The synthase, 3-deoxy-d-manno-octulosonate 8-phosphate (KDO8P), is a key enzyme for the lipopolysaccharide (LPS) biosynthesis of gram-negative bacteria and a potential target for developing new antimicrobial agents. In this study, computational molecular modeling methods were used to determine the complete structure of the KDO8P synthase from Neisseria meningitidis and to investigate the molecular mechanism of its inhibition by three bisphosphate inhibitors: BPH1, BPH2, and BPH3. Our results showed that BPH1 presented a protein-ligand complex with the highest affinity, which is in agreement with experimental data. Furthermore, molecular dynamics (MD) simulations showed that BPH1 is more active due to the many effective interactions, most of which are derived from its phosphoenolpyruvate moiety. Conversely, BPH2 exhibited few hydrogen interactions during the MD simulations with key residues located at the active sites of the KDO8P synthase. In addition, we hydroxylated BPH2 to create the hypothetical molecule named BPH3, to investigate the influence of the hydroxyl groups on the affinity of the bisphosphate inhibitors toward the KDO8P synthase. Overall, we discuss the main interactions between the KDO8P synthase and the bisphosphate inhibitors that are potential starting points for the design of new molecules with significant antibiotic activities.

Experimental study and computational modelling of cruzain cysteine protease inhibition by dipeptidyl nitriles.

Chagas disease affects millions of people in Latin America. This disease is caused by the protozoan parasite Trypanossoma cruzi. The cysteine protease cruzain is a key enzyme for the survival and propagation of this parasite lifecycle. Nitrile-based inhibitors are efficient inhibitors of cruzain that bind by forming a covalent bond with this enzyme. Here, three nitrile-based inhibitors dubbed Neq0409, Neq0410 and Neq0570 were synthesized, and the thermodynamic profile of the bimolecular interaction with cruzain was determined using isothermal titration calorimetry (ITC). The result suggests the inhibition process is enthalpy driven, with a detrimental contribution of entropy. In addition, we have used hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) and Molecular Dynamics (MD) simulations to investigate the reaction mechanism of reversible covalent modification of cruzain by Neq0409, Neq0410 and Neq0570. The computed free energy profile shows that the nucleophilic attack of Cys25 on the carbon C1 of inhibitiors and the proton transfer from His162 to N1 of the dipeptidyl nitrile inhibitor take place in a single step. The calculated free energy of the inhibiton reaction is in agreement with covalent experimental binding. Altogether, the results reported here suggests that nitrile-based inhibitors are good candidates for the development of reversible covalent inhibitors of cruzain and other cysteine proteases.

Structural analysis of viral infectivity factor of HIV type 1 and its interaction with A3G, EloC and EloB.

BACKGROUND: The virion infectivity factor (Vif) is an accessory protein, which is essential for HIV replication in host cells. Vif neutralizes the antiviral host protein APOBEC3 through recruitment of the E3 ubiquitin ligase complex. METHODOLOGY: Fifty thousand Vif models were generated using the ab initio relax protocol of the Rosetta algorithm from sets of three- and nine-residue fragments using the fragment Monte Carlo insertion-simulated annealing strategy, which favors protein-like features, followed by an all-atom refinement. In the protocol, a constraints archive was used to define the spatial relationship between the side chains from Cys/His residues and zinc ions that formed the zinc-finger motif that is essential for Vif function. We also performed centroids analysis and structural analysis with respect to the formation of the zinc-finger, and the residue disposal in the protein binding domains. Additionally, molecular docking was used to explore details of Vif-A3G and Vif-EloBC interactions. Furthermore, molecular dynamics simulation was used to evaluate the stability of the complexes Vif-EloBC-A3G and Vif-EloC. PRINCIPAL FINDINGS: The zinc in the HCCH domain significantly alters the folding of Vif and changes the structural dynamics of the HCCH region. Ab initio modeling indicated that the Vif zinc-finger possibly displays tetrahedral geometry as suggested by Mehle et al. (2006). Our model also showed that the residues L146 and L149 of the BC-box motif bind to EloC by hydrophobic interactions, and the residue P162 of the PPLP motif is important to EloB binding. CONCLUSIONS/SIGNIFICANCE: The model presented here is the first complete three-dimensional structure of the Vif. The interaction of Vif with the A3G protein and the EloBC complex is in agreement with empirical data that is currently available in the literature and could therefore provide valuable structural information for advances in rational drug design.