A rapidly convergent simulation method: Mixed Monte Carlo/stochastic dynamics

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid-base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

High-Mannose Oligosaccharide Hemimimetics that Recapitulate the Conformation and Binding Mode to Concanavalin A, DC-SIGN and Langerin

The "carbohydrate chemical mimicry" exhibited by sp2 -iminosugars has been utilized to develop practical syntheses for analogs of the branched high-mannose-type oligosaccharides (HMOs) Man3 and Man5 . In these compounds, the terminal nonreducing Man residues have been substituted with 5,6-oxomethylidenemannonojirimycin (OMJ) motifs. The resulting oligomannoside hemimimetic accurately reproduce the structure, configuration, and conformational behavior of the original mannooligosaccharides, as confirmed by NMR and computational techniques. Binding studies with mannose binding lectins, including concanavalin A, DC-SIGN, and langerin, by enzyme-linked lectin assay and surface plasmon resonance revealed significant variations in their ability to accommodate the OMJ unit in the mannose binding site. Intriguingly, OMJMan segments demonstrated "in line" heteromultivalent effects during binding to the three lectins. Similar to the mannobiose (Man2 ) branches in HMOs, the binding modes involving the external or internal monosaccharide unit at the carbohydrate binding-domain exist in equilibrium, facilitating sliding and recapture processes. This equilibrium, which influences the multivalent binding of HMOs, can be finely modulated upon incorporation of the OMJ sp2 -iminosugar caps. As a proof of concept, the affinity and selectivity towards DC-SIGN and langerin were adjustable by presenting the OMJMan epitope in platforms with diverse architectures and valencies.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Methyl group configuration on acyclic threoninol nucleic acids (aTNAs) impacts supramolecular properties

We have synthesized acyclic allo-threoninol nucleic acids (allo-aTNAs), artificial xeno-nucleic acids (XNAs) that are diastereomers of acyclic threoninol nucleic acids (aTNAs), and have investigated their supramolecular properties. The allo-aTNAs formed homo-duplexes in an antiparallel manner but with lower thermal stability than DNA, whereas aTNAs formed extremely stable homo-duplexes. The allo-aTNAs formed duplexes with complementary aTNAs and serinol nucleic acid (SNA). The affinities of L-allo-aTNA were the highest for L-aTNA and the lowest for D-aTNA, with SNA being intermediate. The affinities of D-allo-aTNA were the reverse. Circular dichroism measurements revealed that L- and D-allo-aTNAs had weak right-handed and left-handed helicities, respectively. The weak helicity of allo-aTNAs likely explains the poor chiral discrimination of these XNAs, which is in contrast to aTNAs that have strong helical orthogonality. Energy-minimized structures of L-allo-aTNA/RNA and L-allo-aTNA/L-allo-aTNA indicated that the methyl group on the allo-aTNA strand is unfavourable for duplex formation. In contrast, the methyl group on L-aTNA likely stabilizes the duplex structure via hydrophobic effects and van der Waals interactions. Thus, the configuration of the methyl group on the XNA scaffold had an unexpectedly large impact on the hybridization ability and structure.

Synthesis, self-assembly and Langerin recognition studies of a resorcinarene-based glycocluster exposing a hyaluronic acid thiodisaccharide mimetic

Herein, we report the synthesis of an octavalent glycocluster exposing a thiodisaccharide mimetic of the repetitive unit of hyaluronic acid, βSGlcA(1 → 3)βSGlcNAc, constructed on a calix[4]resorcinarene scaffold by CuAAC reaction of suitable precursors. This glycocluster showed a strong tendency toward self-aggregation. DOSY-NMR and DLS experiments demonstrated the formation of spherical micelles of d ≅ 6.2 nm, in good agreement. TEM micrographs showed the presence of particles of different sizes, depending on the pH of the starting solution, thus evidencing that the negative charge on the micelle surface due to ionization of the GlcA residues plays an important role in the aggregation process. STD-NMR and DLS experiments provided evidence of the interaction between the synthetic glycocluster and Langerin, a relevant C-type lectin. This interaction was not observed in the STD-NMR experiments performed with the basic disaccharide, providing evidence of a multivalent effect.

Tinker 8: Software Tools for Molecular Design

The Tinker software, currently released as version 8, is a modular molecular mechanics and dynamics package written primarily in a standard, easily portable dialect of Fortran 95 with OpenMP extensions. It supports a wide variety of force fields, including polarizable models such as the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field. The package runs on Linux, macOS, and Windows systems. In addition to canonical Tinker, there are branches, Tinker-HP and Tinker-OpenMM, designed for use on message passing interface (MPI) parallel distributed memory supercomputers and state-of-the-art graphical processing units (GPUs), respectively. The Tinker suite also includes a tightly integrated Java-based graphical user interface called Force Field Explorer (FFE), which provides molecular visualization capabilities as well as the ability to launch and control Tinker calculations.

Investigating the Interaction of Cyclic RGD Peptidomimetics with αVβ₆ Integrin by Biochemical and Molecular Docking Studies

The interaction of a small library of cyclic RGD (Arg-Gly-Asp) peptidomimetics with α_Vβ₆ integrin has been investigated by means of competitive solid phase binding assays to the isolated receptor and docking calculations in the crystal structure of the α_Vβ₆ binding site. To this aim, a rigid receptor-flexible ligand docking protocol has been set up and then applied to predict the binding mode of the cyclic RGD peptidomimetics to α_Vβ₆ integrin. Although the RGD interaction with α_Vβ₆ recapitulates the RGD binding mode observed in α_Vβ₃, differences between the integrin binding pockets can strongly affect the ligand binding ability. In general, the peptidomimetics exhibited IC₅₀ values for integrin α_Vβ₆ (i.e., the concentration of compound required for 50% inhibition of biotinylated fibronectin binding to isolated α_Vβ₆ integrin) in the nanomolar range (77–345 nM), about 10–100 times higher than those for the related α_Vβ₃ receptor, with a single notable ligand displaying a low nanomolar IC₅₀ value (2.3 nM). Insights from the properties of the binding pocket combined with the analysis of the docking poses provided a rationale for ligand recognition and selectivity.

New potent αvβ3 integrin ligands based on azabicycloalkane (γ,α)-dipeptide mimics

We have designed a new synthetic strategy for the preparation of a new class of cyclic RGD integrin ligands in which the azabicycloalkane scaffold can be envisaged as a (γ,α) dipeptide mimic. The synthesis and in vitro biological evaluation of these RGD derivatives, as well as the computational study of their conformational properties and binding modes to αVβ3 integrin are described. Compound has shown to be a promising candidate as αVβ3 integrin antagonist able to interfere with both cell adhesion and movement on vitronectin with no evidence of cytotoxic effects.

Efficient molecular dynamics using geodesic integration and solvent-solute splitting

We present an approach to Langevin dynamics in the presence of holonomic constraints based on decomposition of the system into components representing geodesic flow, constrained impulse and constrained diffusion. We show that a particular ordering of the components results in an integrator that is an order of magnitude more accurate for configurational averages than existing alternatives. Moreover, by combining the geodesic integration method with a solvent-solute force splitting, we demonstrate that stepsizes of at least 8 fs can be used for solvated biomolecules with high sampling accuracy and without substantially altering diffusion rates, approximately increasing by a factor of two the efficiency of molecular dynamics sampling for such systems. The methods described in this article are easily implemented using the standard apparatus of modern simulation codes.

Protofold II: Enhanced Model and Implementation for Kinetostatic Protein Folding1

A reliable prediction of three-dimensional (3D) protein structures from sequence data remains a big challenge due to both theoretical and computational difficulties. We have previously shown that our kinetostatic compliance method (KCM) implemented into the Protofold package can overcome some of the key difficulties faced by other de novo structure prediction methods, such as the very small time steps required by the molecular dynamics (MD) approaches or the very large number of samples needed by the Monte Carlo (MC) sampling techniques. In this paper, we improve the free energy formulation used in Protofold by including the typically underrated entropic effects, imparted due to differences in hydrophobicity of the chemical groups, which dominate the folding of most water-soluble proteins. In addition to the model enhancement, we revisit the numerical implementation by redesigning the algorithms and introducing efficient data structures that reduce the expected complexity from quadratic to linear. Moreover, we develop and optimize parallel implementations of the algorithms on both central and graphics processing units (CPU/GPU) achieving speed-ups up to two orders of magnitude on the GPU. Our simulations are consistent with the general behavior observed in the folding process in aqueous solvent, confirming the effectiveness of model improvements. We report on the folding process at multiple levels, namely, the formation of secondary structural elements and tertiary interactions between secondary elements or across larger domains. We also observe significant enhancements in running times that make the folding simulation tractable for large molecules.

Impact of ribosomal modification on the binding of the antibiotic telithromycin using a combined grand canonical monte carlo/molecular dynamics simulation approach

Resistance to macrolide antibiotics is conferred by mutation of A2058 to G or methylation by Erm methyltransferases of the exocyclic N6 of A2058 (E. coli numbering) that forms the macrolide binding site in the 50S subunit of the ribosome. Ketolides such as telithromycin mitigate A2058G resistance yet remain susceptible to Erm-based resistance. Molecular details associated with macrolide resistance due to the A2058G mutation and methylation at N6 of A2058 by Erm methyltransferases were investigated using empirical force field-based simulations. To address the buried nature of the macrolide binding site, the number of waters within the pocket was allowed to fluctuate via the use of a Grand Canonical Monte Carlo (GCMC) methodology. The GCMC water insertion/deletion steps were alternated with Molecular Dynamics (MD) simulations to allow for relaxation of the entire system. From this GCMC/MD approach information on the interactions between telithromycin and the 50S ribosome was obtained. In the wild-type (WT) ribosome, the 2'-OH to A2058 N1 hydrogen bond samples short distances with a higher probability, while the effectiveness of telithromycin against the A2058G mutation is explained by a rearrangement of the hydrogen bonding pattern of the 2'-OH to 2058 that maintains the overall antibiotic-ribosome interactions. In both the WT and A2058G mutation there is significant flexibility in telithromycin's imidazole-pyridine side chain (ARM), indicating that entropic effects contribute to the binding affinity. Methylated ribosomes show lower sampling of short 2'-OH to 2058 distances and also demonstrate enhanced G2057-A2058 stacking leading to disrupted A752-U2609 Watson-Crick (WC) interactions as well as hydrogen bonding between telithromycin's ARM and U2609. This information will be of utility in the rational design of novel macrolide analogs with improved activity against methylated A2058 ribosomes.

A fragment-based approach to the SAMPL3 Challenge

The success of molecular fragment-based design depends critically on the ability to make predictions of binding poses and of affinity ranking for compounds assembled by linking fragments. The SAMPL3 Challenge provides a unique opportunity to evaluate the performance of a state-of-the-art fragment-based design methodology with respect to these requirements. In this article, we present results derived from linking fragments to predict affinity and pose in the SAMPL3 Challenge. The goal is to demonstrate how incorporating different aspects of modeling protein-ligand interactions impact the accuracy of the predictions, including protein dielectric models, charged versus neutral ligands, ΔΔGs solvation energies, and induced conformational stress. The core method is based on annealing of chemical potential in a Grand Canonical Monte Carlo (GC/MC) simulation. By imposing an initially very high chemical potential and then automatically running a sequence of simulations at successively decreasing chemical potentials, the GC/MC simulation efficiently discovers statistical distributions of bound fragment locations and orientations not found reliably without the annealing. This method accounts for configurational entropy, the role of bound water molecules, and results in a prediction of all the locations on the protein that have any affinity for the fragment. Disregarding any of these factors in affinity-rank prediction leads to significantly worse correlation with experimentally-determined free energies of binding. We relate three important conclusions from this challenge as applied to GC/MC: (1) modeling neutral ligands--regardless of the charged state in the active site--produced better affinity ranking than using charged ligands, although, in both cases, the poses were almost exactly overlaid; (2) simulating explicit water molecules in the GC/MC gave better affinity and pose predictions; and (3) applying a ΔΔGs solvation correction further improved the ranking of the neutral ligands. Using the GC/MC method under a variety of parameters in the blinded SAMPL3 Challenge provided important insights to the relevant parameters and boundaries in predicting binding affinities using simulated annealing of chemical potential calculations.

Heterobivalent ligands target cell-surface receptor combinations in vivo

A challenge in tumor targeting is to deliver payloads to cancers while sparing normal tissues. A limited number of antibodies appear to meet this challenge as therapeutics themselves or as drug-antibody conjugates. However, antibodies suffer from their large size, which can lead to unfavorable pharmacokinetics for some therapeutic payloads, and that they are targeted against only a single epitope, which can reduce their selectivity and specificity. Here, we propose an alternative targeting approach based on patterns of cell surface proteins to rationally develop small, synthetic heteromultivalent ligands (htMVLs) that target multiple receptors simultaneously. To gain insight into the multivalent ligand strategy in vivo, we have generated synthetic htMVLs that contain melanocortin (MSH) and cholecystokinin (CCK) pharmacophores that are connected via a fluorescent labeled, rationally designed synthetic linker. These ligands were tested in an experimental animal model containing tumors that expressed only one (control) or both (target) MSH and CCK receptors. After systemic injection of the htMVL in tumor-bearing mice, label was highly retained in tumors that expressed both, compared with one, target receptors. Selectivity was quantified by using ex vivo measurement of Europium-labeled htMVL, which had up to 12-fold higher specificity for dual compared with single receptor expressing cells. This proof-of-principle study provides in vivo evidence that small, rationally designed bivalent htMVLs can be used to selectively target cells that express both, compared with single complimentary cell surface targets. These data open the possibility that specific combinations of targets on tumors can be identified and selectively targeted using htMVLs.

Cyclic RGD peptidomimetics containing bifunctional diketopiperazine scaffolds as new potent integrin ligands

The synthesis of eight bifunctional diketopiperazine (DKP) scaffolds is described; these were formally derived from 2,3-diaminopropionic acid and aspartic acid (DKP-1-DKP-7) or glutamic acid (DKP-8) and feature an amine and a carboxylic acid functional group. The scaffolds differ in the configuration at the two stereocenters and the substitution at the diketopiperazinic nitrogen atoms. The bifunctional diketopiperazines were introduced into eight cyclic peptidomimetics containing the Arg-Gly-Asp (RGD) sequence. The resulting RGD peptidomimetics were screened for their ability to inhibit biotinylated vitronectin binding to the purified integrins α(v)β(3) and α(v)β(5), which are involved in tumor angiogenesis. Nanomolar IC(50) values were obtained for the RGD peptidomimetics derived from trans DKP scaffolds (DKP-2-DKP-8). Conformational studies of the cyclic RGD peptidomimetics by (1)H NMR spectroscopy experiments (VT-NMR and NOESY spectroscopy) in aqueous solution and Monte Carlo/Stochastic Dynamics (MC/SD) simulations revealed that the highest affinity ligands display well-defined preferred conformations featuring intramolecular hydrogen-bonded turn motifs and an extended arrangement of the RGD sequence [Cβ(Arg)-Cβ(Asp) average distance ≥8.8 Å]. Docking studies were performed, starting from the representative conformations obtained from the MC/SD simulations and taking as a reference model the crystal structure of the extracellular segment of integrin α(v)β(3) complexed with the cyclic pentapeptide, Cilengitide. The highest affinity ligands produced top-ranked poses conserving all the important interactions of the X-ray complex.

New aryl hydrocarbon receptor homology model targeted to improve docking reliability

The aryl hydrocarbon receptor (AhR) is a ligand-dependent, basic helix-loop-helix Per-ARNT-Sim (PAS) containing transcription factor that can bind and be activated by structurally diverse chemicals, including the toxic environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). As no experimentally determined structures of the AhR ligand binding domain (LBD) are available and previous homology models were only derived from apo template structures, we developed a new model based on holo X-ray structures of the hypoxia-inducible factor 2α (HIF-2α) PAS B domain, targeted to improve the accuracy of the binding site for molecular docking applications. We experimentally confirmed the ability of two HIF-2α crystallographic ligands to bind to the mAhR with relatively high affinity and demonstrated that they are AhR agonists, thus justifying the use of the holo HIF-2α structures as templates. A specific modeling/docking approach was proposed to predict the binding modes of AhR ligands in the modeled LBD. It was validated by comparison of the calculated and the experimental binding affinities of active THS ligands and TCDD for the mAhR and by functional activity analysis using several mAhR mutants generated on the basis of the modeling results. Finally the ability of the proposed approach to reproduce the different affinities of TCDD for AhRs of different species was confirmed, and a first test of its reliability in virtual screening is carried out by analyzing the correlation between the calculated and experimental binding affinities of a set of 14 PCDDs.

Molecular modeling and functional confirmation of a predicted fatty acid binding site of mitochondrial aspartate aminotransferase

Molecular interactions are necessary for proteins to perform their functions. The identification of a putative plasma membrane fatty acid transporter as mitochondrial aspartate aminotransferase (mAsp-AT) indicated that the protein must have a fatty acid binding site. Molecular modeling suggests that such a site exists in the form of a 500-Å(3) hydrophobic cleft on the surface of the molecule and identifies specific amino acid residues that are likely to be important for binding. The modeling and comparison with the cytosolic isoform indicated that two residues (Arg201 and Ala219) were likely to be important to the structure and function of the binding site. These residues were mutated to determine if they were essential to that function. Expression constructs with wild-type or mutated cDNAs were produced for bacteria and eukaryotic cells. Proteins expressed in Escherichia coli were tested for oleate binding affinity, which was decreased in the mutant proteins. 3T3 fibroblasts were transfected with expression constructs for both normal and mutated forms. Plasma membrane expression was documented by indirect immunofluorescence before [(3)H]oleic acid uptake kinetics were assayed. The V(max) for uptake was significantly increased by overexpression of the wild-type protein but changed little after transfection with mutated proteins, despite their presence on the plasma membrane. The hydrophobic cleft in mAsp-AT can serve as a fatty acid binding site. Specific residues are essential for normal fatty acid binding, without which fatty acid uptake is compromised. These results confirm the function of this protein as a fatty acid binding protein.

Triggered release of pharmacophores from [Ni(HAsO₃)]-loaded polymer-caged nanobin enhances pro-apoptotic activity: a combined experimental and theoretical study

Nanoscale drug delivery platforms can provide an attractive therapeutic strategy for cancer treatments, as they can substantially reduce the adverse side effects associated with toxic small-molecule anticancer agents. For enhanced therapeutic efficacy to be achieved with such platforms, a tumor-specific drug-release trigger is a critical requirement. This article reports the use of a pH-sensitive polymer network that surrounds a nanoscale liposome core to trigger the release of both encapsulated hydrophilic, membrane-impermeable Ni(II) cations and amphipathic, membrane-permeable As(III) anticancer agents under acidic conditions commonly encountered in hypoxic tumor tissues and late endosomes. Computational modeling studies provide clear evidence that the acid-triggered drug-release mechanism for this polymer-caged nanobin (PCN) platform arises from a pH- and temperature-responsive conformation change of the cross-linked polymer cage. As a result, the simultaneous release of both of the active agents in this multicomponent therapeutic enhances the pro-apoptotic activity of As(III) while diminishing its acute toxicity, potentially reducing the undesirable side effects commonly associated with this free drug. The ability to engender acid-triggered release of drugs co-encapsulated inside a liposomal template makes drug delivery using PCN an attractive strategy for triggered drug release.

Modelling proteins: conformational sampling and reconstruction of folding kinetics

In the last decades biomolecular simulation has made tremendous inroads to help elucidate biomolecular processes in-silico. Despite enormous advances in molecular dynamics techniques and the available computational power, many problems involve long time scales and large-scale molecular rearrangements that are still difficult to sample adequately. In this review we therefore summarise recent efforts to fundamentally improve this situation by decoupling the sampling of the energy landscape from the description of the kinetics of the process. Recent years have seen the emergence of many advanced sampling techniques, which permit efficient characterisation of the relevant family of molecular conformations by dispensing with the details of the short-term kinetics of the process. Because these methods generate thermodynamic information at best, they must be complemented by techniques to reconstruct the kinetics of the process using the ensemble of relevant conformations. Here we review recent advances for both types of methods and discuss their perspectives to permit efficient and accurate modelling of large-scale conformational changes in biomolecules. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Investigation of protein folding by coarse-grained molecular dynamics with the UNRES force field

Coarse-grained molecular dynamics simulations offer a dramatic extension of the time-scale of simulations compared to all-atom approaches. In this article, we describe the use of the physics-based united-residue (UNRES) force field, developed in our laboratory, in protein-structure simulations. We demonstrate that this force field offers about a 4000-times extension of the simulation time scale; this feature arises both from averaging out the fast-moving degrees of freedom and reduction of the cost of energy and force calculations compared to all-atom approaches with explicit solvent. With massively parallel computers, microsecond folding simulation times of proteins containing about 1000 residues can be obtained in days. A straightforward application of canonical UNRES/MD simulations, demonstrated with the example of the N-terminal part of the B-domain of staphylococcal protein A (PDB code: 1BDD, a three-alpha-helix bundle), discerns the folding mechanism and determines kinetic parameters by parallel simulations of several hundred or more trajectories. Use of generalized-ensemble techniques, of which the multiplexed replica exchange method proved to be the most effective, enables us to compute thermodynamics of folding and carry out fully physics-based prediction of protein structure, in which the predicted structure is determined as a mean over the most populated ensemble below the folding-transition temperature. By using principal component analysis of the UNRES folding trajectories of the formin-binding protein WW domain (PDB code: 1E0L; a three-stranded antiparallel beta-sheet) and 1BDD, we identified representative structures along the folding pathways and demonstrated that only a few (low-indexed) principal components can capture the main structural features of a protein-folding trajectory; the potentials of mean force calculated along these essential modes exhibit multiple minima, as opposed to those along the remaining modes that are unimodal. In addition, a comparison between the structures that are representative of the minima in the free-energy profile along the essential collective coordinates of protein folding (computed by principal component analysis) and the free-energy profile projected along the virtual-bond dihedral angles gamma of the backbone revealed the key residues involved in the transitions between the different basins of the folding free-energy profile, in agreement with existing experimental data for 1E0L .

Effect of structural and conformation modifications, including backbone cyclization, of hydrophilic hexapeptides on their intestinal permeability and enzymatic stability

A library of 18 hexapeptide analogs was synthesized, including sub-libraries of N- or C-methylation of the parent hexapeptide Phe-Gly-Gly-Gly-Gly-Phe, as well as backbone cyclized analogs of each linear analog with various ring sizes. N- or C-methylation as well as cyclization (but not backbone cyclization) have been suggested to improve intestinal permeability and metabolic stability of peptides in general. Here we aimed to assess their applicability to hydrophilic peptides. The intestinal permeability (Papp) of the 18-peptide library was in the range of 0.2-6.8 x 10-6 cm/sec. Based on several tests, we concluded that the absorption mechanism of all tested analogs is paracellular, regardless of the structural or conformational modifications. In all cases, backbone cyclization increased Papp (5-fold) in comparison to the linear analogs due to the smaller 3D size and also dramatically decreased peptide proteolysis by brush border enzymes. N- or C-methylation did not enhance the permeability of the linear analogs in this series.

NMR determination of the bioactive conformation of peloruside A bound to microtubules

We report here on the determination of the conformation of Peloruside A bound to biochemically stabilized microtubules, by using TR-NOESY NMR experiments. As a previous step, the conformation of the free molecule in water solution has also been deduced. Despite the large size of the ring, Peloruside A mainly adopts two conformations in water solution. A conformational selection process takes place, and the microtubules-bound conformer is one of those present in the water solution, different than that existing in chloroform medium. A model of the binding mode to tubulin has also been proposed, by docking the bioactive conformation of peloruside, which involves the alpha-tubulin monomer, in contrast with taxol, which binds to the beta-monomer.

Synthesis and conformational analysis of an α-d-mannopyranosyl-(1→2)-α-d-mannopyranosyl-(1→6)-α-d-mannopyranose mimic

A mimic of a (1-->2),(1-->6)-mannotrioside was synthesized by replacing the central mannose unit with an enantiomerically pure, conformationally stable trans-diaxial cyclohexanediol. The three-dimensional structure of the molecule was investigated by NMR spectroscopy supported by molecular modelling and was compared to the known features of the natural mannotrioside.

Developing a coarse-grained force field for the diblock copolymer poly(styrene-b-butadiene) from atomistic simulation

We have developed a coarse-grained force field for the poly(styrene-b-butadiene) diblock copolymer. We describe the computational methods and discuss how they were applied to develop a coarse-grained force field for this diblock copolymer from the atomistic simulation. The new force field contains three different bonds, four angles, five dihedral angles, and three nonbonded terms. We successfully tested this coarse-grained model against the chain properties, including static and dynamic properties, derived from the atomistic simulation; the results suggest that the coarse-grained force field is an effective model.

Overcoming entropic barrier with coupled sampling at dual resolutions

An enhanced sampling method is proposed for ab initio protein folding simulations. The new method couples a high-resolution model for accuracy and a low-resolution model for efficiency. It aims to overcome the entropic barrier found in the exponentially large protein conformational space when a high-resolution model, such as an all-atom molecular mechanics force field, is used. The proposed method is designed to satisfy the detailed balance condition so that the Boltzmann distribution can be generated in all sampling trajectories in both high and low resolutions. The method was tested on model analytical energy functions and ab initio folding simulations of a beta-hairpin peptide. It was found to be more efficient than replica-exchange method that is used as its building block. Analysis with the analytical energy functions shows that the number of energy calculations required to find global minima and to converge mean potential energies is much fewer with the new method. Ergodic measure shows that the new method explores the conformational space more rapidly. We also studied imperfect low-resolution energy models and found that the introduction of errors in low-resolution models does decrease its sampling efficiency. However, a reasonable increase in efficiency is still observed when the global minima of the low-resolution models are in the vicinity of the global minimum basin of the high-resolution model. Finally, our ab initio folding simulation of the tested peptide shows that the new method is able to fold the peptide in a very short simulation time. The structural distribution generated by the new method at the equilibrium portion of the trajectory resembles that in the equilibrium simulation starting from the crystal structure.

Synthesis and evaluation of bivalent NDP-alpha-MSH(7) peptide ligands for binding to the human melanocortin receptor 4 (hMC4R)

We demonstrate the potential utility of multivalent ligands as targeting agents for cancer imaging or therapy by determining the binding of homobivalent ligands to their corresponding receptors. This manuscript details the synthesis and evaluation of a series of bivalent ligands containing two copies of the truncated heptapeptide version of [Nle4-D-Phe7]-alpha-melanocyte stimulating hormone (NDP-alpha-MSH), referred to as MSH(7). These were connected with various semirigid linkers containing Pro-Gly repeats, with or without flexible poly(ethylene glycol) (PEGO) moieties at their termini. Modeling data suggest a distance of 20-50 A between the ligand binding sites of two adjacent G-protein coupled receptors, GPCRs. These bivalent ligands were observed to bind with higher affinity compared to their monovalent counterparts. Data suggest these ligands may be capable of cross-linking adjacent receptors. An optimal linker length of 25 +/- 10 A, inferred from these ligands, correlated well with the inter-receptor distance estimated through modeling. Although there was no difference in maximal binding affinities between the ligands constructed with the Pro-Gly repeats versus those constructed with the PEGO inserts, the PEGO-containing ligands bound with high affinities over a greater range of linker lengths.

A Combined Computational and Experimental Approach for the Analysis of the Enantioselective Potential of a New Macrocyclic Receptor for N-Protected α-Amino Acids

A new macrocyclic receptor incorporating a thiourea moiety has been synthesised. Crystal structures of the macrocycle showed that the receptor has a rigid backbone but the thiourea moiety can orientate itself to bind to a DMSO solvent molecule. Force-field (MMFFs) calculations were performed to model the macrocycle and its binding properties with respect to N-protected amino acids, which were measured experimentally by NMR titration. Binding free energies were calculated by using the mode integration algorithm (MINTA) or free-energy perturbation (FEP). Excellent qualitative agreement with experiment was obtained. To further exploit the accuracy of the free-energy predictions for this system, the faster free-energy algorithm MINTA was used as a prediction tool to test the binding affinity of the macrocycle towards a series of several other amino acid derivatives, which speeded up considerably the screening process and reduced laboratory costs.

A Simple Model System for the Study of Carbohydrate−Aromatic Interactions

A molecular scaffold was identified which enables the establishment of intramolecular interactions between a monosaccharide and a nearby phenyl ring. A group of molecules containing four different monosaccharides (glucose, galactose, N-acetyl-glucosamine, and N-acetyl-galactosamine) was synthesized and used to investigate the extent and nature of this carbohydrate-arene interaction, as well as the effect on the overall 3D structure of the molecules involved. The sugar-aromatic distance was evaluated by rigorous NMR studies supported by molecular modeling and found to be constant throughout the series, independent of the nature of the sugar and of the conformational behavior of the fragment connecting the two elements. Ab initio calculations at the B3LYP/DZV(2d,p) level of theory enable the analysis of the electronic nature of the interaction. The study shows that, given the opportunity, persistent intramolecular aromatic-sugar interactions can be established and can significantly influence overall molecular shape and energetics. These results have important implications in the design of structural mimics of oligosaccharides.