Biomedical subjects

Daniel R Ripoll

Publications and source records attributed to Daniel R Ripoll.

At least 19 recordsLinked to original sources

PubMed2006

A localized specific interaction alters the unfolding pathways of structural homologues.

Reductive unfolding studies of proteins are designed to provide information about intramolecular interactions that govern the formation (and stabilization) of the native state and about folding/unfolding pathways. By mutating Tyr92 to G, A, or L in the model protein, bovine pancreatic ribonuclease A, and through analysis of temperature factors and molecular dynamics simulations of the crystal structures of these mutants, it is demonstrated that the markedly different reductive unfolding rates and pathways of ribonuclease A and its structural homologue onconase can be attributed to a single, localized, ring-stacking interaction between Tyr92 and Pro93 in the bovine variant. The fortuitous location of this specific stabilizing interaction in a disulfide-bond-containing loop region of ribonuclease A results in the localized modulation of protein dynamics that, in turn, enhances the susceptibility of the disulfide bond to reduction leading to an alteration in the reductive unfolding behavior of the homologues. These results have important implications for folding studies involving topological determinants to obtain folding/unfolding rates and pathways, for protein structure-function prediction through fold recognition, and for predicting proteolytic cleavage sites.

Animals↗

PubMed2005

Coupling between conformation and proton binding in proteins.

Interest centers here on whether the use of a fixed charge distribution of a protein solute, or a treatment that considers proton-binding equilibria by solving the Poisson equation, is a better approach to discriminate native from non-native conformations of proteins. In this analysis of the charge distribution of 7 proteins, we estimate the solvation free energy contribution to the total free energy by exploring the 2(zeta) possible ionization states of the whole molecule, with zeta being the number of ionizable groups in the amino acid sequence, for every conformation in the ensembles of 7 proteins. As an additional consideration of the role of electrostatic interactions in determining the charge distribution of native folds, we carried out a comparison of alternative charge assignment models for the ionizable residues in a set of 21 native-like proteins. The results of this work indicate that (1) for 6 out of 7 proteins, estimation of solvent polarization based on the Generalized Born model with a fixed charge distribution provides the optimal trade-off between accuracy, with respect to the Poisson equation, and speed when compared to the accessible surface area model; for the seventh protein, consideration of all possible ionization states of the whole molecule appears to be crucial to discriminate the native from non-native conformations; (2) significant differences in the degree of ionization and hence the charge distribution for native folds are found between the different charge models examined; (3) the stability of the native state is determined by a delicate balance of all the energy components, and (4) conformational entropy, and hence the dynamics of folding, may play a crucial role for a successful ab initio protein folding prediction.

Ions↗

PubMed2005

On the orientation of the backbone dipoles in native folds.

The role of electrostatic interactions in determining the native fold of proteins has been investigated by analyzing the alignment of peptide bond dipole moments with the local electrostatic field generated by the rest of the molecule with and without solvent effects. This alignment was calculated for a set of 112 native proteins by using charges from a gas phase potential. Most of the peptide dipoles in this set of proteins are on average aligned with the electrostatic field. The dipole moments associated with alpha-helical conformations show the best alignment with the electrostatic field, followed by residues in beta-strand conformations. The dipole moments associated with other secondary structure elements are on average better aligned than in randomly generated conformations. The alignment of a dipole with the local electrostatic field depends on both the topology of the native fold and the charge distribution assumed for all of the residues. The influences of (i) solvent effects, (ii) different sets of charges, and (iii) the charge distribution assumed for the whole molecule were examined with a subset of 22 proteins each of which contains <30 ionizable groups. The results show that alternative charge distribution models lead to significant differences among the associated electrostatic fields, whereas the electrostatic field is less sensitive to the particular set of the adopted charges themselves (empirical conformational energy program for peptides or parameters for solvation energy).

Models, Molecular↗

PubMed2005

Molecular origin of anticooperativity in hydrophobic association.

The structuring of water molecules in the vicinity of nonpolar solutes is responsible for hydrophobic hydration and association thermodynamics in aqueous solutions. Here, we studied the potential of mean force (PMF) for the formation of a dimer and trimers of methane molecules in three specific configurations in explicit water to explain multibody effects in hydrophobic association on a molecular level. We analyzed the packing and orientation of water molecules in the vicinity of the solute to explain the effect of ordering of the water around nonpolar solutes on many-body interactions. Consistent with previous theoretical studies, we observed cooperativity, manifested as a reduction of the height of the desolvation barrier for the trimer in an isosceles triangle geometry, but for linear trimers, we observed only anticooperativity. A simple mechanistic picture of hydrophobic association is drawn. The free energy of hydrophobic association depends primarily on the difference in the number of water molecules in the first solvation shell of a cluster and that in the monomers of a cluster; this can be approximated by the molecular surface area. However, there are unfavorable electrostatic interactions between the water molecules from different parts of the solvation shell of a trimer because of their increased orientation induced by the nonpolar solute. These electrostatic interactions make an anticooperative contribution to the PMF, which is clearly manifested for the linear trimer where the multibody contribution due to changes in the molecular surface area is equal to zero. The information theory model of hydrophobic interactions of Hummer et al. also explains the anticooperativity of hydrophobic association of the linear trimers; however, it predicts anticooperativity with a qualitatively identical distance dependence for nonlinear trimers, which disagrees with the results of simulations.

Hydrophobic and Hydrophilic Interactions↗

PubMed2005

A surface loop of the potato leafroll virus coat protein is involved in virion assembly, systemic movement, and aphid transmission.

Two acidic domains of the Potato leafroll virus (PLRV) coat protein, separated by 55 amino acids and predicted to be adjacent surface features on the virion, were the focus of a mutational analysis. Eleven site-directed mutants were generated from a cloned infectious cDNA of PLRV and delivered to plants by Agrobacterium-mediated mechanical inoculation. Alanine substitutions of any of the three amino acids of the sequence EWH (amino acids 170 to 172) or of D177 disrupted the ability of the coat protein to assemble stable particles and the ability of the viral RNA to move systemically in four host plant species. Alanine substitution of E109, D173, or E176 reduced the accumulation of virus in agrobacterium-infiltrated tissues, the efficiency of systemic infection, and the efficiency of aphid transmission relative to wild-type virus, but the mutations did not affect virion stability. A structural model of the PLRV capsid predicted that the amino acids critical for virion assembly were located within a depression at the center of a coat protein trimer. The other amino acids that affected plant infection and/or aphid transmission were predicted to be located around the perimeter of the depression. PLRV virions play key roles in phloem-limited virus movement in plant hosts as well as in transport and persistence in the aphid vectors. These results identified amino acid residues in a surface-oriented loop of the coat protein that are critical for virus assembly and stability, systemic infection of plants, and movement of virus through aphid vectors.

Animals↗

PubMed2004

Fast and accurate computation of the 13C chemical shifts for an alanine-rich peptide.

The purpose of this work is, first, to present a fast and accurate technique to compute Boltzmann-averaged values of the quantum-chemical 13C chemical shifts for each amino acid in oligopeptides, demonstrated here by an application to the peptide Ac-XXAAAAAAAOO-NH2 (where X denotes diaminobutyric acid, A is alanine, and O is ornithine) [XAO] and, second, to discuss the capability of the 13Calpha and 13Cbeta chemical shifts to distinguish the PP(II) conformation from the alpha-helix and statistical-coil conformations. Use is made of a combination of approaches, summarized as follows: (1) derivation of an ensemble of conformations by using a molecular mechanics technique; (2) use of a clustering procedure to form families and build a reduced set of conformations consisting of the lowest-energy conformations of each family, and (3) computation of the 13C chemical shifts for the lowest-energy conformations of each family, using a quantum-chemical approach that treats a selected residue, or group of residues, with a 6-311+G(2d,p) locally-dense basis set, while the remaining residues in the sequence are treated with a 3-21G basis set. The whole procedure is quite accurate and speeds up the computation of the Boltzmann-averaged values of the 13C-chemical shifts by several orders of magnitude. The present application sheds some light on the conformational preference for alanine and non-alanine residues to occupy the PP(II) helical region of the Ramachandran map.

Alanine↗

PubMed2004

Comparative structural modeling and inference of conserved protein classes in Drosophila seminal fluid.

The constituents of seminal fluid are a complex mixture of proteins and other molecules, most of whose functions have yet to be determined and many of which are rapidly evolving. As a step in elucidating the roles of these proteins and exposing potential functional similarities hidden by their rapid evolution, we performed comparative structural modeling on 28 of 52 predicted seminal proteins produced in the Drosophila melanogaster male accessory gland. Each model was characterized by defining residues likely to be important for structure and function. Comparisons of known protein structures with predicted accessory gland proteins (Acps) revealed similarities undetectable by primary sequence alignments. The structures predict that Acps fall into several categories: regulators of proteolysis, lipid modifiers, immunity/protection, sperm-binding proteins, and peptide hormones. The comparative structural modeling approach indicates that major functional classes of mammalian and Drosophila seminal fluid proteins are conserved, despite differences in reproductive strategies. This is particularly striking in the face of the rapid protein sequence evolution that characterizes many reproductive proteins, including Drosophila and mammalian seminal proteins.

Amino Acid Sequence↗

PubMed2004

Folding of the villin headpiece subdomain from random structures. Analysis of the charge distribution as a function of pH.

The structure of the 36 residue villin headpiece subdomain is investigated with the electrostatically driven Monte Carlo method. The ECEPP/3 (Empirical Conformational Energy Program for Peptides) force field, plus two different continuum solvation models, were used to describe the conformational energy of the chain with both blocked and unblocked N and C termini. A statistical analysis of an ensemble of ab initio generated conformations was carried out, based on a comparison with a set of ten native-like structures derived from published experimental data, by using rigid geometry and NMR-derived constraints obtained at pH 3.7. The ten native-like structures satisfy the NMR-derived constraints. The whole ensemble of conformations of the terminally unblocked villin headpiece sub-domain, generated by using ECEPP/3 with a continuum solvation model, were subsequently evaluated at pH 3.7 with a potential function that includes ECEPP/3 combined with a fast multigrid boundary element method. At pH 3.7, the lowest-energy conformation found during the conformational search satisfies approximately 70% of both the distance and the dihedral-angle constraints, and possesses the characteristic packing of three phenylalanine residues that constitute the main part of the hydrophobic core of the molecule. On the other hand, computations at pH 3.7 and pH 7.0 for the ten native-like structures satisfying the NMR-derived constraints indicate a substantial change in the charge distribution for each type of amino acid residue with the change in pH. The results of this study provide a basis to understand the effect of the interactions, such as hydrophobicity, charge-charge interaction and solvent polarization, on the stability of this small alpha-helical protein.

Carrier Proteins↗

PubMed2004

Dissimilarity in the reductive unfolding pathways of two ribonuclease homologues.

Using DTT(red) as the reducing agent, the kinetics of the reductive unfolding of onconase, a frog ribonuclease, has been examined. An intermediate containing three disulfides, Ir, that is formed rapidly in the reductive pathway, is more resistant to further reduction than the parent molecule, indicating that the remaining disulfides in onconase are less accessible to DTT(red). Disulfide-bond mapping of Ir indicated that it is a single species lacking the (30-75) disulfide bond. The reductive unfolding pattern of onconase is consistent with an analysis of the exposed surface area of the cysteine sulfur atoms in the (30-75) disulfide bond, which reveals that these atoms are about four- and sevenfold, respectively, more exposed than those in the next two maximally exposed disulfides. By contrast, in the reductive unfolding of the homologue, RNase A, there are two intermediates, arising from the reduction of the (40-95) and (65-72) disulfide bonds, which takes place in parallel, and on a much longer time-scale, compared to the initial reduction of onconase; this behavior is consistent with the almost equally exposed surface areas of the cysteine sulfur atoms that form the (40-95) and (65-72) disulfide bonds in RNase A and the fourfold more exposed cysteine sulfur atoms of the (30-75) disulfide bond in onconase. Analysis and in silico mutation of the residues around the (40-95) disulfide bond in RNase A, which is analogous to the (30-75) disulfide bond of onconase, reveal that the side-chain of tyrosine 92 of RNase A, a highly conserved residue among mammalian pancreatic ribonucleases, lies atop the (40-95) disulfide bond, resulting in a shielding of the corresponding sulfur atoms from the solvent; such burial of the (30-75) sulfur atoms is absent from onconase, due to the replacement of Tyr92 by Arg73, which is situated away from the (30-75) disulfide bond and into the solvent, resulting in the large exposed surface-area of the cysteine sulfur atoms forming this bond. Removal of Tyr92 from RNase A resulted in the relatively rapid reduction of the mutant to form a single intermediate (des [40-95] Y92A), i.e. it resulted in an onconase-like reductive unfolding behavior. The reduction of the P93A mutant of RNase A proceeds through a single intermediate, the des [40-95] P93A species, as in onconase. Although mutation of Pro93 to Ala does not increase the exposed surface area of the (40-95) cysteine sulfur atoms, structural analysis of the mutant reveals that there is greater flexibility in the (40-95) disulfide bond compared to the (65-72) disulfide bond that may make the (40-95) disulfide bond much easier to expose, consistent with the reductive unfolding pathway and kinetics of P93A. Mutation of Tyr92 to Phe92 in RNase A has no effect on its reductive unfolding pathway, suggesting that the hydrogen bond between the hydroxyl group of Tyr92 and the carbonyl group of Lys37 has no impact on the local unfolding free energy required to expose the (40-95) disulfide bond. Thus, these data shed light on the differences between the reductive unfolding pathways of the two homologous proteins and provide a structural basis for the origin of this difference.

Animals↗

PubMed2004

Polyproline II helix conformation in a proline-rich environment: a theoretical study.

Interest centers here on whether a polyproline II helix can propagate through adjacent non-proline residues, and on shedding light on recent experimental observations suggesting the presence of significant PP(II) structure in a short alanine-based peptide with no proline in the sequence. For this purpose, we explored the formation of polyproline II helices in proline-rich peptides with the sequences Ac-(Pro)(3)-X-(Pro)(3)-Gly-Tyr-NH(2), with X = Pro (PPP), Ala (PAP), Gln (PQP), Gly (PGP), and Val (PVP), and Ac-(Pro)(3)-Ala-Ala-(Pro)(3)-Gly-Tyr-NH(2) (PAAP), by using a theoretical approach that includes a solvent effect as well as cis <--> trans isomerization of the peptide groups and puckering conformations of the pyrrolidine ring of the proline residues. Since (13)C chemical shifts have proven to be useful for identifying secondary-structure preferences in proteins and peptides, and because values of the dihedral angles (phi,psi) are the main determinants of their magnitudes, we have, therefore, computed the Boltzmann-averaged (13)C chemical shifts for the guest residues in the PXP peptide (X = Pro, Ala, Gln, Gly, and Val) with a combination of approaches, involving molecular mechanics, statistical mechanics, and quantum mechanics. In addition, an improved procedure was used to carry out the conformational searches and to compute the solvent polarization effects faster and more accurately than in previous work. The current theoretical work and additional experimental evidence show that, in short proline-rich peptides, alanine decreases the polyproline II helix content. In particular, the theoretical evidence accumulated in this work calls into question the proposal that alanine has a strong preference to adopt conformations in the polyproline II region of the Ramachandran map.

Amino Acid Sequence↗

PubMed2003

Specificity determinants and diversification of the Brassica self-incompatibility pollen ligand.

Self-incompatibility in crucifers is effected by allele-specific interactions between the highly polymorphic stigmatic S locus receptor kinase (SRK) and its pollen ligand, the S locus cysteine-rich protein (SCR). Here we show that specificity in SCR function is determined by four contiguous amino acids in one variant, indicating that the minimum sequence requirement for gaining a new specificity can be low. We also provide evidence for an extraordinarily high degree of evolutionary flexibility in SCR, whereby SCR can tolerate extensive amino acid changes within the limits of maintaining the same predicted overall structure. This remarkable adaptability suggests a hypothesis for generation of new self-incompatibility specificities by gradual modification of SRK-SCR affinities and, more generally, for functional specialization within families of homologous ligands and receptors.

Amino Acid Sequence↗

PubMed2003

Atomically detailed folding simulation of the B domain of staphylococcal protein A from random structures.

The conformational space of the 10-55 fragment of the B-domain of staphylococcal protein A has been investigated by using the electrostatically driven Monte Carlo (EDMC) method. The ECEPP/3 (empirical conformational energy program for peptides) force-field plus two different continuum solvation models, namely SRFOPT (Solvent Radii Fixed with atomic solvation parameters OPTimized) and OONS (Ooi, Oobatake, Némethy, and Scheraga solvation model), were used to describe the conformational energy of the chain. After an exhaustive search, starting from two different random conformations, three of four runs led to native-like conformations. Boltzmann-averaged root-mean-square deviations (RMSD) for all of the backbone heavy atoms with respect to the native structure of 3.35 A and 4.54 A were obtained with SRFOPT and OONS, respectively. These results show that the protein-folding problem can be solved at the atomic detail level by an ab initio procedure, starting from random conformations, with no knowledge except the amino acid sequence. To our knowledge, the results reported here correspond to the largest protein ever folded from a random conformation by an initial-value formulation with a full atomic potential, without resort to knowledge-based information.

Hydrogen Bonding↗

PubMed2003

Clp protease complexes from photosynthetic and non-photosynthetic plastids and mitochondria of plants, their predicted three-dimensional structures, and functional implications.

Tetradecameric Clp protease core complexes in non-photosynthetic plastids of roots, flower petals, and in chloroplasts of leaves of Arabidopsis thaliana were purified based on native mass and isoelectric point and identified by mass spectrometry. The stoichiometry between the subunits was determined. The protease complex consisted of one to three copies of five different serine-type protease Clp proteins (ClpP1,3-6) and four non-proteolytic ClpR proteins (ClpR1-4). Three-dimensional homology modeling showed that the ClpP/R proteins fit well together in a tetradecameric complex and also indicated unique contributions for each protein. Lateral exit gates for proteolysis products are proposed. In addition, ClpS1,2, unique to land plants, tightly interacted with this core complex, with one copy of each per complex. The three-dimensional modeling show that they do fit well on the axial sites of the ClpPR cores. In contrast to plastids, plant mitochondria contained a single approximately 320-kDa homo-tetradecameric ClpP2 complex, without association of ClpR or ClpS proteins. It is surprising that the Clp core composition appears identical in all three plastid types, despite the remarkable differences in plastid proteome composition. This suggests that regulation of plastid proteolysis by the Clp machinery is not through differential regulation of ClpP/R/S gene expression, but rather through substrate recognition mechanisms and regulated interaction of chaperone-like molecules (ClpS1,2 and others) to the ClpP/R core.

Adenosine Triphosphatases↗

PubMed2003

Packing helices in proteins by global optimization of a potential energy function.

An efficient method has been developed for packing alpha-helices in proteins. It treats alpha-helices as rigid bodies and uses a simplified Lennard-Jones potential with Miyazawa-Jernigan contact-energy parameters to describe the interactions between the alpha-helical elements in this coarse-grained system. Global conformational searches to generate packing arrangements rapidly are carried out with a Monte Carlo-with-minimization type of approach. The results for 42 proteins show that the approach reproduces native-like folds of alpha-helical proteins as low-energy local minima of this highly simplified potential function.

Protein Structure, Secondary↗

PubMed2003

Molecular simulation study of cooperativity in hydrophobic association: clusters of four hydrophobic particles.

The multibody contribution to the potential of mean force (PMF) of hydrophobic association of four methane molecules in water was investigated by means of umbrella-sampling molecular dynamics. Two systems were considered: (i). a trigonal pyramid with three methane molecules at contact distance forming a fixed base, the fourth molecule being placed on the top with variable distance from the base; and (ii). a regular uniformly expanding tetrahedron. Methane-methane distances as far as 12.5 A, i.e. beyond the second solvent-separated minimum of the PMF, were considered to address the baseline problem. In contrast to the small effect in the three-body case studied previously (Protein Sci 9 (2000) 1235), the multibody contribution was found to amount to approximately 0.2 kcal/mol per methane-methane pair, or approximately 25% of the depth of the contact minimum in the PMF. The main effect of the multibody contribution to the PMF is a reduction of the height of the barrier between the contact and solvent separated minima and a narrowing of the region of its maximum, while the region of the contact minimum is affected only weakly. The reduction of the barrier is due to four-body contributions. The cooperative contributions to the PMF agree very well with those computed from the molecular surface of the systems under consideration, which further supports earlier observations that the molecular surface can be used with good accuracy to describe the energetics of hydrophobic association.

Hydrophobic and Hydrophilic Interactions↗

PubMed2003

Unblocked statistical-coil tetrapeptides in aqueous solution: quantum-chemical computation of the carbon-13 NMR chemical shifts.

We recently reported a theoretical characterization of representative ensembles of statistical-coil conformations for tetrapeptides with unblocked termini in aqueous solution, at pH 7. The results showed good agreement between the computed Boltzmann-averaged and experimentally-determined values for both the vicinal coupling constants (3)J(NHalpha) and the alpha-proton chemical shifts. Here, we carry out a cluster analysis of the ensembles of conformations generated in that study, and use them to compute the Boltzmann-averaged values of the quantum-chemical (13)C chemical shifts for different amino acids in the unblocked tetrapeptides GGXA (where X stands for Phe, Arg, His, Glu, Ile, Lys, Gln, Tyr, Leu, Thr, Ala, Gly and Val). The values of the (13)C chemical shifts in these thirteen amino acids (for which experimental data are available) were computed by using Density Functional Theory with a 6-311+G(2d,p) basis set. Good agreement is found in terms of both the correlation coefficient (R) and standard deviations of the difference between the computed Bolztmann-averaged and the NMR-determined values for the (13)C chemical shifts. These results suggest that it may be possible to build a reliable theoretically-derived database of chemical shifts for statistical-coil residues. The results of the current study contribute to our understanding of the relations between chemical shifts, dihedral angles and vicinal coupling constants, (3)J(NHalpha). In addition, they can shed light as to how the statistical-coil conformation is related to the conformational preference of more structured states, such as the alpha-helical conformation.

Amino Acid Sequence↗

PubMed2002

Helix-coil transitions re-visited.

The thermally-induced helix-coil transition in polyamino acids is a good model for determining the helix-forming propensities of amino acids but not for the two-state folding/unfolding transition in globular proteins. The equilibrium and kinetic treatments of the helix-coil transition are summarized here together with a description of applications to various types of homopolymers and copolymers. Attention is then focused on the helix-coil transition in poly-L-alanine as an example of a non-polar polyamino acid. To render such a non-polar polymer water soluble, it is necessary to introduce polar amino acids such as lysines, but care must be taken as to the location of such polar residues. If they are attached as end groups, as in a triblock copolymer, they do not perturb the helix-forming tendency of the central poly-L-alanine block significantly, but if they are introduced within the sequence of alanine residues, then the hydration properties of the lysines dominate the behavior of the resulting copolymer, thereby leading to erroneous values of the parameters characterizing the helix-forming tendency of the alanines. Neutral but polar residues, such as glutamines, also exhibit hydration-dominating properties but less so than charged lysines. Some details of the calculations for an alanine/glutamine copolymer are presented here. It is concluded that random copolymers based on a neutral water-soluble host provide reliable information about the helix-forming tendencies of amino acid residues that are introduced as guests among such neutral host residues.

Amino Acid Sequence↗

PubMed2002

Evolution of physics-based methodology for exploring the conformational energy landscape of proteins.

The evolution of our physics-based computational methods for determining protein conformation without the introduction of secondary-structure predictions, homology modeling, threading, or fragment coupling is described. Initial use of a hard-sphere potential captured much of the structural properties of polypeptide chains, and subsequent more refined force fields, together with efficient methods of global optimization provide indications that progress is being made toward an understanding of the interresidue interactions that underlie protein folding.

Computational Biology↗