help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

Originally published as Biophys J. BioFAST on December 13, 2004.
doi:10.1529/biophysj.104.048611
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.104.048611v1
88/3/2003    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Merlino, A.
Right arrow Articles by Mazzarella, L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Merlino, A.
Right arrow Articles by Mazzarella, L.
Biophysical Journal 88:2003-2012 (2005)
© 2005 The Biophysical Society

Open Interface and Large Quaternary Structure Movements in 3D Domain Swapped Proteins: Insights from Molecular Dynamics Simulations of the C-Terminal Swapped Dimer of Ribonuclease A

Antonello Merlino * {dagger}, Marc Antoine Ceruso {ddagger}, Luigi Vitagliano § and Lelio Mazzarella * {dagger} §

* Centro Interdipartimentale Ricerca e Management, Complesso Ristrutturato S. Andrea delle Dame, 80138, Naples, Italy; {dagger} Dipartimento di Chimica, Università degli Studi di Napoli "Federico II", 80125 Naples, Italy; {ddagger} Department of Physiology and Biophysics, Mount Sinai School of Medicine, New York, New York USA; and § Istituto di Biostrutture e Bioimmagini, CNR, 80134 Naples, Italy

Correspondence: Address reprint requests to Lelio Mazzarella, Dipartimento di Chimica, Università degli Studi di Napoli "Federico II", Complesso Universitario di Monte Sant'Angelo, Via Cynthia, 80126 Napoli, Italy. Tel.: 39-081674279; Fax: 39-081674090; E-mail: lelio.mazzarella{at}unina.it or mazzarella{at}chemistry.unina.it.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Bovine pancreatic ribonuclease (RNase A) forms two three-dimensional (3D) domain swapped dimers. Crystallographic investigations have revealed that these dimers display completely different quaternary structures: one dimer (N-dimer), which presents the swapping of the N-terminal helix, is characterized by a compact structure, whereas the other (C-dimer), which is stabilized by the exchange of the C-terminal end, shows a rather loose assembly of the two subunits. The dynamic properties of monomeric RNase A and of the N-dimer have been extensively characterized. Here, we report a molecular dynamics investigation carried out on the C-dimer. This computational experiment indicates that the quaternary structure of the C-dimer undergoes large fluctuations. These motions do not perturb the proper folding of the two subunits, which retain the dynamic properties of RNase A and the N-dimer. Indeed, the individual subunits of the C-dimer display the breathing motion of the ß-sheet structure, which is important for the enzymatic activity of pancreatic-like ribonucleases. In contrast to what has been observed for the N-dimer, the breathing motion of the two subunits of the C-dimer is not coupled. This finding suggests that the intersubunit communications in a 3D domain swapped dimer strongly rely on the extent of the interchain interface. Furthermore, the observation that the C-dimer is endowed with a high intrinsic flexibility holds interesting implications for the specific properties of 3D domain swapped dimers. Indeed, a survey of the quaternary structures of the other 3D domain swapped dimers shows that large variations are often observed when the structural determinations are conducted in different experimental conditions. The 3D domain swapping phenomenon coupled with the high flexibility of the quaternary structure may be relevant for protein-protein recognition, and in particular for the pathological aggregations.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Many proteins perform their biological function in oligomeric states (Goodsell and Olson, 2000Go). Oligomeric proteins present several advantages when compared to the monomeric counterparts. Multimeric proteins often show increased thermal and/or chemical stability, and modulated activity through cooperative effects. Furthermore, oligomerization also leads to the formation of large protein assembly with minimal amounts of DNA. This aspect minimizes the number of errors that may occur in transcription.

The aggregation of monomers to form an oligomeric protein is a complex molecular process whose details are still obscure. Generally, the transition from monomeric to oligomeric states relies on the generation of novel chain-chain interfaces, which are progressively stabilized by evolution. Recently, a great attention has been devoted to a different mechanism denoted as three-dimensional (3D) domain swapping (Bennett et al., 1995Go; Schlunegger et al., 1997Go). The aggregation via 3D domain swapping is based on the mutual exchange of a domain of one chain with the same domain of a second chain (Bennett et al., 1994Go, 1995Go; Schlunegger et al., 1997Go; Newcomer, 2001Go, 2002Go; Liu and Eisenberg, 2002Go). Accordingly, the structure of the oligomeric protein is stabilized by contacts that preexist in the momomeric form (Closed-interface or C-interface) as well as by novel chain-chain interfaces (Open-interface or O-interface) (Bennett et al., 1995Go; Schlunegger et al., 1997Go).

Among 3D domain swapped proteins, pancreatic-like ribonucleases have played an important role, as they were the subjects of several pioneering studies (see, for example, Crestfield et al., 1962Go). Structural evidences on the occurrence of 3D domain swapping in ribonucleases were derived more than 20 years ago for bovine seminal ribonuclease (BS-RNase) (Capasso et al., 1983Go). This enzyme also exhibits the intriguing characteristic of forming two equilibrium isomers: the swapped (M x M) (Mazzarella et al., 1993Go) and the nonswapped (M = M) (Piccoli et al., 1992Go; Berisio et al., 2003Go) dimer. More recently, structural investigations carried out on the two dimers formed by bovine pancreatic ribonuclease (RNase A) in a variety of experimental conditions (Crestfield et al., 1962Go; Park and Raines, 2000Go; Gotte et al., 2003Go) have demonstrated that they adopt two different types of 3D domain swapping (Liu et al., 1998Go, 2001Go). Indeed, in one dimer (hereafter denoted as C-Dimer) the swapped fragment is the C-terminal ß-strand (Liu et al., 2001Go), whereas in the other one (N-dimer) the exchanged region is the N-terminal {alpha}-helix (Liu et al., 2001Go). As a consequence, the N- and the C-dimer display completely different shapes as the N-dimer presents a very compact structure (Liu et al., 1998Go), whereas C-dimer only shows a very loose O-interface (Liu et al., 2001Go) (Fig. 1).



View larger version (36K):
[in this window]
[in a new window]
 
FIGURE 1  Ribbon diagram of RNase dimers.

 
Despite the increasing interest for 3D domain swapping, the available data on the dynamic properties of 3D domain swapped dimers are rather limited. The first extensive characterizations of the dynamic behavior of a 3D domain swapped dimer compared to its monomeric counterpart have been reported only very recently (Sekijima et al., 2003Go; Japelj et al., 2004Go; Merlino et al., 2004bGo). In previous investigations, we demonstrated that the functional collective motions of the ß-sheet structure (Merlino et al., 2002Go; Vitagliano et al., 1998Go, 2002Go), characteristic of pancreatic-like ribonucleases (Merlino et al., 2003Go), are coupled in the two subunits of the N-dimer, thus suggesting a possible mechanism for the regulation of its catalytic activity (Merlino et al., 2004bGo).

Here, we report the results of molecular dynamics studies of the C-dimer in water. The simulations reveal that the C-Dimer exhibits extensive movements of the two subunits. These findings, together with an examination of the quaternary structure variations of other 3D domain swapped dimers indicate that 3D domain swapping may be an effective mechanism to generate oligomers endowed with an unusual flexibility. Furthermore, the coupling of the breathing motions of the two subunits observed in the N-Dimer, is not detected in the C-dimer. These results underline the role of the compactness of the O-interface in the intersubunit communications.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
System and definitions
Following the definition of Eisenberg and co-workers (Bennett et al., 1995Go; Schlunegger et al., 1997Go), in a 3D domain swapped dimer, the interface that preexists in the monomeric form is termed Closed interface (C-interface) (Fig. 1). On the other hand, the interface that forms upon dimerization is denoted as Open interface (O-interface). The loop connecting the body with its swapped fragment of each subunit, termed hinge peptide, undergoes a dramatic structural transition in the 3D domain swapping process.

The C-dimer is characterized by an elongated shape and by a rather loose O-interface (Fig. 1) (Liu et al., 2001Go). The swapped fragment is represented by the C-terminal ß-strand (residues 116–124). Each subunit is composed by three {alpha}-helices and a V-shaped ß-sheet constituted by two arms that will be hereafter denoted as V1 and V2. V1 consists of residues 61–63, 71–75, and 105–111 of one subunit and of residues 116–124 of the other, whereas V2 embodies residues 42–46, 82–87, and 96–101. The three helices encompass residues 3–13 (helix I), 25–35 (helix II), and 50–60 (helix III). The C-interface consists of two strands (105–111 and 116–124) of the V1 arm. The hinge peptide and the main body of the protein are composed by residues 112–115 and 1–111, respectively. The O-interface is limited to the pairing of the two hinge peptides (Fig. 1). The local structure of this region corresponds to a short ß-sheet interrupted by Pro-114. Because some catalytic residues belong to the C-terminal ß-strand region, in the C-dimer, two different structural units containing all the functional residues, each composed by the body of one subunit and the C-terminal strand of the other, may be identified.

Simulation procedure
The molecular dynamics simulations were performed with the program GROMACS (van der Spoel et al., 1994Go). The starting model for the simulation was the molecule composed of the A/B chains taken from the Protein Data Bank (PDB) entry 1F0V (Liu et al., 2001Go). The dimer was immersed in a rectangular box containing simple-point-charge water molecules (Berendsen et al., 1981Go). The ionization state of charged residues was set to mimic a neutral pH environment (Merlino et al., 2002Go). To neutralize the system, eight water molecules were replaced by chloride counterions. For comparative purposes, the molecular dynamics studies were carried out by using the same protocol applied for the simulation of RNase A (Merlino et al., 2002Go) and the N-dimer (Merlino et al., 2004bGo). Briefly, to allow relaxation of solvent molecules, the energy of the system was preliminarily minimized by keeping fixed the protein atoms. The resulting system, containing 58,264 atoms, was submitted to 20 ps of molecular dynamics at 300 K. The energy of the system was then minimized without restraints, before starting constant temperature molecular dynamics at 300 K. The overall timescale of the simulation was 3 ns. All bond lengths were constrained using the SHAKE algorithm (Ryckaert et al., 1977Go). Nonbonded cutoffs of 10 Å for Lennard-Jones and 13 Å for Coulomb potentials were used. A dielectric constant of {varepsilon} = 1 and an integration time step of 0.002 ps were used. The coordinates were saved at every 0.2 ps.

To verify the reliability of the findings of this MD study, a second independent 3-ns simulation has been carried out at 300 K using a different set of initial velocities. The results of the two simulations are in good agreement. This represents a good indication of the completeness of the conformational sampling of each simulation. Therefore, only the results obtained from one simulation will be described.

Analysis of the trajectories
To assess the quality of the simulations, the trajectories were checked using GROMACS routines and the Metaphor software suite (M. A. Ceruso, unpublished data). To determine coupling between intra- and interstructural unit motions in the C-dimer, atomic correlation matrices were calculated and normalized as

where {Delta}Ri and {Delta}Rj are the displacements from the mean positions of the ith atom and the jth atom, and <> represent the time average over the equilibrated portion of the trajectory. Cij values vary in the interval –1 < Cij < 1. To eliminate the contribution of relative subunit motions, interstructural unit correlation matrices were generated after an independent fit of each structural unit trajectory on the same reference. To compare the motions of the two structural units, a difference atomic correlation matrix was built with elements defined as

where Cij(A) and Cij(B) are the values of Cij calculated for the two structural units. {Delta}Cij values vary in the range –2 < {Delta}Cij < 2.

Essential degrees of freedom were extracted from the trajectories according to the essential dynamics method (Amadei et al., 1993Go) that is based on the construction of the covariance matrix of the coordinate fluctuations. The covariance matrix is diagonalized to obtain the eigenvectors and eigenvalues that provide information about correlated motion throughout the protein. The eigenvectors represent the directions of motion and the eigenvalues the amount of motion along each eigenvector. The eigenvectors are then sorted according to their eigenvalues in descending order. Usually, the first 10–20 eigenvectors suffice to describe almost all conformational substates accessible to the protein (Ceruso et al., 1999Go).

The internal fluctuations of the C-dimer and N-dimer structural units and of the RNase A were compared by evaluating their principal subspaces (first 10 eigenvectors). As already done to compare the concerted motions of several protein systems (Ceruso et al., 1999Go, 2003Go; Grottesi et al., 2002Go; Merlino et al., 2002Go, 2003Go), the overlap of the principal subspaces of each pair of proteins was evaluated by estimating root mean square inner product (RMSIP).

The RMSIP is defined as:

where and are the ith and jth eigenvectors from the set a and set b, respectively.

The RMSIP was also used to ascertain the convergence in the essential space between two halves of the equilibrated trajectories and to evaluate the similarity of the two trajectories.

The analyses of the internal fluctuations were performed omitting the peptide regions 15–22 and 112–115, whose local structure changes in the dimerization process that leads to the formation of the N- and the C-dimer, respectively.

To ensure that the structural and dynamic properties of N-dimer derived from the simulations are free from initial nonequilibrium effects, all of the analyses were performed in the portion of trajectories from 1200 to 3000 ps. In this region the essential subspace converged as revealed by calculating the RMSIP value between the two halves of the trajectory (Table 1).


View this table:
[in this window]
[in a new window]
 
TABLE 1  Root mean square inner products for projection of the first 10 C{alpha} eigenvectors of a set onto the first 10 of another

 
Selection and analysis of domain swapped dimers endowed with high mobility
A survey of the PDB database (Berman et al., 2002Go) was carried out to identify 3D domain swapped dimers that exhibit significantly different quaternary structures as a consequence of mutations, ligand binding, and/or differences in crystal packing. The search was performed on the January, 2004 release of the PDB, which contains ~22,000 protein structures.

To evaluate the rearrangement of the structural units occurring in different structures of a given domain swapped dimer, a stepwise superimposing procedure was carried out using the programs LSQKAB of the CCP4 suite (Bailey, 1994Go) and DynDom (Hayward and Lee, 2002Go). In particular, as already described elsewhere (Vitagliano et al., 1998Go; Merlino et al., 2004cGo), we evaluated the rotation required to overlay the second structural unit, after the optimal superimposition of the first.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
C-Dimer shows extensive movements of its quaternary structure
To evaluate the stability of the MD simulation, several time-dependent properties of the C-dimer were analyzed. As shown in Fig. 2 A (black line), the C-dimer shows high fluctuations of its quaternary structure during the simulation. Indeed, the overall C{alpha} root mean square deviation (RMSD) of the structure during the simulation versus the starting x-ray model fluctuates between 2 and 6 Å.



View larger version (29K):
[in this window]
[in a new window]
 
FIGURE 2  (A) Root mean square positional deviations from the x-ray starting structure computed on C{alpha} atoms of the dimer (black line) and of the bodies of the two subunits (gray lines) as function of time, (B) distances between the centers of mass of the two structural units as function of time, and (C) snapshots from the simulation.

 
It should be noted, however, that the individual structural units (see Methods for the definition) of the protein are very stable, with RMSD close to 2 Å (Fig. 2 A, gray lines). The stability of each unit is confirmed by the analysis of the time evolution of the secondary structure elements and of the total number of hydrogen bonds. The overall stability of the system is also corroborated by the trend of the potential interaction energy, which is virtually constant during the whole trajectory.

These data indicate that the C-dimer undergoes large variations of its quaternary structure that do not perturb the proper folding of each subunit. These variations are also evidenced by the time evolution of the distance between the centers of mass of the two structural units of the dimer (Fig. 2 B). Interestingly, although the system significantly alters its quaternary structure during the simulation the total solvent accessible surface area display a limited fluctuation with a mean value of 725.7 ± 3.5 nm2, which is in close agreement with that (727.2 nm2) derived from the x-ray structure of C-Dimer (PDB code 1F0V). Similarly, limited fluctuations of hydrophobic (mean value 278.1 ± 1.7 nm2) and hydrophilic (mean value 446.2 ± 2.2 nm2) solvent exposed regions are observed. Finally, it is worth mentioning that the C-dimer preserves the overall twofold symmetry during the simulation. For most of the simulation, the two structural units are related by a rotation of 168–180°. Only occasionally (between 2700 and 2800 ps) this rotation is close to 150°. These findings, along with a visual inspection of a number of structures extracted from the trajectory, suggest that the fluctuations of the quaternary structure essentially correspond to a hinge-bending motion of the two structural units (Fig. 2 C). A significant twist of the C-dimer structure is only observed in the region 2700–2800 ps (Fig. 2 C).

The motions of two structural units of the C-dimer are similar
Fig. 3 A shows the C{alpha} root mean square fluctuations (RMSF) of the two structural units of the C-dimer. The two structural units show comparable fluctuations. The atomic correlation matrices are also extremely similar, as shown by the difference atomic correlation matrix, which does not exhibit any significant ({Delta}Cij > 0.7) feature (Fig. 3 B). The concerted motions were also investigated by using the root mean square inner product (RMSIP; see Methods for the definition) between the first 10 eigenvectors derived from the diagonalization of the covariance matrix of the C-dimer structural unit trajectories (Table 1). The high value of the RMSIP (0.69) represents a further indication that the two structural units of the dimer undergo the same motions.



View larger version (35K):
[in this window]
[in a new window]
 
FIGURE 3  (A) Root mean square fluctuations of the C-dimer amino acids calculated from the simulation in comparison with those of RNase A ({diamond}). The symbols • and {triangleup} refer to the A and B subunits, respectively. (B) Difference intrastructural unit correlation matrix; (C) root mean square fluctuations of the C-dimer amino acids in comparison with those of the N-dimer ({diamondsuit} and {circ}, for the A and B subunits, respectively). (D) Comparison between the A subunit B-factor derived from the x-ray structure ({blacktriangleup}) and the root mean square fluctuations of the A subunit C{alpha} atoms of the C-dimer ({blacksquare}). (E) Motions of one of the two structural units along the first principal eigenvector. A grayscale is used to represent the motions in a film-like fashion. Vectors are used as qualitative indicators of the direction of the C{alpha} atom motions.

 
The motions of the C-dimer structural units are similar to those of the N-dimer and to RNase A
The fluctuations of the two structural units of the C-dimer are very similar to those observed for the N-dimer (Fig. 3 C) (Merlino et al., 2004bGo) and for RNase A (Merlino et al., 2002Go), although in the last case the region 50–60 displays a slightly higher mobility (Fig. 3 A). In addition, as observed for the N-dimer (Merlino et al., 2004bGo) and RNase A (Merlino et al., 2002Go), the fluctuations of the C{alpha} atom of the C-dimer structural units derived from the simulations are in good agreement with the experimentally determined B-factors (Fig. 3 D).

The principal subspace derived for the C-dimer structural units largely overlaps that of RNase A and of the N-dimer units, as shown by the RMSIP values reported in Table 1. As expected on the basis of this similarity, the analysis of the first eigenvectors indicates that the ß-sheet regions of the C-Dimer units exhibit the breathing motions (Fig. 3 E), which has been well characterized in experimental and theoretical studies carried out on RNase A (Berisio et al., 2002Go; Merlino et al., 2002Go; Vitagliano et al., 2002Go), N-dimer (Merlino et al., 2004bGo), as well as on other members of the pancreatic-like superfamily (Vitagliano et al., 1998Go; Merlino et al., 2003Go, 2004bGo).

This is in line with the general finding that the dynamic properties strongly depend on the architecture of the proteins (Ceruso et al., 1999Go, 2003Go, 2004Go; Merlino et al., 2003Go, 2004aGo). However, it should be pointed out that, in the case of the C-dimer, the ß-sheet scaffold of each unit of the protein has a composite nature, because one of the ß-strands comes from the partner chain. Thus, the presence of a swapped strand in the ß-sheet region of the protein does not prevent the breathing motion characteristic of pancreatic-like ribonucleases.

Breathing motions of two C-dimer structural units are not coupled
The analysis of the interstructural unit correlated motions of the C-dimer was conducted by evaluating atomic correlation matrices. In Fig. 4, A and B, positive and negative correlations are reported. As can be noted from the figure, the correlation matrices for C-dimer do not show any significant peaks (Cij is always <0.3), indicating that the breathing motions of the two structural units (see above) are not correlated. It is worth mentioning that a completely different picture was derived for the N-dimer (Merlino et al., 2004bGo), where a strong correlation of these motions, likely transmitted through the extended ß-sheet structure of the O-interface, was observed. Therefore, the N- and the C-dimer display different trends for the coupling of the structural unit breathing motions, although the local and global flexibility of individual units is very similar.



View larger version (39K):
[in this window]
[in a new window]
 
FIGURE 4  Positive (A) and negative (B) interstructural unit correlation matrices. As described in the Methods section, in this analysis the residues 16–22 and 112–115 have been omitted.

 
Motions at the C- and O-interface
The C-interface, which corresponds to the interactions between the two ß-strands 105–111 and 116–124, is well preserved during the simulation. Indeed, the hydrogen bonds established among the main-chain atoms of the two strands are thoroughly conserved, as observed for RNase A (Merlino et al., 2002Go). As a representative example of this general trend, the hydrogen bond formed by the O atom of the Val-116 of the first subunit and N atom of Glu-111 of the other is shown in Fig. 5 A. The dynamic behavior of the active site region of the C-dimer is also similar (data not shown) to that reported for RNase A (Vitagliano et al., 2000Go; Merlino et al., 2002Go).



View larger version (40K):
[in this window]
[in a new window]
 
FIGURE 5  Interactions in the C- and O-interface regions are: (A) distance between O atom of Val-116 of the first subunit and N atom of Gln-111 of the other as function of time; (B) distance between O atom of Asn-113 of the first subunit and N atom of Asn-113 of the other as function of time; (C) distance between N atom of Asn-113 of the first subunit and O atom of Asn-113 of the other as function of time; and (D) distance between Nd2 atom of Asn-113 of the first subunit and Od1 atom of Asn-113 of the other as function of time and vice versa.

 
The analysis of the dynamic properties of the O-interface provides interesting information. This region consists of a short antiparallel ß-sheet formed by residues 111–113 of one subunit and 111–113 of the other (Liu et al., 2001Go). Notably, the main-chain hydrogen bonds, which characterize the antiparallel ß-sheet, are rather well preserved during the simulation, despite the large overall motions of the dimer: as an example, the hydrogen bonds between the N and O atoms of the Asn-113 of the two subunits are shown in Fig. 5, B and C. In the analysis of the dynamic properties of the O-interface, a particular attention was devoted to the hydrogen bond between the side chains of Asn-113 of the two subunits. Eisenberg and co-workers (Liu et al., 2001Go) suggested that this Asn-Asn H-bond could be considered as an example of the polar zipper interactions occurring in amyloid fiber formation. As shown in Fig. 5 D, this hydrogen bond is present in several phases of the trajectory. However, the specific geometry of the two Asn side chains changes during the trajectory, with the formation of the hydrogen bond that involves both the O{delta}1 and the N{delta}2 atoms of the two residues.

Extensive movements in 3D domain swapped dimers
The large interstructural unit motions of the C-dimer prompted us to carry out a survey on other 3D domain swapped dimers to verify whether extensive movements of the quaternary structures can be considered a common feature of 3D domain swapping proteins. This hypothesis was checked by analyzing the structures of 3D domain swapped dimers, which were either reported in different crystal forms or solved by NMR. In particular, we evaluated the rotation required to overlay the second structural units, after the optimal superimposition of the first (Table 2).


View this table:
[in this window]
[in a new window]
 
TABLE 2  Quaternary structure variation of 3D domain swapping dimers

 
As shown in the table, many 3D domain swapped dimers, whose structure has been determined in different contexts, present large quaternary structure variations. Extensive variations occur for wild-type cyanovirin between the crystal and the solution state (Yang et al., 1999Go; Bewley and Clore, 2000Go; Barrientos et al., 2002Go). Large quaternary structure variations are also compatible with the NMR data on human cystatin (Staniforth et al., 2001Go). A remarkable reorganization of the quaternary structure are also observed comparing the structures of the transcription factor NusG from Aquifex aeolicus derived from two crystal forms (Knowlton et al., 2003Go), and two dimers of the Grb2-SH2-Ac-pYYNV complex present in the same asymmetric unit (Schiering et al., 2000Go). Smaller but significant reorganization of the quaternary structure are observed for BS-RNase (Vitagliano et al., 1998Go; Merlino et al., 2004cGo), whose quaternary structure is tightly constrained by two intersubunit disulfide bridges (Mazzarella et al., 1993Go). It is also worth noting that single point mutations significantly alter the structure of protein L (O'Neill et al., 2001Go) and cyanovirin (Botos et al., 2002Go), and to a minor extent of the T4 endonuclease VII (Raaijmakers et al., 1999Go, 2001Go). On the other hand, little variations are observed in two structures of the nerve factor binding domain of the TrkA receptor (Ultsch et al., 1999Go; Robertson et al., 2001Go) and minor alterations are produced by single point mutations in rat Cd2 (Murray et al., 1998Go) and by ligand binding in the Fha domain of cell cycle checkpoint protein Chfr (Stavridi et al., 2002Go). Finally, huge variations are observed comparing homologous 3D domain swapped dimers. Indeed, a rotation of 130° is observed when the structures of the N-dimer (Liu et al., 1998Go) and of the human pancreatic ribonuclease dimer (Canals et al., 2001Go) are compared. Similarly, large reorganizations are detected by superimposing the ligand binding domains of TrkA, TrkB, and TrkC receptor structures (Ultsch et al., 1999Go), 41–44% sequence identity, and comparing the structures of human (Janowski et al., 2001Go) and chicken cystatin (Staniforth et al., 2001Go).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
The structural characterization of RNase A dimers carried out by Eisenberg and co-workers (Liu et al., 1998Go, 2001Go) has highlighted several fascinating features: i), a protein can swap both N- and C-terminal ends; ii), the overall quaternary structure of the two dimers is completely different; and iii), the 3D domain swapping may be important for protein misfolding and aggregation. With the aim of gaining dynamic information on these intriguing systems, we have undertaken a thorough investigation of these dimers by molecular dynamics simulations. In previous articles, we studied the dynamic properties of the RNase A monomer (Merlino et al., 2002Go) and of the N-dimer (Merlino et al., 2004bGo). The results of these studies have shown that the monomer and the individual structural units of the N-dimer are endowed with a similar breathing motion of the ß-sheet structure. Interestingly, the motions of the structural units of the N-dimer are highly coupled (Merlino et al., 2004bGo). Here we report a molecular dynamic study on the C-terminal swapped dimer of RNase A. As far as local and global motions of the individual structural units are concerned, the findings of these investigations are in line with the previously reported results (Merlino et al., 2002Go, 2004bGo). Indeed, the active site and the overall breathing motion of the units strictly follow the trends that emerged from the studies on RNase A and the N-dimer. As for the case of the N-terminal end exchange (Merlino et al., 2004bGo), the swapping of the C-terminus does not affect the local dynamic behavior of the enzyme, despite the fact that the swapped strand is part of the ß-sheet structure.

However, the MD investigation on the C-dimer shows distinctive features. In contrast to the data reported for the N-Dimer (Merlino et al., 2004bGo), the ß-sheet motions of the two structural units are not coupled in the C-dimer. In other words, from the dynamic point of view the two structural units of the C-dimer are independent. This finding may be ascribed to the different assembly of the two dimers. In fact, the N- and the C-dimer are characterized by different O-interface areas, being that the N-dimer is more compact. Therefore, the larger interface of the N-dimer is able to transmit interunit communications, whereas the loose connection of the two structural units in the C-dimer does not. This observation suggests that large intersubunit surfaces are essential for signal transmission and could be essential for allosteric modulation.

The results of the MD simulation also show that the loose O-interface allows large fluctuations of the quaternary structure of the C-dimer. Specifically, the RMSD of the whole dimer from the crystallographic structure fluctuates in a wide range (2.0–6.0 Å), despite the fact that main bodies of the two subunits present almost constant RMSD values (~2.0 Å). A different trend was observed for the N-dimer, which presents a rather constant RMSD (~3.5 Å) from the starting x-ray structure in the equilibrated region of the trajectory (Merlino et al., 2004bGo). The large fluctuations of the C-dimer suggest that the x-ray structure is likely conditioned by crystal packing. It should be recalled that the C-dimer crystals were grown in the presence of the substrate analog dCpdG, which links different molecules of the dimer in the crystal. This interaction, important for packing, favors the crystallization process (Liu et al., 2001Go).

An analysis of the structures deposited in PDB shows that large quaternary structure movements are frequent among 3D domain swapped dimers. Particularly evident is the conformational variability of cyanovirin, which displays three distinct quaternary structures (Yang et al., 1999Go; Bewley and Clore, 2000Go; Barrientos et al., 2002Go). Indeed, large variations are not only observed when crystal and NMR models are compared, but also between structures derived from different crystal forms (Yang et al., 1999Go; Bewley and Clore, 2000Go; Barrientos et al., 2002Go). Similarly, the lowest-energy NMR conformers of cystatin display very diverse subunit orientations (Staniforth et al., 2001Go). Very recently, Janowski et al. (2004)Go have reported that the different domain swapping dimers contained in the asymmetric unit of a modified form of cystatin show a significant degree of variation, as a consequence of the high flexibility of the open interface.

In 3D domain swapped dimers, large structural variations may also be induced by single point mutations, mainly occurring in the hinge peptide region. An uncommon quaternary structure variability is also observed among homologous proteins that generate 3D domain swapped dimers.

The unusual motions of the 3D domain swapped dimers reported here are due to the peculiar way in which they are assembled. An intrinsic difference can be identified between nonswapped and swapped dimers. In the former, dimerization requires the formation of a novel interface surface, which should be large enough to prevent dissociation. On the contrary, the stability of a swapped dimer may be ensured by the C-interface, which preexists in the monomer, and does not necessarily depend on the contacts between the main bodies of the two subunits. In other words, large interfaces between the main bodies of the two subunits, which are important for providing stability in nonswapped dimer, are optional in the swapping ones. Notably, the two swapped dimers of RNase A represent two emblematic cases: indeed, the N-dimer presents a rather rigid O-interface with limited structural flexibility, on the other hand the C-dimer, which is characterized by a loose O-interface, undergoes large structural fluctuations.

Therefore, it can be suggested that the swapped dimers endowed with loose O-interfaces may undergo large fluctuations, if the hinge peptide is sufficiently flexible. It should be noted that the inspection of the interacting surfaces observed in 3D domain swapped dimers showed that the surface area of the O-interface is often quite limited (Schlunegger et al., 1997Go) with few hydrogen bonds and salt bridges (Sinha et al., 2001Go). Consequently, the absence of strongly stabilizing interactions between the two units does not constrain their motion, by producing a lower barrier than would need to be overcome in such extensive motions.

The high flexibility may also be an important factor that contributes to the formation of "loose" 3D swapping dimers and higher oligomers. Indeed, it can be speculated that the entropy costs associated with the oligomerization processes are somewhat limited. In other words, in a dimer, the loose O-interface allows a large conformational freedom, which reduces the unfavorable entropic contribution to dimer formation.

Several biological roles have been proposed for the 3D domain swapping. Among these, it has been suggested that it can be important for the modulation of catalytic actions (Vitagliano et al., 1999Go), for immunosuppressive and antitumoral activities (Vescia and Tramontano, 1981Go; Cafaro et al., 1995Go) and for the formation of the amyolid fibrils (Liu et al., 2001Go). Here we show that high flexibility is a distinctive property of several 3D domain swapped dimers. This finding suggests the intriguing idea that the coupling of 3D domain swapping with the high flexibility could be important for a number of biological processes that rely on protein-protein recognition step, including the pathological aggregations.


    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
This work was supported by Italian Ministero dell'Istruzione, Università e Ricerca (PRIN 2002) and Regione Campania ‘Legge 5’.


    FOOTNOTES
 
Abbreviations used: 3D, three-dimensional; MD, molecular dynamics; RMSD, root mean square deviations; RMSF, root mean square fluctuations; RMSIP, root mean square inner product; RNase A, bovine pancreatic ribonuclease; C-dimer, C-terminal swapped dimer of ribonuclease A; N-dimer, N-terminal swapped dimer of ribonuclease A.

Submitted on June 29, 2004; accepted for publication December 2, 2004.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS
 DISCUSSION
 ACKNOWLEDGEMENTS
 REFERENCES
 
Amadei, A., A. B. Linssen, and H. J. C. Berendsen. 1993. Essential dynamics of proteins. Proteins. 17:412–425.[CrossRef][Medline]

Bailey, S. 1994. The CCP4 suite: programs for protein crystallography. Acta Crystallogr. D. 50:760–763.[CrossRef][Medline]

Barrientos, L. G., J. M. Louis, I. Botos, T. Mori, Z. Han, B. R. O'Keefe, M. R. Boyd, A. Wlodawer, and A. M. Gronenborn. 2002. The domain-swapped dimer of cyanovirin-N is in a metastable folded state: reconciliation of X-ray and NMR structures. Structure. 10:673–686.[Medline]

Bennett, M. J., S. Choe, and D. Eisenberg. 1994. Refined structure of dimeric diphtheria toxin at 2.0 A resolution. Protein Sci. 3:1444–1463.[Abstract]

Bennett, M. J., M. P. Schlunegger, and D. Eisenberg. 1995. 3D domain swapping: a mechanism for oligomer assembly. Protein Sci. 4:2455–2468.[Abstract]

Berendsen, H. J. C., J. P. M. Postma, W. F. van Gusteren, and J. Hermans. 1981. Interaction models for water in relation to protein hydration. Dordrecht:Reidel:331–342.

Berisio, R., F. Sica, C. De Lorenzo, A. Di Fiore, R. Piccoli, A. Zagari, and L. Mazzarella. 2003. Crystal structure of the dimeric unswapped form of bovine seminal ribonuclease. FEBS Lett. 554:105–110.[CrossRef][Medline]

Berisio, R., F. Sica, V. S. Lamzin, K. S. Wilson, A. Zagari, and L. Mazzarella. 2002. Atomic resolution structures of ribonuclease A at six pH values. Acta Crystallogr. D. 58:441–450.[CrossRef][Medline]

Berman, H. M., T. Battistuz, T. N. Bhat, W. F. Bluhm, P. E. Bourne, K. Burkhardt, Z. Feng, G. L. Gilliland, L. Iype, S. Jain, P. Fagan, J. Marvin, D. Padilla, V. Ravichandran, B. Schneider, N. Thanki, H. Weissig, J. D. Westbrook, and C. Zardecki. 2002. The Protein Data Bank. Acta Crystallogr. D. 58:899–907.[CrossRef][Medline]

Bewley, C. A., and G. M. Clore. 2000. Determination of the relative orientation of the two halves of the domain-swapped dimer of cyanoviryn-N in solution using dipolar couplings and rigid body minimization. J. Am. Chem. Soc. 122:6009–6016.[CrossRef]

Botos, I., T. Mori, L. K. Cartner, M. R. Boyd, and A. Wlodawer. 2002. Domain-swapped structure of a mutant of cyanovirin-N. Biochem. Biophys. Res. Commun. 294:184–190.[CrossRef][Medline]

Cafaro, V., C. De Lorenzo, R. Piccoli, A. Bracale, M. R. Mastronicola, A. Di Donato, and G. D'Alessio. 1995. The antitumor action of seminal ribonuclease and its quaternary conformations. FEBS Lett. 359:31–34.[CrossRef][Medline]

Canals, A., J. Pous, A. Guasch, A. Benito, M. Ribo, M. Vilanova, and M. Coll. 2001. The structure of an engineered domain-swapped ribonuclease dimer and its implications for the evolution of proteins toward oligomerization. Structure. 9:967–976.[Medline]

Capasso, S., F. Giordano, C. A. Mattia, L. Mazzarella, and A. Zagari. 1983. Refinement of the structure of bovine seminal ribonuclease. Biopolymers. 22:327–332.[CrossRef][Medline]

Ceruso, M. A., A. Amadei, and A. Di Nola. 1999. Mechanics and dynamics of B1 domain of protein G: role of packing and surface hydrophobic residues. Protein Sci. 8:147–160.[Abstract]

Ceruso, M. A., A. Grottesi, and A. Di Nola. 2003. Dynamic effects of mutations within two loops of cytochrome c551 from Pseudomonas aeruginosa. Proteins. 50:222–229.[CrossRef][Medline]

Ceruso, M. A., X. Periole, and H. Weinstein. 2004. Molecular dynamics simulations of transducin: interdomain and front to back communication in activation and nucleotide exchange. J. Mol. Biol. 338:469–481.[CrossRef][Medline]

Crestfield, A. M., W. H. Stein, and S. Moore. 1962. On the aggregation of bovine pancreatic ribonuclease. Arch. Biochem. Biophys. 1:217–222.

Goodsell, D. S., and A. Olson. 2000. Structural symmetry and protein function. Annu. Rev. Biophys. Biomol. Struct. 29:105–153.[CrossRef][Medline]

Gotte, G., F. Vottariello, and M. Libonati. 2003. Thermal aggregation of ribonuclease A. A contribution to understanding the role of the domain swapping in protein aggragation. J. Biol. Chem. 278:10763–10769.[Abstract/Free Full Text]

Grottesi, A., M. A. Ceruso, A. Colosimo, and A. Di Nola. 2002. Molecular dynamics study of a hyperthermophilic and a mesophilic rubredoxin. Proteins. 46:287–294.[CrossRef][Medline]

Hayward, S., and R. A. Lee. 2002. Improvement in the analysis of domain motions in proteins from conformational change: DynDom version 1.50. J. Mol. Graph. Model. 21:181–183.[CrossRef][Medline]

Janowski, R., M. Abrahamson, A. Grubb, and M. Jaskolski. 2004. Domain swapping in N-truncated human cystatin C. J. Mol. Biol. 341:151–160.[CrossRef][Medline]

Janowski, R., M. Kozak, E. Jankowska, Z. Grzonka, A. Grubb, M. Abrahamson, and M. Jaskolski. 2001. Human cystatin C, an amyloidogenic protein, dimerizes through three-dimensional domain swapping. Nat. Struct. Biol. 8:316–320.[CrossRef][Medline]

Japelj, B., J. P. Waltho, and R. Jerala. 2004. Comparison of backbone dynamics of monomeric and domain-swapped stefin A. Proteins. 54:500–512.[CrossRef][Medline]

Knowlton, J. R., M. Bubunenko, M. Andrykovitch, W. Guo, K. M. Routzahn, D. S. Waugh, D. L. Court, and X. Ji. 2003. A spring-loaded state of NusG in its functional cycle is suggested by X-ray crystallography and supported by site-directed mutants. Biochemistry. 42:2275–2281.[CrossRef][Medline]

Liu, Y., and D. Eisenberg. 2002. 3D domain swapping: as domains continue to swap. Protein Sci. 11:1285–1299.[Abstract/Free Full Text]

Liu, Y., G. Gotte, M. Libonati, and D. Eisenberg. 2001. A domain-swapped RNase A dimer with implications for amyloid formation. Nat. Struct. Biol. 8:211–214.[CrossRef][Medline]

Liu, Y., P. J. Hart, M. P. Schlunegger, and D. Eisenberg. 1998. The crystal structure of a 3D domain-swapped dimer of RNase A at a 2.1-A resolution. Proc. Natl. Acad. Sci. USA. 95:3437–3442.[Abstract/Free Full Text]

Mazzarella, L., S. Capasso, D. Demasi, G. Di Lorenzo, C. A. Mattia, and A. Zagari. 1993. Bovine seminal ribonuclease: structure at 1.9 Å resolution. Acta Crystallogr. D. 49:389–402.[CrossRef][Medline]

Merlino, A., G. Graziano, and L. Mazzarella. 2004a. Structural and dynamic effect of {alpha}-helix deletion in Sso7d: implication for protein thermal stability. Proteins. 57:692–701.[CrossRef][Medline]

Merlino, A., L. Vitagliano, M. A. Ceruso, A. Di Nola, and L. Mazzarella. 2002. Global and local motions in ribonuclease A: a molecular dynamics study. Biopolymers. 65:274–283.[CrossRef][Medline]

Merlino, A., L. Vitagliano, M. A. Ceruso, and L. Mazzarella. 2003. Subtle functional collective motions in pancreatic-like ribonucleases:from ribonuclease A to angiogenin. Proteins. 53:101–110.[CrossRef][Medline]

Merlino, A., L. Vitagliano, M. A. Ceruso, and L. Mazzarella. 2004b. Dynamic properties of the N-terminal swapped dimer of ribonuclease A. Biophys. J. 86:2383–2391.[Abstract/Free Full Text]

Merlino, A., L. Vitagliano, F. Sica, A. Zagari, and L. Mazzarella. 2004c. Population shift versus induced fit: the case of bovine seminal ribonuclease swapping dimer. Biopolymers. 73:689–695.[CrossRef][Medline]

Murray, A. J., J. G. Head, J. J. Barker, and R. L. Brady. 1998. Engineering an intertwined form of CD2 for stability and assembly. Nat. Struct. Biol. 5:778–782.[CrossRef][Medline]

Newcomer, M. E. 2001. Trading places. Nat. Struct. Biol. 8:282–284.[CrossRef][Medline]

Newcomer, M. E. 2002. Protein folding and three-dimensional domain swapping: a strained relationship? Curr. Opin. Struct. Biol. 12:48–53.[CrossRef][Medline]

O'Neill, J. W., D. E. Kim, K. Johnsen, D. Baker, and K. Y. Zhang. 2001. Single-site mutations induce 3D domain swapping in the B1 domain of protein L from Peptostreptococcus magnus. Structure. 9:1017–1027.[Medline]

Park, C., and R. T. Raines. 2000. Dimer formation by a monomeric protein. Protein Sci. 9:2026–2033.[Abstract]

Piccoli, R., M. Tamburrini, G. Piccialli, A. Di Donato, A. Parente, and G. D'Alessio. 1992. The dual-mode quaternary structure of seminal RNase. Proc. Natl. Acad. Sci. USA. 89:1870–1874.[Abstract/Free Full Text]

Raaijmakers, H., I. Toro, R. Birkenbihl, B. Kemper, and D. Suck. 2001. Conformational flexibility in T4 endonuclease VII revealed by crystallography: implications for substrate binding and cleavage. J. Mol. Biol. 308:311–323.[CrossRef][Medline]

Raaijmakers, H., O. Vix, I. Toro, S. Golz, B. Kemper, and D. Suck. 1999. X-ray structure of T4 endonuclease VII: a DNA junction resolvase with a novel fold and unusual domain-swapped dimer architecture. EMBO J. 18:1447–1458.[CrossRef][Medline]

Robertson, A. G., M. J. Banfield, S. J. Allen, J. A. Dando, G. G. Mason, S. J. Tyler, G. S. Bennett, S. D. Brain, A. R. Clarke, R. L. Naylor, G. K. Wilcock, R. L. Brady, and D. Dawbarn. 2001. Identification and structure of the nerve growth factor binding site on TrkA. Biochem. Biophys. Res. Commun. 282:131–141.[CrossRef][Medline]

Ryckaert, J. P., G. Ciccotti, and H. J. C. Berendsen. 1977. Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics simulations of n-alkanes. J. Comput. Phys. 23:327–341.[CrossRef]

Schiering, N., E. Casale, P. Caccia, P. Giordano, and C. Battistini. 2000. Dimer formation through domain swapping in the crystal structure of the Grb2–SH2-Ac-pYVNV complex. Biochemistry. 39:13376–13382.[CrossRef][Medline]

Schlunegger, M. P., M. J. Bennett, and D. Eisenberg. 1997. Oligomer formation by 3D domain swapping: a model for protein assembly and misassembly. Adv. Protein Chem. 50:61–122.[Medline]

Sekijima, M., C. Motono, S. Yamasaki, K. Kaneko, and Y. Akiyama. 2003. Molecular dynamics simulation of dimeric and monomeric forms of human prion protein: insight into dynamics and properties. Biophys. J. 85:1176–1185.[Abstract/Free Full Text]

Sinha, N., C. J. Tsai, and R. Nussinov. 2001. A proposed structural model for amyloid fibril elongation: domain swapping forms an interdigitating beta-structure polymer. Protein Eng. 14:93–103.[Abstract/Free Full Text]

Staniforth, R. A., S. Giannini, L. D. Higgins, M. J. Conroy, A. M. Hounslow, R. Jerala, C. J. Craven, and J. P. Waltho. 2001. Three-dimensional domain swapping in the folded and molten-globule states of cystatins, an amyloid-forming structural superfamily. EMBO J. 20:4774–4781.[CrossRef][Medline]

Stavridi, E. S., Y. Huyen, I. R. Loreto, D. M. Scolnick, T. D. Halazonetis, N. P. Pavletich, and P. D. Jeffrey. 2002. Crystal structure of the FHA domain of the Chfr mitotic checkpoint protein and its complex with tungstate. Structure. 10:891–899.[Medline]

Ultsch, M. H., C. Wiesmann, L. C. Simmons, J. Henrich, M. Yang, D. Reilly, S. H. Bass, and A. M. de Vos. 1999. Crystal structures of the neurotrophin-binding domain of TrkA, TrkB and TrkC. J. Mol. Biol. 290:149–159.[CrossRef][Medline]

van der Spoel, D., R. van Druner, and H. J. C. Berendsen. 1994. GROningen MAchine for Chemical Simulation. Department of Biophysical Chemistry, BIOSON Research Institute, Groningen, The Netherlands.

Vescia, S., and D. Tramontano. 1981. Antitumoral action of bovine seminal ribonuclease. Mol. Cell. Biochem. 36:125–128.[CrossRef][Medline]

Vitagliano, L., S. Adinolfi, A. Riccio, F. Sica, A. Zagari, and L. Mazzarella. 1998. Binding of a substrate analog to a domain swapping protein: X-ray structure of the complex of bovine seminal ribonuclease with uridylyl(2',5')adenosine. Protein Sci. 7:1691–1699.[Abstract]

Vitagliano, L., S. Adinolfi, F. Sica, A. Merlino, A. Zagari, and L. Mazzarella. 1999. A potential allosteric subsite generated by domain swapping in bovine seminal ribonuclease. J. Mol. Biol. 293:569–577.[CrossRef][Medline]

Vitagliano, L., A. Merlino, A. Zagari, and L. Mazzarella. 2000. Productive and nonproductive binding to ribonuclease A: X-ray structure of two complexes with uridylyl(2',5')guanosine. Protein Sci. 9:1217–1225.[Abstract]

Vitagliano, L., A. Merlino, A. Zagari, and L. Mazzarella. 2002. Reversible substrate-induced domain motions in ribonuclease A. Proteins. 46:97–104.[CrossRef][Medline]

Yang, F., C. A. Bewley, J. M. Louis, K. R. Gustafson, M. R. Boyd, A. M. Gronenborn, G. M. Clore, and A. Wlodawer. 1999. Crystal structure of cyanovirin-N, a potent HIV-inactivating protein, shows unexpected domain swapping. J. Mol. Biol. 288:403–412.[CrossRef][Medline]




This article has been cited by other articles:


Home page
J. Biol. Chem.Home page
A. Merlino, L. Mazzarella, A. Carannante, A. D. Fiore, A. D. Donato, E. Notomista, and F. Sica
The Importance of Dynamic Effects on the Enzyme Activity: X-RAY STRUCTURE AND MOLECULAR DYNAMICS OF ONCONASE MUTANTS
J. Biol. Chem., May 6, 2005; 280(18): 17953 - 17960.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
biophysj.104.048611v1
88/3/2003    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Merlino, A.
Right arrow Articles by Mazzarella, L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Merlino, A.
Right arrow Articles by Mazzarella, L.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS