Ariela Vergara-Jaque

Computationally efficient methodology for atomic-level characterization of dendrimer-drug complexes: A comparison of amine- and acetyl-terminated PAMAM

PAMAM dendrimers have been widely studied as a novel means for controlled drug delivery; however, computational study of dendrimer-drug complexation is made difficult by the conformational flexibility of dendrimers and the nonspecific nature of the dendrimer-drug interactions. Conventional protocols for studying drug binding have been designed primarily for protein substrates, and, therefore, there is a need to establish new protocols to deal with the unique aspects of dendrimers. In this work, we generate cavities in generation-5 polyamidoamine (PAMAM) dendrimers at selected distances from the center of mass of the dendrimer for the insertion of the model drug: dexamethasone 21-phosphate or Dp21. The complexes are then allowed to equilibrate with distance between centers of mass of the drug and dendrimers confined to selected ranges; the free energy of complexation is estimated by the MM-GBSA (MM, molecular mechanics; GB, generalized Born; SA, surface area) method. For both amine- and modified acetyl-terminated PAMAM at both low and neutral pH, the most favorable free energy of complexation is associated with Dp21 at distance of 15-20 Å from the center of mass of the dendrimer and that smaller or larger distances yield considerably weaker affinity. In agreement with experimental results, we find acetyl-terminated PAMAM at neutral pH to form the least stable complex with Dp21. The greatest affinity is seen in the case of acetyl-terminated PAMAM at low pH, which appears to be due a complex balance of different contributions, which cannot be attributed to electrostatics, van der Waals interactions, hydrogen bonds, or charge-charge interactions alone.

Computational study on the interaction of N1 substituted pyrazole derivatives with B-raf kinase: an unusual water wire hydrogen-bond network and novel interactions at the entrance of the active site.

Docking and molecular dynamics (MD) simulations of N1 substituted pyrazole derivatives complexed with B-Raf kinase were performed to gain insight into the structural and energetic preferences of these inhibitors. First, a comparative study of fully automated docking programs AutoDock, ICM, GLIDE, and Surflex-Dock in closely approximating the X-ray crystal structure of the inhibitor (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime was performed. Afterward, the dynamics of the above-mentioned compound and the less active analogous compounds with 1-methyl-4-piperidinyl and tetrahydro-2H-pyran-4-yl groups at position N1 of pyrazole ring inside the B-Raf active site were analyzed by MD simulations. We found that the most active compound has stable interactions with residues Ile463 and His539 at the entrance of the B-Raf active site. Those interactions were in very good agreement with more reliable quantum mechanics/molecular mechanics calculations performed on the torsional angle phi between the pyrazole ring and the substituents at position N1. In addition, we identified a water wire connecting N2 of the pyrazole ring, Cys532, and Ser536, which is composed of three water molecules for the most active compound. We found some differences in the water wire hydrogen-bond network formed by less active compounds. We suggest that the differences between these structural features are responsible for the differences in activity among the studied compounds.

Mutation Update of the CLCN5 Gene Responsible for Dent Disease 1

Dent disease is a rare X‐linked tubulopathy characterized by low molecular weight proteinuria, hypercalciuria, nephrocalcinosis and/or nephrolithiasis, progressive renal failure, and variable manifestations of other proximal tubule dysfunctions. It often progresses over a few decades to chronic renal insufficiency, and therefore molecular characterization is important to allow appropriate genetic counseling. Two genetic subtypes have been described to date: Dent disease 1 is caused by mutations of the CLCN5 gene, coding for the chloride/proton exchanger ClC‐5; and Dent disease 2 by mutations of the OCRL gene, coding for the inositol polyphosphate 5‐phosphatase OCRL‐1. Herein, we review previously reported mutations (n = 192) and their associated phenotype in 377 male patients with Dent disease 1 and describe phenotype and novel (n = 42) and recurrent mutations (n = 24) in a large cohort of 117 Dent disease 1 patients belonging to 90 families. The novel missense and in‐frame mutations described were mapped onto a three‐dimensional homology model of the ClC‐5 protein. This analysis suggests that these mutations affect the dimerization process, helix stability, or transport. The phenotype of our cohort patients supports and extends the phenotype that has been reported in smaller studies.

Predicting Adsorption Affinities of Small Molecules on Carbon Nanotubes Using Molecular Dynamics Simulation

Computational techniques have the potential to accelerate the design and optimization of nanomaterials for applications such as drug delivery and contaminant removal; however, the success of such techniques requires reliable models of nanomaterial surfaces as well as accurate descriptions of their interactions with relevant solutes. In the present work, we evaluate the ability of selected models of naked and hydroxylated carbon nanotubes to predict adsorption equilibrium constants for about 30 small aromatic compounds with a variety of functional groups. The equilibrium constants determined using molecular dynamics coupled with free-energy calculation techniques are directly compared to those derived from experimental measurements. The calculations are highly predictive of the relative adsorption affinities of the compounds, with excellent correlation (r ≥ 0.9) between calculated and measured values of the logarithm of the adsorption equilibrium constant. Moreover, the agreement in absolute terms is also reasonable, with average errors of less than one decade. We also explore possible effects of surface loading, although we demonstrate that they are negligible for the experimental conditions considered. Given the degree of reliability demonstrated, we move on to employing the in silico techniques in the design of nanomaterials, using the optimization of adsorption affinity for the herbacide atrazine as an example. Our simulations suggest that, compared to other modifications of graphenic carbon, polyvinylpyrrolidone conjugation gives the highest affinity for atrazine-substantially greater than that of graphenic carbon alone-and may be useful as a nanomaterial for delivery or sequestration of atrazine.

Molecular basis of drug resistance in A/H1N1 virus

New mutants of human influenza virus (A/H1N1) exhibit resistance to antiviral drugs. The mechanism whereby they develop insensitivity to these medications is, however, not yet completely understood. A crystallographic structure of A/H1N1 neuraminidase has been published recently. Using molecular dynamic simulations, it is now possible to characterize at the atomic level the mechanism that underlies the loss of binding affinity of the drugs. In this study, free-energy perturbation was used to evaluate the relative binding free energies of Tamiflu and Relenza with H274Y, N294S, and Y252H neuraminidase mutants. Our results demonstrate a remarkable correlation between theoretical and experimental data, which quantitatively confirms that the mutants are resistant to Tamiflu but are still strongly inhibited by Relenza. The simulations further reveal the key interactions that govern the affinity of the two drugs for each mutant. This information is envisioned to prove useful for the design of novel neuraminidase inhibitors and for the characterization of new potential mutants.

Investigation of the differences in activity between hydroxycycloalkyl N1 substituted pyrazole derivatives as inhibitors of B-Raf kinase by using docking, molecular dynamics, QM/MM, and fragment-based de novo design: study of binding mode of diastereomer compounds.

N1 substituted pyrazole derivatives show diverse B-Raf kinase inhibitory activities when different hydroxy-substituted cycloalkyl groups are placed at this position. Docking, molecular dynamics (MD) simulations, and hybrid calculation methods (Quantum Mechanics/Molecular Mechanics (QM/MM)) were performed on the complexes, in order to explain these differences. Docking of the inhibitors showed the same orientation that X-ray crystal structure of the analogous (1E)-5-[1-(4-piperidinyl)-3-(4-pyridinyl)-1H-pyrazol-4-yl]-2,3-dihydro-1H-inden-1-one oxime. MD simulations of the most active diastereomer compounds containing cis- and trans-3-hydroxycyclohexyl substituents showed stable interactions with residue Ile463 at the entrance of the B-Raf active site. On the other hand, the less active diastereomer compounds containing cis- and trans-2-hydroxycyclopentyl substituents showed interactions with inner residues Asn580 and Ser465. We found that the differences in activity can be explained by considering the dynamic interactions between the inhibitors and their surrounding residues within the B-Raf binding site. We also explained the activity trend by using a testing scoring function derived from more reliable QM/MM calculations. In addition, we search for new inhibitors from a virtual screening carried out by fragment-based de novo design. We generated a set of approximately 200 virtual compounds, which interact with Ile463 and fulfill druglikeness properties according to Lipinski, Veber, and Ghose rules.

Study of the affinity between the protein kinase PKA and peptide substrates derived from kemptide using molecular dynamics simulations and MM/GBSA

We have carried out a protocol in computational biochemistry including molecular dynamics (MD) simulations and MM/GBSA free energy calculations on the complex between the protein kinase A (PKA) and the specific peptide substrate Kemptide (LRRASLG). We made the same calculations on other PKA complexes that contain Kemptide derivatives (with mutations of the arginines, and with deletions of N and C-terminal amino acids). We predicted shifts in the free energy changes from the free PKA to PKA-substrate complex (ΔΔGE→ES) when Kemptide structure is modified (we consider that the calculated shifts correlate with the experimental shifts of the free energy changes from the free PKA to the transition states (ΔΔGE→TS) determined by the catalytic efficiency (kcat/KM) changes). Our results demonstrate that it is possible to predict the kinetic properties of protein kinases using simple computational biochemistry methods. As an additional benefit, these methods give detailed molecular information that permit the analysis of the atomic forces that contribute to the affinity between protein kinases and their substrates.

In situ and in silico evaluation of amine- and folate-terminated dendrimers as nanocarriers of anesthetics

The search for new nano-systems for targeted biomedical applications and controlled drug release has attracted significant attention in polymer chemistry, pharmaceutics, and biomaterial science. Controlled drug delivery has many advantages over conventional drug administration, such as reduction of side effects, maintaining a stable plasma level concentration and improving the quality of life of patients. In this study, PAMAM G5 dendrimers and PAMAM G5-folic acid conjugates (PAMAM G5-FA) are synthesized and characterized by mass spectrometry (MALDI-MS). Controlled release studies at different pH values show that PAMAM G5-FA is a good candidate as a carrier for tramadol and morphine, while mathematical modeling is conducted, suggesting that the release process is governed by a diffusion mechanism. In addition, using molecular dynamics simulations, we investigate the structural and energetic properties that facilitate the encapsulation of tramadol and morphine by unmodified and functionalized PAMAM-G5 dendrimers at low, neutral and high pH. Our results correlate well with experimental data, confirming that tramadol and morphine may be encapsulated both by functionalized PAMAM dendrimers and unmodified PAMAM. Moreover, the simulations further reveal that hydrogen-bond and electrostatic interactions govern the affinity the dendrimers for both drugs. This information is envisioned to prove useful for the encapsulation of other drugs and for the design of novel functionalized dendrimers.

K+ Conduction and Mg2+ Blockade in a Shaker Kv-Channel Single Point Mutant with an Unusually High Conductance

Potassium channels exhibit a large diversity of single-channel conductances. Shaker is a low-conductance K-channel in which Pro475→Asp, a single-point mutation near the internal pore entrance, promotes 6- to 8-fold higher unitary current. To assess the mechanism for this higher conductance, we measured Shaker-P475D single-channel current in a wide range of symmetrical K(+) concentrations and voltages. Below 300 mM K(+), the current-to-voltage relations (i-V) showed inward rectification that disappeared at 1000 mM K(+). Single-channel conductance reached a maximum of ∼190 pS at saturating [K(+)], a value 4- to 5-fold larger than that estimated for the native channel. Intracellular Mg(2+) blocked this variant with ∼100-fold higher affinity. Near zero voltage, blockade was competitively antagonized by K(+); however, at voltages >100 mV, it was enhanced by K(+). This result is consistent with a lock-in effect in a single-file diffusion regime of Mg(2+) and K(+) along the pore. Molecular-dynamics simulations revealed higher K(+) density in the pore, especially near the Asp-475 side chains, as in the high-conductance MthK bacterial channel. The molecular dynamics also showed that K(+) ions bound distally can coexist with other K(+) or Mg(2+) in the cavity, supporting a lock-in mechanism. The maximal K(+) transport rate and higher occupancy could be due to a decrease in the electrostatic energy profile for K(+) throughout the pore, reducing the energy wells and barriers differentially by ∼0.7 and ∼2 kT, respectively.

Study of the affinity between the protein kinase PKA and homoarginine-containing peptides derived from kemptide: Free energy perturbation (FEP) calculations

Protein kinases (PKs) discriminate between closely related sequences that contain serine, threonine, and/or tyrosine residues. Such specificity is defined by the amino acid sequence surrounding the phosphorylatable residue, so that it is possible to identify an optimal recognition motif (ORM) for each PK. The ORM for the protein kinase A (PKA), a well‐known member of the PK family, is the sequence RRX(S/T)X, where arginines at the −3 and −2 positions play a key role with respect to the primed phosphorylation site. In this work, differential affinities of PKA for the peptide substrate Kemptide (LRRASLG) and mutants that substitute the arginine residues by the unnatural peptide homoarginine were evaluated through molecular dynamics (MD) and free energy perturbation (FEP) calculations. The FEP study for the homoarginine mutants required previous elaboration of a CHARMM “arginine to homoarginine” (R2B) hybrid topology file which is available in this manuscript as Supporting Information. Mutants substituting the arginine residues by alanine, lysine, and histidine were also considered in the comparison by using the same protocol. FEP calculations allowed estimating the free energy changes from the free PKA to PKA‐substrate complex (ΔΔGE→ES) when Kemptide structure was mutated. Both ΔΔGS→ES values for homoarginine mutants were predicted with a difference below 1 kcal/mol. In addition, FEP correctly predicted that all the studied mutations decrease the catalytic efficiency of Kemptide for PKA.