PubMed Health⌕ Search

Biomedical subjects

S Gnanakaran

Publications and source records attributed to S Gnanakaran.

10 recordsLinked to original sources
PubMed2006

Clade-specific differences between human immunodeficiency virus type 1 clades B and C: diversity and correlations in C3-V4 regions of gp120.

Current knowledge of human immunodeficiency virus type 1 envelope (Env) glycoprotein structure and function is based on studies of clade B viruses. We present evidence of sequence and structural differences in viral glycoprotein gp120 between clades B and C. In clade C, the C3 region alpha2-helix exhibits high sequence entropy at the polar face but maintains its amphipathicity, whereas in clade B it accommodates hydrophobic residues. The V4 hypervariable domain in clade C is shorter than that in clade B. Generally, shorter V4 loops are incompatible with a glycine occurring in the alpha2-helix in clade C, an intriguing association that could be exploited to inform Env immunogen design.

Genetic Variation↗
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↗
PubMed2006

Consistent patterns of change during the divergence of human immunodeficiency virus type 1 envelope from that of the inoculated virus in simian/human immunodeficiency virus-infected macaques.

We have analyzed changes to proviral Env gp120 sequences and the development of neutralizing antibodies (NAbs) during 1 year of simian/human immunodeficiency virus SHIV-89.6P infection in 11 Macaca nemestrina macaques. Seven macaques had significant env divergence from that of the inoculum, and macaques with greater divergence had higher titers of homologous NAbs. Substitutions in sequons encoding potential N-linked glycosylation sites (PNGs) were among the first to be established, although overall the total number of sequons did not increase significantly. The majority (19 of 23) of PNGs present in the inoculum were conserved in the sequences from all macaques. Statistically significant variations in PNGs occurred in multiple macaques within constrained regions we term "hot spots," resulting in the selection of sequences more similar to the B consensus. These included additions on V1, the N-terminal side of V4, and the outer region of C2. Complex mutational patterns resulted in convergent PNG shifts in V2 and V5. Charge changes in Env V1V2, resulting in a net acidic charge, and a proline addition in V5 occurred in several macaques. Molecular modeling of the 89.6P sequence showed that the conserved glycans lie on the silent face of Env and that many are proximal to disulfide bonds, while PNG additions and shifts are proximal to the CD4 binding site. Nonsynonymous-to-synonymous substitution ratios suggest that these changes result from selective pressure. This longitudinal and cross-sectional study of mutations in human immunodeficiency virus (HIV) env in the SHIV background provides evidence that there are more constraints on the configuration of the glycan shield than were previously appreciated.

Amino Acid Sequence↗
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↗
PubMed2005

Simulations of the pressure and temperature unfolding of an alpha-helical peptide.

We study by molecular simulations the reversible folding/unfolding equilibrium as a function of density and temperature of a solvated alpha-helical peptide. We use an extension of the replica exchange molecular dynamics method that allows for density and temperature Monte Carlo exchange moves. We studied 360 thermodynamic states, covering a density range from 0.96 to 1.14 g.cm(-3) and a temperature range from 300 to 547.6 K. We simulated 10 ns per replica for a total simulation time of 3.6 micros. We characterize the structural, thermodynamic, and hydration changes as a function of temperature and pressure. We also calculate the compressibility and expansivity of unfolding. We find that pressure does not affect the helix-coil equilibrium significantly and that the volume change upon pressure unfolding is small and negative (-2.3 ml/mol). However, we find significant changes in the coordination of water molecules to the backbone carbonyls. This finding predicts that changes in the chemical shifts and IR spectra with pressure can be due to changes in coordination and not only changes in the helical content. A simulation of the IR spectrum shows that water coordination effects on frequency shifts are larger than changes due to elastic structural changes in the peptide.

Carbon↗
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 of a highly conserved diverging turn motif from the SH3 domain.

Recent NMR structural characterization studies showed that a seven-residue segment (FKKGERL) from the src SH3 domain adopts the nativelike diverging type II beta-turn in aqueous solution in support of the prediction based on the I-sites library of sequence structural motifs. We study the conformational variability and folding/unfolding thermodynamics of this peptide in explicit solvent using replica-exchange molecular dynamics simulations, which greatly enhances the sampling of the conformational space. This peptide samples three main free energy basins (nativelike, intermediate, and unfolded) separated by small barriers. The nativelike basin is fractionally populated (DeltaG(300K) = 0.4 kcal/mol) with structures that satisfy a subset of the NMR-derived constraints. The intrinsic stability of the diverging turn is examined in relationship to the nature of three specific contacts: a turn-hydrogen bond, a mainchain-to-sidechain hydrogen bond, and an end-to-end hydrophobic contact. We have carried out simulations of mutants at the highly conserved GE positions in the sequence. The mutation E5D destabilizes the isolated diverging turn motif, contrary to the observation that this mutation stabilizes the fyn SH3 domain. The G4T mutation also destabilizes the isolated diverging turn; however, the extent of destabilization is smaller than that of the reverse mutation in the drk SH3.

Amino Acid Motifs↗
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↗
PubMed2001

Conformational preferences and vibrational frequency distributions of short peptides in relation to multidimensional infrared spectroscopy.

Molecular dynamics simulations of the structural distributions and the associated amide-I vibrational modes are carried out for dialanine peptide in water and carbon tetrachloride. The various manifestations in nonlinear-infrared spectroscopic experiments of the distributions of conformations of solvated dialanine are examined. The two-dimensional infrared (2D-IR) spectrum of dialanine exhibits the coupling between the amide oscillators and the correlations of the frequency fluctuations. An internally hydrogen-bonded conformation exists in CCl(4) but not in H(2)O where two externally hydrogen-bonded forms are preferred. Simulations of solvated dialanine show how the 2D-IR spectra expose the underlying structural distributions and dynamics that are not deducible from linear-infrared spectra. In H(2)O the 2D-IR shows cross-peaks from large coupling in the alpha-helical conformer and an elongated higher frequency diagonal peak, reflecting the broader distribution of structures for the more flexible acetyl end. In CCl(4), the computed cross-peak portion of the 2D-IR shows evidence of two amide-I transitions in the high-frequency region which are not apparent from the diagonal peak profile. The vibrational frequency inhomogeneity of the amide-I band arises from fluctuations of the instantaneous normal modes of these conformers rather than the shifts induced by hydrogen bonding. The simulation shows that there are correlations between fluctuations of the acetyl and amino end frequencies in H(2)O that arise from mechanical coupling and not from hydrogen bonding at the two ends of the molecule. The angular relationships between the two amide units which also show up in 2D-IR were computed, and spectral manifestations of them are discussed. The simulations also permit a calculation of the rate of energy transfer from one side of the molecule to the other. From these calculations, 2D-IR spectroscopy in conjunction with simulations is seen to be a promising tool for determining dynamics of structure changes in dipeptides.

Dipeptides↗