Search PubMed⌕ Search

Biomedical subjects

A E García

Publications and source records attributed to A E García.

At least 19 recordsLinked to original sources
PubMed2002

Structure of Met-enkephalin in explicit aqueous solution using replica exchange molecular dynamics.

Replica exchange molecular dynamics (MD) simulations of Met-enkephalin in explicit solvent reveal helical and nonhelical structures. Four predominant structures of Met-enkephalin are sampled with comparable probabilities (two helical and two nonhelical). The energy barriers between these configurations are low, suggesting that Met-enkephalin switches easily between configurations. This is consistent with the requirement that Met-enkephalin be sufficiently flexible to bind to several different receptors. Replica exchange simulations of 32 ns are shown to sample approximately five times more configurational space than constant temperature MD simulations of the same duration. The energy landscape for the replica exchange simulation is presented. A detailed study of replica trajectories demonstrates that the significant increases in temperature provided by the replica exchange technique enable transitions from nonhelical to helical structures that would otherwise be prevented by kinetic trapping. Met-enkephalin (Type Entrez Proteins; Value A61445; Service Entrez Proteins).

Algorithms↗
PubMed2001

Molecular dynamics simulations of pressure effects on hydrophobic interactions.

We report results on the pressure effects on hydrophobic interactions obtained from molecular dynamics simulations of aqueous solutions of methanes in water. A wide range of pressures that is relevant to pressure denaturation of proteins is investigated. The characteristic features of water-mediated interactions between hydrophobic solutes are found to be pressure-dependent. In particular, with increasing pressure we find that (1) the solvent-separated configurations in the solute-solute potential of mean force (PMF) are stabilized with respect to the contact configurations; (2) the desolvation barrier increases monotonically with respect to both contact and solvent-separated configurations; (3) the locations of the minima and the barrier move toward shorter separations; and (4) pressure effects are considerably amplified for larger hydrophobic solutes. Together, these observations lend strong support to the picture of the pressure denaturation process proposed previously by Hummer et al. (Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 1552): with increasing pressure, the transfer of water into protein interior becomes key to the pressure denaturation process, leading to the dissociation of close hydrophobic contacts and subsequent swelling of the hydrophobic protein interior through insertions of water molecules. The pressure dependence of the PMF between larger hydrophobic solutes shows that pressure effects on the interaction between hydrophobic amino acids may be considerably amplified compared to those on the methane-methane PMF.

Computer Simulation↗
PubMed2001

Amino acid substitution at position 99 affects the rate of CRP subunit exchange.

We investigated the characteristics of CRP having amino acid substitutions at position 99. Analysis of amino acid residue proximity to cAMP in molecular dynamics (MD) simulations of the CRP:(cAMP)(2) complex [García, A. E., and Harman, J. G. (1996) Protein Sci. 5, 62-71] showed repositioning of tyrosine 99 (Y99) to interact with the equatorial exocyclic oxygen atom of cAMP. To test the role of Y99 in cAMP-mediated CRP activation, Y99 was substituted with alanine (A) or phenylalanine (F). Cells that contained the WT or mutant forms of CRP induced beta-galactosidase in the presence of cAMP. Purified WT, Y99A, and Y99F CRP showed only a 3- to 4-fold difference in cAMP affinity. There were no apparent differences between the three forms of CRP in cAMP binding cooperativity, in CRP:(cAMP)(1) complex binding to lacP DNA, in the formation of CRP:cAMP:RNAP complexes at lacP, or in CRP efficacy in mediating lacP activity in vitro. The apo-form of Y99A CRP was more sensitive to protease than the apo-form of either WT CRP or Y99F CRP. Whereas the WT or Y99F CRP:(cAMP)(1) complexes were cleaved by protease at hinge-region peptide bonds, the Y99A CRP:(cAMP)(1) complex was cleaved at peptide bonds located at the subunit interface. The rates of subunit exchange for Y99A CRP, both in the apo-form and in a 1:1 complex with cAMP, were significantly greater than that measured for WT CRP. The results of this study show that tyrosine 99 contributes significant structural stability to the CRP dimer, specifically in stabilizing subunit association.

Amino Acid Sequence↗
PubMed2001

Exploring the energy landscape of a beta hairpin in explicit solvent.

We studied the energy landscape of the peptide Ace-GEWTYDDATKTFTVTE-Nme, taken from the C-terminal fragment (41-56) of protein G, in explicit aqueous solution by a highly parallel replica-exchange approach that combines molecular dynamics trajectories with a temperature exchange Monte Carlo process. The combined trajectories in T and configurational space allow a replica to overcome a free energy barrier present at one temperature by increasing T, changing configurations, and cooling in a self-regulated manner, thus allowing sampling of broad regions of configurational space in short (nanoseconds) time scales. The free energy landscape of this system over a wide range of temperatures shows that the system preferentially adopts a beta hairpin structure. However, the peptide also samples other stable ensembles where the peptide adopts helices and helix-turn-helix states, among others. The helical states become increasingly stable at low temperatures, but are slightly less stable than the beta turn ensemble. The energy landscape is rugged at low T, where substates are separated by large energy barriers. These barriers disappear at higher T (approximately 330 K), where the system preferentially adopts a "molten globule" state with structures similar to the beta hairpin.

Amino Acid Sequence↗
PubMed2001

Helix nucleation kinetics from molecular simulations in explicit solvent.

We study the reversible folding/unfolding of short Ala and Gly-based peptides by molecular dynamics simulations of all-atom models in explicit water solvent. A kinetic analysis shows that the formation of a first alpha-helical turn occurs within 0.1-1 ns, in agreement with the analyses of laser temperature jump experiments. The unfolding times exhibit Arrhenius temperature dependence. For a rapidly nucleating all-Ala peptide, the helix nucleation time depends only weakly on temperature. For a peptide with enthalpically competing turn-like structures, helix nucleation exhibits an Arrhenius temperature dependence, corresponding to the unfolding of enthalpic traps in the coil ensemble. An analysis of structures in a "transition-state ensemble" shows that helix-to-coil transitions occur predominantly through breaking of hydrogen bonds at the helix ends, particularly at the C-terminus. The temperature dependence of the transition-state ensemble and the corresponding folding/unfolding pathways illustrate that folding mechanisms can change with temperature, possibly complicating the interpretation of high-temperature unfolding simulations. The timescale of helix formation is an essential factor in molecular models of protein folding. The rapid helix nucleation observed here suggests that transient helices form early in the folding event.

Alanine↗
PubMed2001

Functional dynamics of the hydrophobic cleft in the N-domain of calmodulin.

Molecular dynamics studies of the N-domain (amino acids 1-77; CaM(1-77)) of Ca2+-loaded calmodulin (CaM) show that a solvent exposed hydrophobic cleft in the crystal structure of CaM exhibits transitions from an exposed (open) to a buried (closed) state over a time scale of nanoseconds. As a consequence of burying the hydrophobic cleft, the R(g) of the protein is reduced by 1.5 A. Based on this prediction, x-ray scattering experiments were conducted on this domain over a range of concentrations. Models built from the scattering data show that the R(g) and general shape is consistent with the simulation studies of CaM(1-77). Based on these observations we postulate a model in which the conformation of CaM fluctuates between two different states that expose and bury this hydrophobic cleft. In aqueous solution the closed state dominates the population, while in the presence of peptides, the open state dominates. This inherent flexibility of CaM may be the key to its versatility in recognizing structurally distinct peptide sequences. This model conflicts with the currently accepted hypothesis based on observations in the crystal structure, where upon Ca2+ binding the hydrophobic cleft is exposed to solvent. We postulate that crystal packing forces stabilize the protein conformation toward the open configuration.

Biophysical Phenomena↗
PubMed2000

Conformational diffusion and helix formation kinetics.

The time, temperature, and sequence dependences of helix formation kinetics of fully atomistic peptide models in explicit solvent are described quantitatively by a diffusive search within the coil state with barrierless transitions into the helical state. Conformational diffusion leads to nonexponential kinetics and jump-width dependences in temperature jump experiments.

Kinetics↗
PubMed2000

Water penetration and escape in proteins.

The kinetics of water penetration and escape in cytochrome c (cyt c) is studied by molecular dynamics (MD) simulations at various temperatures. Water molecules that penetrate the protein interior during the course of an MD simulation are identified by monitoring the number of water molecules in the first coordination shell (within 3.5 A) of each water molecule in the system. Water molecules in the interior of cyt c have 0-3 water molecules in their first hydration shell and this coordination number persists for extended periods of time. At T = 300 K we identify over 200 events in which water molecules penetrate the protein and reside inside for at least 5 picoseconds (ps) within a 1.5 nanoseconds (ns) time period. Twenty-seven (27) water molecules reside for at least 300 ps, 17 water molecules reside in the protein interior for times longer than 500 ps, and two interior water molecules do not escape; at T = 360 K one water molecule does not escape; at 430 K all water molecules exchange. Some of the internal water molecules show mean square displacements (MSD) of 1 A2 characteristic of structural waters. Others show MSD as large as 12 A2, suggesting that some of these water molecules occupy transient cavities and diffuse extensively within the protein. Motions of protein-bound water molecules are rotationally hindred, but show large librations. Analysis of the kinetics of water escape in terms of a survival time correlation function shows a power law behavior in time that can be interpreted in terms of a broad distribution of energy barriers, relative to kappa BT, for water exchange. At T = 300 K estimates of the roughness of the activation energy distribution is 4-10 kJ/mol (2-4 kappa BT). Activation enthalpies for water escape are 6-23 kJ/mol. The difference in activation entropies between fast exchanging (0.01 ns) and slow exchanging (0.1-1 ns) water molecules is -27 J/K/mol. Dunitz (Science 1997;264:670.) has estimated the maximum entropy loss of a water molecule due to binding to be 28 J/K/mol. Therefore, our results suggest that the entropy of interior water molecules is similar to entropy of bulk water.

Animals↗
PubMed1999

Pressure-induced protein-folding/unfolding kinetics.

We use an off-lattice minimalist model to describe the effects of pressure in slowing down the folding/unfolding kinetics of proteins when subjected to increasingly larger pressures. The potential energy function used to describe the interactions between beads in the model includes the effects of pressure on the pairwise interaction of hydrophobic groups in water. We show that pressure affects the participation of contacts in the transition state. More significantly, pressure exponentially decreases the chain reconfigurational diffusion coefficient. These results are consistent with experimental results on the kinetics of pressure-denaturation of staphylococcal nuclease.

Computer Simulation↗
PubMed1999

Conformational dynamics of cytochrome c: correlation to hydrogen exchange.

We study the dynamical fluctuations of horse heart cytochrome c by molecular dynamics (MD) simulations in aqueous solution, at four temperatures: 300 K, 360 K, 430 K, and 550 K. Each simulation covers a production time of at least 1.5 nanoseconds (ns). The conformational dynamics of the system is analyzed in terms of collective motions that involve the whole protein, and local motions that involve the formation and breaking of intramolecular hydrogen bonds. The character of the MD trajectories can be described within the framework of rugged energy landscape dynamics. The MD trajectories sample multiple conformational minima, with basins in protein conformational space being sampled for a few hundred picoseconds. The trajectories of the system in configurational space can be described in terms of diffusion of a particle in real space with a waiting time distribution due to partial trapping in shallow minima. As a consequence of the hierarchical nature of the dynamics, the mean square displacement autocorrelation function, <|x(t) - x(0)|2>, exhibits a power law dependence on time, with an exponent of around 0.5 for times shorter than 100 ps, and an exponent of 1.75 for longer times. This power law behavior indicates that the system exhibits suppressed diffusion (sub-diffusion) in sampling of configurational space at time scales shorter than 100 ps, and enhanced (super-diffusion) at longer time scales. The multi-basin feature of the trajectories is present at all temperatures simulated. Structural changes associated with inter-basin displacements correspond to collective motions of the Omega loops and coiled regions and relative motions of the alpha-helices as rigid bodies. Similar motions may be involved in experimentally observed amide hydrogen exchange. However, some groups showing large correlated motions do not expose the amino hydrogens to the solvent. We show that large fluctuations are not necessarily correlated to hydrogen exchange. For example, regions of the proteins forming alpha helices and turns show significant fluctuations, but as rigid bodies, and the hydrogen bonds involved in the formation of these structures do not break in proportion to these fluctuations. Proteins 1999;36:175-191. Published 1999 Wiley-Liss, Inc.

Amides↗
PubMed1998

Folding funnels and frustration in off-lattice minimalist protein landscapes.

A full quantitative understanding of the protein folding problem is now becoming possible with the help of the energy landscape theory and the protein folding funnel concept. Good folding sequences have a landscape that resembles a rough funnel where the energy bias towards the native state is larger than its ruggedness. Such a landscape leads not only to fast folding and stable native conformations but, more importantly, to sequences that are robust to variations in the protein environment and to sequence mutations. In this paper, an off-lattice model of sequences that fold into a beta-barrel native structure is used to describe a framework that can quantitatively distinguish good and bad folders. The two sequences analyzed have the same native structure, but one of them is minimally frustrated whereas the other one exhibits a high degree of frustration.

Animals↗
PubMed1998

The pressure dependence of hydrophobic interactions is consistent with the observed pressure denaturation of proteins.

Proteins can be denatured by pressures of a few hundred MPa. This finding apparently contradicts the most widely used model of protein stability, where the formation of a hydrophobic core drives protein folding. The pressure denaturation puzzle is resolved by focusing on the pressure-dependent transfer of water into the protein interior, in contrast to the transfer of nonpolar residues into water, the approach commonly taken in models of protein unfolding. Pressure denaturation of proteins can then be explained by the pressure destabilization of hydrophobic aggregates by using an information theory model of hydrophobic interactions. Pressure-denatured proteins, unlike heat-denatured proteins, retain a compact structure with water molecules penetrating their core. Activation volumes for hydrophobic contributions to protein folding and unfolding kinetics are positive. Clathrate hydrates are predicted to form by virtually the same mechanism that drives pressure denaturation of proteins.

Hydrostatic Pressure↗
PubMed1997

Hydration of an alpha-helical peptide: comparison of theory and molecular dynamics simulation.

We present a statistical mechanical description of biomolecular hydration that accurately describes the hydrophobic and hydrophilic hydration of a model alpha-helical peptide. The local density of water molecules around a biomolecule is obtained by means of a potential-of-mean-force (PMF) expansion in terms of pair- and triplet-correlation functions of bulk water and dilute solutions of nonpolar atoms. The accuracy of the method is verified by comparing PMF results with the local density and site-site correlation functions obtained by molecular dynamics simulations of a model alpha-helix in solution. The PMF approach quantitatively reproduces all features of the peptide hydration determined from the molecular dynamics simulation. Regions of hydrophobic hydration near the C alpha and C beta atoms along the helix are well reproduced. The hydration of exposed polar groups at the N- and C-termini of the helix are also well described by the theory. A detailed comparison of the local hydration by means of site-site radial distribution functions evaluated with the PMF theory shows agreement with the molecular dynamics simulations. The formulation of this theory is general and can be applied to any biomolecular system. The accuracy, speed of computation, and local character of this theory make it especially suitable for studying large biomolecular systems.

Computer Simulation↗
PubMed1997

Variations on a theme by Debye and Waller: from simple crystals to proteins.

Debye and Waller showed how to adjust scattering intensities in diffraction experiments for harmonic motions of atoms about an average, static reference configuration. However, many motions, particularly in biological molecules as compared to simple crystals, are far from harmonic. We show how, using a variety of simple anharmonic, multiconformational models, it is possible to construct a variety of Generalized Debye-Waller Factors, and understand their meaning. A central result for these cases is that, in principle, intensity factors cannot be obtained from true total mean square displacements of the atoms. We make the distinction between the intensity factors for unimodal quasiharmonic motions and those for the anharmonic, multimodal (valley hopping) motions. Only the former affect the conventional B factors.

Models, Chemical↗
PubMed1996

An information theory model of hydrophobic interactions.

A molecular model of poorly understood hydrophobic effects is heuristically developed using the methods of information theory. Because primitive hydrophobic effects can be tied to the probability of observing a molecular-sized cavity in the solvent, the probability distribution of the number of solvent centers in a cavity volume is modeled on the basis of the two moments available from the density and radial distribution of oxygen atoms in liquid water. The modeled distribution then yields the probability that no solvent centers are found in the cavity volume. This model is shown to account quantitatively for the central hydrophobic phenomena of cavity formation and association of inert gas solutes. The connection of information theory to statistical thermodynamics provides a basis for clarification of hydrophobic effects. The simplicity and flexibility of the approach suggest that it should permit applications to conformational equilibria of nonpolar solutes and hydrophobic residues in biopolymers.

Journal Article↗
PubMed1996

Simulations of CRP:(cAMP)2 in noncrystalline environments show a subunit transition from the open to the closed conformation.

The CRP:cAMP complex functions as a transcription factor that facilitates RNA polymerase recognition of several bacterial promoters. Detailed crystal structure information is available for CRP:(cAMP)2 and for CRP:(cAMP)2 complexed with DNA. In the crystalline environment, CRP:(cAMP)2 subunits are asymmetrically related; one subunit has a closed conformation and the other has an open conformation. The CRP:(cAMP)2 complexed with DNA shows both subunits in a closed conformation. We have studied the molecular dynamics of CRP:(cAMP)2 in noncrystalline environments. CRP:(cAMP)2 was simulated for 625 ps in vacuo and for 140 ps in solution. The crystal structure of CRP:(cAMP)2 in the absence of DNA was used as the initial conformation. Molecule optimal dynamic coordinates (MODCs) (García A, 1992, Phys Rev Lett 68:2696) were used to analyze protein conformations sampled during the course of the simulations. Two MODCs define a transition of the open subunit to a closed subunit conformation during the first 125 ps of simulation in vacuo; the resulting subunit conformation is similar to that observed in CRP:(cAMP)2:DNA crystals. Simulation of CRP:(cAMP)2 in solution showed that a transition from the open to the closed state also occurs when water is explicitly included in the calculations. These calculations suggest that the asymmetric conformation of CRP:(cAMP)2 is stabilized by crystal lattice interactions. The predicted solution conformation is more symmetric, with both subunits in a closed conformation.

Adenosine Monophosphate↗
PubMed1996

Three-dimensional quantitative structure-activity relationships of steroid aromatase inhibitors.

Inhibition of aromatase, a cytochrome P450 that converts androgens to estrogens, is relevant in the therapeutic control of breast cancer. We investigate this inhibition using a three-dimensional quantitative structure-activity relationship (3D QSAR) method known as Comparative Molecular Field Analysis, CoMFA [Cramer III, R.D. et al., J. Am. Chem. Soc., 110 (1988) 5959]. We analyzed the data for 50 steroid inhibitors [Numazawa, M. et al., J. Med. Chem., 37 (1994) 2198, and references cited therein] assayed against androstenedione on human placental microsomes. An initial CoMFA resulted in a three-component model for log(l/Ki), with an explained variance r2 of 0.885, and a cross-validated q2 of 0.673. Chemometric studies were performed using GOLPE [Baroni, M. et al., Quant. Struct.-Act. Relatsh., 12 (1993) 9]. The CoMFA/GOLPE model is discussed in terms of robustness, predictivity, explanatory power and simplicity. After randomized exclusion of 25 or 10 compounds (repeated 25 times), the q2 for one component was 0.62 and 0.61, respectively, while r2 was 0.674. We demonstrate that the predictive r2 based on the mean activity (Ym) of the training set is misleading, while the test set Ym-based predictive r2 index gives a more accurate estimate of external predictivity. Using CoMFA, the observed differences in aromatase inhibition among C6-substituted steroids are rationalized at the atomic level. The CoMFA fields are consistent with known, potent inhibitors of aromatase, not included in the model. When positioned in the same alignment, these compounds have distinct features that overlap with the steric and electrostatic fields obtained in the CoMFA model. The presence of two hydrophobic binding pockets near the aromatase active site is discussed: a steric bulk tolerant one, common for C4, C6-alpha and C7-alpha substituents, and a smaller one at the C6-beta region.

Androstenedione↗