Biomedical subjects

Themis Lazaridis

Publications and source records attributed to Themis Lazaridis.

18 recordsLinked to original sources

PubMed2006

Energetic determinants of oligomeric state specificity in coiled coils.

The coiled coil is one of the simplest and best-studied protein structural motifs, consisting of two to five helices wound around each other. Empirical rules have been established on the tendency of different core sequences to form a certain oligomeric state but the physical forces behind this specificity are unclear. In this work, we model four sequences onto the structures of dimeric, trimeric, tetrameric, and pentameric coiled coils. We first examine the ability of an effective energy function (EEF1.1) to discriminate the correct oligomeric state for a given sequence. We find that inclusion of the translational, rotational, and side-chain conformational entropy is necessary for discriminating the native structures from their misassembled counterparts. The decomposition of the effective energy into residue contributions yields theoretical values for the oligomeric propensity of different residue types at different heptad positions. We find that certain calculated residue propensities are general and consistent with existing rules. For example, leucine at d favors dimers, leucine at a favors tetramers or pentamers, and isoleucine at a favors trimers. Other residue propensities are sequence context dependent. For example, glutamine at d favors trimers in one context and pentamers in another. Charged residues at e and g positions usually destabilize higher oligomers due to higher desolvation. Nonpolar residues at these positions confer pentamer specificity when combined with certain residues at positions a and d. Specifically, the pair Leua-Alag' or the inverse was found to stabilize the pentamer. The small energy gap between the native and misfolded counterparts explains why a few mutations at the core sites are sufficient to induce a change in the oligomeric state of these peptides. A large number of possible experiments are suggested by these results.

Amino Acid Sequence↗

PubMed2006

Modeling a spin-labeled fusion peptide in a membrane: implications for the interpretation of EPR experiments.

Site-directed spin-labeling and electron paramagnetic resonance are powerful tools for studying structure and conformational dynamics of proteins, especially in membranes. The position of the spin label is used as an indicator of the position of the site to which it is attached. The interpretation of these experiments is based on the assumptions that the spin label does not affect the peptide configuration and that it has a fixed orientation and distance with respect to the protein backbone. Here, the validity of these assumptions is examined through implicit membrane molecular dynamics simulations of the influenza hemagglutinin fusion peptide that has been labeled with methanethiosulfonate spin label. We find that the methanethiosulfonate spin label can occasionally induce peptide orientations that differ from those adopted by the wild-type peptide. Furthermore, the spin-label resides, on average, several Angstroms deeper in the membrane than the corresponding backbone C(alpha)-atom even at sites pointing toward the solvent. The nitroxide spin label exhibits flexibility and adopts various configurations depending on the surrounding residues.

Amino Acid Sequence↗

PubMed2006

Calculating the free energy of association of transmembrane helices.

A large number of experimental studies have been devoted to quantifying the interaction between transmembrane (TM) helices in detergent micelles and, more recently, in bilayers. Theoretical calculation of association free energy of TM helices would be useful for predicting the propensity of given sequences to oligomerize and for understanding the difference between association in micelles and in bilayers. In this article, the theoretical foundation for calculating the standard association free energy of TM helices is laid out and is applied to glycophorin A in both micelles and bilayers. The standard association free energy is decomposed into the effective energy, translational, rotational, and conformational entropy terms. The effective energy of association is obtained by molecular dynamics simulations in an implicit membrane model. The translational and rotational entropy of association is calculated from the probability distribution of the translational and rotational degrees of freedom obtained from the molecular dynamics simulations. The side-chain conformational entropy of association is estimated from the probability distribution obtained by rigid rotation of all side-chain dihedral angles. The calculated standard association free energy of glycophorin A in N-dodecylphosphocholine micelles is in good agreement with the experimental value. The translational entropy cost is larger, whereas the rotational entropy cost is smaller in bilayers than in micelles. The standard association free energy in 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayers is calculated to be approximately 1.3 kcal/mol more favorable than in N-dodecylphosphocholine micelles, consistent with available experimental data.

Cell Membrane↗

PubMed2006

Voltage-dependent energetics of alamethicin monomers in the membrane.

The implicit membrane model IMM1 is extended to include the effect of transmembrane potential and used to investigate the optimal membrane binding configurations and energies for alamethicin helices. In the absence of voltage, the lowest energy configuration is on the membrane surface with a tilt allowing the N terminus to be fully buried. Slightly higher in energy is an also tilted configuration with the N terminus deeper in the membrane and almost crossing the membrane. In 26A membranes and in the presence of 0.1V voltage, the TM orientation becomes lower in energy. This is consistent with the assumption that voltage induces a transition from the interfacial to the inserted (TM) orientation. This effect of voltage is smaller in thicker membranes. The results are compared to previous experimental and theoretical studies and the findings are discussed in relation to the mechanism of channel formation by alamethicin.

Alamethicin↗

PubMed2006

Energetics of the native and non-native states of the glycophorin transmembrane helix dimer.

Using an implicit membrane model (IMM1), we examine whether the structure of the transmembrane domain of Glycophorin A (GpA) could be predicted based on energetic considerations alone. The energetics of native GpA shows that van der Waals interactions make the largest contribution to stability. Although specific electrostatic interactions are stabilizing, the overall electrostatic contribution is close to zero. The GXXXG motif contributes significantly to stability, but residues outside this motif contribute almost twice as much. To generate non-native states a global conformational search was done on two segments of GpA: an 18-residue peptide (GpA74-91) that is embedded in the membrane and a 29-residue peptide (GpA70-98) that has additional polar residues flanking the transmembrane region. Simulated annealing was done on a large number of conformations generated from parallel, antiparallel, left- and right-handed starting structures by rotating each helix at 20 degrees intervals around its helical axis. Several crossing points along the helix dimer were considered. For 18-residue parallel topology, an ensemble of native-like structures was found at the lowest effective energy region; the effective energy is lowest for a right-handed structure with an RMSD of 1.0 A from the solid-state NMR structure with correct orientation of the helices. For the 29-residue peptide, the effective energies of several left-handed structures were lower than that of the native, right-handed structure. This could be due to deficiencies in modeling the interactions between charged sidechains and/or omission of the sidechain entropy contribution to the free energy. For 18-residue antiparallel topology, both IMM1 and a Generalized Born model give effective energies that are lower than that of the native structure. In contrast, the Poisson-Boltzmann solvation model gives lower effective energy for the parallel topology, largely because the electrostatic solvation energy is more favorable for the parallel structure. IMM1 seems to underestimate the solvation free energy advantage when the CO and NH dipoles just outside the membrane are parallel. This highlights the importance of electrostatic interactions even when these are not obvious by looking at the structures.

Amino Acid Sequence↗

PubMed2006

Calculations of pH-dependent binding of proteins to biological membranes.

Binding of proteins to membranes is often accompanied by titration of ionizable residues and is, therefore, dependent on pH. We present a theoretical treatment and computational approach for predicting absolute, pH-dependent membrane binding free energies. The standard free energy of binding, DeltaG, is defined as -RTln(P(b)/P(f)), where P(b) and P(f) are the amounts of bound and free protein. The apparent pK(a) of binding is the pH value at which P(b) and P(f) are equal. Proteins bind to the membrane in the pH range where DeltaG is negative. The components of the binding free energy are (a) the free energy cost of ionization state changes (DeltaG(ion)), (b) the effective energy of transfer from solvent to the membrane surface, (c) the translational/rotational entropy cost of binding, and (d) an ideal entropy term that depends on the relative volume of the bound and free state and therefore depends on lipid concentration. Calculation of the first term requires determination of pK(a) values in solvent and on the membrane surface. All energies required by the method are obtained from molecular dynamics trajectories on an implicit membrane (IMM1-GC). The method is tested on pentalysine and the helical peptide VEEKS, derived from the membrane-binding domain of phosphocholine cytidylyltransferase. The agreement between the measured and the calculated free energies of binding of pentalysine is good. The extent of membrane binding of VEEKS is, however, underestimated compared to experiment. Calculations of the interaction energy between two VEEKS helices on the membrane suggest that the discrepancy is mainly due to the neglect of protein-protein interactions on the membrane surface.

Aspartic Acid↗

PubMed2006

Thermodynamics of buried water clusters at a protein-ligand binding interface.

The structure of the complex of cyclophilin A (CypA) with cyclosporin A (CsA, 1) shows a cluster of four water molecules buried at the binding interface, which is rearranged when CsA is replaced by (5-hydroxynorvaline)-2-cyclosporin (2). The thermodynamic contributions of each bound water molecule in the two complexes are explored with the inhomogeneous fluid solvation theory and molecular dynamics simulations. Water (WTR) 133 in complex 1 contributes little to the binding affinity, while WTR6 and 7 in complex 2 play an essential role in mediating protein-ligand binding with a hydrogen bond network. The calculations reveal that the rearrangement of the water molecules contributes favorably to the binding affinity, even though one of them is displaced going from ligand 1 to 2. Another favorable contribution comes from the larger protein-ligand interactions of ligand 2. However, these favorable contributions are not sufficient to overcome the unfavorable desolvation free energy change and the conformational entropy of the hydroxylpropyl group of ligand 2 in the complex, leading to a lower binding affinity of ligand 2. These physical insights may be useful in the development of improved scoring functions for binding affinity prediction.

Binding Sites↗

PubMed2005

Implicit solvent simulations of DPC micelle formation.

The formation of micelles by dodecylphosphocholine (DPC) is modeled by treating the surfactants in atomic detail and the solvent implicitly, in the spirit of the EEF1 solvation model for proteins. The solvation parameters of the DPC atoms are carried over from those of similar atoms in proteins. A slight adjustment of the parameters for the headgroup was found necessary for obtaining an aggregation number consistent with experiment. Molecular dynamics simulations of 960 DPC molecules at different concentrations are used to obtain the aggregation number, the micelle size distribution, and the CMC. At 20 mM concentration we obtain an aggregation number of 53-56 and a CMC of 1.25 mM, values close to the experimental ones. At 100 mM the aggregation number increases to 90. Simulations of individual micelles of varying size show that the effective energy per surfactant molecule is initially a decreasing function of aggregation number but stabilizes at about 60 molecules. The van der Waals term and the desolvation of nonpolar groups contribute to micellization, whereas the desolvation of polar groups opposes it. From the difference between the effective energy and the free energy (calculated from the CMC), the translational and rotational entropy contributions to the free energy are estimated at about 7 kcal/mol per monomer. The micelles obtained here are more irregular than those obtained in explicit water simulations. This modeling approach allows the study of larger surfactant aggregates for longer times and the extraction of thermodynamic in addition to structural information.

Computer Simulation↗

PubMed2005

Implicit solvent simulations of peptide interactions with anionic lipid membranes.

A recently developed implicit membrane model (IMM1) is supplemented with a Gouy-Chapman term describing counterion-screened electrostatic interactions of a solute with negatively charged membrane lipids. The new model is tested on peptides that bind to anionic membranes. Pentalysine binds just outside the plane of negative charge, whereas Lys-Phe peptides insert their aromatic rings into the hydrophobic core. Melittin and magainin 2 bind more strongly to anionic than to neutral membranes and in both cases insert their hydrophobic residues into the hydrocarbon core. The third domain of Antennapedia homeodomain (penetratin) binds as an alpha-helix in the headgroup region. Cardiotoxin II binds strongly to anionic membranes but marginally to neutral ones. In all cases, the location and configuration of the peptides are consistent with experimental data, and the effective energy changes upon binding compare favorably with experimental binding free energies. The model opens the way to exploring the effect of membrane charge on the location, conformation, and dynamics of a large variety of biologically active peptides on membranes.

Animals↗

PubMed2005

The contribution of C alpha-H...O hydrogen bonds to membrane protein stability depends on the position of the amide.

Structural analyses of membrane proteins reveal a large number of C(alpha)-H...O contacts between transmembrane helices, presumed to be hydrogen bonds. Recent experiments produced conflicting results for the contribution of such hydrogen bonds to membrane protein stability. An FTIR study estimated an energy of -0.88 kcal/mol for the G79-C(alpha)-H...I76-O hydrogen bond in glycophorin A, whereas a mutagenesis study showed that the A51-C(alpha)-H...T24-O(gamma) hydrogen bond does not stabilize bacteriorhodopsin. Here, we reconcile these results using molecular mechanics calculations and an implicit membrane model (IMM1). With explicit hydrogen atoms, the potential energy of the G79-C(alpha)-H...I76-O interaction in GpA ranges from -0.54 to -0.9 kcal/mol and its contribution to stability (effective energy) from -0.49 to -0.83 kcal/mol, depending on the structural model used. The average values of these quantities in GpA-like motifs are similar. In bR, the corresponding numbers for the A51-C(alpha)-H...T24-O(gamma) interaction are +0.15 and +0.32 kcal/mol. The difference results from the different arrangement of the interacting groups and specifically the position of the acceptor with respect to the C(alpha) and N atoms. This conclusion likely applies to soluble proteins as well.

Alanine↗

PubMed2005

The effect of water displacement on binding thermodynamics: concanavalin A.

Interactions at the binding interface of biomolecular complexes are often mediated by ordered water molecules. In this work, we considered two concanavalin A-carbohydrate complexes. In the first, a water molecule is buried at the binding interface. In the second, this water molecule is displaced by a modification of the ligand (Clarke, C.; Woods, R. J.; Gluska, J.; Cooper, A.; Nutley, M. A.; Boons, G. J. J. Am. Chem. Soc. 2001, 123, 12238-12247). We computed the contribution of this water molecule to the thermodynamic properties using statistical mechanical formulas for the energy and entropy and molecular dynamics simulations. Other contributions to the binding affinity, including desolvation, entropy of conformational restriction, and interaction between the ligand and protein, were also computed. The thermodynamic consequences of displacement of the ordered water molecule by ligand modification were in qualitative agreement with experimental data. The free energy contribution of the water molecule (-17.2 kcal/mol; -19.2 enthalpic and +2 entropic) was nearly equivalent to the additional protein-ligand interactions in trimannoside 2 (-18.9 kcal/mol). The two structural ions interact more strongly with the water than with the hydroxyl of trimannoside 2, thus favoring trimannoside 1. The contributions from desolvation and conformational entropy are much smaller but significant, compared to the binding free energy difference. The picture that emerges is that the final outcome of water displacement is sensitive to the details of the binding site and cannot be predicted by simple empirical rules.

Binding Sites↗

PubMed2003

Effective energy function for proteins in lipid membranes.

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.

Amino Acids↗

PubMed2003

Thermodynamic contributions of the ordered water molecule in HIV-1 protease.

Binding between biomolecules is usually accompanied by the formation of direct interactions with displacement of water from the binding sites. In some cases, however, the interactions are mediated by ordered water molecules, whose effect on binding affinity and the other thermodynamic functions is unclear. In this work, we compute the contribution of one such water molecule, the strongly bound water molecule at the binding site of HIV-1 protease, to the thermodynamic properties using statistical mechanical formulas for the energy and entropy. The requisite correlation functions are obtained by molecular dynamics simulations. We find that the entropic penalty of ordering is large but is outweighed by the favorable water-protein interactions. We also find a large negative contribution from this water molecule to the heat capacity. This approach could be useful in rational drug design by estimating which bound water molecules would be most favorable to displace.

Entropy↗

PubMed2003

Potentials of mean force between ionizable amino acid side chains in water.

Potentials of mean force (PMF) between all possible ionizable amino acid side chain pairs in various protonation states were calculated using explicit solvent molecular dynamics simulations with umbrella sampling and the weighted histogram analysis method. The side chains were constrained in various orientations inside a spherical cluster of 200 water molecules. Beglov and Roux's Spherical Solvent Boundary Potential was used to account for the solvent outside this sphere. This approach was first validated by calculating PMFs between monatomic ions (K(+), Na(+), Cl(-)) and comparing them to results from the literature and results obtained using Ewald summation. The strongest interaction (-4.5 kcal/mol) was found for the coaxial Arg(+).Glu(-) pair. Many like-charged side chains display a remarkable lack of repulsion, and occasionally a weak attraction. The PMFs are compared to effective energy curves obtained with common implicit solvation models, namely Generalized Born (GB), EEF1, and uniform dielectric of 80. Overall, the EEF1 curves are too attractive, whereas the GB curves in most cases match the minima of the PMF curves quite well. The uniform dielectric model, despite some fortuitous successes, is grossly inadequate.

Amino Acids↗

PubMed2003

Investigation of pathways for the low-pH conformational transition in influenza hemagglutinin.

Targeted molecular dynamics simulations were used to study the conformational transition of influenza hemagglutinin (HA) from the native conformation to putative fusogenic or postfusion conformations populated at low pH. Three pathways for this conformational change were considered. Complete dissociation of the globular domains of HA was observed in one pathway, whereas smaller rearrangements were observed in the other two. The fusion peptides became exposed and moved toward the target membrane, although occasional movement toward the viral membrane was also observed. The effective energy profiles along the paths show multiple barriers. The final low-pH structures, which are consistent with available experimental data, are comparable in effective energy to native HA. As a control, the uncleaved precursor HA0 was also forced along the same pathway. In this case both the final energy and the energy barrier were much higher than in the cleaved protein. This study suggests that 1) as proposed, the native conformation is the global minimum energy conformation for the uncleaved precursor but a metastable state for cleaved HA; 2) the spring-loaded conformational change is energetically plausible in full-length HA; and 3) complete globular domain dissociation is not necessary for extension of the coiled coil and fusion peptide exposure, but the model with complete dissociation has lower energy.

Computer Simulation↗

PubMed2003

Thermodynamics of protein folding: a microscopic view.

Statistical thermodynamics provides a powerful theoretical framework for analyzing, understanding and predicting the conformational properties of biomolecules. The central quantity is the potential of mean force or effective energy as a function of conformation, which consists of the intramolecular energy and the solvation free energy. The intramolecular energy can be reasonably described by molecular mechanics-type functions. While the solvation free energy is more difficult to model, useful results can be obtained with simple approximations. Such functions have been used to estimate the intramolecular energy contribution to protein stability and obtain insights into the origin of thermodynamic functions of protein folding, such as the heat capacity. With reasonable decompositions of the various energy terms, one can obtain meaningful values for the contribution of one type of interaction or one chemical group to stability. Future developments will allow the thermodynamic characterization of ever more complex biological processes.

Algorithms↗

PubMed2002

Contributions to the binding free energy of ligands to avidin and streptavidin.

The free energy of binding of a ligand to a macromolecule is here formally decomposed into the (effective) energy of interaction, reorganization energy of the ligand and the macromolecule, conformational entropy change of the ligand and the macromolecule, and translational and rotational entropy loss of the ligand. Molecular dynamics simulations with implicit solvation are used to evaluate these contributions in the binding of biotin, biotin analogs, and two peptides to avidin and streptavidin. We find that the largest contribution opposing binding is the protein reorganization energy, which is calculated to be from 10 to 30 kcal/mol for the ligands considered here. The ligand reorganization energy is also significant for flexible ligands. The translational/rotational entropy is 4.5-6 kcal/mol at 1 M standard state and room temperature. The calculated binding free energies are in the correct range, but the large statistical uncertainty in the protein reorganization energy precludes precise predictions. For some complexes, the simulations show multiple binding modes, different from the one observed in the crystal structure. This finding is probably due to deficiencies in the force field but may also reflect considerable ligand flexibility.

Avidin↗

PubMed2002

Distance and exposure dependent effective dielectric function.

In an effort to develop a dielectric screening function for molecular dynamics simulations of biomolecules in implicit solvent, effective dielectric constants (D(eff)) for a large number of atom pairs in a typical globular protein are calculated by continuum electrostatics. Plots of D(eff) versus the intercharge distance are in general sigmoidal with the characteristics of the curve depending on the distance of the two charges from the dielectric boundary and, secondarily, on the extent to which the area surrounding each charge is occupied by solvent (the "exposure"). The D(eff) values were fitted to an empirical, analytical function of these parameters that reproduces the data reasonably well, although considerable scatter exists in the range of D(eff) from 30 to 80. In the system used for parameterization, the mean square deviation of electrostatic interaction energies with this function is 0.48 kcal/mol, compared to 1.45 for an analytical Generalized Born model and 1.52 for the linear distance-dependent dielectric model. When tested in other proteins of varying size and compactness, the present function is superior to both of the above models, except for a fully unfolded polypeptide chain, where the Generalized Born model is superior.

Journal Article↗