GROMACS: Fast, flexible, and free

Electrostatic similarity of proteins: application of three dimensional spherical harmonic decomposition

We present a method for describing and comparing global electrostatic properties of biomolecules based on the spherical harmonic decomposition of electrostatic potential data. Unlike other approaches our method does not require any prior three dimensional structural alignment. The electrostatic potential, given as a volumetric data set from a numerical solution of the Poisson or Poisson-Boltzmann equation, is represented with descriptors that are rotation invariant. The method can be applied to large and structurally diverse sets of biomolecules enabling to cluster them according to their electrostatic features.

Understanding structural/functional properties of immunoconjugates for cancer therapy by computational approaches

Monoclonal antibodies coupled to highly toxic molecules (immunoconjugates) are currently being developed for cancer therapy. We have used an in silico procedure for evaluating some physicochemical properties of two tumor-targeting anti-HER2 immunoconjugates: (a) the single-chain antibody scFv(FRP5) linked to a bacterial toxin, that has been recently progressed to phase I clinical trial in human cancer; (b) the putative molecule formed by the intrinsically stable scFv(800E6), which has been proposed as toxin carrier to cancer cells in human therapy, joined to the same toxin of (a). Structural models of the immunoconjugates have been built by homology modeling and assessed by molecular dynamics simulations. The trajectories have been analyzed to extract some biochemical properties and to assess the potential effects of the toxin on the structure and dynamics of the anti-HER2 antibodies. The results of the computational approach indicate that the antibodies maintain their correct folding even in presence of the toxin, whereas a certain stiffness in correspondence of some structural regions is observed. Furthermore, the toxin does not seem to affect the antibody solubility, whereas it enhances the structural stability. The proposed computational approach represent a promising tool for analyzing some physicochemical properties of immunoconjugates and for predicting the effects of the linked toxin on structure, dynamics, and functionality of the antibodies.

Protein structure prediction using basin-hopping

Associative memory Hamiltonian structure prediction potentials are not overly rugged, thereby suggesting their landscapes are like those of actual proteins. In the present contribution we show how basin-hopping global optimization can identify low-lying minima for the corresponding mildly frustrated energy landscapes. For small systems the basin-hopping algorithm succeeds in locating both lower minima and conformations closer to the experimental structure than does molecular dynamics with simulated annealing. For large systems the efficiency of basin-hopping decreases for our initial implementation, where the steps consist of random perturbations to the Cartesian coordinates. We implemented umbrella sampling using basin-hopping to further confirm when the global minima are reached. We have also improved the energy surface by employing bioinformatic techniques for reducing the roughness or variance of the energy surface. Finally, the basin-hopping calculations have guided improvements in the excluded volume of the Hamiltonian, producing better structures. These results suggest a novel and transferable optimization scheme for future energy function development.

Binless estimation of the potential of mean force

For a system in thermal equilibrium, described by classical statistical mechanics, we derive an unbiased estimator for the marginal probability distribution of a coordinate of interest, rho( x). This result provides a "binless" method for estimating the potential of mean force, Phi = -beta (-1) ln rho, eliminating the need to construct histograms or perform numerical thermodynamic integration. In our method, the distribution that we seek to compute is expressed as the sum of a reference distribution, rho 0(x)essentially an initial guess or estimate of rho( x)and a correction term. While the method is valid for arbitrary rho 0, we speculate that an accurate choice of the reference distribution improves the convergence of the method. Using a model molecule, simulated both in vacuum and in solvent, we validate our proposed approach and compare its performance with the histogram and thermodynamic integration methods. We also discuss and validate an extension in which our approach is used in combination with a biasing force, meant to improve uniform sampling of the coordinate of interest.

The selectivity of K+ ion channels: testing the hypotheses

How K(+) channels are able to conduct certain cations yet not others remains an important but unresolved question. The recent elucidation of the structure of NaK, an ion channel that conducts both Na(+) and K(+) ions, offers an opportunity to test the various hypotheses that have been put forward to explain the selectivity of K(+) ion channels. We test the snug-fit, field-strength, and over-coordination hypotheses by comparing their predictions to the results of classical molecular dynamics simulations of the K(+) selective channel KcsA and the less selective channel NaK embedded in lipid bilayers. Our results are incompatible with the so-called strong variant of the snug-fit hypothesis but are consistent with the over-coordination hypothesis and neither confirm nor refute the field-strength hypothesis. We also find that the ions and waters in the NaK selectivity filter unexpectedly move to a new conformation in seven K(+) simulations: the two K(+) ions rapidly move from site S4 to S2 and from the cavity to S4. At the same time, the selectivity filter narrows around sites S1 and S2 and the carbonyl oxygen atoms rotate 20 degrees -40 degrees inwards toward the ion. These motions diminish the large structural differences between the crystallographic structures of the selectivity filters of NaK and KcsA and appear to allow the binding of ions to S2 of NaK at physiological temperature.

Molecular dynamics study on the biophysical interactions of seven green tea catechins with lipid bilayers of cell membranes

Molecular dynamics simulations were performed to study the interactions of bioactive catechins (flavonoids) commonly found in green tea with lipid bilayers, as a model for cell membranes. Previously, multiple experimental studies rationalized catechin's anticarcinogenic, antibacterial, and other beneficial effects in terms of physicochemical molecular interactions with the cell membranes. To contribute toward understanding the molecular role of catechins on the structure of cell membranes, we present simulation results for seven green tea catechins in lipid bilayer systems representative of HepG2 cancer cells. Our simulations show that the seven tea catechins evaluated have a strong affinity for the lipid bilayer via hydrogen bonding to the bilayer surface, with some of the smaller catechins able to penetrate underneath the surface. Epigallocatechin-gallate (EGCG) showed the strongest interaction with the lipid bilayer based on the number of hydrogen bonds formed with lipid headgroups. The simulations also provide insight into the functional characteristics of the catechins that distinguish them as effective compounds to potentially alter the lipid bilayer properties. The results on the hydrogen-bonding effects, described here for the first time, may contribute to a better understanding of proposed multiple molecular mechanisms of the action of catechins in microorganisms, cancer cells, and tissues.

S-nitrosylation of syntaxin 1 at Cys(145) is a regulatory switch controlling Munc18-1 binding

Exocytosis is regulated by NO in many cell types, including neurons. In the present study we show that syntaxin 1a is a substrate for S-nitrosylation and that NO disrupts the binding of Munc18-1 to the closed conformation of syntaxin 1a in vitro. In contrast, NO does not inhibit SNARE {SNAP [soluble NSF (N-ethylmaleimide-sensitive fusion protein) attachment protein] receptor} complex formation or binding of Munc18-1 to the SNARE complex. Cys(145) of syntaxin 1a is the target of NO, as a non-nitrosylatable C145S mutant is resistant to NO and novel nitrosomimetic Cys(145) mutants mimic the effect of NO on Munc18-1 binding in vitro. Furthermore, expression of nitrosomimetic syntaxin 1a in living cells affects Munc18-1 localization and alters exocytosis release kinetics and quantal size. Molecular dynamic simulations suggest that NO regulates the syntaxin-Munc18 interaction by local rearrangement of the syntaxin linker and H3c regions. Thus S-nitrosylation of Cys(145) may be a molecular switch to disrupt Munc18-1 binding to the closed conformation of syntaxin 1a, thereby facilitating its engagement with the membrane fusion machinery.

The negatively charged amino acids in the lumenal loop influence the pigment binding and conformation of the major light-harvesting chlorophyll a/b complex of photosystem II

The major chlorophyll (Chl) a/b complexes of photosystem II (LHCIIb), in addition to their primary light-harvesting function, play key roles in the organization of the granal ultrastructure of the thylakoid membranes and in various regulatory processes. These functions depend on the structural stability and flexibility of the complexes. The lumenal side of LHCIIb is exposed to broadly variable pH environments, due to the build-up and decay of the pH gradient during photosynthesis. Therefore, the negatively charged amino acids in the lumenal loop might be of paramount importance for adjusting the structure and functions of LHCIIb. In order to clarify the structural roles of these residues, we investigated the pigment stoichiometries, absorption, linear and circular dichroism spectra of the reconstituted LHCIIb complexes, in which the negatively charged amino acids in the lumenal loop were exchanged to neutral ones (E94G, E107V and D111V). The mutations influenced the pigment binding and the molecular architecture of the complexes. Exchanging E94 to G destabilized the 3(10) helix in the lumenal loop structure and led to an acquired pH sensitivity of the LHCIIb structure. We conclude that these amino acids are important not only for pigment binding in the complexes, but also in stabilizing the conformation of LHCIIb at different pHs.

Effects of oxidation, pH and lipids on amyloidogenic peptide structure: implications for fibril formation?

We have performed experimental and computational studies to investigate the influences of phospholipids, methionine oxidation and acidic pH on amyloid fibril formation by a peptide derived from human apolipoprotein C-II (apoC-II), a known component of proteinaceous atherosclerotic plaques. Fibril growth monitored by thioflavin T fluorescence revealed inhibition under lipid-rich and oxidising conditions. We subsequently performed fully-solvated atomistic molecular dynamics (MD) simulations of the peptide monomer to study its conformations under both fibril favouring (neutral and low pH) and inhibiting (lipid-rich and oxidising) conditions. Examination of the chain topology, backbone hydrogen-bonding patterns and aromatic sidechain orientations of the peptide under different conditions reveals that, while the peptide adopts similar structures under the fibril-favouring conditions, significantly different structures are obtained under fibril-disruptive conditions. Based on our results, we advance hypotheses for the roles of peptide conformation on aggregation and fibrillisation propensities.

Thr729 in human topoisomerase I modulates anti-cancer drug resistance by altering protein domain communications as suggested by molecular dynamics simulations

The role of Thr729 in modulating the enzymatic function of human topoisomerase I has been characterized by molecular dynamics (MD) simulation. In detail, the structural–dynamical behaviour of the Thr729Lys and the Thr729Pro mutants have been characterized because of their in vivo and in vitro functional properties evidenced in the accompanying paper. Both mutants can bind to the DNA substrate and are enzymatically active, but while Thr729Lys is resistant even at high concentration of the camptothecin (CPT) anti-cancer drug, Thr729Pro shows only a mild reduction in drug sensitivity and in DNA binding. MD simulations show that the Thr729Lys mutation provokes a structural perturbation of the CPT-binding pocket. On the other hand, the Thr729Pro mutant maintains the wild-type structural scaffold, only increasing its rigidity. The simulations also show the complete abolishment, in the Thr729Lys mutant, of the protein communications between the C-terminal domain (where the active Tyr723 is located) and the linker domain, that plays an essential role in the control of the DNA rotation, thus explaining the distributive mode of action displayed by this mutant.

Not only enthalpy: large entropy contribution to ion permeation barriers in single-file channels

The effect of channel length on the barrier for potassium ion permeation through single-file channels has been studied by means of all-atom molecular dynamics simulations. Using series of peptidic gramicidin-like and simplified ring-structured channels, both embedded in model membranes, we obtained two distinct types of behavior: saturation of the central free energy barriers for peptidic channels and a linear increase in simplified ring-structured channels with increasing channel length. The saturation of the central free energy barrier for the peptidic channels occurs at relatively short lengths, and it is correlated with the desolvation from the bulk water. Remarkably, decomposition of free energy barriers into enthalpic and entropic terms reveals an entropic cost for ion permeation. Furthermore, this entropic cost dominates the ion permeation free energy barrier, since the corresponding free energy contribution is higher than the enthalpic barrier. We conclude that the length dependence of the free energy is enthalpy-dominated, but the entropy is the major contribution to the permeation barrier. The decrease in rotational water motion and the reduction of channel mobility are putative origins for the overall entropic penalty.

The catalytic site of glutathione peroxidases

In GPxs, the redox-active Se or S, is at hydrogen bonding distance from Gln and Trp residues that contribute to catalysis. From sequence homology of >400 sequences and modeling of the DmGPx as a paradigm, Asn136 emerged as a fourth essential component of the active site. Mutational substitution of Asn136 by His, Ala, or Asp results in a dramatic decline of specific activity. Kinetic analysis indicates that k(+1), the rate constant for the oxidation of the enzyme, decreases by two to three orders of magnitude, whereas the reductive steps characterized by k'(+2) are less affected. Accordingly, MS/MS analysis shows that in Asn136 mutants, the peroxidatic Cys45 stays largely reduced also in the presence of a hydroperoxide, whereas in the wild-type enzyme, it is oxidized, forming a disulfide with the resolving Cys. Computational calculation of pK(a) values indicates that the residues facing the catalytic thiol, Asn136, Gln80, and, to a lesser extent Trp135, contribute to the dissociation of the thiol group, Asn136 being most relevant. These data disclose that the catalytic site of GPxs has to be redrawn as a tetrad, including Asn136, and suggest a mechanism accounting for the extraordinary catalytic efficiency of GPxs.

Construction of the free energy landscape of biomolecules via dihedral angle principal component analysis

A systematic approach to construct a low-dimensional free energy landscape from a classical molecular dynamics (MD) simulation is presented. The approach is based on the recently proposed dihedral angle principal component analysis (dPCA), which avoids artifacts due to the mixing of internal and overall motions in Cartesian coordinates and circumvents problems associated with the circularity of angular variables. Requiring that the energy landscape reproduces the correct number, energy, and location of the system's metastable states and barriers, the dimensionality of the free energy landscape (i.e., the number of essential components) is obtained. This dimensionality can be determined from the distribution and autocorrelation of the principal components. By performing an 800 ns MD simulation of the folding of hepta-alanine in explicit water and using geometric and kinetic clustering techniques, it is shown that a five-dimensional dPCA energy landscape is a suitable and accurate representation of the full-dimensional landscape. In the second step, the dPCA energy landscape can be employed (e.g., in a Langevin simulation) to facilitate a detailed investigation of biomolecular dynamics in low dimensions. Finally, several ways to visualize the multidimensional energy landscape are discussed.

A theoretical study of repeating sequence in HRP II: a combination of molecular dynamics simulations and (17)O quadrupole coupling tensors

Histidine rich protein II derived peptide (HRP II 169-182) was investigated by molecular dynamics, MD, simulation and (17)O electric field gradient, EFG, tensor calculations. MD simulation was performed in water at 300 K with alpha-helix initial structure. It was found that peptide loses its initial alpha-helix structure rapidly and is converted to random coil and bent secondary structures. To understand how peptide structure affects EFG tensors, initial structure and final conformations resulting from MD simulations were used to calculate (17)O EFG tensors of backbone carbonyl oxygens. Calculations were performed using B3LYP method and 6-31+G basis set. Calculated (17)O EFG tensors were used to evaluate quadrupole coupling constants, QCC, and asymmetry parameters, eta(Q). Difference between the calculated QCC and eta(Q) values revealed how hydrogen-bonding interactions affect EFG tensors at the sites of each oxygen nucleus.

Importin-beta: structural and dynamic determinants of a molecular spring

The beta-karyopherin/RanGTP system constitutes the largest known family of cellular cargo transporters. The flexibility of the karyopherin transport receptors is the key to their versatility in binding cargoes of different shape and size. Despite strong binding of the Ran complex, the comparably low energy associated with GTP hydrolysis suffices to drive dissociation and fuel the transport cycle. Here, we elucidate the drastic structural dynamics of the prototypic karyopherin, importin-beta, and show that its flexibility also solves this energetic puzzle. Our nonequilibrium atomistic simulations reveal fast conformational changes, validated by small-angle X-ray scattering data, and unusually large structural fluctuations. The characteristic dynamic patterns of importin-beta and the observed unfolding pathway of the IBB domain suggest a cooperative mechanism of importin-beta function in the nucleus. We propose a molecular model in which the stored energy and structural dynamics account for an exchange pathway that explains the high observed rates of nucleocytoplasmic transport. Karyopherins utilize a mechanism of entropy/enthalpy control that might be a general feature of highly flexible proteins involved in protein-protein interactions.

The investigation of interactions of kappa-Hefutoxin1 with the voltage-gated potassium channels: a computational simulation

Kappa-Hefutoxin1 is a K(+) channel-blocking toxin from the scorpion Heterometrus fluvipes. It is a 22-residue protein that adapts a novel fold of two parallel helices linked by two disulfide bridges without beta-sheets. Recognition of interactions of kappa-Hefutoxin1 with the voltage-gated potassium channels, Kv1.1, Kv1.2, and Kv1.3, was studied by 3D-Dock software package. All structures of kappa-Hefutoxin1 were considered during the simulations, which indicated that even small changes in the structure of kappa-Hefutoxin1 considerably affected both the recognition and the binding between kappa-Hefutoxin1 and the Kv1 channels. kappa-Hefutoxin1 is located around the extracellular part of the Kv1 channels, making contacts with its helices. Lys 19, Tyr 5, Arg 6, Trp 9, or Arg 10 in the toxin and residues Asp 402, His 404, Thr 407,Gly 401, and Asp 386 in each subunit of the Kv potassium channel are the key residues for the toxin-channel recognition. Moreover, the simulation result demonstrates that the hydrophobic interactions are important in interaction of negatively charged toxins with potassium channels. The results of our docking/molecular dynamics simulations indicate that our 3D model structure of the kappa-Hefutoxin1-complex is both reasonable and acceptable and could be helpful for smarter drug design and the blocking agents of Kv1 channels.

CHARMM-GUI: a web-based graphical user interface for CHARMM

CHARMM is an academic research program used widely for macromolecular mechanics and dynamics with versatile analysis and manipulation tools of atomic coordinates and dynamics trajectories. CHARMM-GUI, http://www.charmm-gui.org, has been developed to provide a web-based graphical user interface to generate various input files and molecular systems to facilitate and standardize the usage of common and advanced simulation techniques in CHARMM. The web environment provides an ideal platform to build and validate a molecular model system in an interactive fashion such that, if a problem is found through visual inspection, one can go back to the previous setup and regenerate the whole system again. In this article, we describe the currently available functional modules of CHARMM-GUI Input Generator that form a basis for the advanced simulation techniques. Future directions of the CHARMM-GUI development project are also discussed briefly together with other features in the CHARMM-GUI website, such as Archive and Movie Gallery.

Protein structure analysis online

Computational biology/chemistry tools are used in most areas of life/health science research. These methods are continually being developed and their use can present difficulties for both experienced and novice investigators. To facilitate the use of these applications, many packages have been implemented online during these last 5 years. This unit focuses on online computational methods with a special emphasis on structural refinement/atomic simulations, protein electrostatic calculations, searches for functional sites, searches for druggable pockets, protein docking and small molecule docking, and prediction of potential impact of amino acid variations on the structure and function of the protein molecules.

A dihydropyridine receptor alpha1s loop region critical for skeletal muscle contraction is intrinsically unstructured and binds to a SPRY domain of the type 1 ryanodine receptor

The II-III loop of the dihydropyridine receptor (DHPR) alpha(1s) subunit is a modulator of the ryanodine receptor (RyR1) Ca(2+) release channel in vitro and is essential for skeletal muscle contraction in vivo. Despite its importance, the structure of this loop has not been reported. We have investigated its structure using a suite of NMR techniques which revealed that the DHPR II-III loop is an intrinsically unstructured protein (IUP) and as such belongs to a burgeoning structural class of functionally important proteins. The loop does not possess a stable tertiary fold: it is highly flexible, with a strong N-terminal helix followed by nascent helical/turn elements and unstructured segments. Its residual structure is loosely globular with the N and C termini in close proximity. The unstructured nature of the II-III loop may allow it to easily modify its interaction with RyR1 following a surface action potential and thus initiate rapid Ca(2+) release and contraction. The in vitro binding partner for the II-III was investigated. The II-III loop interacts with the second of three structurally distinct SPRY domains in RyR1, whose function is unknown. This interaction occurs through two preformed N-terminal alpha-helical regions and a C-terminal hydrophobic element. The A peptide corresponding to the helical N-terminal region is a common probe of RyR function and binds to the same SPRY domain as the full II-III loop. Thus the second SPRY domain is an in vitro binding site for the II-III loop. The possible in vivo role of this region is discussed.

Functional and computational studies of the ligand-associated metal binding site of beta3 integrins

A combination of experimental and computational approaches was used to provide a structural context for the role of the beta3 integrin subunit ligand-associated metal binding site (LIMBS) in the binding of physiological ligands to beta3 integrins. Specifically, we have carried out (1) adhesion assays on cells expressing normal alphaIIbeta3, normal alphaVbeta3, or the corresponding beta3 D217A LIMBS mutants; and (2) equilibrium and nonequilibrium (steered) molecular dynamics (MD) simulations of eptifibatide in complex with either a fully hydrated normal alphaIIbeta3 integrin fragment (alphaIIb beta-propeller and the beta3 betaA (I-like), hybrid, and PSI domains) or the equivalent beta3 D217A mutant. Normal alphaIIbeta3 expressing cells adhered to immobilized fibrinogen and echistatin, whereas cells expressing the alphaIIbeta3 D217A LIMBS mutant failed to adhere to either ligand. Similarly, the equivalent alphaVbeta3 mutant was unable to support adhesion to vitronectin or fibrinogen. The alphaIIbeta3 D217A mutation increased the binding of mAb AP5, which recognizes a ligand-induced binding site (LIBS) in the beta3 PSI domain, indicating that this mutation induced allosteric changes in the protein. Steered MD simulating the unbinding of eptifibatide from either normal alphaIIbeta3 or the equivalent beta3 D217A mutant suggested that the reduction in ligand binding caused by the LIMBS mutant required the loss of both the LIMBS and the metal ion-dependent adhesion site (MIDAS) metal ions. Our computational results indicate that the LIMBS plays a crucial role in ligand binding to alphaIIbeta3 by virtue of its effects on the coordination of the MIDAS.

Thermodynamics of peptide insertion and aggregation in a lipid bilayer

A variety of biomolecules mediate physiological processes by inserting and reorganizing in cell membranes, and the thermodynamic forces responsible for their partitioning are of great interest. Recent experiments provided valuable data on the free energy changes associated with the transfer of individual amino acids from water to membrane. However, a complete picture of the pathways and the associated changes in energy of peptide insertion into a membrane remains elusive. To this end, computational techniques supplement the experimental data with atomic-level details and shed light on the energetics of insertion. Here, we employed the technique of umbrella sampling in a total 850 ns of all-atom molecular dynamics simulations to study the free energy profile and the pathway of insertion of a model hexapeptide consisting of a tryptophan and five leucines (WL5). The computed free energy profile of the peptide as it travels from bulk solvent through the membrane core exhibits two minima: a local minimum at the water−membrane interface or the headgroup region and a global minimum at the hydrophobic−hydrophilic interface close to the lipid glycerol region. A rather small barrier of roughly 1 kcal mol⁻¹ exists at the membrane core, which is explained by the enhanced flexibility of the peptide when deeply inserted. Combining our results with those in the literature, we present a thermodynamic model for peptide insertion and aggregation which involves peptide aggregation upon contact with the membrane at the solvent−lipid headgroup interface.

Modeling and molecular dynamics simulation of the human gonadotropin-releasing hormone receptor in a lipid bilayer

In the present study, a model for the human gonadotropin-releasing hormone receptor embedded in an explicit lipid bilayer was developed. The final conformation was obtained by extensive molecular dynamics simulations of a homology model based on the bovine rhodopsin crystal structure. The analysis of the receptor structure allowed us to detect a number of specific contacts between different amino acid residues, as well as water- and lipid-mediated interactions. These interactions were stable in six additional independent 35 ns long simulations at 310 and 323 K, which used the refined model as the starting structure. All loops, particularly the extracellular loop 2 and the intracellular loop 3, exhibited high fluctuations, whereas the transmembrane helices were more static. Although other models of this receptor have been previously developed, none of them have been subjected to extensive molecular dynamics simulations, and no other three-dimensional structure is publicly available. Our results suggest that the presence of ions as well as explicit solvent and lipid molecules are critical for the structure of membrane protein models, and that molecular dynamics simulations are certainly useful for their refinement.

Peptide aggregation and pore formation in a lipid bilayer: a combined coarse-grained and all atom molecular dynamics study

We present a simulation study where different resolutions, namely coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, are used sequentially to combine the long timescale reachable by CG simulations with the high resolution of AA simulations, to describe the complete processes of peptide aggregation and pore formation by alamethicin peptides in a hydrated lipid bilayer. In the 1-micros CG simulations the peptides spontaneously aggregate in the lipid bilayer and exhibit occasional transitions between the membrane-spanning and the surface-bound configurations. One of the CG systems at t = 1 micros is reverted to an AA representation and subjected to AA simulation for 50 ns, during which water molecules penetrate the lipid bilayer through interactions with the peptide aggregates, and the membrane starts leaking water. During the AA simulation significant deviations from the alpha-helical structure of the peptides are observed, however, the size and arrangement of the clusters are not affected within the studied time frame. Solid-state NMR experiments designed to match closely the setup used in the molecular dynamics simulations provide strong support for our finding that alamethicin peptides adopt a diverse set of configurations in a lipid bilayer, which is in sharp contrast to the prevailing view of alamethicin oligomers formed by perfectly aligned helical alamethicin peptides in a lipid bilayer.

HIV-1 capsid assembly inhibitor (CAI) peptide: structural preferences and delivery into human embryonic lung cells and lymphocytes

The Human immunodeficiency virus 1 derived capsid assembly inhibitor peptide (HIV-1 CAI-peptide) is a promising lead candidate for anti-HIV drug development. Its drawback, however, is that it cannot permeate cells directly. Here we report the transport of the pharmacologically active CAI-peptide into human lymphocytes and Human Embryonic Lung cells (HEL) using the BioShuttle platform. Generally, the transfer of pharmacologically active substances across membranes, demonstrated by confocal laser scanning microscopy (CLSM), could lead to a loss of function by changing the molecule's structure. Molecular dynamics (MD) simulations and circular dichroism (CD) studies suggest that the CAI-peptide has an intrinsic capacity to form a helical structure, which seems to be critical for the pharmacological effect as revealed by intensive docking calculations and comparison with control peptides. This coupling of the CAI-peptide to a BioShuttle-molecule additionally improved its solubility. Under the conditions described, the HIV-1 CAI peptide was transported into living cells and could be localized in the vicinity of the mitochondria.

Evidence for proton shuffling in a thioredoxin-like protein during catalysis

Proteins of the thioredoxin (Trx) superfamily catalyze disulfide-bond formation, reduction and isomerization in substrate proteins both in prokaryotic and in eukaryotic cells. All members of the Trx family with thiol-disulfide oxidoreductase activity contain the characteristic Cys-X-X-Cys motif in their active site. Here, using Poisson-Boltzmann-based protonation-state calculations based on 100-ns molecular dynamics simulations, we investigate the catalytic mechanism of DsbL, the most oxidizing Trx-like protein known to date. We observed several correlated transitions in the protonation states of the buried active-site cysteine and a neighboring lysine coupled to the exposure of the active-site thiolate. These results support the view of an internal proton shuffling mechanism during oxidation crucial for the uptake of two electrons from the substrate protein. Intramolecular disulfide-bond formation is probably steered by the conformational switch facilitating interaction with the active-site thiolate. A consistent catalytic mechanism for DsbL, probably conferrable to other proteins of the same class, is presented. Our results suggest a functional role of hydration entropy of active-site groups.