PubMed Health⌕ Search

Biomedical subjects

Angel E García

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

14 recordsLinked to original sources
PubMed2006

Local structure formation in simulations of two small proteins.

Massively parallel all-atom, explicit solvent molecular dynamics simulations were used to explore the formation and existence of local structure in two small alpha-helical proteins, the villin headpiece and the helical fragment B of protein A. We report on the existence of transient helices and combinations of helices in the unfolded ensemble, and on the order of formation of helices, which appears to largely agree with previous experimental results. Transient local structure is observed even in the absence of overall native structure. We also calculate sets of residue-residue pairs that are statistically predictive of the formation of given local structures in our simulations.

Computer Simulation↗
PubMed2006

Atomic-level description of amyloid beta-dimer formation.

All-atom simulations have been carried out on a monomer and dimer of the aggregation-prone fragment (16-22) of amyloid beta peptide, which is implicated in Alzheimer's disease. The replica exchange molecular dynamics method, which has been successfully applied to peptide folding, is utilized as a means to sample the configurational space with proper Boltzmann weighting so that the structural, motional, and thermodynamic description of self-assembly can be obtained. The free energy landscape showing the delicate balance between different monomer and dimer conformations is mapped along carefully chosen reaction coordinates. The canonical ensembles at 38 different temperatures are used to describe the thermodynamics and the relative stabilities of at least six different dimer conformations including that of parallel and antiparallel orientations. We also delineate the nature of the molecular forces that activate and stabilize these different dimer conformations as a function of temperature, especially as related to secondary structural propensity of monomer. We identify parallel loop dimer conformations that are stabilized due to specific interactions with water molecules.

Amyloid beta-Peptides↗
PubMed2005

Molecular dynamics study of water penetration in staphylococcal nuclease.

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.

Algorithms↗
PubMed2005

Helix-coil transition of alanine peptides in water: force field dependence on the folded and unfolded structures.

The force fields used in classical modeling studies are semiempirical in nature and rely on their validation by comparison of simulations with experimental data. The all-atom replica-exchange molecular dynamics (REMD) methodology allows us to calculate the thermodynamics of folding/unfolding of peptides and small proteins, and provides a way of evaluating the reliability of force fields. We apply the REMD to obtain equilibrium folding/unfolding thermodynamics of a 21-residue peptide containing only alanine residues in explicit aqueous solution. The thermodynamics of this peptide is modeled with both the OPLS/AA/L and the A94/MOD force fields. We find that the helical content and the values for the helix propagation and nucleation parameters for this alanine peptide are consistent with measurements on similar peptides and with calculations using the modified AMBER force field (A94/MOD). The nature of conformations, both folded and unfolded, that contributes to the helix-coil transition profile, however, is quite different between these two force fields.

Alanine↗
PubMed2004

Reversible temperature and pressure denaturation of a protein fragment: a replica exchange molecular dynamics simulation study.

We determine the reversible folding-unfolding of the C-terminal (41-56) fragment of protein G as a function of density and temperature using replica exchange molecular dynamics simulations. We employ a total of 253 replicas, covering the temperature range between 320 and 515 K and the density range between 0.96 and 1.16 g cm(-3). Using the root mean square deviation from the folded structure as a quantitative measure, we are able to obtain the fraction of folded states, and can thus establish the free energy difference between the folded and the unfolded states of the protein fragment as a function of temperature and pressure. For the pressure denaturation the weakening of the hydrophobic interaction between the bulky side chains is found to be crucial at lower temperatures, leading to an apparent destabilization of the folded backbone structure at elevated pressures.

Biophysics↗
PubMed2004

Nature of structural inhomogeneities on folding a helix and their influence on spectral measurements.

Extensive conformational sampling and calculations of vibrational coupling provide a quantitative basis for the structurally inhomogeneous spectra of the amide unit in aqueous solutions containing folded and unfolded state distributions of helices. Replica exchange molecular dynamics simulations of the capped helical peptide, AA(AAKAA)(3)AAY, is carried out over a range of temperatures, where the system populates the folded and unfolded states. This sampling defines a set of ensembles that characterizes the conformational variability for configurations identified by their fraction of helical content. The effects of hydrogen bonding, both internal and external (with water), and the coupling between amide-I modes are computed as a function of temperature and helical content. End-to-end distance and coupling distributions are also computed. The solvent H-bonding, which is present at all temperatures, shifts the amide-I band toward lower frequency compared with the unsolvated band. Upon thermal denaturation of the peptide, the amide-I band shifts to higher frequency because the increase in solvent H-bonding fails to compensate for the loss in internal (helical) H-bonds. The loss of uniformity of the mode coupling along the helix at higher temperatures accounts for the well-known thermal broadening of the amide IR spectrum. The calculated inhomogeneities of segments of the peptide predict experimental properties of isotope-edited helices.

Amides↗
PubMed2003

Folding a protein in a computer: an atomic description of the folding/unfolding of protein A.

We study the folding mechanism of a three-helix bundle protein at atomic resolution, including effects of explicit water. Using replica exchange molecular dynamics we perform enough sampling over a wide range of temperatures to obtain the free energy, entropy, and enthalpy surfaces as a function of structural reaction coordinates. Simulations were started from different configurations covering the folded and unfolded states. Because many transitions between all minima at the free energy surface are observed, a quantitative determination of the free energy barriers and the ensemble of configurations associated with them is now possible. The kinetic bottlenecks for folding can be determined from the thermal ensembles of structures on the free energy barriers, provided the kinetically determined transition-state ensembles are similar to those determined from free energy barriers. A mechanism incorporating the interplay among backbone ordering, sidechain packing, and desolvation arises from these calculations. Large Phi values arise not only from native contacts, which mostly form at the transition state, but also from contacts already present in the unfolded state that are partially destroyed at the transition.

Biophysical Phenomena↗
PubMed2003

Simulation of the folding equilibrium of alpha-helical peptides: a comparison of the generalized Born approximation with explicit solvent.

We compare simulations using the generalized Born/surface area (GB/SA) implicit solvent model with simulations using explicit solvent (transferable intermolecular potential 3 point, TIP3P) to test the GB/SA algorithm. We use the replica exchange molecular dynamics method to sample the conformational phase space of two alpha-helical peptides, A21 and the Fs, by using two different classical potentials and both water models. We find that when using GB/SA: (i) A21 is predicted to be more helical than the Fs peptide at all temperatures; (ii) the native structure of the Fs peptide is predicted to be a helical bundle instead of a single helix; and (iii) the persistence length and most probable end-to-end distance are too large in the unfolded state when compared against the explicit solvent simulations. We find that the potential of mean force in the phi(psi) plane is markedly different in the two solvents, making the two simulated peptides respond differently when the backbone torsions are perturbed. A fit of the temperature melting curves obtained in these simulations to a Lifson-Roig model finds that the GB/SA model has an unphysically large nucleation parameter, whereas the explicit solvent model produces values similar to experiment.

Algorithms↗
PubMed2003

Interplay among tertiary contacts, secondary structure formation and side-chain packing in the protein folding mechanism: all-atom representation study of protein L.

Experimental and theoretical results suggest that, since proteins are energetically minimally frustrated, the native fold, or topology, plays a primary role in determining the structure of the transition state ensemble and on-pathway intermediate states in protein folding. Although the central role of native state topology in determining the folding mechanism is thought to be a quite general result-at least for small two-state folding proteins-there are remarkable exceptions. Recent experimental findings have shown that topology alone cannot always determine the folding mechanism, and demonstrated that the balance between topology and energetics is very delicate. This balance seems to be particularly critical in proteins with a highly symmetrical native structure, such as proteins L and G, which have similar native structure topology but fold by different mechanisms. Simplified, C(alpha)-atom only protein models have shown not be sufficient to differentiate these mechanisms. An all-atom Gō model provides a valuable intermediate model between structurally simplified protein representations and all-atom protein simulations with explicit/implicit solvent descriptions. We present here a detailed study of an all-atom Gō-like representation of protein L, in close comparison with the experimental results and with the results obtained from a simple C(alpha)-atom representation of the same protein. We also perform simulations for protein G, where we obtain a folding mechanism in which the protein symmetry is broken exactly in the opposite way to protein L as has been observed experimentally. A detailed analysis for protein L also shows that the role of specific residues is correctly and quantitatively reproduced by the all-atom Gō model over almost the entire protein.

Bacterial Proteins↗
PubMed2003

Role of backbone hydration and salt-bridge formation in stability of alpha-helix in solution.

We test molecular level hypotheses for the high thermal stability of alpha-helical conformations of alanine-based peptides by performing detailed atomistic simulations of a 20-amino-acid peptide with explicit treatment of water. To assess the contribution of large side chains to alpha-helix stability through backbone desolvation and salt-bridge formation, we simulate the alanine-rich peptide, Ac-YAEAAKAAEAAKAAEAAKAF-Nme, referred to as the EK peptide, that has three pairs of "i, i + 3" glutamic acid(-) and lysine(+) substitutions. Efficient configurational sampling of the EK peptide over a wide temperature range enabled by the replica exchange molecular dynamics technique allows characterization of the stability of alpha-helix with respect to heat-induced unfolding. We find that near ambient temperatures, the EK peptide predominately samples alpha-helical configurations with 80% fractional helicity at 300 K. The helix melts over a broad range of temperatures with melting temperature, T(m), equal to 350 K, that is significantly higher than the T(m) of a 21-residue polyalanine peptide, A(21). Salt-bridges between oppositely charged Glu(-) and Lys(+) side chains can, in principle, provide thermal stability to alpha-helical conformers. For the specific EK peptide sequence, we observe infrequent formation of Glu-Lys salt-bridges (with approximately 10-20% probability) and therefore we conclude that salt-bridge formation does not contribute significantly to the EK peptide's helical stability. However, lysine side chains are found to shield specific "i, i + 4" backbone hydrogen bonds from water, indicating that large side-chain substituents can play an important role in stabilizing alpha-helical configurations of short peptides in aqueous solution through mediation of water access to backbone hydrogen bonds. These observations have implications on molecular engineering of peptides and biomolecules in the design of their thermostable variants where the shielding mechanism can act in concert with other factors such as salt-bridge formation, thereby increasing thermal stability considerably.

Alanine↗
PubMed2003

Peptide folding simulations.

Developments in the design of small peptides that mimic proteins in complexity, recent advances in nanosecond time-resolved spectroscopy methods to study peptides and the development of modern, highly parallel simulation algorithms have come together to give us a detailed picture of peptide folding dynamics. Two newly implemented simulation techniques, parallel replica dynamics and replica exchange molecular dynamics, can now describe directly from simulations the kinetics and thermodynamics of peptide formation, respectively. Given these developments, the simulation community now has the tools to verify and validate simulation protocols and models (forcefields).

Computer Simulation↗
PubMed2002

Alpha-helical stabilization by side chain shielding of backbone hydrogen bonds.

We study atomic models of the thermodynamics of the structural transition of peptides that form alpha-helices. The effect of sequence variation on alpha-helix formation for alanine-rich peptides, Ac-Ala21-methyl amide (A21) and Ac-A5 (AAARA)3A-methyl amide (Fs peptide), is investigated by atomic simulation studies of the thermodynamics of the helix-coil transition in explicit water. The simulations show that the guanidinium group in the Arg side chains in the Fs peptide interacts with the carbonyl group four amino acids upstream in the chain and desolvates backbone hydrogen bonds. This desolvation can be directly correlated with a higher probability of hydrogen bond formation. We find that Fs has higher helical content than A21 at all temperatures. A small modification in the amber force field reproduces the experimental helical content and helix-coil transition temperatures for the Fs peptide.

Algorithms↗
PubMed2002

Protein folding mediated by solvation: water expulsion and formation of the hydrophobic core occur after the structural collapse.

The interplay between structure-search of the native structure and desolvation in protein folding has been explored using a minimalist model. These results support a folding mechanism where most of the structural formation of the protein is achieved before water is expelled from the hydrophobic core. This view integrates water expulsion effects into the funnel energy landscape theory of protein folding. Comparisons to experimental results are shown for the SH3 protein. After the folding transition, a near-native intermediate with partially solvated hydrophobic core is found. This transition is followed by a final step that cooperatively squeezes out water molecules from the partially hydrated protein core.

Animals↗