PubMed Health⌕ Search

Biomedical subjects

R W Pastor

Publications and source records attributed to R W Pastor.

At least 19 recordsLinked to original sources
PubMed1998

Solution structure and dynamics of linked cell attachment modules of mouse fibronectin containing the RGD and synergy regions: comparison with the human fibronectin crystal structure.

We report the three-dimensional solution structure of the mouse fibronectin cell attachment domain consisting of the linked ninth and tenth type III modules, mFnFn3(9,10). Because the tenth module contains the RGD cell attachment sequence while the ninth contains the synergy region, mFnFn3(9,10) has the cell attachment activity of intact fibronectin. Essentially complete signal assignments and approximately 1800 distance and angle restraints were derived from multidimensional heteronuclear NMR spectra. These restraints were used with a hybrid distance geometry/simulated annealing protocol to generate an ensemble of 20 NMR structures having no distance or angle violations greater than 0.3 A or 3 degrees. Although the beta-sheet core domains of the individual modules are well-ordered structures, having backbone atom rmsd values from the mean structure of 0.51(+/-0.12) and 0.40(+/-0.07) A, respectively, the rmsd of the core atom coordinates increases to 3.63(+/-1.41) A when the core domains of both modules are used to align the coordinates. The latter result is a consequence of the fact that the relative orientation of the two modules is not highly constrained by the NMR restraints. Hence, while structures of the beta-sheet core domains of the NMR structures are very similar to the core domains of the crystal structure of hFnFn3(9,10), the ensemble of NMR structures suggests that the two modules form a less extended and more flexible structure than the fully extended rod-like crystal structure. The radius of gyration, Rg, of mFnFn3(9,10) derived from small-angle neutron scattering measurements, 20.5(+/-0.5) A, agrees with the average Rg calculated for the NMR structures, 20.4 A, and is ca 1 A less than the value of Rg calculated for the X-ray structure. The values of the rotational anisotropy, D ||/D perpendicular, derived from an analysis of 15N relaxation data, range from 1.7 to 2.1, and are significantly less than the anisotropy of 2.67 predicted by hydrodynamic modeling of the crystal coordinates. In contrast, hydrodynamic modeling of the NMR coordinates yields anisotropies in the range of 1.9 to 2.7 (average 2.4(+/-0.2)), with NMR structures bent by more than 20 degrees relative the crystal structure having calculated anisotropies in best agreement with experiment. In addition, the relaxation parameters indicate that several loops in mFnFn3(9,10), including the RGD loop, are flexible on the nanosecond to picosecond time-scale. Taken together, our results suggest that, in solution, the limited set of interactions between the mFnFn3(9,10) modules position the RGD and synergy regions to interact specifically with cell surface integrins, and at the same time permit sufficient flexibility that allows mFnFn3(9,10) to adjust for some variation in integrin structure or environment.

Amino Acid Sequence↗
PubMed1997

Molecular dynamics simulation of unsaturated lipid bilayers at low hydration: parameterization and comparison with diffraction studies.

A potential energy function for unsaturated hydrocarbons is proposed and is shown to agree well with experiment, using molecular dynamics simulations of a water/octene interface and a dioleoyl phosphatidylcholine (DOPC) bilayer. The simulation results verify most of the assumptions used in interpreting the DOPC experiments, but suggest a few that should be reconsidered. Comparisons with recent results of a simulation of a dipalmitoyl phosphatidylcholine (DPPC) lipid bilayer show that disorder is comparable, even though the temperature, hydration level, and surface area/lipid for DOPC are lower. These observations highlight the dramatic effects of unsaturation on bilayer structure.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed1997

Length scales of lipid dynamics and molecular dynamics.

Following a brief overview of length scales and system size in computer simulation, it is demonstrated that a simulation sized lipid bilayer (typically a 50 x 50 A2 patch) is in the regime where stretching dominates undulation, while the reverse holds for a flaccid macroscopic membrane. Then it is estimated that current system sizes of membrane simulations must be increased by at least a factor of 10 before thermodynamic limits are approached for quantities such as surface tension.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed1996

On simulating lipid bilayers with an applied surface tension: periodic boundary conditions and undulations.

As sketched in Fig. 1, a current molecular dynamics computer simulation of a lipid bilayer fails to capture significant features of the macroscopic system, including long wavelength undulations. Such fluctuations are intrinsically connected to the value of the macroscopic (or thermodynamic) surface tension (cf. Eqs. 1 and 9; for a related treatment, see Brochard et al., 1975, 1976). Consequently, the surface tension that might be evaluated in an MD simulation should not be expected to equal the surface tension obtained from macroscopic measurements. Put another way, the largest of the three simulations presented here contained over 16,000 atoms and required substantial computer time to complete, but modeled a system of only 36 lipids per side. From this perspective it is not surprising that the system is not at the thermodynamic limit. An important practical consequence of this effect is that simulations with fluctuating area should be carried out with a nonzero applied surface tension (gamma 0 of Fig. 2) even when the macroscopic tension is zero, or close to zero. Computer simulations at fixed surface area, which can explicitly determine pressure anisotropy at the molecular level, should ultimately lend insight into the value of gamma 0, including its dependence on lipid composition and other membrane components. As we have noted and will describe further in separate publications (Feller et al., 1996; Feller et al., manuscript in preparation), surface tensions obtained from simulations can be distorted by inadequate initial conditions and convergence, and are sensitive to potential energy functions, force truncation methods, and system size; it is not difficult, in fact, to tune terms in the potential energy function so as to yield surface tensions close to zero. This is why parameters should be tested extensively on simpler systems, for example, monolayers. The estimates of gamma 0 that we have presented here should be regarded as qualitative, and primarily underscore the assertion that the surface tension of a microscopically flat, simulation-sized patch is significantly greater than zero. As the simulation cell length increases, the surface tension that would be evaluated (or should be applied) decreases; in the limit of micrometer-sized simulation cells, gamma would approach zero or its appropriate thermodynamic value. The theories presented here also imply that the estimation of bilayer surface tension from monolayer data should take the degree of flatness into account. These conclusions are independent of the precise values of parameters such as bending constants. In conclusion, from the simulator's perspective, the question "What is the surface tension of a bilayer?" is better phrased as "What is the value of the applied surface tension necessary to simulate a particular experimental system with a given number of lipids?". As we have shown, the answer to the second question varies, but it should not be assumed a priori to equal zero.

Biophysical Phenomena↗
PubMed1995

A method for characterizing transition concertedness from polymer dynamics computer simulations.

A statistical method based on classifying the transitions among a set of dihedral angles within an "energy transfer window" is developed, and used to analyze Brownian (BD) and molecular dynamics (MD) simulations of the acyl chains in a lipid bilayer, and MD of neat hexadecane. It is shown for the BD simulation that when a transition of the dihedral angle in the center of the chain occurs, a transition of a particular next nearest neighbor (or angle 2-apart) will follow concertedly with a probability of approximately 0.10 within a time window of approximately 3 ps. The MD bilayer simulations, which are based on a more flexible model of the hydrocarbon chains, yield corresponding concerted transition probabilities of approximately 0.083 and window sizes of 1-2 ps. An analysis of angles 4-apart yields concerted transition probabilities of 0.03 and 0.04 for the BD and MD bilayer simulations, respectively, and window sizes close to those of the corresponding 2-apart cases. Statistical hypothesis testing very strongly rejects the assertion that these follower transitions are occurring at random. Similar analysis reveals marginal or no evidence of concertedness between 1-apart (nearest neighbor) and between 3-apart dihedral angle transitions. The pattern of concertedness for hexadecane is qualitatively similar to that of the lipid chains, although concertedness is somewhat stronger for the 3-apart transitions and somewhat weaker for those 4-apart. Finally, it is suggested that the diffusion of small solute molecules in membranes is better facilitated by nonconcerted transitions, which are associated with relatively large displacements of the chains, than by concerted transitions, which do little to change the chain shape.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed1993

Molecular dynamics simulations of a lipid bilayer and of hexadecane: an investigation of membrane fluidity.

Molecular dynamics simulations of a fluid-phase dipalmitoyl phosphatidylcholine lipid bilayer in water and of neat hexadecane are reported and compared with nuclear magnetic resonance spin-lattice relaxation and quasi-elastic neutron scattering data. On the 100-picosecond time scale of the present simulations, there is effectively no difference in the reorientational dynamics of the carbons in the membrane interior and in pure hexadecane. Given that the calculated fast reorientational correlation times and the "microscopic" lateral diffusion of the lipids show excellent agreement with the experimental results, it is concluded that the apparently high viscosity of the membrane is more closely related to molecular interactions on the surface rather than in the interior.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed1992

Backbone dynamics of calmodulin studied by 15N relaxation using inverse detected two-dimensional NMR spectroscopy: the central helix is flexible.

The backbone dynamics of Ca(2+)-saturated recombinant Drosophila calmodulin has been studied by 15N longitudinal and transverse relaxation experiments, combined with 15N(1H) NOE measurements. Results indicate a high degree of mobility near the middle of the central helix of calmodulin, from residue K77 through S81, with order parameters (S2) in the 0.5-0.6 range. The anisotropy observed in the motion of the two globular calmodulin domains is much smaller than expected on the basis of hydrodynamic calculations for a rigid dumbbell type structure. This indicates that, for the purposes of 15N relaxation, the tumbling of the N-terminal (L4-K77) and C-terminal (E82-S147) lobes of calmodulin is effectively independent. A slightly shorter motional correlation time (tau c approximately 6.3 ns) is obtained for the C-terminal domain compared to the N-terminal domain (tau c approximately 7.1 ns), in agreement with the smaller size of the C-terminal domain. A high degree of mobility, with order parameters of approximately 0.5, is also observed in the loop that connects the first with the second EF-hand type calcium binding domain and in the loop connecting the third and fourth calcium binding domain.

Animals↗
PubMed1992

Langevin dynamics of peptides: the frictional dependence of isomerization rates of N-acetylalanyl-N'-methylamide.

The rate constant for the transition between the equatorial and axial conformations of N-acetylalanyl-N'-methylamide has been determined from Langevin dynamics (LD) simulations with no explicit solvent. The isomerization rate is maximum at collision frequency gamma = 2 ps-1, shows diffusive character for gamma greater than or equal to 10 ps-1, but does not approach zero even at gamma = 0.01 ps-1. This behavior differs from that found for a one-dimensional bistable potential and indicates that both collisional energy transfer with solvent and vibrational energy transfer between internal modes are important in the dynamics of barrier crossing for this system. It is suggested that conformational searches of peptides be carried out using LD with a collision frequency that maximizes the isomerization rate (i.e., gamma approximately 2 ps-1). This method is expected to be more efficient than either molecular dynamics in vacuo (which corresponds to LD with gamma = 0) or molecular dynamics in solvent (where dynamics is largely diffusive).

Alanine↗
PubMed1992

Conformational states of a TT mismatch from molecular dynamics simulation of duplex d (CGCGATTCGCG).

The TT mismatch region in duplex d (CGCGATTCGCG) was studied using a 500-ps molecular dynamics (MD) simulation in water, and a series of 1-ps MD simulations and energy minimizations in vacuum. The DNA maintained its duplex structure, although the mismatch region showed significantly higher flexibility than the GC regions. The predominant conformation in the 500-ps MD simulation involved an average -42 degrees propeller twist between T6 and T'6, and a -22 degree buckle between A5 and T'7. One hydrogen bond was formed between T6 and T'6, and another between T6 and the O2 of T'7, with both Watson-Crick hydrogen bonds between A5 and T'7 remaining intact. The minimizations resulted in conformations with the equivalent hydrogen-bonding pattern, as well as ones with "wobble pair" hydrogen bonds between T6 and T'6. However, the wobble pair conformation was found to be unstable in the water simulation.

Base Sequence↗
PubMed1991

Mean field stochastic boundary molecular dynamics simulation of a phospholipid in a membrane.

Computer simulations of phospholipid membranes have been carried out by using a combined approach of molecular and stochastic dynamics and a mean field based on the Marcelja model. First, the single-chain mean field simulations of Pastor et al. [(1988) J. Chem. Phys. 89, 1112-1127] were extended to a complete dipalmitoylphosphatidylcholine molecule; a 102-ns Langevin dynamics simulation is presented and compared with experiment. Subsequently, a hexagonally packed seven-lipid array was simulated with Langevin dynamics and a mean field at the boundary and with molecular dynamics (and no mean field) in the center. This hybrid method, mean field stochastic boundary molecular dynamics, reduces bias introduced by the mean field and eliminates the need for periodic boundary conditions. As a result, simulations extending to tens of nanoseconds may be carried out by using a relatively small number of molecules to model the membrane environment. Preliminary results of a 20-ns simulation are reported here. A wide range of motions, including overall reorientation with a nanosecond decay time, is observed in both simulations, and good agreement with NMR, IR, and neutron diffraction data is found.

Computer Simulation↗
PubMed1991

Model for the structure of the lipid bilayer.

A detailed model for the structure and dynamics of the interior of the lipid bilayer in the liquid crystal phase is presented. The model includes two classes of motion: (i) the internal dynamics of the chains, determined from Brownian dynamics simulations with a continuous version of the Marcelja mean-field potential, and (ii) noncollective reorientation (axial rotation and wobble) of the entire molecule, introduced by a cone model. The basic unit of the model is a single lipid chain with field parameters adjusted to fit the 2H order parameters and the frequency-dependent 13C NMR T1 relaxation times of dipalmitoyl phosphatidylcholine bilayers. The chain configurations obtained from the trajectory are used to construct a representation of the bilayer. The resulting lipid assembly is consistent with NMR, neutron diffraction, surface area, and density data. It indicates that a high degree of chain disorder and entanglement exists in biological membranes.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed1989

Theoretically determined three-dimensional structures for amphipathic segments of the HIV-1 gp41 envelope protein.

Three-dimensional computer models for two segments of the C terminus of gp41, the transmembrane AIDS envelope protein, which may form amphipathic alpha-helices, have been generated using structure prediction techniques combined with energy minimization and molecular dynamics simulations. Regions gp41(772-790) and gp41(828-848) of the HXB2 strain of HIV-1 display extraordinarily high hydrophobic moment maxima as alpha-helices and when in an antiparallel conformation exhibit charge complementarity, implying that they may bind with each other and associate with the membrane. The feasibility of this hypothesis was tested in a series of computer simulations of these peptides, extended by several residues to include additional charge pairing. Beginning with a trial structure in the form of antiparallel alpha-helices of segments 770-794 and 824-856, systematic axial rotations and displacements were used to generate alternative initial states. Molecular dynamics simulations with alpha-helical torsional restraints yielded several approximately cylindrical dimeric structures highly stabilized by numerous salt links and other hydrogen bonds. This suggests that these two regions may fold back on each other in antiparallel fashion to form a loop in the tertiary structure over residues 770-856, with the loop closed by membrane-associated amphipathic alpha-helices with charged sides facing each other. We speculate that such structures could aggregate to form channels or otherwise destabilize the membrane, thereby contributing to the cytopathic effects of the gp120-gp41 complex.

Amino Acid Sequence↗
PubMed1987

Theoretically determined three-dimensional structure for the repeating tetrapeptide unit of the circumsporozoite coat protein of the malaria parasite Plasmodium falciparum.

A model for the three-dimensional structure of the repeating Asn-Pro-Asn-Ala tetrapeptide of the immunodominant circumsporozoite protein of Plasmodium falciparum has been developed. A trial structure in the form of a type I beta turn with asparagine side chains hydrogen-bonded to the backbone peptide linkages was used as a starting point. A repeating oligomer of this trial structure was modeled using energy minimization and molecular dynamics computer simulations in conjunction with image boundary conditions. The most stable structure generated is a right-handed 12(38) helix, which is unlike any previously identified protein secondary structure. The helix has 12 residues per turn, corresponding to an angle of twist of 120 degrees per tetrapeptide unit, and a pitch of 4.95 A, corresponding to a rise of 1.65 A per tetrapeptide unit. It is highly stabilized by extensive hydrogen bonding, with each tetrapeptide unit acting as an acceptor for five hydrogen bonds and as a donor for five hydrogen bonds to residues in adjacent turns as well as having four weak internal hydrogen bonds. A number of nearly isoenergetic variations on the most stable structure that still retained the basic 12(38) helical motif were also discovered. The implications of these structures for vaccine development are discussed.

Animals↗
PubMed1986

A recognition site on synthetic helical oligonucleotides for monoclonal anti-native DNA autoantibody.

The binding site in native DNA for a murine monoclonal anti-DNA autoantibody was investigated by measurements of competitive binding of a series of synthetic helical oligonucleotides. The antibody bound to a (dG-dC)3 or (dG-dC)4 core in the center of a base-paired octadecanucleotide. Reactions of analogues containing modifications or substitutions at specific sites indicated that the antibody bound to portions of cytosine and guanine in the major groove, a limited region of the backbone, and the 2-amino group of one guanine in the minor groove. For these interactions to occur, the antibody combining site would straddle the backbone of one of the helical strands of DNA.

Animals↗