Development and Evaluation of Topical Zinc Oxide Nanogels Formulation Using Dendrobium anosmum and Its Effect on Acne Vulgaris.

Zinc oxide nanoparticles have high levels of biocompatibility, a low impact on environmental contamination, and suitable to be used as an ingredient for environmentally friendly skincare products. In this study, biogenically synthesized zinc oxide nanoparticles using Dendrobium anosum are used as a reducing and capping agent for topical anti-acne nanogels, and the antimicrobial effect of the nanogel is assessed on Cutibacterium acne and Staphylococcus aureus. Dendrobium anosmum leaf extract was examined for the presence of secondary metabolites and its total amount of phenolic and flavonoid content was determined. Both the biogenically and chemogenic-synthesized zinc oxide nanoparticles were compared using UV-Visible spectrophotometer, FE-SEM, XRD, and FTIR. To produce the topical nanogel, the biogenic and chemogenic zinc oxide nanoparticles were mixed with a carbomer and hydroxypropyl-methyl cellulose (HPMC) polymer. The mixtures were then tested for physical and chemical characteristics. To assess their anti-acne effectiveness, the mixtures were tested against C. acne and S. aureus. The biogenic zinc oxide nanoparticles have particle sizes of 20 nm and a high-phase purity. In comparison to chemogenic nanoparticles, the hydrogels with biogenically synthesized nanoparticles was more effective against Gram-positive bacteria. Through this study, the hybrid nanogels was proven to be effective against the microbes that cause acne and to be potentially used as a green product against skin infections.

Discrimination of Native-like States of Membrane Proteins with Implicit Membrane-based Scoring Functions.

A scoring protocol based on implicit membrane-based scoring functions and a new protocol for optimizing the positioning of proteins inside the membrane was evaluated for its capacity to discriminate native-like states from misfolded decoys. A decoy set previously established by the Baker lab (Proteins: Struct., Funct., Genet. 2006, 62, 1010-1025) was used along with a second set that was generated to cover higher resolution models. The Implicit Membrane Model 1 (IMM1), IMM1 model with CHARMM 36 parameters (IMM1-p36), generalized Born with simple switching (GBSW), and heterogeneous dielectric generalized Born versions 2 (HDGBv2) and 3 (HDGBv3) were tested along with the new HDGB van der Waals (HDGBvdW) model that adds implicit van der Waals contributions to the solvation free energy. For comparison, scores were also calculated with the distance-scaled finite ideal-gas reference (DFIRE) scoring function. Z-scores for native state discrimination, energy vs root-mean-square deviation (RMSD) correlations, and the ability to select the most native-like structures as top-scoring decoys were evaluated to assess the performance of the scoring functions. Ranking of the decoys in the Baker set that were relatively far from the native state was challenging and dominated largely by packing interactions that were captured best by DFIRE with less benefit of the implicit membrane-based models. Accounting for the membrane environment was much more important in the second decoy set where especially the HDGB-based scoring functions performed very well in ranking decoys and providing significant correlations between scores and RMSD, which shows promise for improving membrane protein structure prediction and refinement applications. The new membrane structure scoring protocol was implemented in the MEMScore web server ( http://feiglab.org/memscore ).

Doing the methylene shuffle--further insights into the inhibition of mitotic kinesin Eg5 with S-trityl L-cysteine.

S-Trityl L-cysteine (STLC) is an inhibitor of the mitotic kinesin Eg5 with potential as an antimitotic chemotherapeutic agent. We previously reported the crystal structure of the ligand-protein complex, and now for the first time, have quantified the interactions using a molecular dynamics based approach. Based on these data, we have explored the SAR of the trityl head group using the methylene shuffle strategy to expand the occupation of one of the hydrophobic pockets. The most potent compounds exhibit strong (<100 nM) inhibition of Eg5 in the basal ATPase assay and inhibit growth in a variety of tumour-derived cell lines.

How do carbon nanotubes serve as carriers for gemcitabine transport in a drug delivery system?

Aiming at understanding the molecular properties of the encapsulation of the anticancer drug gemcitabine in the single-walled carbon nanotube (SWCNT), molecular dynamics (MD) simulations were applied to the two scenarios; that of gemcitabine filling inside the SWCNT, and that of the drug in the free state. Inside the SWCNT, the cytosine ring of gemcitabine was found to form a π-π stacking conformation with the SWCNT surface, and this movement is not along the centerline of the tube from one end to the other of the tube where the distance from the center of gravity of the molecule to the surface is 4.7 Å. A tilted angle of 19° was detected between the cytosine ring of gemcitabine and the inner surface of SWCNT. In comparison to its conformation in the free form, no significant difference was observed on the torsion angle between the five- (ribose) and the six- (cytosine) membered rings. However, gemcitabine inside the SWCNT was found to have a lower number of solvating water molecules but with a stronger net solvation than the drug in the free state. This is due to the collaborative interactions between gemcitabine and the surface of the SWCNT. In addition, the steered molecular dynamics simulation (SMD) approach was employed to investigate the binding free energy for gemcitabine moving from one end to another end throughout the SWCNT. In excellent agreement with that yielded from the classical MD, the SMD energy profile confirms that the drug molecule prefers to locate inside the SWCNT.

Ranking ligand affinity for the DNA minor groove by experiment and simulation.

The structural and thermodynamic basis for the strength and selectivity of the interactions of minor groove binders (MGBs) with DNA is not fully understood. In 2003, we reported the first example of a thiazole-containing MGB that bound in a phase-shifted pattern that spanned six base pairs rather than the usual four (for tricyclic distamycin-like compounds). Since then, using DNA footprinting, NMR spectroscopy, isothermal titration calorimetry, and molecular dynamics, we have established that the flanking bases around the central four being read by the ligand have subtle effects on recognition. We have investigated the effect of these flanking sequences on binding and the reasons for the differences and established a computational method to rank ligand affinity against varying DNA sequences.

A detailed binding free energy study of 2:1 ligand-DNA complex formation by experiment and simulation.

In 2004, we used NMR to solve the structure of the minor groove binder thiazotropsin A bound in a 2:1 complex to the DNA duplex, d(CGACTAGTCG)2. In this current work, we have combined theory and experiment to confirm the binding thermodynamics of this system. Molecular dynamics simulations that use polarizable or non-polarizable force fields with single and separate trajectory approaches have been used to explore complexation at the molecular level. We have shown that the binding process invokes large conformational changes in both the receptor and ligand, which is reflected by large adaptation energies. This is compensated for by the net binding free energy, which is enthalpy driven and entropically opposed. Such a conformational change upon binding directly impacts on how the process must be simulated in order to yield accurate results. Our MM-PBSA binding calculations from snapshots obtained from MD simulations of the polarizable force field using separate trajectories yield an absolute binding free energy (-15.4 kcal mol(-1)) very close to that determined by isothermal titration calorimetry (-10.2 kcal mol(-1)). Analysis of the major energy components reveals that favorable non-bonded van der Waals and electrostatic interactions contribute predominantly to the enthalpy term, whilst the unfavorable entropy appears to be driven by stabilization of the complex and the associated loss of conformational freedom. Our results have led to a deeper understanding of the nature of side-by-side minor groove ligand binding, which has significant implications for structure-based ligand development.

Accurate prediction of protonation state as a prerequisite for reliable MM-PB(GB)SA binding free energy calculations of HIV-1 protease inhibitors.

Binding free energies were calculated for the inhibitors lopinavir, ritonavir, saquinavir, indinavir, amprenavir, and nelfinavir bound to HIV-1 protease. An MMPB/SA-type analysis was applied to conformational samples from 3 ns explicit solvent molecular dynamics simulations of the enzyme-inhibitor complexes. Binding affinities and the sampled conformations of the inhibitor and enzyme were compared between different HIV-1 protease protonation states to find the most likely protonation state of the enzyme in the complex with each of the inhibitors. The resulting set of protonation states leads to good agreement between calculated and experimental binding affinities. Results from the MMPB/SA analysis are compared with an explicit/implicit hybrid scheme and with MMGB/SA methods. It is found that the inclusion of explicit water molecules may offer a slight advantage in reproducing absolute binding free energies while the use of the Generalized Born approximation significantly affects the accuracy of the calculated binding affinities.

Structural and dynamical properties of different protonated states of mutant HIV-1 protease complexed with the saquinavir inhibitor studied by molecular dynamics simulations.

To understand the basis of drug resistance, particularly of the HIV-1 PR, three molecular dynamics (MD) simulations of HIV-1 PR mutant species, G48V, complexed with saquinavir (SQV) in explicit aqueous solution with three protonation states, diprotonation on Asp25 and Asp25' (Di-pro) and monoprotonation on each Asp residue (Mono-25 and Mono-25'). For all three states, H-bonds between saquinavir and HIV-1 PR were formed only in the two regions, flap and active site. It was found that conformation of P2 subsite of SQV in the Mono-25 state differs substantially from the other two states. The rotation about 177 degrees from the optimal structure of the wild type was observed, the hydrogen bond between P2 and the flap residue (Val48) was broken and indirect hydrogen bonds with the three residues (Asp29, Gly27, and Asp30) were found instead. In terms of complexation energies, interaction energy of -37.3 kcal/mol for the Mono-25 state is significantly lower than those of -30.7 and -10.7kcal/mol for the Mono-25' and Di-pro states, respectively. It was found also that protonation at the Asp25 leads to a better arrangement in the catalytic dyad, i.e., the Asp25-Asp25' interaction energy of -8.8 kcal/mol of the Mono-25 is significantly lower than that of -2.6kcal/mol for the Mono-25' state. The above data suggest us to conclude that interaction in the catalytic area should be used as criteria to enhance capability in drug designing and drug screening instead of using the total inhibitor/enzyme interaction.

Structure, dynamics and solvation of HIV-1 protease/saquinavir complex in aqueous solution and their contributions to drug resistance: molecular dynamic simulations.

As it is known that the understanding of the basic properties of the enzyme/inhibitor complex leads directly to enhancing the capability in drug designing and drug discovery. Molecular dynamics simulations have been performed to examine detailed information on the structure and dynamical properties of the HIV-1 PR complexed with saquinavir in the three protonated states, monoprotonates at Asp25 (Mono-25) and Asp25'(Mono-25') and diprotonate (Di-Pro) at both Asp25 and Asp25'. The obtained results support clinical data which reveal that Ile84 and Gly48 are two of the most frequent residues where mutation toward a protease inhibitor takes place. In contrast to the Ile84 mutation due to high displacement of Ile84 in the presence of saquinavir, source of the Gly48 mutation was observed to be due to the limited space in the HIV-1 PR pocket. The Gly48 was, on one side, found to form strong hydrogen bonds with saquinavir, while on the other side this residue was repelled by the hydrophobic Phe53 residue. In terms of inhibitor/enzyme binding, interactions between saquinavir and a catalytic triad of the HIV-1 PR were calculated using the ab initio method. The results show an order of the binding energy of Mono25

Insights into saquinavir resistance in the G48V HIV-1 protease: quantum calculations and molecular dynamic simulations.

The spread of acquired immune deficiency syndrome has increasingly become a great concern owing largely to the failure of chemotherapies. The G48V is considered the key signature residue mutation of HIV-1 protease developing with saquinavir therapy. Molecular dynamics simulations of the wild-type and the G48V HIV-1 protease complexed with saquinavir were carried out to explore structure and interactions of the drug resistance. The molecular dynamics results combined with the quantum-based and molecular mechanics Poisson-Boltzmann surface area calculations indicated a monoprotonation took place on D25, one of the triad active site residues. The inhibitor binding of the triad residues and its interaction energy in the mutant were similar to those in the wild-type. The overall structure of both complexes is almost identical. However, the steric conflict of the substituted valine results in the conformational change of the P2 subsite and the disruption of hydrogen bonding between the -NH of the P2 subsite and the backbone -CO of the mutated residue. The magnitude of interaction energy changes was comparable to the experimental K(i) data. The designing for a new drug should consider a reduction of steric repulsion on P2 to enhance the activity toward this mutant strain.