Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), a method to estimate interaction free energies, has been increasingly used in the study of biomolecular interactions. Recently, this method has also been applied as a scoring function in computational drug design. Here a new tool g_mmpbsa, which implements the MM-PBSA approach using subroutines written in-house or sourced from the GROMACS and APBS packages is described. g_mmpbsa was developed as part of the Open Source Drug Discovery (OSDD) consortium. Its aim is to integrate high-throughput molecular dynamics (MD) simulations with binding energy calculations. The tool provides options to select alternative atomic radii and different nonpolar solvation models including models based on the solvent accessible surface area (SASA), solvent accessible volume (SAV), and a model which contains both repulsive (SASA-SAV) and attractive components (described using a Weeks-Chandler-Andersen like integral method). We showcase the effectiveness of the tool by comparing the calculated interaction energy of 37 structurally diverse HIV-1 protease inhibitor complexes with their experimental binding free energies. The effect of varying several combinations of input parameters such as atomic radii, dielectric constant, grid resolution, solute-solvent dielectric boundary definition, and nonpolar models was investigated. g_mmpbsa can also be used to estimate the energy contribution per residue to the binding energy. It has been used to identify those residues in HIV-1 protease that are most critical for binding a range of inhibitors.

We sought to clarify on the hitherto unresolved role of N-terminal transmembrane segments (TMS) of cytochrome P450 (CYP) and its' reductase (CPR) in protein interaction/catalysis. TMS analyses show little evolutionary conservation in CYPs. The conserved CPR's TMS poses limited scope for predictable/consistent hetero-recognition with the wide bevy of CYPs' TMS, as evident from preliminary analyses and TMhit server predictions for inter-helical binding. Further, experimentations with four different CPR preparations (preps) and two liver microsomal CYPs (2C9 and 2E1) shows that the hydroxylated product formation rate is not quantitatively correlated to the extent of integrity of the CPR N-terms. Incorporation of cytochrome b (5) in some reactions afforded similar rates while employing either fully intact or partially intact CPR. A survey of literature shows that liver microsomal CYPs function quite well even without the TMS or with significantly altered TMS. These observations negate the hypothesis that N-term TMS of CPR or CYP is obligatory for CYP-CPR interaction and catalysis. Also, in CYP2E1-mediated hydroxylation of para-nitrophenol, the extent of intactness or truncation did not significantly affect the CPR preps' catalytic role at very low or high substrate concentrations. To interpret these results, we draw support from recently published research on reduced nicotinamide adenide dinucleotide phosphate oxidase (Takac et al., J Biol Chem, 286:13304-13313, 2011) and from our pertinent earlier works. We infer that CPR' free TMS segment could alter the diffusible reactive oxygen species' dynamics in the microenvironment, thereby altering the reaction outcome. Based on the evidence, we conclude that TMS merely facilitates "interaction/catalysis" by anchoring the CYP and CPR in the lipid interface.

No abstract

No abstract

Resistance development in insects significantly threatens the important benefits obtained by insecticide usage in vector control of disease-transmitting insects. Discovery of new chemical entities with insecticidal activity is highly desired in order to develop new insecticide candidates. Here, we present the design, synthesis, and biological evaluation of phenoxyacetamide-based inhibitors of the essential enzyme acetylcholinesterase 1 (AChE1). AChE1 is a validated insecticide target to control mosquito vectors of, e.g., malaria, dengue, and Zika virus infections. The inhibitors combine a mosquito versus human AChE selectivity with a high potency also for the resistance-conferring mutation G122S; two properties that have proven challenging to combine in a single compound. Structure−activity relationship analyses and molecular dynamics simulations of inhibitor−protein complexes have provided insights that elucidate the molecular basis for these properties. We also show that the inhibitors demonstrate in vivo insecticidal activity on disease-transmitting mosquitoes. Our findings support the concept of noncovalent, selective, and resistance-breaking inhibitors of AChE1 as a promising approach for future insecticide development.